Researcher profile

Toshiki Saito

Toshiki Saito contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
29works
0followers
2topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

29 published item(s)

preprint2024arXiv

The ALCHEMI atlas: principal component analysis reveals starburst evolution in NGC 253

Molecular lines are powerful diagnostics of the physical and chemical properties of the interstellar medium (ISM). These ISM properties, which affect future star formation, are expected to differ in starburst galaxies from those of more quiescent galaxies. We investigate the ISM properties in the central molecular zone of the nearby starburst galaxy NGC 253 using the ultra-wide millimeter spectral scan survey from the ALMA Large Program ALCHEMI. We present an atlas of velocity-integrated images at a 1".6 resolution of 148 unblended transitions from 44 species, including the first extragalactic detection of HCNH$^+$ and the first interferometric images of C$_3$H$^+$, NO, HCS$^+$. We conduct a principal component analysis (PCA) on these images to extract correlated chemical species and to identify key groups of diagnostic transitions. To the best of our knowledge, our dataset is currently the largest astronomical set of molecular lines to which PCA has been applied. The PCA can categorize transitions coming from different physical components in NGC 253 such as i) young starburst tracers characterized by high-excitation transitions of HC$_3$N and complex organic molecules (COMs) versus tracers of on-going star formation (radio recombination lines) and high-excitation transitions of CCH and CN tracing PDRs, ii) tracers of cloud-collision-induced shocks (low-excitation transitions of CH$_3$OH, HNCO, HOCO$^+$, and OCS) versus shocks from star-formation-induced outflows (high-excitation transitions of SiO), as well as iii) outflows showing emission from HOC$^+$, CCH, H$_3$O$^+$, CO isotopologues, HCN, HCO$^+$, CS, and CN. Our findings show these intensities vary with galactic dynamics, star formation activities, and stellar feedback.

preprint2022arXiv

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

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

preprint2022arXiv

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

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

preprint2022arXiv

AGN-driven Cold Gas Outflow of NGC 1068 Characterized by Dissociation-Sensitive Molecules

Recent developments in (sub-)millimeter facilities have drastically changed the amount of information obtained from extragalactic spectral scans. In this paper, we present a feature extraction technique using principal component analysis (PCA) applied to arcsecond-resolution (1.0-2.0 arcsec = 72-144 pc) spectral scan datasets for the nearby type-2 Seyfert galaxy, NGC 1068, using Band 3 of the Atacama Large Millimeter/submillimeter Array. We apply PCA to 16 well-detected molecular line intensity maps convolved to a common 150 pc resolution. In addition, we include the [SIII]/[SII] line ratio and [CI] $^3P_1$-$^3P_0$ maps in the literature, both of whose distributions show remarkable resemblance with that of a kpc-scale biconical outflow from the central AGN. We identify two prominent features: (1) central concentration at the circumnuclear disk (CND) and (2) two peaks across the center that coincide with the biconical outflow peaks. The concentrated molecular lines in the CND are mostly high-dipole molecules (e.g., H$^{13}$CN, HC$_3$N, and HCN). Line emissions from molecules known to be enhanced in irradiated interstellar medium, CN, C$_2$H, and HNC, show similar concentrations and extended components along the bicone, suggesting that molecule dissociation is a dominant chemical effect of the cold molecular outflow of this galaxy. Although further investigation should be made, this scenario is consistent with the faintness or absence of the emission lines from CO isotopologues, CH$_3$OH, and N$_2$H$^+$, in the outflow, which are easily destroyed by dissociating photons and electrons.

preprint2022arXiv

Detection of extended millimeter emission in the host galaxy of 3C273 and its implications for QSO feedback via high dynamic range ALMA imaging

We estimate the amount of negative feedback energy injected into the ISM of the host galaxy of 3C273, a prototypical radio loud quasar. We obtained 93, 233 and 343 GHz continuum images with the Atacama Large Millimeter/Sub-millimeter Array (ALMA). After self calibration and point source subtraction, we reach an image dynamic range of $\sim 85000$ at 93\ GHz, $\sim 39000$ at 233\ GHz and $\sim 2500$ at 343\ GHz. These are currently the highest image dynamic range obtained using ALMA. We detect spatially extended millimeter emission associated with the host galaxy, cospatial with the Extended Emission Line Region (EELR) observed in the optical. The millimeter spectral energy distribution and comparison with centimeter data show that the extended emission cannot be explained by dust thermal emission, synchrotron or thermal bremsstrahlung arising from massive star formation. We interpret the extended millimeter emission as thermal bremsstrahlung from gas directly ionized by the central source. The extended flux indicates that at least $\sim 7\%$ of the bolometric flux of the nuclear source was used to ionize atomic hydrogen in the host galaxy. The ionized gas is estimated to be as massive as $10^{10}$ to $10^{11}\ \mathrm{M_\odot}$, but the molecular gas fraction with respect to the stellar mass is consistent with other ellipticals, suggesting that direct ionization ISM by the QSO may not be sufficient to suppress star formation, or we are witnessing a short timescale before negative feedback becomes observable. The discovery of a radio counterpart to EELRs provides a new pathway to studying the QSO-host ISM interaction.

