Search

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., Comment: 29 pages, 22 figures, accepted for publication in A&A

We present PHANGS-ALMA, the first survey to map CO J=2-1 line emission at ~1" ~ 100pc spatial resolution from a representative sample of 90 nearby (d<~20 Mpc) galaxies that lie on or near the z=0 "main sequence" of star-forming galaxies. CO line emission traces the bulk distribution of molecular gas, which is the cold, star-forming phase of the interstellar medium. At the resolution achieved by PHANGS-ALMA, each beam reaches the size of a typical individual giant molecular cloud (GMC), so that these data can be used to measure the demographics, life-cycle, and physical state of molecular clouds across the population of galaxies where the majority of stars form at z=0. This paper describes the scientific motivation and background for the survey, sample selection, global properties of the targets, ALMA observations, and characteristics of the delivered ALMA data and derived data products. As the ALMA sample serves as the parent sample for parallel surveys with VLT/MUSE, HST, AstroSat, VLA, and other facilities, we include a detailed discussion of the sample selection. We detail the estimation of galaxy mass, size, star formation rate, CO luminosity, and other properties, compare estimates using different systems and provide best-estimate integrated measurements for each target. We also report the design and execution of the ALMA observations, which combine a Cycle~5 Large Program, a series of smaller programs, and archival observations. Finally, we present the first 1" resolution atlas of CO emission from nearby galaxies and describe the properties and contents of the first PHANGS-ALMA public data release., Comment: 76 pages, 33 figures. Accepted for publication in the Astrophysical Journal Supplement series. Full resolution version and the image atlas to appear as a figure set in the published version can be found https://sites.google.com/view/phangs/publications . Data release coming soon to the ALMA archive and CADC temporarily available at http://phangs.org/data

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

