Showing 1–50 of 198 results for author: Bolatto, A D

JWST Observations of Starbursts: Cold Clouds and Plumes Launching in the M82 Outflow

Authors: Deanne B. Fisher, Alberto D. Bolatto, John Chisholm, Drummond Fielding, Rebecca C. Levy, Elizabeth Tarantino, Martha L. Boyer, Serena A. Cronin, Laura A. Lopez, J. D. Smith, Danielle A. Berg, Sebastian Lopez, Sylvain Veilleux, Paul P. van der Werf, Torsten Böker, Leindert A. Boogaard, Laura Lenkić, Simon C. O. Glover, Vicente Villanueva, Divakara Mayya, Thomas S. -Y. Lai, Daniel A. Dale, Kimberly L. Emig, Fabian Walter, Monica Relaño , et al. (6 additional authors not shown)

Abstract: In this paper we study the filamentary substructure of 3.3 $μ$m PAH emission from JWST/NIRCam observations in the base of the M82 star-burst driven wind. We identify plume-like substructure within the PAH emission with widths of $\sim$50 pc. Several of the plumes extend to the edge of the field-of-view, and thus are at least 200-300 pc in length. In this region of the outflow, the vast majority (… ▽ More In this paper we study the filamentary substructure of 3.3 $μ$m PAH emission from JWST/NIRCam observations in the base of the M82 star-burst driven wind. We identify plume-like substructure within the PAH emission with widths of $\sim$50 pc. Several of the plumes extend to the edge of the field-of-view, and thus are at least 200-300 pc in length. In this region of the outflow, the vast majority ($\sim$70\%) of PAH emission is associated with the plumes. We show that those structures contain smaller scale "clouds" with widths that are $\sim$5-15 pc, and they are morphologically similar to the results of "cloud-crushing" simulations. We estimate the cloud-crushing time-scales of $\sim$0.5-3 Myr, depending on assumptions. We show this time scale is consistent with a picture in which these observed PAH clouds survived break-out from the disk rather than being destroyed by the hot wind. The PAH emission in both the midplane and the outflow is shown to tightly correlate with that of Pa$α$ emission (from HST/NICMOS data), at the scale of both plumes and clouds, though the ratio of PAH-to-Pa$α$ increases at further distances from the midplane. Finally, we show that the outflow PAH emission is suppressed in regions of the M82 wind that are bright in X-ray emission. Overall, our results are broadly consistent with a picture in which cold gas in galactic outflows is launched via hierarchically structured plumes, and those small scale clouds are more likely to survive the wind environment when collected into the larger plume structure. △ Less

Submitted 6 May, 2024; originally announced May 2024.

Comments: Submitted to MNRAS

Modeling Ionized Gas in the Small Magellanic Cloud: The Wolf-Rayet Nebula N76

Authors: Elizabeth Tarantino, Alberto D. Bolatto, Rémy Indebetouw, Mónica Rubio, Karin M. Sandstrom, J. -D T. Smith, Daniel Stapleton, Mark Wolfire

Abstract: We present Cloudy modeling of infrared emission lines in the Wolf-Rayet (WR) nebula N76 caused by one of the most luminous and hottest WR stars in the low metallicity Small Magellanic Cloud. We use spatially resolved mid-infrared Spitzer/IRS and far-infrared Herschel/PACS spectroscopy to establish the physical conditions of the ionized gas. The spatially resolved distribution of the emission allow… ▽ More We present Cloudy modeling of infrared emission lines in the Wolf-Rayet (WR) nebula N76 caused by one of the most luminous and hottest WR stars in the low metallicity Small Magellanic Cloud. We use spatially resolved mid-infrared Spitzer/IRS and far-infrared Herschel/PACS spectroscopy to establish the physical conditions of the ionized gas. The spatially resolved distribution of the emission allows us to constrain properties much more accurately than using spatially integrated quantities. We construct models with a range of constant hydrogen densities between n$_H$ = 4 - 10 cm$^{-3}$ and a stellar wind-blown cavity of 10 pc which reproduces the intensity and shape of most ionized gas emission lines, including the high ionization lines [OIV] and [NeV], as well as [SIII], [SIV], [OIII], and [NeIII]. Our models suggest that the majority of [SiII] emission (91%) is produced at the edge of the HII region around the transition between ionized and atomic gas while very little of the [CII] emission (<5%) is associated with the ionized gas. The physical conditions of N76 are characterized by a hot HII region with a maximum electron temperature of T$_e$ ~ 24,000 K, electron densities that range from n$_e$ ~ 4 to 12 cm$^{-3}$, and high ionization parameters of log(U) ~ -1.15 to -1.77. By analyzing a low metallicty WR nebula with a single ionization source, this work gives valuable insights on the impact WR stars have to the galaxy-integrated ionized gas properties in nearby dwarf galaxies. △ Less

Submitted 11 April, 2024; originally announced April 2024.

Comments: 27 pages, 10 figures, 2 tables. Accepted in the Astrophysical Journal

DUVET: Resolved direct metallicity measurements in the outflow of starburst galaxy NGC 1569

Authors: Magdalena J. Hamel-Bravo, Deanne B. Fisher, Danielle Berg, Bjarki Björgvinsson, Alberto D. Bolatto, Alex J. Cameron, John Chisholm, Drummond B. Fielding, Rodrigo Herrera-Camus, Glenn G. Kacprzak, Miao Li, Barbara Mazzilli Ciraulo, Anna F. McLeod, Daniel K. McPherson, Nikole M. Nielsen, Bronwyn Reichardt Chu, Ryan J. Rickards Vaught, Karin Sandstrom

Abstract: We present the results of direct-method metallicity measurements in the disk and outflow of the low-metallicity starburst galaxy NGC 1569. We use Keck Cosmic Web Imager observations to map the galaxy across 54$\arcsec$ (800 pc) along the major axis and 48$\arcsec$ (700 pc) along the minor axis with a spatial resolution of 1$\arcsec$ ($\sim$15 pc). We detect common strong emission lines ([\ion{O}{I… ▽ More We present the results of direct-method metallicity measurements in the disk and outflow of the low-metallicity starburst galaxy NGC 1569. We use Keck Cosmic Web Imager observations to map the galaxy across 54$\arcsec$ (800 pc) along the major axis and 48$\arcsec$ (700 pc) along the minor axis with a spatial resolution of 1$\arcsec$ ($\sim$15 pc). We detect common strong emission lines ([\ion{O}{III}] $λ$5007, H$β$, [\ion{O}{II}] $λ$3727) and the fainter [\ion{O}{III}] $λ$4363 auroral line, which allows us to measure electron temperature ($T_e$) and metallicity. Theory suggests that outflows drive metals out of the disk driving observed trends between stellar mass and gas-phase metallicity. Our main result is that the metallicity in the outflow is similar to that of the disk, $Z_{\rm out} / Z_{\rm ISM} \approx 1$. This is consistent with previous absorption line studies in higher mass galaxies. Assumption of a mass-loading factor of $\dot{M}_{\rm out}/{\rm SFR}\sim3$ makes the metal-loading of NGC 1569 consistent with expectations derived from the mass-metallicity relationship. Our high spatial resolution metallicity maps reveal a region around a supermassive star cluster (SSC-B) with distinctly higher metallicity and higher electron density, compared to the disk. Given the known properties of SSC-B the higher metallicity and density of this region are likely the result of star formation-driven feedback acting on the local scale. Overall, our results are consistent with the picture in which metal-enriched winds pollute the circumgalactic medium surrounding galaxies, and thus connect the small-scale feedback processes to large-scale properties of galaxy halos. △ Less

Submitted 6 April, 2024; originally announced April 2024.

Comments: 15 pages, 11 figures, accepted by MNRAS

Confusion of extragalactic sources in the far infrared: a baseline assessment of the performance of PRIMAger in intensity and polarization

Authors: Matthieu Béthermin, Alberto D. Bolatto, François Boulanger, Charles M. Bradford, Denis Burgarella, Laure Ciesla, James Donnellan, Brandon S. Hensley, Jason Glenn, Guilaine Lagache, Enrique Lopez-Rodriguez, Seb Oliver, Alexandra Pope, Marc Sauvage

Abstract: Because of their limited angular resolution, far-infrared telescopes are usually affected by confusion phenomenon. Since several galaxies can be located in the same instrumental beam, only the brightest objects emerge from the fluctuations caused by fainter sources. The probe far-infrared mission for astrophysics imager (PRIMAger) will observe the mid- and far-infrared (25-235 $μ$m) sky both in in… ▽ More Because of their limited angular resolution, far-infrared telescopes are usually affected by confusion phenomenon. Since several galaxies can be located in the same instrumental beam, only the brightest objects emerge from the fluctuations caused by fainter sources. The probe far-infrared mission for astrophysics imager (PRIMAger) will observe the mid- and far-infrared (25-235 $μ$m) sky both in intensity and polarization. We aim to provide predictions of the confusion level and its consequences for future surveys. We produced simulated PRIMAger maps affected only by the confusion noise using the simulated infrared extragalactic sky (SIDES) semi-empirical simulation. We then estimated the confusion limit in these maps and extracted the sources using a basic blind extractor. By comparing the input galaxy catalog and the extracted source catalog, we derived various performance metrics as completeness, purity, and the accuracy of various measurements. In intensity, we predict that the confusion limit increases rapidly with increasing wavelength. The confusion limit in polarization is more than 100x lower. The measured flux density is dominated by the brightest galaxy in the beam, but other objects also contribute at longer wavelength (~30% at 235 $μ$m). We also show that galaxy clustering has a mild impact on confusion in intensity (up to 25%), while it is negligible in polarization. In intensity, a basic blind extraction will be sufficient to detect galaxies at the knee of the luminosity function up to z~3 and 10$^{11}$ M$_\odot$ main-sequence galaxies up to z~5. In polarization for a conservative sensitivity, we expect ~8 000 detections up to z=2.5 opening a totally new window on the high-z dust polarization. Finally, we show that intensity surveys at short wavelength and polarization surveys at long wavelength tend to reach confusion at similar depth. There is thus a strong synergy. △ Less

Submitted 28 April, 2024; v1 submitted 5 April, 2024; originally announced April 2024.

Comments: 21 pages, 15 figures, 4 tables, submitted to A&A, minor update on the names of the sensitivity scenarii (for consistency with other papers)

JWST Observations of Starbursts: Polycyclic Aromatic Hydrocarbon Emission at the Base of the M 82 Galactic Wind

Authors: Alberto D. Bolatto, Rebecca C. Levy, Elizabeth Tarantino, Martha L. Boyer, Deanne B. Fisher, Adam K. Leroy, Serena A. Cronin, Ralf S. Klessen, J. D. Smith, Dannielle A. Berg, Torsten Boeker, Leindert A. Boogaard, Eve C. Ostriker, Todd A. Thompson, Juergen Ott, Laura Lenkic, Laura A. Lopez, Daniel A. Dale, Sylvain Veilleux, Paul P. van der Werf, Simon C. O. Glover, Karin M. Sandstrom, Evan D. Skillman, John Chisholm, Vicente Villanueva , et al. (15 additional authors not shown)

Abstract: We present new observations of the central 1 kpc of the M 82 starburst obtained with the James Webb Space Telescope (JWST) near-infrared camera (NIRCam) instrument at a resolution ~0.05"-0.1" (~1-2 pc). The data comprises images in three mostly continuum filters (F140M, F250M, and F360M), and filters that contain [FeII] (F164N), H2 v=1-0 (F212N), and the 3.3 um PAH feature (F335M). We find promine… ▽ More We present new observations of the central 1 kpc of the M 82 starburst obtained with the James Webb Space Telescope (JWST) near-infrared camera (NIRCam) instrument at a resolution ~0.05"-0.1" (~1-2 pc). The data comprises images in three mostly continuum filters (F140M, F250M, and F360M), and filters that contain [FeII] (F164N), H2 v=1-0 (F212N), and the 3.3 um PAH feature (F335M). We find prominent plumes of PAH emission extending outward from the central starburst region, together with a network of complex filamentary substructure and edge-brightened bubble-like features. The structure of the PAH emission closely resembles that of the ionized gas, as revealed in Paschen alpha and free-free radio emission. We discuss the origin of the structure, and suggest the PAHs are embedded in a combination of neutral, molecular, and photoionized gas. △ Less

Submitted 21 April, 2024; v1 submitted 29 January, 2024; originally announced January 2024.

Comments: Submitted to The Astrophysical Journal

Hidden Gems on a Ring: Infant Massive Clusters and Their Formation Timeline Unveiled by ALMA, HST, and JWST in NGC 3351

Authors: Jiayi Sun, Hao He, Kyle Batschkun, Rebecca C. Levy, Kimberly Emig, M. Jimena Rodriguez, Hamid Hassani, Adam K. Leroy, Eva Schinnerer, Eve C. Ostriker, Christine D. Wilson, Alberto D. Bolatto, Elisabeth A. C. Mills, Erik Rosolowsky, Janice C. Lee, Daniel A. Dale, Kirsten L. Larson, David A. Thilker, Leonardo Ubeda, Bradley C. Whitmore, Thomas G. Williams, Ashley. T. Barnes, Frank Bigiel, Melanie Chevance, Simon C. O. Glover , et al. (16 additional authors not shown)

Abstract: We study young massive clusters (YMCs) in their embedded "infant" phase with $\sim0.\!^{\prime\prime}1$ ALMA, HST, and JWST observations targeting the central starburst ring in NGC 3351, a nearby Milky Way analog galaxy. Our new ALMA data reveal 18 bright and compact (sub-)millimeter continuum sources, of which 8 have counterparts in JWST images and only 6 have counterparts in HST images. Based on… ▽ More We study young massive clusters (YMCs) in their embedded "infant" phase with $\sim0.\!^{\prime\prime}1$ ALMA, HST, and JWST observations targeting the central starburst ring in NGC 3351, a nearby Milky Way analog galaxy. Our new ALMA data reveal 18 bright and compact (sub-)millimeter continuum sources, of which 8 have counterparts in JWST images and only 6 have counterparts in HST images. Based on the ALMA continuum and molecular line data, as well as ancillary measurements for the HST and JWST counterparts, we identify 14 sources as infant star clusters with high stellar and/or gas masses (${\sim}10^5\;\mathrm{M_\odot}$), small radii (${\lesssim}\,5\;\mathrm{pc}$), large escape velocities ($6{-}10\;\mathrm{km/s}$), and short free-fall times ($0.5{-}1\;\mathrm{Myr}$). Their multiwavelength properties motivate us to divide them into four categories, likely corresponding to four evolutionary stages from starless clumps to exposed HII region-cluster complexes. Leveraging age estimates for HST-identified clusters in the same region, we infer an evolutionary timeline going from $\sim$1-2 Myr before cluster formation as starless clumps, to $\sim$4-6 Myr after as exposed HII region-cluster complexes. Finally, we show that the YMCs make up a substantial fraction of recent star formation across the ring, exhibit an non-uniform azimuthal distribution without a very coherent evolutionary trend along the ring, and are capable of driving large-scale gas outflows. △ Less

Submitted 10 April, 2024; v1 submitted 25 January, 2024; originally announced January 2024.

Comments: 27 pages, 12 figures; ApJ accepted

The EDGE-CALIFA Survey: An Extragalactic Database for Galaxy Evolution Studies

Authors: Tony Wong, Yixian Cao, Yufeng Luo, Alberto D. Bolatto, Sebastián F. Sánchez, Jorge K. Barrera-Ballesteros, Leo Blitz, Dario Colombo, Helmut Dannerbauer, Alex Green, Veselina Kalinova, Ferzem Khan, Andrew Kim, Eduardo A. D. Lacerda, Adam K. Leroy, Rebecca C. Levy, Xincheng Lin, Yuanze Luo, Erik W. Rosolowsky, Mónica Rubio, Peter Teuben, Dyas Utomo, Vicente Villanueva, Stuart N. Vogel, Xinyu Wang

Abstract: The EDGE-CALIFA survey provides spatially resolved optical integral field unit (IFU) and CO spectroscopy for 125 galaxies selected from the CALIFA Data Release 3 sample. The Extragalactic Database for Galaxy Evolution (EDGE) presents the spatially resolved products of the survey as pixel tables that reduce the oversampling in the original images and facilitate comparison of pixels from different i… ▽ More The EDGE-CALIFA survey provides spatially resolved optical integral field unit (IFU) and CO spectroscopy for 125 galaxies selected from the CALIFA Data Release 3 sample. The Extragalactic Database for Galaxy Evolution (EDGE) presents the spatially resolved products of the survey as pixel tables that reduce the oversampling in the original images and facilitate comparison of pixels from different images. By joining these pixel tables to lower dimensional tables that provide radial profiles, integrated spectra, or global properties, it is possible to investigate the dependence of local conditions on large-scale properties. The database is freely accessible and has been utilized in several publications. We illustrate the use of this database and highlight the effects of CO upper limits on the inferred slopes of the local scaling relations between stellar mass, star formation rate (SFR), and H$_2$ surface densities. We find that the correlation between H$_2$ and SFR surface density is the tightest among the three relations. △ Less

Submitted 23 January, 2024; originally announced January 2024.

Comments: 21 pages, accepted for publication in ApJS, see DOIs below for code and data access

Hot Gas Outflow Properties of the Starburst Galaxy NGC 4945

Authors: Natalia Porraz Barrera, Sebastian Lopez, Laura A. Lopez, Adi Foord, Dustin D. Nguyen, Todd A. Thompson, Smita Mathur, Alberto D. Bolatto

Abstract: We analyze 330 ks of {\it Chandra} X-ray imaging and spectra of the nearby, edge-on starburst and Seyfert Type 2 galaxy NGC 4945 to measure the hot gas properties along the galactic outflows. We extract and model spectra from 15 regions extending from $-$0.55 kpc to $+$0.85 kpc above and below the galactic disk to determine the best-fit parameters and metal abundances. We find that the hot gas tem… ▽ More We analyze 330 ks of {\it Chandra} X-ray imaging and spectra of the nearby, edge-on starburst and Seyfert Type 2 galaxy NGC 4945 to measure the hot gas properties along the galactic outflows. We extract and model spectra from 15 regions extending from $-$0.55 kpc to $+$0.85 kpc above and below the galactic disk to determine the best-fit parameters and metal abundances. We find that the hot gas temperatures and number densities peak in the central regions and decrease along the outflows. These profiles are inconsistent with a spherical, adiabatically-expanding wind model, suggesting the need to include mass loading and/or a non-spherical outflow geometry. We estimate the mass outflow rate of the hot wind to be $1.6\:M_{\odot}~\rm{yr}^{-1}$. Emission from charge exchange is detected in the northern outflow, and we estimate it contributes 12\% to the emitted, broad-band ($0.5-7$ keV) X-ray flux. △ Less

Submitted 11 April, 2024; v1 submitted 13 December, 2023; originally announced December 2023.

Comments: 11 pages, 7 figures; submitted (13 December 2023), accepted ApJ (11 April 2023)

The EDGE-CALIFA survey: Molecular Gas and Star Formation Activity Across the Green Valley

Authors: Vicente Villanueva, Alberto D. Bolatto, Stuart N. Vogel, Tony Wong, Adam K. Leroy, Sebastian F. Sanchez, Rebecca C. Levy, Erik Rosolowsky, Dario Colombo, Veselina Kalinova, Serena Cronin, Peter Teuben, Monica Rubio, Zein Bazzi

Abstract: We present a $^{12}$CO($J$=2-1) survey of 60 local galaxies using data from the Atacama Large Millimeter/submillimeter Compact Array as part of the Extragalactic Database for Galaxy Evolution: the ACA EDGE survey. These galaxies all have integral field spectroscopy from the CALIFA survey. Compared to other local galaxy surveys, ACA EDGE is designed to mitigate selection effects based on CO brightn… ▽ More We present a $^{12}$CO($J$=2-1) survey of 60 local galaxies using data from the Atacama Large Millimeter/submillimeter Compact Array as part of the Extragalactic Database for Galaxy Evolution: the ACA EDGE survey. These galaxies all have integral field spectroscopy from the CALIFA survey. Compared to other local galaxy surveys, ACA EDGE is designed to mitigate selection effects based on CO brightness and morphological type. Of the 60 galaxies in ACA EDGE, 36 are on the star-formation main sequence, 13 are on the red sequence, and 11 lie in the ``green valley" transition between these sequences. We test how star formation quenching processes affect the star formation rate (SFR) per unit molecular gas mass, SFE$_{\rm mol}=$SFR/$M_{\rm mol}$, and related quantities in galaxies with stellar masses $10\leq$log[$M_\star/$M$_\odot$]$\leq11.5$ covering the full range of morphological types. We observe a systematic decrease of the molecular-to-stellar mass fraction ($R^{\rm mol}_{\star}$) with decreasing level of star formation activity, with green valley galaxies having also lower SFE$_{\rm mol}$ than galaxies on the main sequence. On average, we find that the spatially resolved SFE$_{\rm mol}$ within the bulge region of green valley galaxies is lower than in the bulges of main sequence galaxies if we adopt a constant CO-to-H$_2$ conversion factor, $α_{\rm CO}$. While efficiencies in main sequence galaxies remain almost constant with galactocentric radius, in green valley galaxies we note a systematic increase of SFE$_{\rm mol}$, $R^{\rm mol}_{\star}$, and specific star formation rate, sSFR, with increasing radius. Our results suggest that although gas depletion (or removal) seems to be the most important driver of the star-formation quenching in galaxies transiting through the green valley, a reduction in star formation efficiency is also required during this stage. △ Less

Submitted 6 December, 2023; originally announced December 2023.

Resolved Measurements of the CO-to-H$_2$ Conversion Factor in 37 Nearby Galaxies

Authors: I-Da Chiang, Karin M. Sandstrom, Jeremy Chastenet, Alberto D. Bolatto, Eric W. Koch, Adam K. Leroy, Jiayi Sun, Yu-Hsuan Teng, Thomas G. Williams

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

Submitted 26 December, 2023; v1 submitted 1 November, 2023; originally announced November 2023.

Comments: 27 pages, 13 figures, 4 tables, submitted to ApJ

SuperCAM CO(3-2) APEX survey at 6 pc resolution in the Small Magellanic Clouds

Authors: H. P. Saldaño, M. Rubio, A. D. Bolatto, K. Sandstrom, B. J. Swift, C. Verdugo, K. Jameson, C. K. Walker, C. Kulesa, J. Spilker, P. Bergman, G. A. Salazar

Abstract: We present the CO(3-2) APEX survey at 6 pc resolution of the bar of the SMC. We aboard the CO analysis in the SMC-Bar comparing the CO(3-2) survey with that of the CO(2-1) of similar resolution. We study the CO(3-2)-to-CO(2-1) ratio (R32) that is very sensitive to the environment properties (e.g., star-forming regions). We analyzed the correlation of this ratio with observational quantities that t… ▽ More We present the CO(3-2) APEX survey at 6 pc resolution of the bar of the SMC. We aboard the CO analysis in the SMC-Bar comparing the CO(3-2) survey with that of the CO(2-1) of similar resolution. We study the CO(3-2)-to-CO(2-1) ratio (R32) that is very sensitive to the environment properties (e.g., star-forming regions). We analyzed the correlation of this ratio with observational quantities that trace the star formation as the local CO emission, the Spitzer color [70/160], and the total IR surface brightness measured from the Spitzer and Herschel bands. For the identification of the CO(3-2) clouds, we used the CPROPS algorithm, which allowed us to measure the physical properties of the clouds. We analyzed the scaling relationships of such physical properties. We obtained an R32 of 0.65 as a median value for the SMC, with a standard deviation of 0.3. We found that R32 varies from region to region, depending on the star formation activity. In regions dominated by HII and photo-dissociated regions (e.g., N22, N66), R32 tends to be higher than the median values. Meanwhile, lower values were found toward quiescent clouds. We also found that R32 correlates positively with the IR color [70/160] and the total IR surface brightness. This finding indicates that R32 increases with environmental properties like the dust temperature, the total gas density, and the radiation field. We have identified 225 molecular clouds with sizes R > 1.5 pc and signal-to-noise (S/N) ratio > 3 and only 17 well-resolved CO(3-2) clouds increasing the S/N ratio to > 5. These 17 clouds follow consistent scaling relationships to the inner Milky Way clouds but with some departure. The CO(3-2) tends to be less turbulent and less luminous than the inner Milky Way clouds of similar size. Finally, we estimated a median virial-based CO-to-H2 conversion factor of 12.6_{-7}^{+10} Msun/(K km s^{-1} pc^{2}) for the total sample. △ Less

Submitted 31 October, 2023; originally announced November 2023.

Journal ref: A&A 687, A26 (2024)

Star Formation Efficiency in Nearby Galaxies Revealed with a New CO-to-H2 Conversion Factor Prescription

Authors: Yu-Hsuan Teng, I-Da Chiang, Karin M. Sandstrom, Jiayi Sun, Adam K. Leroy, Alberto D. Bolatto, Antonio Usero, Eve C. Ostriker, Miguel Querejeta, Jeremy Chastenet, Frank Bigiel, Mederic Boquien, Jakob den Brok, Yixian Cao, Melanie Chevance, Ryan Chown, Dario Colombo, Cosima Eibensteiner, Simon C. O. Glover, Kathryn Grasha, Jonathan D. Henshaw, Maria J. Jimenez-Donaire, Daizhong Liu, Eric J. Murphy, Hsi-An Pan , et al. (2 additional authors not shown)

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

Submitted 24 November, 2023; v1 submitted 24 October, 2023; originally announced October 2023.

Comments: 12 pages of main text, 4 figures, accepted for publication in ApJ

[CII] Spectral Mapping of the Galactic Wind and Starbursting Disk of M82 with SOFIA

Authors: Rebecca C. Levy, Alberto D. Bolatto, Elizabeth Tarantino, Adam K. Leroy, Lee Armus, Kimberly L. Emig, Rodrigo Herrera-Camus, Daniel P. Marrone, Elisabeth Mills, Oliver Ricken, Juergen Stutzki, Sylvain Veilleux, Fabian Walter

Abstract: M82 is an archetypal starburst galaxy in the local Universe. The central burst of star formation, thought to be triggered by M82's interaction with other members in the M81 group, is driving a multiphase galaxy-scale wind away from the plane of the disk that has been studied across the electromagnetic spectrum. Here, we present new velocity-resolved observations of the [CII] 158$μ$m line in the ce… ▽ More M82 is an archetypal starburst galaxy in the local Universe. The central burst of star formation, thought to be triggered by M82's interaction with other members in the M81 group, is driving a multiphase galaxy-scale wind away from the plane of the disk that has been studied across the electromagnetic spectrum. Here, we present new velocity-resolved observations of the [CII] 158$μ$m line in the central disk and the southern outflow of M82 using the upGREAT instrument onboard SOFIA. We also report the first detections of velocity-resolved ($ΔV = 10$ km s$^{-1}$) [CII] emission in the outflow of M82 at projected distances of $\approx1-2$ kpc south of the galaxy center. We compare the [CII] line profiles to observations of CO and HI and find that likely the majority ($>55$%) of the [CII] emission in the outflow is associated with the neutral atomic medium. We find that the fraction of [CII] actually outflowing from M82 is small compared to the bulk gas outside the midplane (which may be in a halo or tidal streamers), which has important implications for observations of [CII] outflows at higher redshift. Finally, by comparing the observed ratio of the [CII] and CO intensities to models of photodissociation regions, we estimate that the far-ultraviolet (FUV) radiation field in the disk is $\sim10^{3.5}~G_0$, in agreement with previous estimates. In the outflow, however, the FUV radiation field is 2-3 orders of magnitudes lower, which may explain the high fraction of [CII] arising from the neutral medium in the wind. △ Less

Submitted 27 September, 2023; originally announced September 2023.

Comments: 22 pages, 9 figures, Accepted for publication in ApJ

VERTICO VII: Environmental quenching caused by suppression of molecular gas content and star formation efficiency in Virgo Cluster galaxies

Authors: Toby Brown, Ian D. Roberts, Mallory Thorp, Sara L. Ellison, Nikki Zabel, Christine D. Wilson, Yannick M. Bahé, Dhruv Bisaria, Alberto D. Bolatto, Alessandro Boselli, Aeree Chung, Luca Cortese, Barbara Catinella, Timothy A. Davis, María J. Jiménez-Donaire, Claudia D. P. Lagos, Bumhyun Lee, Laura C. Parker, Rory Smith, Kristine Spekkens, Adam R. H. Stevens, Vicente Villanueva, Adam B. Watts

Abstract: We study how environment regulates the star formation cycle of 33 Virgo Cluster satellite galaxies on 720 parsec scales. We present the first resolved star-forming main sequence for cluster galaxies, dividing the sample based on their global HI properties and comparing to a control sample of field galaxies. HI-poor cluster galaxies have reduced star formation rate (SFR) surface densities with resp… ▽ More We study how environment regulates the star formation cycle of 33 Virgo Cluster satellite galaxies on 720 parsec scales. We present the first resolved star-forming main sequence for cluster galaxies, dividing the sample based on their global HI properties and comparing to a control sample of field galaxies. HI-poor cluster galaxies have reduced star formation rate (SFR) surface densities with respect to both HI-normal cluster and field galaxies (0.5 dex), suggesting that mechanisms regulating the global HI content are responsible for quenching local star formation. We demonstrate that the observed quenching in HI-poor galaxies is caused by environmental processes such as ram pressure stripping (RPS) simultaneously reducing molecular gas surface density and star formation efficiency (SFE), compared to regions in HI-normal systems (by 0.38 and 0.22 dex, respectively). We observe systematically elevated SFRs that are driven by increased molecular gas surface densities at fixed stellar mass surface density in the outskirts of early-stage RPS galaxies, while SFE remains unchanged with respect to the field sample. We quantify how RPS and starvation affect the star formation cycle of inner and outer galaxy discs as they are processed by the cluster. We show both are effective quenching mechanisms with the key difference being that RPS acts upon the galaxy outskirts while starvation regulates the star formation cycle throughout disc, including within the truncation radius. For both processes, the quenching is caused by a simultaneous reduction in molecular gas surface densities and SFE at fixed stellar mass surface density. △ Less

Submitted 21 August, 2023; originally announced August 2023.

Comments: 17 pages, 1 table, 5 figures, accepted for publication in ApJ

DUVET Survey: Mapping Outflows in the Metal-Poor Starburst Mrk 1486

Authors: Daniel K. McPherson, Deanne B. Fisher, Nikole M. Nielsen, Glenn G. Kacprzak, Bronwyn Reichardt Chu, Alex J. Cameron, Alberto D. Bolatto, John Chisholm, Drummond B. Fielding, Danielle Berg, Rodrigo Herrera-Camus, Miao Li, Ryan J. Rickards Vaught, Karin Sandstrom

Abstract: We present a method to characterize star-formation driven outflows from edge-on galaxies and apply this method to the metal-poor starburst galaxy, Mrk 1486. Our method uses the distribution of emission line flux (from H$β$ and [OIII] 5007) to identify the location of the outflow and measure the extent above the disk, the opening angle, and the transverse kinematics. We show that this simple techni… ▽ More We present a method to characterize star-formation driven outflows from edge-on galaxies and apply this method to the metal-poor starburst galaxy, Mrk 1486. Our method uses the distribution of emission line flux (from H$β$ and [OIII] 5007) to identify the location of the outflow and measure the extent above the disk, the opening angle, and the transverse kinematics. We show that this simple technique recovers a similar distribution of the outflow without requiring complex modelling of line-splitting or multi-Gaussian components, and is therefore applicable to lower spectral resolution data. In Mrk 1486 we observe an asymmetric outflow in both the location of the peak flux and total flux from each lobe. We estimate an opening angle of $17-37^{\circ}$ depending on the method and assumptions adopted. Within the minor axis outflows, we estimate a total mass outflow rate of $\sim2.5$ M$_{\odot}$ yr$^{-1}$, which corresponds to a mass loading factor of $η=0.7$. We observe a non-negligible amount of flux from ionized gas outflowing along the edge of the disk (perpendicular to the biconical components), with a mass outflow rate $\sim0.9$ M$_{\odot}$ yr$^{-1}$. Our results are intended to demonstrate a method that can be applied to high-throughput, low spectral resolution observations, such as narrow band filters or low spectral resolution IFS that may be more able to recover the faint emission from outflows. △ Less

Submitted 13 August, 2023; originally announced August 2023.

Comments: 12 Pages, 6 Figures

The EDGE-CALIFA Survey: Spatially Resolved 13CO(1-0) Observations and Variations in 12CO(1-0)/13CO(1-0) in Nearby Galaxies on kpc Scales

Authors: Yixian Cao, Tony Wong, Alberto D. Bolatto, Adam Leroy, Erik W. Rosolowsky, Dyas Utomo, Sebastian Sanchez, Jorge Barrera-Ballesteros, Rebecca Levy, Dario Colombo, Leo Blitz, Stuart Vogel, Johannes Puschnig, Vicente Villanueva, Monica Rubio

Abstract: We present 13CO(1-0) observations for the EDGE-CALIFA survey, which is a mapping survey of 126 nearby galaxies at a typical spatial resolution of 1.5 kpc. Using detected 12CO(1-0) emission as a prior, we detect 13CO(1-0) in 41 galaxies via integrated line flux over the entire galaxy, and in 30 galaxies via integrated line intensity in resolved synthesized beams. Incorporating our CO observations a… ▽ More We present 13CO(1-0) observations for the EDGE-CALIFA survey, which is a mapping survey of 126 nearby galaxies at a typical spatial resolution of 1.5 kpc. Using detected 12CO(1-0) emission as a prior, we detect 13CO(1-0) in 41 galaxies via integrated line flux over the entire galaxy, and in 30 galaxies via integrated line intensity in resolved synthesized beams. Incorporating our CO observations and optical IFU spectroscopy, we perform a systematic comparison between the line ratio R12/13 and the properties of the stars and ionized gas. Higher R12/13 values are found in interacting galaxies than in non-interacting galaxies. The global R12/13 slightly increases with infrared color F60/F100, but appears insensitive to other host galaxy properties such as morphology, stellar mass, or galaxy size. We also present annulus-averaged R12/13 profiles for our sample up to a galactocentric radius of 0.4r25 (~6 kpc), taking into account the 13CO(1-0) non-detections by spectral stacking. The radial profiles of R12/13 are quite flat across our sample. Within galactocentric distances of 0.2r25, azimuthally-averaged R12/13 increases with star formation rate. However, the Spearman rank correlation tests show the azimuthally-averaged R12/13 does not strongly correlate with any other gas or stellar properties in general, especially beyond 0.2r25 from the galaxy centers. Our findings suggest that in the complex environments in galaxy disks, R12/13 is not a sensitive tracer for ISM properties. Dynamical disturbances, like galaxy interactions or the presence of a bar, also have an overall impact on R12/13, which further complicate the interpretations of R12/13 variations. △ Less

Submitted 13 June, 2023; originally announced June 2023.

Comments: 35pages, 11 figure, 6 tables, accepted for publication in ApJS

Exploring the Impact of Galactic Interactions and Mergers on the Central Star Formation of APEX/EDGE-CALIFA Galaxies

Authors: Y. Garay-Solis, J. K. Barrera-Ballesteros, D. Colombo, S. F. Sánchez, A. Z. Lugo-Aranda, V. Villanueva, T. Wong, A. D. Bolatto

Abstract: Galactic interactions and subsequent mergers are a paramount channel for galaxy evolution. In this work, we use the data from 236 star forming CALIFA galaxies with integrated molecular gas observations in their central region (approximately within an effective radius) -- from the APEX millimeter telescope and the CARMA millimeter telescope array. This sample includes isolated (126 galaxies) and in… ▽ More Galactic interactions and subsequent mergers are a paramount channel for galaxy evolution. In this work, we use the data from 236 star forming CALIFA galaxies with integrated molecular gas observations in their central region (approximately within an effective radius) -- from the APEX millimeter telescope and the CARMA millimeter telescope array. This sample includes isolated (126 galaxies) and interacting galaxies in different merging stages (110 galaxies; from pairs, merging and post-merger galaxies). We show that the impact of interactions and mergers in the center of galaxies is revealed as an increase in the fraction of molecular gas (compared to isolated galaxies). Furthermore, our results suggest that the change in star formation efficiency is the main driver for both an enhancement and/or suppression of the central star formation -- except in merging galaxies where the enhanced star formation appears to be driven by an increase of molecular gas. We suggest that gravitational torques due to the interaction and subsequent merger transport cold molecular gas inwards, increasing the gas fraction without necessarily increasing star formation. △ Less

Submitted 5 June, 2023; originally announced June 2023.

Comments: 3 Figures, 13 Pages. Accepted for publication in ApJ

The Physical Drivers and Observational Tracers of CO-to-H2 Conversion Factor Variations in Nearby Barred Galaxy Centers

Authors: Yu-Hsuan Teng, Karin M. Sandstrom, Jiayi Sun, Munan Gong, Alberto D. Bolatto, I-Da Chiang, Adam K. Leroy, Antonio Usero, Simon C. O. Glover, Ralf S. Klessen, Daizhong Liu, Miguel Querejeta, Eva Schinnerer, Frank Bigiel, Yixian Cao, Melanie Chevance, Cosima Eibensteiner, Kathryn Grasha, Frank P. Israel, Eric J. Murphy, Lukas Neumann, Hsi-An Pan, Francesca Pinna, Mattia C. Sormani, J. D. T. Smith , et al. (2 additional authors not shown)

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

Submitted 10 April, 2023; originally announced April 2023.

Comments: Accepted for publication in ApJ; 30 pages of main text + 3 appendices

VERTICO V: The environmentally driven evolution of the inner cold gas discs of Virgo cluster galaxies

Authors: Adam B. Watts, Luca Cortese, Barbara Catinella, Toby Brown, Christine D. Wilson, Nikki Zabel, Ian D. Roberts, Timothy A. Davis, Mallory Thorp, Aeree Chung, Adam R. H. Stevens, Sara L. Ellison, Kristine Spekkens, Laura C. Parker, Yannick M. Bahé, Vicente Villanueva, María Jiménez-Donaire, Dhruv Bisaria, Alessandro Boselli, Alberto D. Bolatto, Bumhyun Lee

Abstract: The quenching of cluster satellite galaxies is inextricably linked to the suppression of their cold interstellar medium (ISM) by environmental mechanisms. While the removal of neutral atomic hydrogen (HI) at large radii is well studied, how the environment impacts the remaining gas in the centres of galaxies, which are dominated by molecular gas, is less clear. Using new observations from the Virg… ▽ More The quenching of cluster satellite galaxies is inextricably linked to the suppression of their cold interstellar medium (ISM) by environmental mechanisms. While the removal of neutral atomic hydrogen (HI) at large radii is well studied, how the environment impacts the remaining gas in the centres of galaxies, which are dominated by molecular gas, is less clear. Using new observations from the Virgo Environment traced in CO survey (VERTICO) and archival HI data, we study the HI and molecular gas within the optical discs of Virgo cluster galaxies on 1.2-kpc scales with spatially resolved scaling relations between stellar ($Σ_{\star}$), HI ($Σ_\mathrm{HI}$), and molecular gas ($Σ_\mathrm{mol}$) surface densities. Adopting HI deficiency as a measure of environmental impact, we find evidence that, in addition to removing the HI at large radii, the cluster processes also lower the average $Σ_\mathrm{HI}$ of the remaining gas even in the central 1.2 kpc. The impact on molecular gas is comparatively weaker than on the HI, and we show that the lower $Σ_\mathrm{mol}$ gas is removed first. In the most HI-deficient galaxies, however, we find evidence that environmental processes reduce the typical $Σ_\mathrm{mol}$ of the remaining gas by nearly a factor of 3. We find no evidence for environment-driven elevation of $Σ_\mathrm{HI}$ or $Σ_\mathrm{mol}$ in HI-deficient galaxies. Using the ratio of $Σ_\mathrm{mol}$-to-$Σ_\mathrm{HI}$ in individual regions, we show that changes in the ISM physical conditions, estimated using the total gas surface density and midplane hydrostatic pressure, cannot explain the observed reduction in molecular gas content. Instead, we suggest that direct stripping of the molecular gas is required to explain our results. △ Less

Submitted 13 March, 2023; originally announced March 2023.

Comments: Accepted for publication in PASA. Main text is 19 pages including 12 figures and 3 tables, plus 1 appendix. A 2.5 min, high-level summary can be found at https://youtu.be/7djMmVEpDVc

Star Formation Laws and Efficiencies across 80 Nearby Galaxies

Authors: Jiayi Sun, Adam K. Leroy, Eve C. Ostriker, Sharon Meidt, Erik Rosolowsky, Eva Schinnerer, Christine D. Wilson, Dyas Utomo, Francesco Belfiore, Guillermo A. Blanc, Eric Emsellem, Christopher Faesi, Brent Groves, Annie Hughes, Eric W. Koch, Kathryn Kreckel, Daizhong Liu, Hsi-An Pan, Jerome Pety, Miguel Querejeta, Alessandro Razza, Toshiki Saito, Amy Sardone, Antonio Usero, Thomas G. Williams , et al. (15 additional authors not shown)

Abstract: We measure empirical relationships between the local star formation rate (SFR) and properties of the star-forming molecular gas on 1.5 kpc scales across 80 nearby galaxies. These relationships, commonly referred to as "star formation laws," aim at predicting the local SFR surface density from various combinations of molecular gas surface density, galactic orbital time, molecular cloud free-fall ti… ▽ More We measure empirical relationships between the local star formation rate (SFR) and properties of the star-forming molecular gas on 1.5 kpc scales across 80 nearby galaxies. These relationships, commonly referred to as "star formation laws," aim at predicting the local SFR surface density from various combinations of molecular gas surface density, galactic orbital time, molecular cloud free-fall time, and the interstellar medium dynamical equilibrium pressure. Leveraging a multiwavelength database built for the PHANGS survey, we measure these quantities consistently across all galaxies and quantify systematic uncertainties stemming from choices of SFR calibrations and the CO-to-H$_2$ conversion factors. The star formation laws we examine show 0.3-0.4 dex of intrinsic scatter, among which the molecular Kennicutt-Schmidt relation shows a $\sim$10% larger scatter than the other three. The slope of this relation ranges $β\approx0.9{-}1.2$, implying that the molecular gas depletion time remains roughly constant across the environments probed in our sample. The other relations have shallower slopes ($β\approx0.6{-}1.0$), suggesting that the star formation efficiency (SFE) per orbital time, the SFE per free-fall time, and the pressure-to-SFR surface density ratio (i.e., the feedback yield) may vary systematically with local molecular gas and SFR surface densities. Last but not least, the shapes of the star formation laws depend sensitively on methodological choices. Different choices of SFR calibrations can introduce systematic uncertainties of at least 10-15% in the star formation law slopes and 0.15-0.25 dex in their normalization, while the CO-to-H$_2$ conversion factors can additionally produce uncertainties of 20-25% for the slope and 0.10-0.20 dex for the normalization. △ Less

Submitted 23 February, 2023; originally announced February 2023.

Comments: 10 pages main text + 2 appendices. ApJL in press. Data products available at https://www.canfar.net/storage/list/phangs/RELEASES/Sun_etal_2023 . Slides summarizing key results can be found at https://www.dropbox.com/s/5gsegexeo9n0t05/Sun_et_PHANGS_2023.pptx?dl=0

Journal ref: ApJL, 945, L19 (2023)

CO Excitation in High-z Main Sequence Analogues: Resolved CO(4-3)/CO(3-2) Line Ratios in DYNAMO Galaxies

Authors: Laura Lenkić, Alberto D. Bolatto, Deanne B. Fisher, Roberto Abraham, Karl Glazebrook, Rodrigo Herrera-Camus, Rebecca C. Levy, Danail Obreschkow, Carolyn G. Volpert

Abstract: The spectral line energy distribution of carbon monoxide contains information about the physical conditions of the star forming molecular hydrogen gas; however, the relation to local radiation field properties is poorly constrained. Using ~ 1-2 kpc scale ALMA observations of CO(3-2) and CO(4-3), we characterize the CO(4-3)/CO(3-2) line ratios of local analogues of main sequence galaxies at z ~ 1-2… ▽ More The spectral line energy distribution of carbon monoxide contains information about the physical conditions of the star forming molecular hydrogen gas; however, the relation to local radiation field properties is poorly constrained. Using ~ 1-2 kpc scale ALMA observations of CO(3-2) and CO(4-3), we characterize the CO(4-3)/CO(3-2) line ratios of local analogues of main sequence galaxies at z ~ 1-2, drawn from the DYNAMO sample. We measure CO(4-3)/CO(3-2) across the disk of each galaxy and find a median line ratio of $R_{43} = 0.54^{+0.16}_{-0.15}$ for the sample. This is higher than literature estimates of local star-forming galaxies and is consistent with multiple lines of evidence that indicate DYNAMO galaxies, despite residing in the local Universe, resemble main-sequence galaxies at z ~ 1-2. Comparing to existing lower resolution CO(1-0) observations, we find $R_{41}$ and $R_{31}$ values in the range $\sim 0.2-0.3$ and $\sim 0.4-0.8$ respectively. We combine our kpc-scale resolved line ratio measurements with HST observations of H$α$ to investigate the relation to star formation rate surface density and compare this relation to expectations from models. We find increasing CO(4-3)/CO(3-2) with increasing star formation rate surface density; however, models over-predict the line ratios across the range of star formation rate surface densities we probe, particularly at the lower range. Finally, SOFIA observations with HAWC+ and FIFI-LS reveal low dust temperatures and no deficit of [CII] emission with respect to the total infrared luminosity. △ Less

Submitted 12 January, 2023; originally announced January 2023.

Comments: 22 pages, 9 figures, 4 tables, ApJ accepted

PHANGS-JWST First Results: Mapping the 3.3 micron Polycyclic Aromatic Hydrocarbon Vibrational Band in Nearby Galaxies with NIRCam Medium Bands

Authors: Karin Sandstrom, Jérémy Chastenet, Jessica Sutter, Adam K. Leroy, Oleg V. Egorov, Thomas G. Williams, Alberto D. Bolatto, Médéric Boquien, Yixian Cao, Daniel A. Dale, Janice C. Lee, Erik Rosolowsky, Eva Schinnerer, Ashley. T. Barnes, Francesco Belfiore, F. Bigiel, Mélanie Chevance, Kathryn Grasha, Brent Groves, Hamid Hassani, Annie Hughes, Ralf S. Klessen, J. M. Diederik Kruijssen, Kirsten L. Larson, Daizhong Liu , et al. (6 additional authors not shown)

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

Submitted 2 January, 2023; originally announced January 2023.

Comments: submitted to AAS journals and revised according to referee comments, part of a PHANGS-JWST Focus Issue to appear in ApJ

PHANGS-JWST First Results: Tracing the Diffuse ISM with JWST Imaging of Polycyclic Aromatic Hydrocarbon Emission in Nearby Galaxies

Authors: Karin M. Sandstrom, Eric W. Koch, Adam K. Leroy, Erik Rosolowsky, Eric Emsellem, Rowan J. Smith, Oleg V. Egorov, Thomas G. Williams, Kirsten L. Larson, Janice C. Lee, Eva Schinnerer, David A. Thilker, Ashley. T. Barnes, Francesco Belfiore, F. Bigiel, Guillermo A. Blanc, Alberto D. Bolatto, Médéric Boquien, Yixian Cao, Jérémy Chastenet, Mélanie Chevance, I-Da Chiang, Daniel A. Dale, Christopher M. Faesi, Simon C. O. Glover , et al. (21 additional authors not shown)

Abstract: JWST observations of polycyclic aromatic hydrocarbon (PAH) emission provide some of the deepest and highest resolution views of the cold interstellar medium (ISM) in nearby galaxies. If PAHs are well mixed with the atomic and molecular gas and illuminated by the average diffuse interstellar radiation field, PAH emission may provide an approximately linear, high resolution, high sensitivity tracer… ▽ More JWST observations of polycyclic aromatic hydrocarbon (PAH) emission provide some of the deepest and highest resolution views of the cold interstellar medium (ISM) in nearby galaxies. If PAHs are well mixed with the atomic and molecular gas and illuminated by the average diffuse interstellar radiation field, PAH emission may provide an approximately linear, high resolution, high sensitivity tracer of diffuse gas surface density. We present a pilot study that explores using PAH emission in this way based on MIRI observations of IC 5332, NGC 628, NGC 1365, and NGC 7496 from the PHANGS-JWST Treasury. Using scaling relationships calibrated in Leroy et al. (2022), scaled F1130W provides 10--40 pc resolution and 3$σ$ sensitivity of $Σ_{\rm gas} \sim 2$ M$_\odot$ pc$^{-2}$. We characterize the surface densities of structures seen at $< 7$ M$_\odot$ pc$^{-2}$ in our targets, where we expect the gas to be HI-dominated. We highlight the existence of filaments, inter-arm emission, and holes in the diffuse ISM at these low surface densities. Below $\sim 10$ M$_\odot$ pc$^{-2}$ for NGC 628, NGC 1365, and NGC 7496 the gas distribution shows a ``Swiss cheese''-like topology due to holes and bubbles pervading the relatively smooth distribution of diffuse ISM. Comparing to recent galaxy simulations, we observe similar topology for the low surface density gas, though with notable variations between simulations with different setups and resolution. Such a comparison of high resolution, low surface density gas with simulations is not possible with existing atomic and molecular gas maps, highlighting the unique power of JWST maps of PAH emission. △ Less

Submitted 21 December, 2022; originally announced December 2022.

Comments: 18 pages, 6 figures, accepted as part of a PHANGS-JWST Focus Issue to appear in ApJ

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

PHANGS-JWST First Results: A Global and Moderately Resolved View of Mid-Infrared and CO Line Emission from Galaxies at the Start of the JWST Era

Authors: Adam K. Leroy, Alberto D. Bolatto, Karin Sandstrom, Erik Rosolowsky, Ashley. T. Barnes, F. Bigiel, Médéric Boquien, Jakob S. den Brok, Yixian Cao, Jérémy Chastenet, Mélanie Chevance, I-Da Chiang, Ryan Chown, Dario Colombo, Sara L. Ellison, Eric Emsellem, Kathryn Grasha, Jonathan D. Henshaw, Annie Hughes, Ralf S. Klessen, Eric W. Koch, Jaeyeon Kim, Kathryn Kreckel, J. M. Diederik Kruijssen, Kirsten L. Larson , et al. (19 additional authors not shown)

Abstract: We explore the relationship between mid-infrared (mid-IR) and CO rotational line emission from massive star-forming galaxies, which is one of the tightest scalings in the local universe. We assemble a large set of unresolved and moderately ($\sim 1$ kpc) spatially resolved measurements of CO (1-0) and CO (2-1) intensity, $I_{\rm CO}$, and mid-IR intensity, $I_{\rm MIR}$, at 8, 12, 22, and 24$μ$m.… ▽ More We explore the relationship between mid-infrared (mid-IR) and CO rotational line emission from massive star-forming galaxies, which is one of the tightest scalings in the local universe. We assemble a large set of unresolved and moderately ($\sim 1$ kpc) spatially resolved measurements of CO (1-0) and CO (2-1) intensity, $I_{\rm CO}$, and mid-IR intensity, $I_{\rm MIR}$, at 8, 12, 22, and 24$μ$m. The $I_{\rm CO}$ vs. $I_{\rm MIR}$ relationship is reasonably described by a power law with slopes $0.7{-}1.2$ and normalization $I_{\rm CO} \sim 1$ K km s$^{-1}$ at $I_{\rm MIR} \sim 1$ MJy sr$^{-1}$. Both the slopes and intercepts vary systematically with choice of line and band. The comparison between the relations measured for CO~(1-0) and CO (2-1) allow us to infer that $R_{21} \propto I_{\rm MIR}^{0.2}$, in good agreement with other work. The $8μ$m and $12μ$m bands, with strong PAH features, show steeper CO vs. mid-IR slopes than the $22μ$m and $24μ$m, consistent with PAH emission arising not just from CO-bright gas but also from atomic or CO-dark gas. The CO-to-mid-IR ratio correlates with global galaxy stellar mass ($M_\star$) and anti-correlates with SFR/$M_\star$. At $\sim 1$ kpc resolution, the first four PHANGS-JWST targets show CO to mid-IR relationships that are quantitatively similar to our larger literature sample, including showing the steep CO-to-mid-IR slopes for the JWST PAH-tracing bands, although we caution that these initial data have a small sample size and span a limited range of intensities. △ Less

Submitted 27 December, 2022; v1 submitted 19 December, 2022; originally announced December 2022.

Comments: 29 pages, 13 figures, key quantitative results in Table 3, Accepted as part of a PHANGS-JWST Focus Issue to appear in ApJ

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

VERTICO III: The Kennicutt-Schmidt relation in Virgo cluster galaxies

Authors: M. J. Jiménez-Donaire, T. Brown, C. D. Wilson, I. D. Roberts, N. Zabel, S. L. Ellison, M. Thorp, V. Villanueva, R. Chown, D. Bisaria, A. D. Bolatto, A. Boselli, B. Catinella, A. Chung, L. Cortese, T. A. Davis, C. D. P. Lagos, B. Lee, L. C. Parker, K. Spekkens, A. R. H. Stevens, J. Sun

Abstract: In this VERTICO science paper we aim to study how the star formation process depends on galactic environment and gravitational interactions in the context of galaxy evolution. We explore the scaling relation between the star formation rate (SFR) surface density and the molecular gas surface density, also known as the Kennicutt-Schmidt (KS) relation, in a subsample of Virgo cluster spiral galaxies.… ▽ More In this VERTICO science paper we aim to study how the star formation process depends on galactic environment and gravitational interactions in the context of galaxy evolution. We explore the scaling relation between the star formation rate (SFR) surface density and the molecular gas surface density, also known as the Kennicutt-Schmidt (KS) relation, in a subsample of Virgo cluster spiral galaxies. We use new ACA and TP observations from the VERTICO-ALMA Large Program at 720pc resolution to resolve the molecular gas content, as traced by the 12CO(2-1) transition, across the disks of 37 spiral galaxies in the Virgo cluster. In combination with archival observations, we estimate the parameters of the KS relation for the entire ensemble of galaxies, and within individual galaxies. We find the KS slope for the entire population to be N=0.97+/-0.07, with a characteristic molecular gas depletion time of 1.86Gyr for our full sample, in agreement with previous work in isolated star-forming galaxies. In individual galaxies, we find KS slopes ranging between 0.69 and 1.40, and typical star formation efficiencies (SFE) that can vary from galaxy to galaxy by a factor of ~4. These galaxy-to-galaxy variations account for ~0.20dex in scatter in the ensemble KS relation, which is characterized by a 0.42dex scatter. We find that the HI-deficient galaxies in the Virgo cluster show a steeper resolved KS relation and lower molecular gas efficiencies than HI-normal cluster galaxies. While the molecular gas content in Virgo cluster galaxies appears to behave similarly to that in isolated galaxies, our VERTICO sample shows that cluster environments play a key role in regulating star formation. The environmental mechanisms affecting the HI galaxy content also have a direct impact in the SFE of molecular gas in cluster galaxies, leading to longer depletion times in HI-deficient members. △ Less

Submitted 29 November, 2022; originally announced November 2022.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 671, A3 (2023)

CO(2-1) Survey at 9 parsec resolution in the SMC

Authors: H. P. Saldaño, M. Rubio, A. D. Bolatto, C. Verdugo, K. E. Jameson, A. K. Leroy

Abstract: The Small Magellanic Cloud (SMC) is the closest low-metallicity galaxy to the Milky Way where the dynamical state of molecular clouds can be analyzed. We present a CO(2-1) survey at 9 pc resolution obtained with the APEX telescope in an extensive region of the SMC and characterize the properties of the molecular clouds. We study the dynamical state and stability of these clouds uniformly. We ident… ▽ More The Small Magellanic Cloud (SMC) is the closest low-metallicity galaxy to the Milky Way where the dynamical state of molecular clouds can be analyzed. We present a CO(2-1) survey at 9 pc resolution obtained with the APEX telescope in an extensive region of the SMC and characterize the properties of the molecular clouds. We study the dynamical state and stability of these clouds uniformly. We identify 177 molecular clouds within the SMC by using CPROPS, of which 124 clouds are fully resolved with signal-to-noise ratio $>$ 5. The scaling relationships show that the SMC clouds are (on average) less turbulent and less luminous than their inner Milky Way counterparts of similar size by a factor of 2 and 3, respectively, while for a fixed linewidth, the SMC clouds are over-luminous by a factor of 3.5. Using the virial masses, we derive a CO-to-H2 conversion factor for the SMC CO clouds of 10.5 M$_{\odot}(K km s^{-1} pc^{2})^{-1}$, measured at 9 pc resolution. We also determine a dust-based conversion factor of 28 M$_{\odot}(K km s^{-1} pc^{2})^{-1}$, obtained at 12 pc resolution. We find that the SMC clouds appear to be in approximate gravitational virial equilibrium. We find that the cumulative mass functions based on both the luminous mass and the virial mass are steeper than $\frac{dN}{dM} \propto M^{-2}$, suggesting that most of the molecular mass of the SMC is contained in low-mass clouds. △ Less

Submitted 14 November, 2022; originally announced November 2022.

Comments: 22 pages, 18 figures. Submitted to the Astronomy & Astrophysics (A&A) journal

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

DUVET: Spatially Resolved Observations of Star Formation Regulation via Galactic Outflows in a Starbursting Disk Galaxy

Authors: Bronwyn Reichardt Chu, Deanne B. Fisher, Alberto D. Bolatto, John Chisholm, Drummond Fielding, Danielle Berg, Alex J. Cameron, Karl Glazebrook, Rodrigo Herrera-Camus, Glenn G. Kacprzak, Laura Lenkić, Miao Li, Daniel K. McPherson, Nikole M. Nielsen, Danail Obreschkow, Ryan J. Rickards Vaught, Karin Sandstrom

Abstract: We compare 500~pc scale, resolved observations of ionised and molecular gas for the $z\sim0.02$ starbursting disk galaxy IRAS08339+6517, using measurements from KCWI and NOEMA. We explore the relationship of the star formation driven ionised gas outflows with colocated galaxy properties. We find a roughly linear relationship between the outflow mass flux ($\dotΣ_{\rm out}$) and star formation rate… ▽ More We compare 500~pc scale, resolved observations of ionised and molecular gas for the $z\sim0.02$ starbursting disk galaxy IRAS08339+6517, using measurements from KCWI and NOEMA. We explore the relationship of the star formation driven ionised gas outflows with colocated galaxy properties. We find a roughly linear relationship between the outflow mass flux ($\dotΣ_{\rm out}$) and star formation rate surface density ($Σ_{\rm SFR}$), $\dotΣ_{\rm out}\proptoΣ_{\rm SFR}^{1.06\pm0.10}$, and a strong correlation between $\dotΣ_{\rm out}$ and the gas depletion time, such that $\dotΣ_{\rm out} \propto t_{dep}^{-1.1\pm0.06}$. Moreover, we find these outflows are so-called ``breakout" outflows, according to the relationship between the gas fraction and disk kinematics. Assuming that ionised outflow mass scales with total outflow mass, our observations suggest that the regions of highest $Σ_{\rm SFR}$ in IRAS08 are removing more gas via the outflow than through the conversion of gas into stars. Our results are consistent with a picture in which the outflow limits the ability for a region of a disk to maintain short depletion times. Our results underline the need for resolved observations of outflows in more galaxies. △ Less

Submitted 3 November, 2022; originally announced November 2022.

Comments: 16 pages, 7 figures, Submitted to ApJ

VERTICO IV: Environmental Effects on the Gas Distribution and Star Formation Efficiency of Virgo Cluster Spirals

Authors: Vicente Villanueva, Alberto D. Bolatto, Stuart Vogel, Tobias Brown, Christine D. Wilson, Nikki Zabel, Sara Ellison, Adam R. H. Stevens, Maria Jesus Jimenez Donaire, Kristine Spekkens, Mallory Thorp, Timothy A. Davis, Laura C. Parker, Ian D. Roberts, Dhruv Bisaria, Alessandro Boselli, Barbara Catinella, Aeree Chung, Luca Cortese, Bumhyun Lee, Adam Watts

Abstract: We measure the molecular-to-atomic gas ratio, $R_{\rm mol}$, and the star formation rate (SFR) per unit molecular gas mass, SFE$_{\rm mol}$, in 38 nearby galaxies selected from the Virgo Environment Traced in CO (VERTICO) survey. We determine their scale-lengths for the molecular and stellar components and find a roughly 3:5 ratio between them compared to $\sim$1:1 in field galaxies, indicating th… ▽ More We measure the molecular-to-atomic gas ratio, $R_{\rm mol}$, and the star formation rate (SFR) per unit molecular gas mass, SFE$_{\rm mol}$, in 38 nearby galaxies selected from the Virgo Environment Traced in CO (VERTICO) survey. We determine their scale-lengths for the molecular and stellar components and find a roughly 3:5 ratio between them compared to $\sim$1:1 in field galaxies, indicating that the CO emission is more centrally concentrated than the stars. We compute $R_{\rm mol}$ as a function of different physical quantities. While the spatially-resolved $R_{\rm mol}$ on average decreases with increasing radius, we find that the mean molecular-to-atomic gas ratio within the stellar effective radius $R_{\rm e}$, $R_{\rm mol}(r<R_{\rm e})$, shows a systematic increase with the level of H$_{\rm I}$, truncation and/or asymmetry (H$_{\rm I}$ perturbation). Analysis of the molecular- and the atomic-to-stellar mass ratios within $R_{\rm e}$, $R^{\rm mol}_{\star}(r<R_{\rm e})$ and $R^{\rm atom}_{\star}(r<R_{\rm e})$, shows that VERTICO galaxies have increasingly lower $R^{\rm atom}_{\star}(r<R_{\rm e})$ for larger levels of H$_{\rm I}$perturbation (compared to field galaxies matched in stellar mass), but no significant change in $R^{\rm mol}_{\star}(r<R_{\rm e})$. We also measure a clear systematic decrease of the SFE$_{\rm mol}$ within $R_{\rm e}$, SFE$_{\rm mol}(r<R_{\rm e})$, with increasingly perturbed H$_{\rm I}$. Therefore, compared to galaxies from the field, VERTICO galaxies are more compact in CO emission in relation to their stellar distribution, but increasingly perturbed atomic gas increases their $R_{\rm mol}$ and decreases the efficiency with which their molecular gas forms stars. (abridged) △ Less

Submitted 1 November, 2022; v1 submitted 11 October, 2022; originally announced October 2022.

X-ray Properties of NGC 253's Starburst-Driven Outflow

Authors: Sebastian Lopez, Laura A. Lopez, Dustin D. Nguyen, Todd A. Thompson, Smita Mathur, Alberto D. Bolatto, Neven Vulic, Amy Sardone

Abstract: We analyze image and spectral data from $\approx$365~ks of observations from the {\it Chandra} X-ray Observatory of the nearby, edge-on starburst galaxy NGC 253 to constrain properties of the hot phase of the outflow. We focus our analysis on the $-$1.1 to $+$0.63 kpc region of the outflow and define several regions for spectral extraction where we determine best-fit temperatures and metal abundan… ▽ More We analyze image and spectral data from $\approx$365~ks of observations from the {\it Chandra} X-ray Observatory of the nearby, edge-on starburst galaxy NGC 253 to constrain properties of the hot phase of the outflow. We focus our analysis on the $-$1.1 to $+$0.63 kpc region of the outflow and define several regions for spectral extraction where we determine best-fit temperatures and metal abundances. We find that the temperatures and electron densities peak in the central $\sim$250 pc region of the outflow and decrease with distance. These temperature and density profiles are in disagreement with an adiabatic spherically expanding starburst wind model and suggest the presence of additional physics such as mass loading and non-spherical outflow geometry. Our derived temperatures and densities yield few-Myr cooling times in the nuclear region, which may imply that the hot gas can undergo bulk radiative cooling as it escapes along the minor axis. Our metal abundances of O, Ne, Mg, Si, S, and Fe all peak in the central region and decrease with distance along the outflow, with the exception of Ne which maintains a flat distribution. The metal abundances indicate significant dilution outside of the starburst region. We also find estimates on the mass outflow rates which are $2.8\:M_{\odot}/\rm{yr}$ in the northern outflow and $3.2\:M_{\odot}/\rm{yr}$ in the southern outflow. Additionally, we detect emission from charge exchange and find it has a significant contribution ($20-42$\%) to the total broad-band ($0.5-7$ keV) X-ray emission in the central and southern regions of the outflow. △ Less

Submitted 19 January, 2023; v1 submitted 19 September, 2022; originally announced September 2022.

Comments: 15 pages, 9 figures; submitted (19 September 2022), accepted ApJ (22 November 2022)

Extragalactic Science with the Experiment for Cryogenic Large-aperture Intensity Mapping

Authors: Anthony R. Pullen, Patrick C. Breysse, Trevor Oxholm, Eric R. Switzer, Christopher J. Anderson, Emily Barrentine, Alberto D. Bolatto, Giuseppe Cataldo, Thomas Essinger-Hileman, Abhishek Maniyar, Thomas Stevenson, Rachel S. Somerville, Carrie Volpert, Edward Wollack, Shengqi Yang, L. Y. Aaron Yung, Zilu Zhou

Abstract: The EXperiment for Cryogenic Large-Aperture Intensity Mapping (EXCLAIM) is a balloon-borne cryogenic telescope that will survey the spectrum of diffuse emission from both the Milky Way and the cosmic web to probe star formation, the interstellar medium, and galaxy evolution across cosmic time. EXCLAIM's primary extragalactic science survey maps 305 deg$^2$ along the celestial equator with an R=512… ▽ More The EXperiment for Cryogenic Large-Aperture Intensity Mapping (EXCLAIM) is a balloon-borne cryogenic telescope that will survey the spectrum of diffuse emission from both the Milky Way and the cosmic web to probe star formation, the interstellar medium, and galaxy evolution across cosmic time. EXCLAIM's primary extragalactic science survey maps 305 deg$^2$ along the celestial equator with an R=512 spectrometer over the frequency range ν=420-540 GHz, targeting emission of the [CII] line over redshifts 2.5<z<3.5 and several CO lines for z<1. Cross-correlation with galaxy redshift catalogs isolates line emission from the large-scale structure at target redshifts. In this paper, we forecast the sensitivity for both the two-point and conditional one-point cross-correlation. We predict that EXCLAIM will detect both the [CII]-QSO cross-power spectrum and the conditional voxel intensity distribution (CVID) at various redshifts under a broad range of [CII] intensity models, allowing it to differentiate among these models in the literature. These forecasts for the power spectra include the effects of line interlopers and continuum foreground contamination. We then convert the joint [CII] constraints from both the cross-power spectrum and the CVID into constraints on the [CII] halo luminosity-mass relation $L_\mathrm{[CII]}(M)$ model parameters and the star formation rate density (SFRD) from [CII] emission. We also develop sensitivity estimates for CO, showing the ability to differentiate between models. △ Less

Submitted 6 September, 2022; originally announced September 2022.

Comments: 19 pages, 14 figures, will submit to MNRAS

The EDGE-CALIFA survey: The role of spiral arms and bars in driving central molecular gas concentrations

Authors: Si-Yue Yu, Veselina Kalinova, Dario Colombo, Alberto D. Bolatto, Tony Wong, Rebecca C. Levy, Vicente Villanueva, Sebastián F. Sánchez, Luis C. Ho, Stuart N. Vogel, Peter Teuben, Mónica Rubio

Abstract: Shocks and torques produced by non-axisymmetric structures such as spiral arms and bars may transport gas to galaxy central regions. We test this hypothesis by studying the dependence of concentration of CO luminosity ($C_{CO}$), molecular gas ($C_{mol}$), and star formation rate ($C_{SFR}$) in central $\sim$ 2 kpc on the $strength$ of non-axisymmetric disk structure using a sample of 57 disk gala… ▽ More Shocks and torques produced by non-axisymmetric structures such as spiral arms and bars may transport gas to galaxy central regions. We test this hypothesis by studying the dependence of concentration of CO luminosity ($C_{CO}$), molecular gas ($C_{mol}$), and star formation rate ($C_{SFR}$) in central $\sim$ 2 kpc on the $strength$ of non-axisymmetric disk structure using a sample of 57 disk galaxies selected from the EDGE-CALIFA survey. $C_{mol}$ is calculated using a CO-to-H$_2$ conversion factor that decreases with higher metallicity and higher stellar surface density. We find that $C_{mol}$ is systematically 0.22 dex lower than $C_{CO}$. We confirm that high $C_{mol}$ and strong non-axisymmetric disk structure are more common in barred galaxies than in unbarred galaxies. However, we find that spiral arms also increase $C_{mol}$. We show that there is a good correlation between $C_{mol}$ and the $strength$ of non-axisymmetric structure (which can be due to a bar, spiral arms, or both). This suggests that the stronger the bars and spirals, the more efficient the galaxy is at transporting cold gas to its center. Despite the small subsample size, $C_{mol}$ of the four Seyferts are not significantly reduced compared to inactive galaxies of similar disk structure, implying that the AGN feedback in Seyferts may not notably affect the molecular gas distribution in the central $\sim$2kpc. We find that $C_{SFR}$ tightly correlates with $C_{mol}$ in both unbarred and barred galaxies. Likewise, elevated $C_{SFR}$ is found in galaxies with strong disk structure. Our results suggest that the disk structure, either spirals or bars, can transport gas to the central regions, with higher inflow rates corresponding to stronger structure, and consequently boost central star formation. Both spirals and bars play, therefore, an essential role in the secular evolution of disk galaxies. △ Less

Submitted 9 September, 2022; v1 submitted 31 August, 2022; originally announced August 2022.

Comments: 24 pages; accepted for publication in Astronomy & Astrophysics; typos corrected

Journal ref: A&A 666, A175 (2022)

Sequential Star Formation in the Young SMC Region NGC 602: Insights from ALMA

Authors: Theo J. O'Neill, Remy Indebetouw, Karin Sandstrom, Alberto D. Bolatto, Katherine E. Jameson, Lynn R. Carlson, Molly K. Finn, Margaret Meixner, Elena Sabbi, Marta Sewilo

Abstract: NGC 602 is a young, low-metallicity star cluster in the "Wing" of the Small Magellanic Cloud. We reveal the recent evolutionary past of the cluster through analysis of high-resolution ($\sim$0.4 pc) Atacama Large Millimeter/submillimeter Array observations of molecular gas in the associated $\textrm{H}\scriptstyle\mathrm{II}$ region N90. We identify 110 molecular clumps ($R <$ 0.8 pc) traced by CO… ▽ More NGC 602 is a young, low-metallicity star cluster in the "Wing" of the Small Magellanic Cloud. We reveal the recent evolutionary past of the cluster through analysis of high-resolution ($\sim$0.4 pc) Atacama Large Millimeter/submillimeter Array observations of molecular gas in the associated $\textrm{H}\scriptstyle\mathrm{II}$ region N90. We identify 110 molecular clumps ($R <$ 0.8 pc) traced by CO emission, and study the relationship between the clumps and associated young stellar objects (YSOs) and pre-main-sequence (PMS) stars. The clumps have high virial parameters (typical $α_{\rm{vir}} = $ 4-11) and may retain signatures of a collision in the last $\lesssim$8 Myr between $\textrm{H}\scriptstyle\mathrm{I}$ components of the adjacent supergiant shell SMC-SGS 1. We obtain a CO-bright-to-H$_2$ gas conversion factor of $X_{CO,B} = (3.4 \pm 0.2) \times 10^{20}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$, and correct observed clump properties for CO-dark H$_2$ gas to derive a total molecular gas mass in N90 of $16,600 \pm 2,400 \ M_\odot$. We derive a recent ($\lesssim 1$ Myr) star formation rate of $130 \pm 30 \ M_{\odot}$ Myr$^{-1}$ with an efficiency of 8 $ \pm$ 3\% assessed through comparing total YSO mass to total molecular gas mass. Very few significant radial trends exist between clump properties or PMS star ages and distance from NGC 602. We do not find evidence for a triggered star formation scenario among the youngest ($\lesssim$2 Myr) stellar generations, and instead conclude that a sequential star formation process in which NGC 602 did not directly cause recent star formation in the region is likely. △ Less

Submitted 29 August, 2022; originally announced August 2022.

Comments: 31 pages, 13 figures, accepted for publication in ApJ

Developing a New Generation of Integrated Micro-Spec Far Infrared Spectrometers for the EXperiment for Cryogenic Large-Aperture Intensity Mapping (EXCLAIM)

Authors: Carolyn G. Volpert, Emily M. Barrentine, Mona Mirzaei, Alyssa Barlis, Alberto D. Bolatto, Berhanu Bulcha, Giuseppe Cataldo, Jake A. Connors, Nicholas Costen, Negar Ehsan, Thomas Essinger-Hileman, Jason Glenn, James P. Hays-Wehle, Larry A. Hess, Alan J. Kogut, Harvey Moseley, Jonas Mugge-Durum, Omid Noroozian, Trevor M. Oxholm, Maryam Rahmani, Thomas Stevenson, Eric R. Switzer, Joseph Watson, Edward J. Wollack

Abstract: The current state of far-infrared astronomy drives the need to develop compact, sensitive spectrometers for future space and ground-based instruments. Here we present details of the $\rm μ$-Spec spectrometers currently in development for the far-infrared balloon mission EXCLAIM. The spectrometers are designed to cover the $\rm 555 - 714\ μ$m range with a resolution of $\rm R\ =\ λ/ Δλ =\ 512$ at t… ▽ More The current state of far-infrared astronomy drives the need to develop compact, sensitive spectrometers for future space and ground-based instruments. Here we present details of the $\rm μ$-Spec spectrometers currently in development for the far-infrared balloon mission EXCLAIM. The spectrometers are designed to cover the $\rm 555 - 714\ μ$m range with a resolution of $\rm R\ =\ λ/ Δλ =\ 512$ at the $\rm 638\ μ$m band center. The spectrometer design incorporates a Rowland grating spectrometer implemented in a parallel plate waveguide on a low-loss single-crystal Si chip, employing Nb microstrip planar transmission lines and thin-film Al kinetic inductance detectors (KIDs). The EXCLAIM $\rm μ$-Spec design is an advancement upon a successful $\rm R = 64\ μ$-Spec prototype, and can be considered a sub-mm superconducting photonic integrated circuit (PIC) that combines spectral dispersion and detection. The design operates in a single $M{=}2$ grating order, allowing one spectrometer to cover the full EXCLAIM band without requiring a multi-order focal plane. The EXCLAIM instrument will fly six spectrometers, which are fabricated on a single 150 mm diameter Si wafer. Fabrication involves a flip-wafer-bonding process with patterning of the superconducting layers on both sides of the Si dielectric. The spectrometers are designed to operate at 100 mK, and will include 355 Al KID detectors targeting a goal of NEP ${\sim}8\times10^{-19}$ $\rm W/\sqrt{Hz}$. We summarize the design, fabrication, and ongoing development of these $\rm μ$-Spec spectrometers for EXCLAIM. △ Less

Submitted 4 August, 2022; originally announced August 2022.

Comments: 9 pages, 5 figures, to appear in the Proceedings of the SPIE Astronomical Telescopes + Instrumentation (2022)

Molecular Cloud Populations in the Context of Their Host Galaxy Environments: A Multiwavelength Perspective

Authors: Jiayi Sun, Adam K. Leroy, Erik Rosolowsky, Annie Hughes, Eva Schinnerer, Andreas Schruba, Eric W. Koch, Guillermo A. Blanc, I-Da Chiang, Brent Groves, Daizhong Liu, Sharon Meidt, Hsi-An Pan, Jerome Pety, Miguel Querejeta, Toshiki Saito, Karin Sandstrom, Amy Sardone, Antonio Usero, Dyas Utomo, Thomas G. Williams, Ashley T. Barnes, Samantha M. Benincasa, Frank Bigiel, Alberto D. Bolatto , et al. (13 additional authors not shown)

Abstract: We present a rich, multiwavelength, multiscale database built around the PHANGS-ALMA CO$\,$(2-1) survey and ancillary data. We use this database to present the distributions of molecular cloud populations and sub-galactic environments in 80 PHANGS galaxies, to characterize the relationship between population-averaged cloud properties and host galaxy properties, and to assess key timescales relevan… ▽ More We present a rich, multiwavelength, multiscale database built around the PHANGS-ALMA CO$\,$(2-1) survey and ancillary data. We use this database to present the distributions of molecular cloud populations and sub-galactic environments in 80 PHANGS galaxies, to characterize the relationship between population-averaged cloud properties and host galaxy properties, and to assess key timescales relevant to molecular cloud evolution and star formation. We show that PHANGS probes a wide range of kpc-scale gas, stellar, and star formation rate (SFR) surface densities, as well as orbital velocities and shear. The population-averaged cloud properties in each aperture correlate strongly with both local environmental properties and host galaxy global properties. Leveraging a variable selection analysis, we find that the kpc-scale surface densities of molecular gas and SFR tend to possess the most predictive power for the population-averaged cloud properties. Once their variations are controlled for, galaxy global properties contain little additional information, which implies that the apparent galaxy-to-galaxy variations in cloud populations are likely mediated by kpc-scale environmental conditions. We further estimate a suite of important timescales from our multiwavelength measurements. The cloud-scale free-fall time and turbulence crossing time are ${\sim}5{-}20$ Myr, comparable to previous cloud lifetime estimates. The timescales for orbital motion, shearing, and cloud-cloud collisions are longer, ${\sim}100$ Myr. The molecular gas depletion time is $1{-}3$ Gyr and shows weak to no correlations with the other timescales in our data. We publish our measurements online and expect them to have broad utility to future studies of molecular clouds and star formation. △ Less

Submitted 14 June, 2022; originally announced June 2022.

Comments: 32 pages + 6 appendices. AJ in press. Data products available at https://www.canfar.net/storage/list/phangs/RELEASES/Sun_etal_2022 . Associated software package available at https://doi.org/10.5281/zenodo.6584842 . Slides summarizing the key results can be found at https://www.dropbox.com/s/7llnrob9wi0isfk/Sun_et_PHANGS_2022.pptx?dl=0

Journal ref: AJ, 164, 43 (2022)

The 30 Doradus Molecular Cloud at 0.4 pc Resolution with the Atacama Large Millimeter/submillimeter Array: Physical Properties and the Boundedness of CO-emitting Structures

Authors: Tony Wong, Luuk Oudshoorn, Eliyahu Sofovich, Alex Green, Charmi Shah, Rémy Indebetouw, Margaret Meixner, Alvaro Hacar, Omnarayani Nayak, Kazuki Tokuda, Alberto D. Bolatto, Mélanie Chevance, Guido De Marchi, Yasuo Fukui, Alec S. Hirschauer, K. E. Jameson, Venu Kalari, Vianney Lebouteiller, Leslie W. Looney, Suzanne C. Madden, Toshikazu Onishi, Julia Roman-Duval, Mónica Rubio, A. G. G. M. Tielens

Abstract: We present results of a wide-field (approximately 60 x 90 pc) ALMA mosaic of CO(2-1) and $^{13}$CO(2-1) emission from the molecular cloud associated with the 30 Doradus star-forming region. Three main emission complexes, including two forming a bowtie-shaped structure extending northeast and southwest from the central R136 cluster, are resolved into complex filamentary networks. Consistent with pr… ▽ More We present results of a wide-field (approximately 60 x 90 pc) ALMA mosaic of CO(2-1) and $^{13}$CO(2-1) emission from the molecular cloud associated with the 30 Doradus star-forming region. Three main emission complexes, including two forming a bowtie-shaped structure extending northeast and southwest from the central R136 cluster, are resolved into complex filamentary networks. Consistent with previous studies, we find that the central region of the cloud has higher line widths at fixed size relative to the rest of the molecular cloud and to other LMC clouds, indicating an enhanced level of turbulent motions. However, there is no clear trend in gravitational boundedness (as measured by the virial parameter) with distance from R136. Structures observed in $^{13}$CO are spatially coincident with filaments and are close to a state of virial equilibrium. In contrast, CO structures vary greatly in virialization, with low CO surface brightness structures outside of the main filamentary network being predominantly unbound. The low surface brightness structures constitute ~10% of the measured CO luminosity; they may be shredded remnants of previously star-forming gas clumps, or alternatively the CO-emitting parts of more massive, CO-dark structures. △ Less

Submitted 13 June, 2022; originally announced June 2022.

Comments: 24 pages; published in ApJ on 15 Jun 2022; associated data can be found at https://doi.org/10.13012/B2IDB-1671495_V1

Journal ref: ApJ, 932, 47 (2022)

The Morpho-Kinematic Architecture of Super Star Clusters in the Center of NGC253

Authors: Rebecca C. Levy, Alberto D. Bolatto, Adam K. Leroy, Mattia C. Sormani, Kimberly L. Emig, Mark Gorski, Laura Lenkić, Elisabeth A. C. Mills, Elizabeth Tarantino, Peter Teuben, Sylvain Veilleux, Fabian Walter

Abstract: The center of the nearby galaxy NGC\,253 hosts a population of more than a dozen super star clusters (SSCs) which are still in the process of forming. The majority of the star formation of the burst is concentrated in these SSCs, and the starburst is powering a multiphase outflow from the galaxy. In this work, we measure the 350~GHz dust continuum emission towards the center of NGC\,253 at 47~mill… ▽ More The center of the nearby galaxy NGC\,253 hosts a population of more than a dozen super star clusters (SSCs) which are still in the process of forming. The majority of the star formation of the burst is concentrated in these SSCs, and the starburst is powering a multiphase outflow from the galaxy. In this work, we measure the 350~GHz dust continuum emission towards the center of NGC\,253 at 47~milliarcsecond (0.8 pc) resolution using data from the Atacama Large Millimeter/submillimeter Array (ALMA). We report the detection of 350~GHz (dust) continuum emission in the outflow for the first time, associated with the prominent South-West streamer. In this feature, the dust emission has a width of $\approx$~8~pc, is located at the outer edge of the CO emission, and corresponds to a molecular gas mass of $\sim~(8-17)\times10^6$~M$_\odot$. In the starburst nucleus, we measure the resolved radial profiles, sizes, and molecular gas masses of the SSCs. Compared to previous work at somewhat lower spatial resolution, the SSCs here break apart into smaller substructures with radii $0.4-0.7$~pc. In projection, the SSCs, dust, and dense molecular gas appear to be arranged as a thin, almost linear, structure roughly 155~pc in length. The morphology and kinematics of this structure can be well explained as gas following $x_2$ orbits at the center of a barred potential. We constrain the morpho-kinematic arrangement of the SSCs themselves, finding that an elliptical, angular momentum-conserving ring is a good description of the both morphology and kinematics of the SSCs. △ Less

Submitted 22 June, 2022; v1 submitted 9 June, 2022; originally announced June 2022.

Comments: 26 pages, 13 figures, accepted for publication in ApJ

Effects of CO-dark Gas on Measurements of Molecular Cloud Stability and the Size-Linewidth Relationship

Authors: Theo J. O'Neill, Remy Indebetouw, Alberto D. Bolatto, Suzanne C. Madden, Tony Wong

Abstract: Stars form within molecular clouds, so characterizing the physical states of molecular clouds is key in understanding the process of star formation. Cloud structure and stability is frequently assessed using metrics including the virial parameter and Larson (1981) scaling relationships between cloud radius, velocity dispersion, and surface density. Departures from the typical Galactic relationship… ▽ More Stars form within molecular clouds, so characterizing the physical states of molecular clouds is key in understanding the process of star formation. Cloud structure and stability is frequently assessed using metrics including the virial parameter and Larson (1981) scaling relationships between cloud radius, velocity dispersion, and surface density. Departures from the typical Galactic relationships between these quantities have been observed in low-metallicity environments. The amount of H$_2$ gas in cloud envelopes without corresponding CO emission is expected to be high under these conditions; therefore, this "CO-dark" gas could plausibly be responsible for the observed variations in cloud properties. We derive simple corrections that can be applied to empirical clump properties (mass, radius, velocity dispersion, surface density, and virial parameter) to account for CO-dark gas in clumps following power-law and Plummer mass density profiles. We find that CO-dark gas is not likely to be the cause of departures from Larson's relationships in low-metallicity regions, but that virial parameters may be systematically overestimated. We demonstrate that correcting for CO-dark gas is critical for accurately comparing the dynamical state and evolution of molecular clouds across diverse environments. △ Less

Submitted 26 May, 2022; originally announced May 2022.

Comments: 19 pages, 5 figures, accepted for publication in ApJ

Snowmass2021 Cosmic Frontier White Paper: Observational Facilities to Study Dark Matter

Authors: Sukanya Chakrabarti, Alex Drlica-Wagner, Ting S. Li, Neelima Sehgal, Joshua D. Simon, Simon Birrer, Duncan A. Brown, Rebecca Bernstein, Alberto D. Bolatto, Philip Chang, Kyle Dawson, Paul Demorest, Daniel Grin, David L. Kaplan, Joseph Lazio, Jennifer Marshall, Eric J. Murphy, Scott Ransom, Brant E. Robertson, Rajeev Singh, Anže Slosar, Tommaso Treu, Yu-Dai Tsai, Benjamin F. Williams

Abstract: We present an overview of future observational facilities that will significantly enhance our understanding of the fundamental nature of dark matter. These facilities span a range of observational techniques including optical/near-infrared imaging and spectroscopy, measurements of the cosmic microwave background, pulsar timing, 21-cm observations of neutral hydrogen at high redshift, and the measu… ▽ More We present an overview of future observational facilities that will significantly enhance our understanding of the fundamental nature of dark matter. These facilities span a range of observational techniques including optical/near-infrared imaging and spectroscopy, measurements of the cosmic microwave background, pulsar timing, 21-cm observations of neutral hydrogen at high redshift, and the measurement of gravitational waves. Such facilities are a critical component of a multi-pronged experimental program to uncover the nature of dark matter, while often providing complementary measurements of dark energy, neutrino physics, and inflation. △ Less

Submitted 11 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021. 31 pages, 6 figures, 2 tables

Cuspy dark matter density profiles in massive dwarf galaxies

Authors: Lauren H. Cooke, Rebecca C. Levy, Alberto D. Bolatto, Joshua D. Simon, Andrew B. Newman, Peter Teuben, Brandon D. Davey, Melvyn Wright, Elizabeth Tarantino, Laura Lenkić, Vicente Villanueva

Abstract: Rotation curves of galaxies probe their total mass distributions, including dark matter. Dwarf galaxies are excellent systems to investigate the dark matter density distribution, as they tend to have larger fractions of dark matter compared to higher mass systems. The core-cusp problem describes the discrepancy found in the slope of the dark matter density profile in the centres of galaxies (… ▽ More Rotation curves of galaxies probe their total mass distributions, including dark matter. Dwarf galaxies are excellent systems to investigate the dark matter density distribution, as they tend to have larger fractions of dark matter compared to higher mass systems. The core-cusp problem describes the discrepancy found in the slope of the dark matter density profile in the centres of galaxies ($β^*$) between observations of dwarf galaxies (shallower cores) and dark matter-only simulations (steeper cusps). We investigate $β^*$ in six nearby spiral dwarf galaxies for which high-resolution CO $J=1-0$ data were obtained with ALMA. We derive rotation curves and decompose the mass profile of the dark matter using our CO rotation curves as a tracer of the total potential and 4.5$μ$m photometry to define the stellar mass distribution. We find $\langleβ^*\rangle = 0.6$ with a standard deviation of $\pm0.1$ among the galaxies in this sample, in agreement with previous measurements in this mass range. The galaxies studied are on the high stellar mass end of dwarf galaxies and have cuspier profiles than lower mass dwarfs, in agreement with other observations. When the same definition of the slope is used, we observe steeper slopes than predicted by the FIRE and NIHAO simulations. This may signal that these relatively massive dwarfs underwent stronger gas inflows toward their centres than predicted by these simulations, that these simulations over-predict the frequency of accretion or feedback events, or that a combination of these or other effects are at work. △ Less

Submitted 1 March, 2022; originally announced March 2022.

Comments: Accepted for publication in MNRAS, 21 pages (4 pages of appendix material), 11 figures (3 appendix figures), code and data products available here: https://github.com/rclevy/RotationCurveTiltedRings

Kiloparsec view of a typical star-forming galaxy when the Universe was $\sim$1 Gyr old II. Regular rotating disk and evidence for baryon dominance on galactic scales

Authors: R. Herrera-Camus, N. M. Förster Schreiber, S. H. Price, H. Übler, A. D. Bolatto, R. L. Davies, D. Fisher, R. Genzel, D. Lutz, T. Naab, A. Nestor, T. Shimizu, A. Sternberg, L. Tacconi, K. Tadaki

Abstract: We present a kinematic analysis of the main-sequence galaxy HZ4 at $z=5.5$. Our study is based on deep, spatially resolved observations of the [CII] 158 $μ$m transition obtained with the Atacama Large Millimeter/Submillimeter Array (ALMA). From the combined analysis of the disk morphology, the two-dimensional velocity structure, and forward-modeling of the one-dimensional velocity and velocity dis… ▽ More We present a kinematic analysis of the main-sequence galaxy HZ4 at $z=5.5$. Our study is based on deep, spatially resolved observations of the [CII] 158 $μ$m transition obtained with the Atacama Large Millimeter/Submillimeter Array (ALMA). From the combined analysis of the disk morphology, the two-dimensional velocity structure, and forward-modeling of the one-dimensional velocity and velocity dispersion profiles, we conclude that HZ4 has a regular rotating disk in place. The intrinsic velocity dispersion in HZ4 is high ($σ_{0}=65.8^{+2.9}_{-3.3}$ km s$^{-1}$), and the ratio between the rotational velocity and the intrinsic velocity dispersion is $V_{\rm rot}/σ_{0}=2.2$. These values are consistent with the expectations from the trends of increasing $σ_{0}$ and decreasing $V_{\rm rot}/σ_{0}$ as a function of redshift observed in main-sequence galaxies up to $z\approx4$. Galaxy evolution models suggest that the high level of turbulence observed in HZ4 can only be achieved if, in addition to stellar feedback, there is radial transport of gas within the disk. Finally, we find that HZ4 is baryon dominated on galactic scales ($\lesssim2\times R_{\rm e}$), with a dark matter fraction at one effective radius of $f_{\rm DM}(R_{\rm e})=0.41^{+0.25}_{-0.22}$. This value is comparable to the dark matter fractions found in lower redshift galaxies that could be the descendants of HZ4: massive ($M_{\star}\approx10^{11}~M_{\odot}$), star-forming galaxies at $z\sim2$, and passive, early type galaxies at $z\approx0$. △ Less

Submitted 1 March, 2022; originally announced March 2022.

Comments: Accepted for publication in A&A Letters, 10 pages, 6 figures

Journal ref: A&A 665, L8 (2022)

Extreme Variation in Star Formation Efficiency Across a Compact, Starbursting Disk Galaxy

Authors: D. B. Fisher, A. D. Bolatto, K. Glazebrook, D. Obreschkow, R. G. Obreschkow, G. G. Kacprzak, N. M. Nielsen

Abstract: We report on the internal distribution of star formation efficiency in IRAS 08339+6517 (hereafter IRAS08), using $\sim$200~pc resolution CO(2-1) observations from NOEMA. The molecular gas depletion time changes by 2 orders-of-magnitude from disk-like values in the outer parts to less than 10$^8$~yr inside the half-light radius. This translates to a star formation efficiency per free-fall time that… ▽ More We report on the internal distribution of star formation efficiency in IRAS 08339+6517 (hereafter IRAS08), using $\sim$200~pc resolution CO(2-1) observations from NOEMA. The molecular gas depletion time changes by 2 orders-of-magnitude from disk-like values in the outer parts to less than 10$^8$~yr inside the half-light radius. This translates to a star formation efficiency per free-fall time that also changes by 2 orders-of-magnitude, reaching 50-100\%, different than local spiral galaxies and typical assumption of constant, low star formation efficiencies. Our target is a compact, massive disk galaxy that has SFR 10$\times$ above the $z=0$ main-sequence; Toomre $Q\approx0.5-0.7$ and high gas velocity dispersion ($σ_{mol}\approx 25$~km~s$^{-1}$). We find that IRAS08 is similar to other rotating, starburst galaxies from the literature in the resolved $Σ_{SFR}\proptoΣ_{mol}^N$ relation. By combining resolved literature studies we find that distance from the main-sequence is a strong indicator of the Kennicutt-Schmidt powerlaw slope, with slopes of $N\approx1.6$ for starbursts from 100-10$^4$~M$_{\odot}$~pc$^{-2}$. Our target is consistent with a scenario in which violent disk instabilities drive rapid inflows of gas. It has low values of Toomre-$Q$, and also at all radii the inflow timescale of the gas is less than the depletion time, which is consistent with the flat metallicity gradients in IRAS08. We consider these results in light of popular star formation theories, in general observations of IRAS08 find the most tension with theories in which star formation efficiency is a constant. Our results argue for the need of high spatial resolution CO observations are a larger number of similar targets. △ Less

Submitted 31 January, 2022; originally announced February 2022.

Comments: Accepted to ApJ

Molecular Gas Properties and CO-to-H2 Conversion Factors in the Central Kiloparsec of NGC 3351

Authors: Yu-Hsuan Teng, Karin M. Sandstrom, Jiayi Sun, Adam K. Leroy, L. Clifton Johnson, Alberto D. Bolatto, J. M. Diederik Kruijssen, Andreas Schruba, Antonio Usero, Ashley T. Barnes, Frank Bigiel, Guillermo A. Blanc, Brent Groves, Frank P. Israel, Daizhong Liu, Erik Rosolowsky, Eva Schinnerer, J. D. Smith, Fabian Walter

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

Submitted 28 December, 2021; v1 submitted 10 November, 2021; originally announced November 2021.

Comments: 31 pages, 17 figures; accepted for publication in ApJ; font size of text and figures adjusted

VERTICO: The Virgo Environment Traced In CO Survey

Authors: Toby Brown, Christine D. Wilson, Nikki Zabel, Timothy A. Davis, Alessandro Boselli, Aeree Chung, Sara L. Ellison, Claudia D. P. Lagos, Adam R. H. Stevens, Luca Cortese, Yannick M. Bahé, Dhruv Bisaria, Alberto D. Bolatto, Claire R. Cashmore, Barbara Catinella, Ryan Chown, Benedikt Diemer, Pascal J. Elahi, Maan H. Hani, María J. Jiménez-Donaire, Bumhyun Lee, Katya Leidig, Angus Mok, Karen Pardos Olsen, Laura C. Parker , et al. (11 additional authors not shown)

Abstract: We present the Virgo Environment Traced in CO (VERTICO) survey, a new effort to map $^{12}$CO($2-1$), $^{13}$CO($2-1$), and C$^{18}$O($2-1$) in 51 Virgo Cluster galaxies with the Atacama Compact Array, part of the Atacama Large Millimeter/submillimeter Array (ALMA). The primary motivation of VERTICO is to understand the physical mechanisms that perturb molecular gas disks, and therefore star forma… ▽ More We present the Virgo Environment Traced in CO (VERTICO) survey, a new effort to map $^{12}$CO($2-1$), $^{13}$CO($2-1$), and C$^{18}$O($2-1$) in 51 Virgo Cluster galaxies with the Atacama Compact Array, part of the Atacama Large Millimeter/submillimeter Array (ALMA). The primary motivation of VERTICO is to understand the physical mechanisms that perturb molecular gas disks, and therefore star formation and galaxy evolution, in dense environments. This first paper contains an overview of VERTICO's design and sample selection, $^{12}$CO($2-1$) observations, and data reduction procedures. We characterize global $^{12}$CO($2-1$) fluxes and molecular gas masses for the 49 detected VERTICO galaxies, provide upper limits for the two non-detections, and produce resolved $^{12}$CO($2-1$) data products (median resolution $= 8^{\prime\prime} \approx 640~{\rm pc}$). Azimuthally averaged $^{12}$CO($2-1$) radial intensity profiles are presented along with derived molecular gas radii. We demonstrate the scientific power of VERTICO by comparing the molecular gas size--mass scaling relation for our galaxies with a control sample of field galaxies, highlighting the strong effect that radius definition has on this correlation. We discuss the drivers of the form and scatter in the size--mass relation and highlight areas for future work. VERTICO is an ideal resource for studying the fate of molecular gas in cluster galaxies and the physics of environment-driven processes that perturb the star formation cycle. Upon public release, the survey will provide a homogeneous legacy dataset for studying galaxy evolution in our closest cluster. △ Less

Submitted 1 November, 2021; originally announced November 2021.

Comments: 68 pages, 13 Figures, 2 Figure Sets, Accepted for publication in ApJS, Online FITS versions of Tables 1, 2, and 3 are available with the journal publication

HI mapping of the Leo Triplet: Morphologies and kinematics of tails and bridges

Authors: Gang Wu, David Martínez-Delgado, Christian Henkel, Pavel Kroupa, Fabian Walter, Nico Krieger, Alberto D. Bolatto, Timothy Robishaw, Joshua D. Simon, Álvaro Ibáñez Pérez, Karl M. Menten, Jarken Esimbek

Abstract: A fully-sampled and hitherto highest resolution and sensitivity observation of neutral hydrogen (HI) in the Leo Triplet (NGC 3628, M 65/NGC 3623, and M 66/NGC 3627) reveals six HI structures beyond the three galaxies. We present detailed results of the morphologies and kinematics of these structures, which can be used for future simulations. In particular, we detect a two-arm structure in the plum… ▽ More A fully-sampled and hitherto highest resolution and sensitivity observation of neutral hydrogen (HI) in the Leo Triplet (NGC 3628, M 65/NGC 3623, and M 66/NGC 3627) reveals six HI structures beyond the three galaxies. We present detailed results of the morphologies and kinematics of these structures, which can be used for future simulations. In particular, we detect a two-arm structure in the plume of NGC 3628 for the first time, which can be explained by a tidal interaction model. The optical counterpart of the plume is mainly associated with the southern arm. The connecting part (base) of the plume (directed eastwards) with NGC 3628 is located at the blueshifted (western) side of NGC 3628. Two bases appear to be associated with the two arms of the plume. A clump with reversed velocity gradient (relative to the velocity gradient of M 66) and a newly detected tail, i.e. M 66SE, is found in the southeast of M 66. We suspect that M 66SE represents gas from NGC 3628 which was captured by M 66 in the recent interaction between the two galaxies. Meanwhile gas is falling toward M 66, resulting in features already previously observed in the southeastern part of M 66, e.g. large line widths and double peaks. An upside-down `Y'-shaped HI gas component (M 65S) is detected in the south of M 65 which suggests that M 65 may also have been involved in the interaction. We strongly encourage modern hydrodynamical simulations of this interacting group of galaxies to reveal the origin of the gaseous debris surrounding all three galaxies. △ Less

Submitted 25 October, 2021; originally announced October 2021.

Comments: 21 pages, 12 figures, accepted for publication in A&A

Journal ref: A&A 658, A25 (2022)

Low-J CO Line Ratios From Single Dish CO Mapping Surveys and PHANGS-ALMA

Authors: Adam K. Leroy, Erik Rosolowsky, Antonio Usero, Karin Sandstrom, Eva Schinnerer, Andreas Schruba, Alberto D. Bolatto, Jiayi Sun, Ashley. T. Barnes, Francesco Belfiore, Frank Bigiel, Jakob S. den Brok, Yixian Cao, I-Da Chiang, Mélanie Chevance, Daniel A. Dale, Cosima Eibensteiner, Christopher M. Faesi, Simon C. O. Glover, Annie Hughes, Maria J. Jiménez Donaire, Ralf S. Klessen, Eric W. Koch, J. M. Diederik Kruijssen, Daizhong Liu , et al. (10 additional authors not shown)

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

Submitted 23 September, 2021; originally announced September 2021.

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

ALMA Imaging of a Galactic Molecular Outflow in NGC4945

Authors: Alberto D. Bolatto, Adam K. Leroy, Rebecca C. Levy, David S. Meier, Elisabeth A. C. Mills, Todd A. Thompson, Kimberly L. Emig, Sylvain Veilleux, Juergen Ott, Mark Gorski, Fabian Walter, Laura A. Lopez, Laura Lenkic

Abstract: We present the ALMA detection of molecular outflowing gas in the central regions of NGC4945, one of the nearest starbursts and also one of the nearest hosts of an active galactic nucleus (AGN). We detect four outflow plumes in CO (3-2) at ~0.3" resolution that appear to correspond to molecular gas located near the edges of the known ionized outflow cone and its (unobserved) counterpart behind the… ▽ More We present the ALMA detection of molecular outflowing gas in the central regions of NGC4945, one of the nearest starbursts and also one of the nearest hosts of an active galactic nucleus (AGN). We detect four outflow plumes in CO (3-2) at ~0.3" resolution that appear to correspond to molecular gas located near the edges of the known ionized outflow cone and its (unobserved) counterpart behind the disk. The fastest and brightest of these plumes has emission reaching observed line-of-sight projected velocities of over 450 km/s beyond systemic, equivalent to an estimated physical outflow velocity v>600 km/s for the fastest emission. Most of these plumes have corresponding emission in HCN or HCO+ (4-3). We discuss a kinematic model for the outflow emission where the molecular gas has the geometry of the ionized gas cone and shares the rotation velocity of the galaxy when ejected. We use this model to explain the velocities we observe, constrain the physical speed of the ejected material, and account for the fraction of outflowing gas that is not detected due to confusion with the galaxy disk. We estimate a total molecular mass outflow rate dMmol/dt~20 Msun/yr flowing through a surface within 100 pc of the disk midplane, likely driven by a combination of the central starburst and AGN. △ Less

Submitted 1 October, 2021; v1 submitted 21 September, 2021; originally announced September 2021.

Comments: 20 pages, 12 figures, accepted by The Astrophysical Journal on September 28, 2021

Giant Star Forming Complexes in High-z Main Sequence Galaxy Analogues: The Internal Structure of Clumps in DYNAMO Galaxies

Authors: Laura Lenkić, Alberto D. Bolatto, Deanne B. Fisher, Karl Glazebrook, Danail Obreschkow, Roberto Abraham, Liyualem Ambachew

Abstract: To indirectly study the internal structure of giant clumps in main sequence galaxies at $z \sim 1-3$, we target very turbulent and gas-rich local analogues from the DYNAMO sample with the Hubble Space Telescope, over a wavelength range of $\sim 200-480$ nm. We present a catalog of 58 clumps identified in six DYNAMO galaxies, including the WFC3/UVIS F225W, F336W, and F467M photometry where the (… ▽ More To indirectly study the internal structure of giant clumps in main sequence galaxies at $z \sim 1-3$, we target very turbulent and gas-rich local analogues from the DYNAMO sample with the Hubble Space Telescope, over a wavelength range of $\sim 200-480$ nm. We present a catalog of 58 clumps identified in six DYNAMO galaxies, including the WFC3/UVIS F225W, F336W, and F467M photometry where the ($225-336$) and ($336-467$) colours are sensitive to extinction and stellar population age respectively. We measure the internal colour gradients of clumps themselves to study their age and extinction properties. We find a marked colour trend within individual clumps, where the resolved colour distributions show that clumps generally have bluer ($336-467$) colours (denoting very young ages) in their centers than at their edges, with little variation in the ($225-336$) colour associated with extinction. Furthermore, we find that clumps whose colours suggest they are older, are preferentially located closer toward the centers of their galaxies, and we find no young clumps at small galactocentric distances. Both results are consistent with simulations of high-redshift star forming systems that show clumps form via violent disk instability, and through dynamic processes migrate to the centers of their galaxies to contribute to bulge growth on timescales of a few 100 Myr, while continually forming stars in their centers. When we compare the DYNAMO clumps to those in these simulations, we find the best agreement with the long-lived clumps. △ Less

Submitted 2 July, 2021; originally announced July 2021.

Comments: 23 pages, 15 figures, resubmitted to MNRAS after first review