preprint2022arXiv

Detection of nitrogen and oxygen in a galaxy at the end of reionization

We present observations of [NII] 205 $μ$m, [OIII] 88 $μ$m and dust emission in a strongly-lensed, submillimeter galaxy (SMG) at $z=6.0$, G09.83808, with the Atacama Large Millimeter/submillimeter Array (ALMA). Both [NII] and [OIII] line emissions are detected at $>12σ$ in the 0.8$"$-resolution maps. Lens modeling indicates that the spatial distribution of the dust continuum emission is well characterized by a compact disk with an effective radius of 0.64$\pm$0.02 kpc and a high infrared surface brightness of $Σ_\mathrm{IR}=(1.8\pm0.3)\times10^{12}~L_\odot$ kpc$^{-2}$. This result supports that G09.83808 is the progenitors of compact quiescent galaxies at $z\sim4$, where the majority of its stars are expected to be formed through a strong and short burst of star formation. G09.83808 and other lensed SMGs show a decreasing trend of the [NII] line to infrared luminosity ratio with increasing continuum flux density ratio between 63 $μ$m and 158 $μ$m, as seen in local luminous infrared galaxies (LIRGs). The decreasing trend can be reproduced by photoionization models with increasing ionization parameters. Furthermore, by combining the [NII]/[OIII] luminosity ratio with far-infrared continuum flux density ratio in G09.83808, we infer that the gas phase metallicity is already $Z\approx 0.5-0.7~Z_\odot$. G09.83808 is likely one of the earliest galaxies that has been chemically enriched at the end of reionization.

preprint2022arXiv

Environmental dependence of the molecular cloud lifecycle in 54 main sequence galaxies

The processes of star formation and feedback, regulating the cycle of matter between gas and stars on the scales of giant molecular clouds (GMCs; $\sim$100pc), play a major role in governing galaxy evolution. Measuring the time-scales of GMC evolution is important to identify and characterise the specific physical mechanisms that drive this transition. By applying a robust statistical method to high-resolution CO and narrow-band H$α$ imaging from the PHANGS survey, we systematically measure the evolutionary timeline from molecular clouds to exposed young stellar regions on GMC scales, across the discs of an unprecedented sample of 54 star-forming main-sequence galaxies (excluding their unresolved centres). We find that clouds live for about $1{-}3$ GMC turbulence crossing times ($5{-}30$Myr) and are efficiently dispersed by stellar feedback within $1{-}5$Myr once the star-forming region becomes partially exposed, resulting in integrated star formation efficiencies of $1{-}8$%. These ranges reflect physical galaxy-to-galaxy variation. In order to evaluate whether galactic environment influences GMC evolution, we correlate our measurements with average properties of the GMCs and their local galactic environment. We find several strong correlations that can be physically understood, revealing a quantitative link between galactic-scale environmental properties and the small-scale GMC evolution. Notably, the measured CO-visible cloud lifetimes become shorter with decreasing galaxy mass, mostly due to the increasing presence of CO-dark molecular gas in such environment. Our results represent a first step towards a comprehensive picture of cloud assembly and dispersal, which requires further extension and refinement with tracers of the atomic gas, dust, and deeply-embedded stars.

preprint2022arXiv

Multi-wavelength and Multi-CO View of The Minor Merger Driven Star Formation in the Nearby LIRG NGC 3110

We present Atacama Large Millimeter/submillimeter Array observations of multiple CO(1-0), $^{13}$CO(1-0), and C$^{18}$O(1-0) lines and 2.9 mm and 1.3 mm continuum emission toward the nearby interacting luminous infrared galaxy NGC 3110, supplemented with similar spatial resolution H$α$, 1.4GHz continuum, and $K$-band data. We estimate the typical CO-to-H$_2$ conversion factor of 1.7 $M_{\odot}$ (K km s$^{-1}$ pc$^2$)$^{-1}$ within the disk using LTE-based and dust-based H$_2$ column densities, and measure the 1-kpc scale surface densities of star formation rate ($Σ_{\rm SFR}$), super star clusters ($Σ_{\rm SSC}$), molecular gas mass, and star formation efficiency (SFE) toward the entire gas disk. These parameters show a peak at the southern part of the southern spiral arm (SFE $\sim$ 10$^{-8.2}$ yr$^{-1}$, $Σ_{\rm SFR}$ $\sim$ 10$^{-0.6}$ $M_{\odot}$ kpc$^{-2}$ yr$^{-1}$, $Σ_{\rm SSC}$ $\sim$ 6.0 kpc$^{-2}$), which is likely attributed to the on-going tidal interaction with the companion galaxy MCG-01-26-013, as well as toward the circumnuclear region. We also find that thermal free-free emission contributes to a significant fraction of the millimeter continuum emission at the southern peak position. Those measurements imply that the peak of the southern arm is an active and young star-forming region, whereas the central part of NGC 3110 is a site of long-continued star formation. We suggest that, during the early stage of the galaxy-galaxy interaction with large mass ratio that in NGC 3110, fragmentation along the main galaxy's arms is an important driver of merger-induced star formation and massive gas inflow results in dusty nuclear starbursts.