We present the highest resolution CO (2-1) observations obtained to date (0.25") of NGC 3256 and use them to determine the detailed properties of the molecular interstellar medium in the central 6 kpc of this merger. Distributions of physical quantities are reported from pixel-by-pixel measurements at 55 and 120 pc scales and compared to disc galaxies observed by PHANGS-ALMA. Mass surface densities range from 8 to 5500 M$_{\odot}$ pc$^{-2}$ and velocity dispersions from 10 to 200 km s$^{-1}$. Peak brightness temperatures as large as 37 K are measured, indicating the gas in NGC 3256 may be hotter than all regions in nearby disc galaxies measured by PHANGS-ALMA. Brightness temperatures even surpass those in the overlap region of NGC 4038/9 at the same scales. The majority of the gas appears unbound with median virial parameters of 7 to 19, although external pressure may bind some of the gas. High internal turbulent pressures of 10$^{5}$ to 10$^{10}$ K cm$^{-3}$ are found. Given the lack of significant trends in surface density, brightness temperature, and velocity dispersion with physical scale we argue the molecular gas is made up of a smooth medium down to 55 pc scales, unlike the more structured medium found in the PHANGS-ALMA disc galaxies., Comment: accepted to MNRAS, 21 pages, 8 figures

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 $\alpha_{\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., Comment: 20 pages, 14 figures, accepted for publication in MNRAS

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' 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's large mass, this phase of the cloud'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's mass growth., Comment: A&A accepted

We report ~2" resolution Atacama Large Millimeter/submillimeter Array observations of the HCN(1-0), HCO+(1-0), CO(1-0), CO(2-1), and CO(3-2) lines towards the nearby merging double-nucleus galaxy NGC 3256. We find that the high density gas outflow traced in HCN(1-0) and HCO+(1-0) emission is co-located with the diffuse molecular outflow emanating from the southern nucleus, where a low-luminosity active galactic nucleus (AGN) is believed to be the dominant source of the far-infrared luminosity. On the other hand, the same lines were undetected in the outflow region associated with the northern nucleus, whose primary heating source is likely related to starburst activity without obvious signs of AGN. Both HCO+(1-0)/CO(1-0) line ratio (i.e. dense gas fraction) and the CO(3-2)/CO(1-0) line ratio are larger in the southern outflow (0.20$\pm$0.04 and 1.3$\pm$0.2, respectively) than in the southern nucleus (0.08$\pm$0.01, 0.7$\pm$0.1, respectively). By investigating these line ratios for each velocity component in the southern outflow, we find that the dense gas fraction increases and the CO(3-2)/CO(1-0) line ratio decreases towards the largest velocity offset. This suggests the existence of a two-phase (diffuse and clumpy) outflow. One possible scenario to produce such a two-phase outflow is an interaction between the jet and the interstellar medium, which possibly triggers shocks and/or star formation associated with the outflow., Comment: 22 pages, 12 figures, 5 tables, accepted for publication in APJ

In external galaxies, molecular composition may be influenced by extreme environments such as starbursts and galaxy mergers. To study such molecular chemistry, we observed the luminous-infrared galaxy and merger NGC 3256 using the Atacama Large Millimeter/sub-millimeter Array. We covered most of the 3-mm and 1.3-mm bands for a multi-species, multi-transition analysis. We first analyzed intensity ratio maps of selected lines such as HCN/HCO$^+$, which shows no enhancement at an AGN. We then compared the chemical compositions within NGC 3256 at the two nuclei, tidal arms, and positions with influence from galactic outflows. We found the largest variation in SiO and CH$_3$OH, species that are likely to be enhanced by shocks. Next, we compared the chemical compositions in the nuclei of NGC 3256, NGC 253, and Arp 220; these galactic nuclei have varying star formation efficiencies. Arp 220 shows higher abundances of SiO and HC$_3$N than NGC 3256 and NGC 253. Abundances of most species do not show strong correlation with the star formation efficiencies, although the CH$_3$CCH abundance seems to have a weak positive correlation with the star formation efficiency. Lastly, the chemistry of spiral arm positions in NGC 3256 is compared with that of W 51, a Galactic molecular cloud complex in a spiral arm. We found higher fractional abundances of shock tracers, and possibly also higher dense gas fraction in NGC 3256 compared with W 51., Comment: 54 pages, 19 figures. Accepted for publication in ApJ. Some extensive figures will only be shown online-only in the published version

We present ALMA $^{12}$CO (J=1-0, 3-2 and 6-5), $^{13}$CO (J=1-0) and C$^{18}$O (J=1-0) observations of the local Ultra Luminous Infrared Galaxy, IRAS 13120-5453 (dubbed "The Yo-yo"). The morphologies of the three isotopic species differ, where $^{13}$CO shows a hole in emission towards the center. We measure integrated brightness temperature line ratios of $^{12}$CO/$^{13}$CO $\geq$ 60 (exceeding 200) and $^{13}$CO/C$^{18}$O $\leq$ 1 in the central region. Assuming optical thin emission, C$^{18}$O is more abundant than $^{13}$CO in several regions. The abundances within the central 500 pc are consistent with enrichment of the ISM via a young starburst ($<$7Myr), a top-heavy initial mass function or a combination of both., Comment: 7 pages, 3 figures, accepted for publication in ApJL

Aims. We present new IRAM Plateau de Bure Interferometer observations of Arp 220 in HCN, HCO$^{+}$, HN$^{13}$C J=1-0, C$_{2}$H N=1-0, SiO J = 2-1, HNCO J$_{k,k'}$ = 5$_{0,4}$ - 4$_{0,4}$, CH$_{3}$CN(6-5), CS J=2-1 and 5-4 and $^{13}$CO J=1-0 and 2-1 and of NGC 6240 in HCN, HCO$^{+}$ J = 1-0 and C$_{2}$H N = 1-0. In addition, we present Atacama Large Millimeter/submillmeter Array science verification observations of Arp 220 in CS J = 4-3 and CH$_{3}$CN(10-9). Various lines are used to analyse the physical conditions of the molecular gas including the [$^{12}$CO]/[$^{13}$CO] and [$^{12}$CO]/[C$^{18}$O] abundance ratios. These observations will be made available to the public. Methods. We create brightness temperature line ratio maps to present the different physical conditions across Arp 220 and NGC 6240. In addition, we use the radiative transfer code RADEX and a Monte Carlo Markov Chain likelihood code to model the $^{12}$CO, $^{13}$CO and C$^{18}$O lines of Arp 220 at ~2" (~700 pc) scales, where the $^{12}$CO and C$^{18}$O measurements were obtained from literature. Results. Line ratios of optically thick lines such as $^{12}$CO show smoothly varying ratios while the line ratios of optically thin lines such as $^{13}$CO show a east-west gradient across Arp 220. The HCN/HCO$^{+}$ line ratio differs between Arp 220 and NGC 6240, where Arp 220 has line ratios above 2 and NGC 6240 below 1. The radiative transfer analysis solution is consistent with a warm (~40 K), moderately dense (~10$^{3.4}$ cm$^{-3}$) molecular gas component averaged over the two nuclei. We find [$^{12}$CO]/[$^{13}$CO] and [$^{12}$CO]/[C$^{18}$O] abundance ratios of ~90 for both. The abundance enhancement of C$^{18}$O can be explained by stellar nucleosynthesis enrichment of the interstellar medium., Comment: 17 pages, 9 figures, accepted for publication in A&A

We present new high-resolution ALMA (13CO J=1-0 and J= 2-1) and CARMA (12CO and 13CO J=1-0) observations of two Luminous Infrared Galaxies (LIRGs): Arp 55 and NGC 2623. The new data are complementary to published and archival Submillimeter Array observations of 12CO J=2-1 and J=3-2. We perform a Bayesian likelihood non-local thermodynamic equilibrium analysis to constrain the molecular gas physical conditions such as temperature, column and volume densities and the [12CO]/[13CO] abundance ratio. For Arp 55, an early/intermediate staged merger, the line measurements are consistent with cold (~10-20 K), dense (>10$^{3.5}$ cm$^{-3}$) molecular gas. For NGC 2623, the molecular gas is warmer (~110 K) and less dense (~10$^{2.7}$ cm$^{-3}$). Since Arp 55 is an early/intermediate stage merger while NGC 2623 is a merger remnant, the difference in physical conditions may be an indicator of merger stage. Comparing the temperature and volume density of several LIRGs shows that the molecular gas, averaged over ~kpc scale, of advanced mergers is in general warmer and less dense than early/intermediate stage mergers. We also find that the [12CO]/[13CO] abundance ratio of NGC 2623 is unusually high (>250) when compared to the Milky Way; however, it follows a trend seen with other LIRGs in literature. This high [12CO]/[13CO] value is very likely due to stellar nucleosynthesis enrichment of the interstellar medium. On the other hand, Arp 55 has a more Galactic [12CO]/[13CO] value with the most probable [12CO]/[13CO] value being 20-30. We measure the CO-to-H2 conversion factor, $\alpha_{CO}$, to be ~0.1 and ~0.7 (3x10$^{-4}$/x$_{CO}$) M$_{\odot}$ (K km s$^{-1}$ pc$^{2}$)$^{-1}$ for Arp 55 and NGC 2623, respectively. Since Arp 55 is an early/intermediate , this suggests that the transition from a Galactic conversion factor to a LIRG value happens at an even earlier merger stage., Comment: 19 pages, 9 figures, Accepted for publication in ApJ

We present high-resolution (1".0) Atacama Large Millimeter/submillimeter Array (ALMA) observations of CO (1-0) and CO (2- 1) rotational transitions toward the nearby IR-luminous merger NGC 1614 supplemented with ALMA archival data of CO (3-2), and CO (6-5) transitions. The CO (6-5) emission arises from the starburst ring (central 590 pc in radius), while the lower-$J$ CO lines are distributed over the outer disk ($\sim$ 3.3 kpc in radius). Radiative transfer and photon dominated region (PDR) modeling reveal that the starburst ring has a single warmer gas component with more intense far-ultraviolet radiation field ($n_{\rm{H_2}}$ $\sim$ 10$^{4.6}$ cm$^{-3}$, $T_{\rm{kin}}$ $\sim$ 42 K, and $G_{\rm{0}}$ $\sim$ 10$^{2.7}$) relative to the outer disk ($n_{\rm{H_2}}$ $\sim$ 10$^{5.1}$ cm$^{-3}$, $T_{\rm{kin}}$ $\sim$ 22 K, and $G_{\rm{0}}$ $\sim$ 10$^{0.9}$). A two-phase molecular interstellar medium with a warm and cold ($>$ 70 K and $\sim$ 19 K) component is also an applicable model for the starburst ring. A possible source for heating the warm gas component is mechanical heating due to stellar feedback rather than PDR. Furthermore, we find evidence for non-circular motions along the north-south optical bar in the lower-$J$ CO images, suggesting a cold gas inflow. We suggest that star formation in the starburst ring is sustained by the bar-driven cold gas inflow, and starburst activities radiatively and mechanically power the CO excitation. The absence of a bright active galactic nucleus can be explained by a scenario that cold gas accumulating on the starburst ring is exhausted as the fuel for star formation, or is launched as an outflow before being able to feed to the nucleus., Comment: 20 pages, 19 figures, 2 tables, accepted for publication in ApJ

The cloud-scale density, velocity dispersion, and gravitational boundedness of the interstellar medium (ISM) vary within and among galaxies. In turbulent models, these properties play key roles in the ability of gas to form stars. New high fidelity, high resolution surveys offer the prospect to measure these quantities across galaxies. We present a simple approach to make such measurements and to test hypotheses that link small-scale gas structure to star formation and galactic environment. Our calculations capture the key physics of the Larson scaling relations, and we show good correspondence between our approach and a traditional "cloud properties" treatment. However, we argue that our method is preferable in many cases because of its simple, reproducible characterization of all emission. Using, low-J 12CO data from recent surveys, we characterize the molecular ISM at 60pc resolution in the Antennae, the Large Magellanic Cloud, M31, M33, M51, and M74. We report the distributions of surface density, velocity dispersion, and gravitational boundedness at 60pc scales and show galaxy-to-galaxy and intra-galaxy variations in each. The distribution of flux as a function of surface density appears roughly lognormal with a 1sigma width of ~0.3 dex, though the center of this distribution varies from galaxy to galaxy. The 60pc resolution line width and molecular gas surface density correlate well, which is a fundamental behavior expected for virialized or free-falling gas. Varying the measurement scale for the LMC and M31, we show that the molecular ISM has higher surface densities, lower line widths, and more self-gravity at smaller scales., Comment: Accepted for publication in the Astrophysical Journal, 36 pages (24+appendix), 21 figures (12+appendix), until publication high resolution version at http://www.astronomy.ohio-state.edu/~leroy.42/cloudscale.pdf

We present high-resolution archival Atacama Large Millimeter/submillimeter Array (ALMA) CO J=3-2 and J=6-5 and HCO+ J=4-3 observations and new CARMA CO and 13CO J=1-0 observations of the luminous infrared galaxy NGC 1614. The high-resolution maps show the previously identified ring-like structure while the CO J=3-2 map shows extended emission that traces the extended dusty features. We combined these new observations with previously published Submillimeter Array CO and 13CO J=2-1 observations to constrain the physical conditions of the molecular gas at a resolution of 230 pc using a radiative transfer code and a Bayesian likelihood analysis. At several positions around the central ring-like structure, the molecular gas is cold (20-40 K) and dense (> 10^{3.0} cm^{-3}) . The only region that shows evidence of a second molecular gas component is the "hole" in the ring. The CO-to-13CO abundance ratio is found to be greater than 130, more than twice the local interstellar medium value. We also measure the CO-to-H_{2} conversion factor, alpha_{CO}, to range from 0.9 to 1.5 M_sol (K km/s pc^{2})^{-1}., Comment: ApJL accepted

We present high-resolution (~2.5") observations of 12CO J=6-5 towards the luminous infrared galaxy VV 114 using the Submillimeter Array. We detect 12CO J=6-5 emission from the eastern nucleus of VV 114 but do not detect the western nucleus or the central region. We combine the new 12CO J=6-5 observations with previously published or archival low-J CO observations, that include 13CO J=1-0 Atacama Large Millimeter/submillimeter Array cycle 0 observations, to analyze the beam-averaged physical conditions of the molecular gas in the eastern nucleus. We use the radiative transfer code RADEX and a Bayesian likelihood code to constrain the temperature (T_kin), density (n(H2)) and column density (N(12CO)) of the molecular gas. We find that the most probable scenario for the eastern nucleus is a cold (T_kin = 38 K), moderately dense (n(H2) = 10^2.89 cm^-3) molecular gas component. We find the most probable 12CO to 13CO abundance ratio ([12CO]/[13CO]) is 229, roughly three times higher than the Milky Way value. This high abundance ratio may explain the observed high 12CO/ 13CO line ratio (> 25). The unusual 13CO J=2-1/J=1-0 line ratio of 0.6 is produced by a combination of moderate 13CO optical depths (tau = 0.4 - 1.1) and extremely subthermal excitation temperatures. We measure the CO-to-H2 conversion factor, alpha_co to be 0.5 M_sol (K km s^-1 pc^2)^-1, which agrees with the widely used factor for ultra luminous infrared galaxies of Downes & Solomon (1998; alpha_co =0.8 M_sol (K km s^-1 pc^2)^-1)., Comment: 33 pages, 7 figures. Accepted for publication in ApJ

We have used high resolution (~2.3") observations of the local (D = 46 Mpc) luminous infrared galaxy Arp 299 to map out the physical properties of the molecular gas which provides the fuel for its extreme star formation activity. The 12CO J=3-2, 12CO J=2-1 and 13CO J=2-1 lines were observed with the Submillimeter Array and the short spacings of the 12CO J=2-1 and J=3-2 observations have been recovered using James Clerk Maxwell Telescope single dish observations. We use the radiative transfer code RADEX to estimate the physical properties (density, column density and temperature) of the different regions in this system. The RADEX solutions of the two galaxy nuclei, IC 694 and NGC 3690, are consistent with a wide range of gas components, from warm moderately dense gas with T_{kin} > 30 K and n(H_{2}) ~ 0.3 - 3 x 10^{3} cm^{-3} to cold dense gas with T_{kin} ~ 10-30 K and n(H_{2}) > 3 x 10^{3} cm^{-3}. The overlap region is shown to have a better constrained solution with T_{\rm{kin}}$ ~ 10-50 K and n(H_{2}) ~ 1-30 x 10^{3} cm^{-3}. We estimate the gas masses and star formation rates of each region in order to derive molecular gas depletion times. The depletion times of all regions (20-60 Myr) are found to be about 2 orders of magnitude lower than those of normal spiral galaxies. This rapid depletion time can probably be explained by a high fraction of dense gas on kiloparsec scales in Arp 299. We estimate the CO-to-H_{2} factor, \alpha_{co} to be 0.4 \pm 0.3 (3 x 10^{-4}/ x_{CO}) M_{sol} (K km s^{-1} pc^{2})^{-1} for the overlap region. This value agrees well with values determined previously for more advanced merger systems., Comment: 24 pages, 4 figures, ApJ accepted