preprint2022arXiv

Spatially-resolved relation between [CI] $^{3}P_{1}$-$^{3}P_{0}$ and $^{12}$CO (1-0) in Arp 220

We present $\sim$0."3 (114 pc) resolution maps of [CI] $^{3}P_{1}$-$^{3}P_{0}$ (hereafter [CI] (1-0)) and $^{12}$CO (1-0) obtained toward Arp 220 with the Atacama Large Millimeter/submillimeter Array. The overall distribution of the [CI] (1-0) emission is consistent with the CO (1-0). While the [CI] (1-0) and CO (1-0) luminosities of the system follow the empirical linear relation for the unresolved ULIRG sample, we find a sublinear relation between [CI] (1-0) and CO (1-0) using the spatially-resolved data. We measure the [CI] (1-0)/CO (1-0) luminosity ratio per pixel in star-forming environments of Arp 220 and investigate its dependence on the CO (3-2)/CO (1-0) ratio ($R_{\rm CO}$). On average, the [CI] (1-0)/CO (1-0) luminosity ratio is almost constant up to $R_{\rm CO} \simeq 1$ and then increases with $R_{\rm CO}$. According to the radiative transfer analysis, a high CI/CO abundance ratio is required in regions with high [CI] (1-0)/CO (1-0) luminosity ratios and $R_{\rm CO} > 1$, suggesting that the CI/CO abundance ratio varies at $\sim$100 pc scale in Arp 220. The [CI] (1-0)/CO (1-0) luminosity ratio depends on multiple factors and may not be straightforward to interpret. We also find the high-velocity components traced by [CI] (1-0) in the western nucleus, likely associated with the molecular outflow. The [CI] (1-0)/CO (1-0) luminosity ratio in the putative outflow is 0.87 $\pm$ 0.28, which is four times higher than the average ratio of Arp 220. While there is a possibility that the [CI] (1-0) and CO (1-0) emission traces different components, we suggest that the high line ratios are likely because of elevated CI/CO abundance ratios based on our radiative transfer analysis. A CI-rich and CO-poor gas phase in outflows could be caused by the irradiation of the cosmic rays, the shock heating, and the intense radiation field.

preprint2022arXiv

The Gas-Star Formation Cycle in Nearby Star-forming Galaxies II. Resolved Distributions of CO and H$α$ Emission for 49 PHANGS Galaxies

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

preprint2022arXiv

The Kiloparsec-scale Neutral Atomic Carbon Outflow in the Nearby Type-2 Seyfert Galaxy NGC 1068: Evidence for Negative AGN Feedback

Active galactic nucleus (AGN) feedback is postulated as a key mechanism for regulating star formation within galaxies. Studying the physical properties of the outflowing gas from AGN is thus crucial for understanding the co-evolution of galaxies and supermassive black holes. Here we report 55 pc resolution ALMA neutral atomic carbon [CI] $^3P_1\text{-}^3P_0$ observations toward the central 1 kpc of the nearby type-2 Seyfert galaxy NGC 1068, supplemented by 55 pc resolution CO($J=1\text{-}0$) observations. We find that [CI] emission within the central kpc is strongly enhanced by a factor of $>$5 compared to the typical [CI]/CO intensity ratio of $\sim$0.2 for nearby starburst galaxies (in units of brightness temperature). The most [CI]-enhanced gas (ratio $>$ 1) exhibits a kpc-scale elongated structure centered at the AGN that matches the known biconical ionized gas outflow entraining molecular gas in the disk. A truncated, decelerating bicone model explains well the kinematics of the elongated structure, indicating that the [CI] enhancement is predominantly driven by the interaction between the ISM in the disk and the highly inclined ionized gas outflow (which is likely driven by the radio jet). Our results strongly favor the "CO dissociation scenario" rather than the "in-situ C formation" one which prefers a perfect bicone geometry. We suggest that the high [CI]/CO intensity ratio gas in NGC 1068 directly traces ISM in the disk that is currently dissociated and entrained by the jet and the outflow, i.e., the "negative" effect of the AGN feedback.

preprint2022arXiv

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

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

preprint2022arXiv

WISDOM Project -- X. The morphology of the molecular ISM in galaxy centres and its dependence on galaxy structure

We use high-resolution maps of the molecular interstellar medium (ISM) in the centres of eighty-six nearby galaxies from the millimetre-Wave Interferometric Survey of Dark Object Masses (WISDOM) and Physics at High Angular Resolution in Nearby GalaxieS (PHANGS) surveys to investigate the physical mechanisms setting the morphology of the ISM at molecular cloud scales. We show that early-type galaxies tend to have smooth, regular molecular gas morphologies, while the ISM in spiral galaxy bulges is much more asymmetric and clumpy when observed at the same spatial scales. We quantify these differences using non-parametric morphology measures (Asymmetry, Smoothness and Gini), and compare these measurements with those extracted from idealised galaxy simulations. We show that the morphology of the molecular ISM changes systematically as a function of various large-scale galaxy parameters, including galaxy morphological type, stellar mass, stellar velocity dispersion, effective stellar mass surface density, molecular gas surface density, star formation efficiency and the presence of a bar. We perform a statistical analysis to determine which of these correlated parameters best predicts the morphology of the ISM. We find the effective stellar mass surface (or volume) density to be the strongest predictor of the morphology of the molecular gas, while star formation and bars maybe be important secondary drivers. We find that gas self-gravity is not the dominant process shaping the morphology of the molecular gas in galaxy centres. Instead effects caused by the depth of the potential well such as shear, suppression of stellar spiral density waves and/or inflow affect the ability of the gas to fragment.

preprint2021arXiv

CO excitation, molecular gas density and interstellar radiation field in local and high-redshift galaxies

We study the Carbon Monoxide (CO) excitation, mean molecular gas density and interstellar radiation field (ISRF) intensity in a comprehensive sample of 76 galaxies from local to high redshift (z~0-6), selected based on detections of their CO transitions J=2-1 and 5-4 and their optical/infrared/(sub-)millimeter spectral energy distributions (SEDs). We confirm the existence of a tight correlation between CO excitation as traced by the CO(5-4)/(2-1) line ratio (R52), and the mean ISRF intensity U as derived from infrared SED fitting using dust SED templates. By modeling the molecular gas density probability distribution function (PDF) in galaxies and predicting CO line ratios with large velocity gradient radiative transfer calculations, we present a framework linking global CO line ratios to the mean molecular hydrogen gas density nH2 and kinetic temperature Tkin. Mapping in this way observed R52 ratios to nH2 and Tkin probability distributions, we obtain positive U-nH2 and U-Tkin correlations, which imply a scenario in which the ISRF in galaxies is mainly regulated by Tkin and (non-linearly) by nH2. A small fraction of starburst galaxies showing enhanced nH2 could be due to merger-driven compaction. Our work demonstrates that ISRF and CO excitation are tightly coupled, and that density-PDF modeling is a promising tool for probing detailed ISM properties inside galaxies.

preprint2021arXiv

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

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

preprint2021arXiv

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

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&#34; \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.

preprint2021arXiv

PHANGS-MUSE: the HII region luminosity function of local star-forming galaxies

We use an unprecedented sample of about 23,000 HII regions detected at an average physical resolution of 67pc in the PHANGS-MUSE sample to study the extragalactic HII region Ha luminosity function (LF). Our observations probe the star-forming disk of 19 nearby spiral galaxies with low inclination and located close to the star formation main sequence at z=0. The mean LF slope $α$ in our sample is =1.73 with a $σ$ of 0.15. We find that $α$ decreases with the galaxy&#39;s star formation rate surface density and argue that this is driven by an enhanced clustering of young stars at high gas surface densities. Looking at the HII regions within single galaxies we find that no significant variations occur between the LF of the inner and outer part of the star-forming disk, whereas the LF in the spiral arm areas is shallower than in the inter-arm areas for six out of the 13 galaxies with clearly visible spiral arms. We attribute these variations to the spiral arms increasing the molecular clouds arm--inter-arm mass contrast and find suggestive evidence that they are more evident for galaxies with stronger spiral arms. Furthermore, we find systematic variations in $α$ between samples of HII regions with high and low ionization parameter q and argue that they are driven by the aging of HII regions.

preprint2021arXiv

The PHANGS-HST Survey: Physics at High Angular resolution in Nearby GalaxieS with the Hubble Space Telescope

The PHANGS program is building the first dataset to enable the multi-phase, multi-scale study of star formation across the nearby spiral galaxy population. This effort is enabled by large survey programs with ALMA, VLT/MUSE, and HST, with which we have obtained CO(2-1) imaging, optical spectroscopic mapping, and high resolution UV-optical imaging, respectively. Here, we present PHANGS-HST, which has obtained five band NUV-U-B-V-I imaging of the disks of 38 spiral galaxies at distances of 4-23 Mpc, and parallel V and I band imaging of their halos, to provide a census of tens of thousands of compact star clusters, and multi-scale stellar associations. The combination of HST, ALMA, and VLT/MUSE observations will yield an unprecedented joint catalog of the observed and physical properties of ~100,000 star clusters, associations, HII regions, and molecular clouds. With these basic units of star formation, PHANGS will systematically chart the evolutionary cycling between gas and stars, across a diversity of galactic environments found in nearby galaxies. We discuss the design of the PHANGS-HST survey, and provide an overview of the HST data processing pipeline and first results, highlighting new methods for selecting star cluster candidates, morphological classification of candidates with convolutional neural networks, and identification of stellar associations over a range of physical scales with a watershed algorithm. We describe the cross-observatory imaging, catalogs, and software products to be released. These high-level science products will seed a broad range of investigations, in particular, the study of embedded stellar populations and dust with JWST, for which a PHANGS Cycle 1 Treasury program to obtain eight band 2-21 $μ$m imaging has been approved.

preprint2020arXiv

A non-corotating gas component in an extreme starburst at z=4.3

We report the detection of a non-corotating gas component in a bright unlensed submillimeter galaxy at z=4.3, COSMOS-AzTEC-1, hosting a compact starburst. ALMA 0.17 and 0.09 arcsec resolution observations of [CII] emission clearly demonstrate that the gas kinematics is characterized by an ordered rotation. After subtracting the best-fit model of a rotating disk, we kinematically identify two residual components in the channel maps. Both observing simulations and analysis of dirty images confirm that these two subcomponents are not artificially created by noise fluctuations and beam deconvolution. One of the two has a velocity offset of 200 km/s and a physical separation of 2 kpc from the primary disk and is located along the kinematic minor axis of disk rotation. We conclude that this gas component is falling into the galaxy from a direction perpendicular to the disk rotation. The accretion of such small non-corotating gas components could stimulate violent disk instability, driving radial gas inflows into the center of galaxies and leading to formation of in-situ clumps such as identified in dust continuum and CO. We require more theoretical studies on high gas fraction mergers with mass ratio of 1:>10 to verify this process.

preprint2020arXiv

Discovery of a [CI]-faint, CO-bright Galaxy: ALMA Observations of the Merging Galaxy NGC 6052

We report sensitive [\ion{C}{1}]~$^3P_1$--$^3P_0$ and $^{12}$CO~$J$=4--3 observations of the nearby merging galaxy NGC 6052 using the Morita (Atacama Compact) Array of ALMA. We detect $^{12}$CO~$J$=4--3 toward the northern part of NGC 6052, but [\ion{C}{1}]~$^3P_1$--$^3P_0$ is not detected with a [\ion{C}{1}]~$^3P_1$--$^3P_0$ to $^{12}$CO~$J$=4--3 line luminosity ratio of$~\lesssim0.07$. According to models of photodissociation regions, the unusual weakness of [\ion{C}{1}]~$^3P_1$--$^3P_0$ relative to $^{12}$CO~$J$=4--3 can be explained if the interstellar medium has a hydrogen density larger than $10^5\,{\rm cm}^{-3}$, conditions that might arise naturally in the ongoing merging process in NGC 6052. Its [\ion{C}{1}]~$^3P_1$--$^3P_0$ emission is also weaker than expected given the molecular gas mass inferred from previous measurements of $^{12}$CO~$J$=1--0 and $^{12}$CO~$J$=2--1. This suggests that [\ion{C}{1}]~$^3P_1$--$^3P_0$ may not be a reliable tracer of molecular gas mass in this galaxy. NGC 6052 is a unique laboratory to investigate how the merger process impacts the molecular gas distribution.

preprint2020arXiv

Dynamical Equilibrium in the Molecular ISM in 28 Nearby Star-Forming Galaxies

We compare the observed turbulent pressure in molecular gas, $P_\mathrm{turb}$, to the required pressure for the interstellar gas to stay in equilibrium in the gravitational potential of a galaxy, $P_\mathrm{DE}$. To do this, we combine arcsecond resolution CO data from PHANGS-ALMA with multi-wavelength data that traces the atomic gas, stellar structure, and star formation rate (SFR) for 28 nearby star-forming galaxies. We find that $P_\mathrm{turb}$ correlates with, but almost always exceeds the estimated $P_\mathrm{DE}$ on kiloparsec scales. This indicates that the molecular gas is over-pressurized relative to the large-scale environment. We show that this over-pressurization can be explained by the clumpy nature of molecular gas; a revised estimate of $P_\mathrm{DE}$ on cloud scales, which accounts for molecular gas self-gravity, external gravity, and ambient pressure, agrees well with the observed $P_\mathrm{turb}$ in galaxy disks. We also find that molecular gas with cloud-scale ${P_\mathrm{turb}}\approx{P_\mathrm{DE}}\gtrsim{10^5\,k_\mathrm{B}\,\mathrm{K\,cm^{-3}}}$ in our sample is more likely to be self-gravitating, whereas gas at lower pressure appears more influenced by ambient pressure and/or external gravity. Furthermore, we show that the ratio between $P_\mathrm{turb}$ and the observed SFR surface density, $Σ_\mathrm{SFR}$, is compatible with stellar feedback-driven momentum injection in most cases, while a subset of the regions may show evidence of turbulence driven by additional sources. The correlation between $Σ_\mathrm{SFR}$ and kpc-scale $P_\mathrm{DE}$ in galaxy disks is consistent with the expectation from self-regulated star formation models. Finally, we confirm the empirical correlation between molecular-to-atomic gas ratio and kpc-scale $P_\mathrm{DE}$ reported in previous works.

preprint2020arXiv

Is this an Early Stage Merger? A Case Study on Molecular Gas and Star Formation Properties of Arp 240

We present new high resolution $^{12}$CO $J$=1-0, $J$=2-1, and $^{13}$CO $J$=1-0 maps of the early stage merger Arp 240 (NGC5257/8) obtained with the Atacama Large Millimeter/submillimeter Array (ALMA). Simulations in the literature suggest that the merger has just completed its first passage; however, we find that this system has a lower global gas fraction but a higher star formation efficiency compared to typical close galaxy pairs, which suggests that this system may already be in an advanced merger stage. We combine the ALMA data with $^{12}$CO $J$=3-2 observations from the Submillimeter Array and carry out RADEX modeling on several different regions. Both the RADEX modeling and a local thermal equilibrium (LTE) analysis show that the regions are most likely to have a CO-to-H$_2$ conversion factor $α_{\mathrm{CO}}$ close to or perhaps even smaller than the typical value for (ultra-)luminous infrared galaxies. Using 33 GHz data from the Very Large Array to measure the star formation rate, we find that most star forming regions have molecular gas depletion times of less than 100 Myr. We calculated the star formation efficiency (SFE) per free-fall time for different regions and find some regions appear to have values greater than 100%. We find these regions generally show evidence for young massive clusters (YMCs). After exploring various factors, we argue that this is mainly due to the fact that radio continuum emission in those regions is dominated by that from YMCs, which results in an overestimate of the SFE per free-fall time.

preprint2020arXiv

PHANGS CO kinematics: disk orientations and rotation curves at 150 pc resolution

We present kinematic orientations and high resolution (150 pc) rotation curves for 67 main sequence star-forming galaxies surveyed in CO (2-1) emission by PHANGS-ALMA. Our measurements are based on the application of a new fitting method tailored to CO velocity fields. Our approach identifies an optimal global orientation as a way to reduce the impact of non-axisymmetric (bar and spiral) features and the uneven spatial sampling characteristic of CO emission in the inner regions of nearby galaxies. The method performs especially well when applied to the large number of independent lines-of-sight contained in the PHANGS CO velocity fields mapped at 1&#39;&#39; resolution. The high resolution rotation curves fitted to these data are sensitive probes of mass distribution in the inner regions of these galaxies. We use the inner slope as well as the amplitude of our fitted rotation curves to demonstrate that CO is a reliable global dynamical mass tracer. From the consistency between photometric orientations from the literature and kinematic orientations determined with our method, we infer that the shapes of stellar disks in the mass range of log($\rm M_{\star}(M_{\odot})$)=9.0-10.9 probed by our sample are very close to circular and have uniform thickness.

preprint2020arXiv

Star formation traced by optical and millimeter hydrogen recombination lines and free-free emissions in the dusty merging galaxy NGC 3256 -- MUSE/VLT and ALMA synergy

A galaxy-galaxy merger and the subsequent triggering of starburst activity are fundamental processes linked to the morphological transformation of galaxies and the evolution of star formation across the history of the Universe. Both nuclear and disk-wide starbursts are assumed to occur during the merger process. However, quantifying both nuclear and disk-wide star formation activity is non-trivial because the nuclear starburst is dusty in the most active merging starburst galaxies. This paper presents a new approach to this problem: combining hydrogen recombination lines in optical, millimeter, and free-free emission. Using NGC~3256 as a case study, H$β$, H40$α$, and free-free emissions are investigated using the Multi Unit Spectroscopic Explorer at the Very Large Telescope of the European Southern Observatory (MUSE/VLT) and the Atacama Large Millimeter/submillimeter Array (ALMA). The H$β$ image obtained by MUSE identifies star-forming regions outside the nuclear regions, suggesting a disk-wide starburst. In contrast, the H40$α$ image obtained by ALMA identifies a nuclear starburst where optical lines are undetected due to dust extinction ($A_{\rm V}\sim25$). Combining both MUSE and ALMA observations, we conclude that the total SFR is $49\pm2~M_{\odot}$~yr$^{-1}$ and the contributions from nuclear and disk-wide starbursts are $\sim34~\%$ and $\sim66~\%$, respectively. This suggests the dominance of disk-wide star formation in NGC~3256. In addition, pixel-by-pixel analyses for disk-wide star-forming regions suggest that shock gas tracers (e.g., CH$_3$OH) are enhanced where gas depletion time ($τ_{\rm gas}$=$M_{\rm gas}/SFR$) is long. This possibly means that merger-induced shocks regulate disk-wide star formation activities.

preprint2020arXiv

The 300-pc Scale ALMA View of [CI] $^3P_1$-$^3P_0$, CO $J$=1-0, and 609 $μ$m Dust Continuum in A Luminous Infrared Galaxy

We present high-quality ALMA Band 8 observations of the [CI] $^3P_1$-$^3P_0$ line and 609 $μ$m dust continuum emission toward the nearby luminous infrared galaxy (LIRG) IRAS F18293-3413, as well as matched resolution (300-pc scale) Band 3 CO $J=$1-0 data, which allow us to assess the use of the [CI] $^3P_1$-$^3P_0$ line as a total gas mass estimator. We find that the [CI] line basically traces structures detected in CO (and dust), and a mean (median) [CI]/CO luminosity ($L&#39;_{\rm [CI]}$/$L&#39;_{\rm CO}$) ratio of 0.17 (0.16) with a scatter of 0.04. However, a pixel-by-pixel comparison revealed that there is a radial $L&#39;_{\rm [CI]}$/$L&#39;_{\rm CO}$ gradient and a superlinear $L&#39;_{\rm CO}$ vs. $L&#39;_{\rm [CI]}$ relation (slope = 1.54 $\pm$ 0.02) at this spatial scale, which can be explained by radial excitation and/or line opacity gradients. Based on the molecular gas masses converted from the dust continuum emission, we found that the CO-to-H$_2$ and [CI]-to-H$_2$ conversion factors are relatively flat across the molecular gas disk with a median value of 3.5$^{+1.9}_{-1.3}$ and 20.7$^{+9.2}_{-4.9}$ $M_{\odot}$ (K km s$^{-1}$ pc$^2$)$^{-1}$, respectively. A non-LTE calculation yields that typical molecular gas properties seen in nearby (U)LIRGs ($n_{\rm H_2}$ = 10$^{3-4}$ cm$^{-3}$, $T_{\rm kin}$ $\sim$ 50 K, and $X_{\rm CI}$ = (0.8-2.3) $\times$ 10$^{-5}$) can naturally reproduce the derived [CI]-to-H$_2$ conversion factor. However, we caution that a careful treatment of the physical gas properties is required in order to measure H$_2$ gas mass distributions in galaxies using a single [CI] line. Otherwise, a single [CI] line is not a good molecular gas estimator in a spatially resolved manner.

preprint2019arXiv

SOFIA/HAWC+ View of an Extremely Luminous Infrared Galaxy, WISE1013+6112

We present far-infrared (FIR) properties of an extremely luminous infrared galaxy (ELIRG) at $z_{\rm spec}$ = 3.703, WISE J101326.25+611220.1 (WISE1013+6112). This ELIRG is selected as an IR-bright dust-obscured galaxy (DOG) based on the photometry from the Sloan digital sky survey (SDSS) and wide-field infrared survey explorer (WISE). In order to derive its accurate IR luminosity, we perform follow-up observations at 89 and 154 $μ$m using the high-resolution airborne wideband camera-plus (HAWC+) on board the 2.7-m stratospheric observatory for infrared astronomy (SOFIA) telescope. We conduct spectral energy distribution (SED) fitting with CIGALE using 15 photometric data (0.4-1300 $μ$m). We successfully pin down FIR SED of WISE1013+6112 and its IR luminosity is estimated to be $L_{\rm IR}$ = (1.62 $\pm$ 0.08) $\times 10^{14}$ $L_{\odot}$, making it one of the most luminous IR galaxies in the universe. We determine the dust temperature of WISE1013+6112 is $T_{\rm dust}$ = 89 $\pm$ 3 K, which is significantly higher than that of other populations such as SMGs and FIR-selected galaxies at similar IR luminosities. The resultant dust mass is $M_{\rm dust}$ = (2.2 $\pm$ 0.1) $\times 10^{8}$ $M_{\odot}$. This indicates that WISE1013+6112 has a significant active galactic nucleus (AGN) and star-forming activity behind a large amount of dust.

preprint2019arXiv

The Gas-Star Formation Cycle in Nearby Star-Forming Galaxies I. Assessment of Multi-scale Variations

The processes regulating star formation in galaxies are thought to act across a hierarchy of spatial scales. To connect extragalactic star formation relations from global and kpc-scale measurements to recent cloud-scale resolution studies, we have developed a simple, robust method that quantifies the scale dependence of the relative spatial distributions of molecular gas and recent star formation. In this paper, we apply this method to eight galaxies with roughly 1 arcsec resolution molecular gas imaging from the PHANGS-ALMA and PAWS surveys that have matched resolution, high quality narrowband Halpha imaging. At a common scale of 140pc, our massive (log(Mstar/Msun)=9.3-10.7), normally star-forming (SFR/Msun/yr=0.3-5.9) galaxies exhibit a significant reservoir of quiescent molecular gas not associated with star formation as traced by Halpha emission. Galactic structures act as backbones for both molecular and HII region distributions. As we degrade the spatial resolution, the quiescent molecular gas disappears, with the most rapid changes occurring for resolutions up to about 0.5kpc. As the resolution becomes poorer, the morphological features become indistinct for spatial scales larger than about 1kpc. The method is a promising tool to search for relationships between the quiescent or star-forming molecular reservoir and galaxy properties, but requires a larger sample size to identify robust correlations between the star-forming molecular gas fraction and global galaxy parameters.

preprint2019arXiv

The headlight cloud in NGC 628: An extreme giant molecular cloud in a typical galaxy disk

Cloud-scale surveys of molecular gas reveal the link between molecular clouds properties and star formation (SF) across a range of galactic environments. Cloud populations in galaxy disks are considered to be representative of the `normal&#39; SF. At high resolution, however, clouds with exceptional gas properties and SF activity may also be observed in normal disk environments. In this paper, we study the brightest cloud traced in CO emission in the disk of NGC628. The cloud is spatially coincident with an extremely bright HII region. We characterize its molecular gas properties and investigate how feedback and large-scale processes influence the properties of the molecular gas. High resolution CO ALMA observations are used to characterize its mass and dynamical state, which are compared to other clouds in NGC628. A LVG analysis is used to constrain the beam-diluted density and temperature of the molecular gas. We analyze the MUSE spectrum using Starburst99 to characterize the young stellar population associated with the HII region. The cloud is massive ($1-2\times10^7$M$_{\odot}$), with a beam-diluted density of $n_{\rm H_2}=5\times10^4$ cm$^{-3}$. It has a low virial parameter, suggesting that its CO emission may be overluminous due to heating by the HII region. A young ($2-4$ Myr), massive $3\times10^{5}$ M$_{\odot}$ stellar population is associated. We argue that the cloud is currently being destroyed by feedback from young massive stars. Due to the cloud&#39;s large mass, this phase of the cloud&#39;s evolution is long enough for the impact of feedback on the excitation of the gas to be observed. Its high mass may be related to its location at a spiral co-rotation radius, where gas experiences reduced galactic shear compared to other regions of the disk, and receives a sustained inflow of gas that can promote the cloud&#39;s mass growth.

preprint2019arXiv

The lifecycle of molecular clouds in nearby star-forming disc galaxies

It remains a major challenge to derive a theory of cloud-scale ($\lesssim100$ pc) star formation and feedback, describing how galaxies convert gas into stars as a function of the galactic environment. Progress has been hampered by a lack of robust empirical constraints on the giant molecular cloud (GMC) lifecycle. We address this problem by systematically applying a new statistical method for measuring the evolutionary timeline of the GMC lifecycle, star formation, and feedback to a sample of nine nearby disc galaxies, observed as part of the PHANGS-ALMA survey. We measure the spatially-resolved ($\sim100$ pc) CO-to-H$α$ flux ratio and find a universal de-correlation between molecular gas and young stars on GMC scales, allowing us to quantify the underlying evolutionary timeline. GMC lifetimes are short, typically 10-30 Myr, and exhibit environmental variation, between and within galaxies. At kpc-scale molecular gas surface densities $Σ_{\rm H_2}\geqslant8$M$_{\odot}$pc$^{-2}$, the GMC lifetime correlates with time-scales for galactic dynamical processes, whereas at $Σ_{\rm H_2}\leqslant8$M$_{\odot}$pc$^{-2}$ GMCs decouple from galactic dynamics and live for an internal dynamical time-scale. After a long inert phase without massive star formation traced by H$α$ (75-90% of the cloud lifetime), GMCs disperse within just 1-5 Myr once massive stars emerge. The dispersal is most likely due to early stellar feedback, causing GMCs to achieve integrated star formation efficiencies of 4-10% These results show that galactic star formation is governed by cloud-scale, environmentally-dependent, dynamical processes driving rapid evolutionary cycling. GMCs and HII regions are the fundamental units undergoing these lifecycles, with mean separations of 100-300 pc in star-forming discs. Future work should characterise the multi-scale physics and mass flows driving these lifecycles.