Source author record

Annie Hughes

Annie Hughes appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

44works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

44 published item(s)

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

Estimate of the carbon footprint of astronomical research infrastructures

The carbon footprint of astronomical research is an increasingly topical issue with first estimates of research institute and national community footprints having recently been published. As these assessments have typically excluded the contribution of astronomical research infrastructures, we complement these studies by providing an estimate of the contribution of astronomical space missions and ground-based observatories using greenhouse gas emission factors that relates cost and payload mass to carbon footprint. We find that worldwide active astronomical research infrastructures currently have a carbon footprint of 20.3$\pm$3.3 MtCO$_2$ equivalent (CO$_2$e) and an annual emission of 1,169$\pm$249 ktCO$_2$e yr$^{-1}$ corresponding to a footprint of 36.6$\pm$14.0 tCO$_2$e per year per astronomer. Compared with contributions from other aspects of astronomy research activity, our results suggest that research infrastructures make the single largest contribution to the carbon footprint of an astronomer. We discuss the limitations and uncertainties of our method and explore measures that can bring greenhouse gas emissions from astronomical research infrastructures towards a sustainable level.

preprint2022arXiv

The carbon footprint of IRAP

We present an assessment of the greenhouse gases emissions of the Institute for Research in Astrophysics and Planetology (IRAP), located in Toulouse (France). It was performed following the established "Bilan Carbone" methodology, over a large scope compared to similar previous studies, including in particular the contribution from the purchase of goods and services as well as IRAP's use of external research infrastructures, such as ground-based observatories and space-borne facilities. The carbon footprint of the institute for the reference year 2019 is 7400 +/- 900 tCO2e. If we exclude the contribution from external research infrastructures to focus on a restricted perimeter over which the institute has some operational control, IRAP's emissions in 2019 amounted to 3300 +/- 400 tCO2e. Over the restricted perimeter, the contribution from purchasing goods and services is dominant, about 40% of the total, slightly exceeding the contribution from professional travel including hotel stays, which accounts for 38%. Local infrastructures make a smaller contribution to IRAP's carbon footprint, about 25% over the restricted perimeter. We note that this repartition may be specific to IRAP, since the energy used to produce the electricity and heating has a relatively low carbon footprint. Over the full perimeter, the large share from the use of ground-based observatories and space-borne facilities and the fact that the majority of IRAP purchases are related to instrument development indicate that research infrastructures represent the most significant challenge for reducing the carbon footprint of research at our institute. With ~260 staff members employed, our results imply that performing research in astronomy and astrophysics at IRAP according to the standards of 2019 produces average GHG emissions of 28 tCO2e/yr per person involved in that activity (Abridged).

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

preprint2021arXiv

Extragalactic magnetism with SOFIA (Legacy Program) -- II: A magnetically-driven flow in the starburst ring of NGC 1097

Galactic bars are frequent in disk galaxies and they may support the transfer of matter towards the central engine of active nuclei. The barred galaxy NGC 1097 has magnetic forces controlling the gas flow at several kpc scales, which suggest that magnetic fields (B-fields) are dynamically important along the bar and nuclear ring. However, the effect of the B-field on the gas flows in the central kpc scale has not been characterized. Using thermal polarized emission at $89$ $μ$m with HAWC+/SOFIA, here, we measure that the polarized flux is spatially located at the contact regions of the outer-bar with the starburst ring. The linear polarization decomposition analysis shows that the $89$ $μ$m and radio ($3.5$ and $6.2$ cm) polarization traces two different modes, $m$, of the B-field: a constant B-field orientation and dominated by $m=0$ at $89$ $μ$m, and a spiral B-field dominated by $m=2$ at radio. We show that the B-field at 89 $μ$m is concentrated in the warmest region of a shock driven by the galactic-bar dynamics in the contact regions between the outer-bar with the starburst ring. Radio polarization traces a superposition of the spiral B-field outside and within the starburst ring. According to Faraday rotation measures between $3.5$ and $6.2$ cm, the radial component of the B-field along the contact regions points toward the galaxy's center on both sides. We conclude that gas streams outside and within the starburst ring follow the B-field, which feeds the black hole with matter from the host galaxy.

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" \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

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 model for the onset of self-gravitation and star formation in molecular gas governed by galactic forces: II. the bottleneck to collapse set by cloud-environment decoupling

In Meidt et al. (2018), we showed that gas kinematics on the scale of individual molecular clouds are not dominated by self-gravity but also track a component that originates with orbital motion in the potential of the host galaxy. This agrees with observed cloud line widths, which show systematic variations from virial motions with environment, pointing at the influence of the galaxy potential. In this paper, we hypothesize that these motions act to slow down the collapse of gas and so help regulate star formation. Extending the results of Meidt et al. (2018), we derive a dynamical collapse timescale that approaches the free-fall time only once the gas has fully decoupled from the galactic potential. Using this timescale we make predictions for how the fraction of free-falling, strongly self-gravitating gas varies throughout the disks of star-forming galaxies. We also use this collapse timescale to predict variations in the molecular gas star formation efficiency, which is lowered from a maximum, feedback-regulated level in the presence of strong coupling to the galactic potential. Our model implies that gas can only decouple from the galaxy to collapse and efficiently form stars deep within clouds. We show that this naturally explains the observed drop in star formation rate per unit gas mass in the Milky Way's CMZ and other galaxy centers. The model for a galactic bottleneck to star formation also agrees well with resolved observations of dense gas and star formation in galaxy disks and the properties of local clouds.

preprint2020arXiv

C18O, 13CO, and 12CO abundances and excitation temperatures in the Orion B molecular cloud: An analysis of the precision achievable when modeling spectral line within the Local Thermodynamic Equilibrium approximation

CO isotopologue transitions are routinely observed in molecular clouds to probe the column density of the gas, the elemental ratios of carbon and oxygen, and to trace the kinematics of the environment. We aim at estimating the abundances, excitation temperatures, velocity field and velocity dispersions of the three main CO isotopologues towards a subset of the Orion B molecular cloud. We use the Cramer Rao Bound (CRB) technique to analyze and estimate the precision of the physical parameters in the framework of local-thermodynamic-equilibrium excitation and radiative transfer with an additive white Gaussian noise. We propose a maximum likelihood estimator to infer the physical conditions from the 1-0 and 2-1 transitions of CO isotopologues. Simulations show that this estimator is unbiased and efficient for a common range of excitation temperatures and column densities (Tex > 6 K, N > 1e14 - 1e15 cm-2). Contrary to the general assumptions, the different CO isotopologues have distinct excitation temperatures, and the line intensity ratios between different isotopologues do not accurately reflect the column density ratios. We find mean fractional abundances that are consistent with previous determinations towards other molecular clouds. However, significant local deviations are inferred, not only in regions exposed to UV radiation field but also in shielded regions. These deviations result from the competition between selective photodissociation, chemical fractionation, and depletion on grain surfaces. We observe that the velocity dispersion of the C18O emission is 10% smaller than that of 13CO. The substantial gain resulting from the simultaneous analysis of two different rotational transitions of the same species is rigorously quantified. The CRB technique is a promising avenue for analyzing the estimation of physical parameters from the fit of spectral lines.

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

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'' 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

Quantitative inference of the $H_2$ column densities from 3 mm molecular emission: A case study towards Orion B

Molecular hydrogen being unobservable in cold molecular clouds, the column density measurements of molecular gas currently rely either on dust emission observation in the far-IR or on star counting. (Sub-)millimeter observations of numerous trace molecules are effective from ground based telescopes, but the relationships between the emission of one molecular line and the H2 column density (NH2) is non-linear and sensitive to excitation conditions, optical depths, abundance variations due to the underlying physico-chemistry. We aim to use multi-molecule line emission to infer NH2 from radio observations. We propose a data-driven approach to determine NH2 from radio molecular line observations. We use supervised machine learning methods (Random Forests) on wide-field hyperspectral IRAM-30m observations of the Orion B molecular cloud to train a predictor of NH2, using a limited set of molecular lines as input, and the Herschel-based dust-derived NH2 as ground truth output. For conditions similar to the Orion B molecular cloud, we obtain predictions of NH2 within a typical factor of 1.2 from the Herschel-based estimates. An analysis of the contributions of the different lines to the predictions show that the most important lines are $^{13}$CO(1-0), $^{12}$CO(1-0), C$^{18}$O(1-0), and HCO$^+$(1-0). A detailed analysis distinguishing between diffuse, translucent, filamentary, and dense core conditions show that the importance of these four lines depends on the regime, and that it is recommended to add the N$_2$H$^+$(1-0) and CH$_3$OH(20-10) lines for the prediction of NH2 in dense core conditions. This article opens a promising avenue to directly infer important physical parameters from the molecular line emission in the millimeter domain. The next step will be to try to infer several parameters simultaneously (e.g., NH2 and far-UV illumination field) to further test the method. [Abridged]

preprint2020arXiv

Tracers of the ionization fraction in dense and translucent gas: I. Automated exploitation of massive astrochemical model grids

The ionization fraction plays a key role in the physics and chemistry of the neutral interstellar medium, from controlling the coupling of the gas to the magnetic field to allowing fast ion-neutral reactions that drive interstellar chemistry. Most estimations of the ionization fraction have relied on deuterated species such as DCO+, whose detection is limited to dense cores representing an extremely small fraction of the volume of the giant molecular clouds they are part of. As large field-of-view hyperspectral maps become available, new tracers may be found. We search for the best observable tracers of the ionization fraction based on a grid of astrochemical models. We build grids of models that sample randomly a large space of physical conditions (unobservable quantities such as gas density, temperature, etc.) and compute the corresponding observables (line intensities, column densities) and the ionization fraction. We estimate the predictive power of each potential tracer by training a Random Forest model to predict the ionization fraction from that tracer, based on these model grids. In both translucent medium and cold dense medium conditions, several observable tracers with very good predictive power for the ionization fraction are found. Several tracers in cold dense medium conditions are found to be better and more widely applicable than the traditional DCO+/HCO+ ratio. We also provide simpler analytical fits for estimating the ionization fraction from the best tracers, and for estimating the associated uncertainties. We discuss the limitations of the present study and select a few recommended tracers in both types of conditions. The method presented here is very general and can be applied to the measurement of any other quantity of interest (cosmic ray flux, elemental abundances, etc.) from any type of model (PDR models, time-dependent chemical models, etc.). (abridged)

preprint2019arXiv

Relations Between Molecular Cloud Structure Sizes and Line Widths in the Large Magellanic Cloud

We present a comparative study of the size-line width relation for substructures within six molecular clouds in the Large Magellanic Cloud (LMC) mapped with the Atacama Large Millimeter/submillimeter Array (ALMA). Our sample extends our previous study, which compared a Planck detected cold cloud in the outskirts of the LMC with the 30 Doradus molecular cloud and found the typical line width for 1 pc radius structures to be 5 times larger in 30 Doradus. By observing clouds with intermediate levels of star formation activity, we find evidence that line width at a given size increases with increasing local and cloud-scale 8$μ$m intensity. At the same time, line width at a given size appears to independently correlate with measures of mass surface density. Our results suggest that both virial-like motions due to gravity and local energy injection by star formation feedback play important roles in determining intracloud dynamics.

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' 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.

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.

preprint2018arXiv

Cloud-Scale Molecular Gas Properties in 15 Nearby Galaxies

We measure the velocity dispersion, $σ$, and surface density, $Σ$, of the molecular gas in nearby galaxies from CO spectral line cubes with spatial resolution $45$-$120$ pc, matched to the size of individual giant molecular clouds. Combining $11$ galaxies from the PHANGS-ALMA survey with $4$ targets from the literature, we characterize ${\sim}30,000$ independent sightlines where CO is detected at good significance. $Σ$ and $σ$ show a strong positive correlation, with the best-fit power law slope close to the expected value for resolved, self-gravitating clouds. This indicates only weak variation in the virial parameter $α_\mathrm{vir}\proptoσ^2/Σ$, which is ${\sim}1.5$-$3.0$ for most galaxies. We do, however, observe enormous variation in the internal turbulent pressure $P_\mathrm{turb}\proptoΣ\,σ^2$, which spans ${\sim}5\rm\;dex$ across our sample. We find $Σ$, $σ$, and $P_\mathrm{turb}$ to be systematically larger in more massive galaxies. The same quantities appear enhanced in the central kpc of strongly barred galaxies relative to their disks. Based on sensitive maps of M31 and M33, the slope of the $σ$-$Σ$ relation flattens at $Σ\lesssim10\rm\;M_\odot\,pc^{-2}$, leading to high $σ$ for a given $Σ$ and high apparent $α_\mathrm{vir}$. This echoes results found in the Milky Way, and likely originates from a combination of lower beam filling factors and a stronger influence of local environment on the dynamical state of molecular gas in the low density regime.

preprint2016arXiv

A Portrait of Cold Gas in Galaxies at 60pc Resolution and a Simple Method to Test Hypotheses That Link Small-Scale ISM Structure to Galaxy-Scale Processes

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.

preprint2016arXiv

Radiative and mechanical feedback into the molecular gas in the Large Magellanic Cloud. I. N159W

We present Herschel SPIRE Fourier Transform Spectrometer (FTS) observations of N159W, an active star-forming region in the Large Magellanic Cloud (LMC). In our observations, a number of far-infrared cooling lines including CO(4-3) to CO(12-11), [CI] 609 and 370 micron, and [NII] 205 micron are clearly detected. With an aim of investigating the physical conditions and excitation processes of molecular gas, we first construct CO spectral line energy distributions (SLEDs) on 10 pc scales by combining the FTS CO transitions with ground-based low-J CO data and analyze the observed CO SLEDs using non-LTE radiative transfer models. We find that the CO-traced molecular gas in N159W is warm (kinetic temperature of 153-754 K) and moderately dense (H2 number density of (1.1-4.5)e3 cm-3). To assess the impact of the energetic processes in the interstellar medium on the physical conditions of the CO-emitting gas, we then compare the observed CO line intensities with the models of photodissociation regions (PDRs) and shocks. We first constrain the properties of PDRs by modelling Herschel observations of [OI] 145, [CII] 158, and [CI] 370 micron fine-structure lines and find that the constrained PDR components emit very weak CO emission. X-rays and cosmic-rays are also found to provide a negligible contribution to the CO emission, essentially ruling out ionizing sources (ultraviolet photons, X-rays, and cosmic-rays) as the dominant heating source for CO in N159W. On the other hand, mechanical heating by low-velocity C-type shocks with ~10 km/s appears sufficient enough to reproduce the observed warm CO.

preprint2016arXiv

The Relationship Between Molecular Gas, HI, and Star Formation in the Low-Mass, Low-Metallicity Magellanic Clouds

The Magellanic Clouds provide the only laboratory to study the effect of metallicity and galaxy mass on molecular gas and star formation at high (~20 pc) resolution. We use the dust emission from HERITAGE Herschel data to map the molecular gas in the Magellanic Clouds, avoiding the known biases of CO emission as a tracer of H$_{2}$. Using our dust-based molecular gas estimates, we find molecular gas depletion times of ~0.4 Gyr in the LMC and ~0.6 SMC at 1 kpc scales. These depletion times fall within the range found for normal disk galaxies, but are shorter than the average value, which could be due to recent bursts in star formation. We find no evidence for a strong intrinsic dependence of the molecular gas depletion time on metallicity. We study the relationship between gas and star formation rate across a range in size scales from 20 pc to ~1 kpc, including how the scatter in molecular gas depletion time changes with size scale, and discuss the physical mechanisms driving the relationships. We compare the metallicity-dependent star formation models of Ostriker, McKee, and Leroy (2010) and Krumholz (2013) to our observations and find that they both predict the trend in the data, suggesting that the inclusion of a diffuse neutral medium is important at lower metallicity.

preprint2015arXiv

Graph-based interpretation of the Molecular Interstellar Medium Segmentation

We present a generalization of the Giant Molecular Cloud (GMC) identification problem based on cluster analysis. The method we designed, SCIMES (Spectral Clustering for Interstellar Molecular Emission Segmentation) considers the dendrogram of emission in the broader framework of graph theory and utilizes spectral clustering to find discrete regions with similar emission properties. For Galactic molecular cloud structures, we show that the characteristic volume and/or integrated CO luminosity are useful criteria to define the clustering, yielding emission structures that closely reproduce "by-eye" identification results. SCIMES performs best on well-resolved, high-resolution data, making it complementary to other available algorithms. Using 12CO(1-0) data for the Orion-Monoceros complex, we demonstrate that SCIMES provides robust results against changes of the dendrogram-construction parameters, noise realizations and degraded resolution. By comparing SCIMES with other cloud decomposition approaches, we show that our method is able to identify all canonical clouds of the Orion-Monoceros region, avoiding the over-division within high resolution survey data that represents a common limitation of several decomposition algorithms. The Orion-Monoceros objects exhibit hierarchies and size-line width relationships typical to the turbulent gas in molecular clouds, although "the Scissors" region deviates from this common description. SCIMES represents a significant step forward in moving away from pixel-based cloud segmentation towards a more physical-oriented approach, where virtually all properties of the ISM can be used for the segmentation of discrete objects.

preprint2015arXiv

Short GMC lifetimes: an observational estimate with the PdBI Arcsecond Whirlpool Survey (PAWS)

We describe and execute a novel approach to observationally estimate the lifetimes of giant molecular clouds (GMCs). We focus on the cloud population between the two main spiral arms in M51 (the inter-arm region) where cloud destruction via shear and star formation feedback dominates over formation processes. By monitoring the change in GMC number densities and properties from one side of the inter-arm to the other, we estimate the lifetime as a fraction of the inter-arm travel time. We find that GMC lifetimes in M51's inter-arm are finite and short, 20 to 30 Myr. Such short lifetimes suggest that cloud evolution is influenced by environment, in which processes can disrupt GMCs after a few free-fall times. Over most of the region under investigation shear appears to regulate the lifetime. As the shear timescale increases with galactocentric radius, we expect cloud destruction to switch primarily to star formation feedback at larger radii. We identify a transition from shear- to feedback-dominated disruption through a change in the behavior of the GMC number density. The signature suggests that shear is more efficient at completely dispersing clouds, whereas feedback transforms the population, e.g. by fragmenting high mass clouds into lower mass pieces. Compared to the characteristic timescale for molecular hydrogen in M51, our short lifetimes suggest that gas can remain molecular while clouds disperse and reassemble. We propose that galaxy dynamics regulates the cycling of molecular material from diffuse to bound (and ultimately star-forming) objects, contributing to long observed molecular depletion times in normal disk galaxies. We also speculate that, in more extreme environments such as elliptical galaxies and concentrated galaxy centers, star formation can be suppressed when the shear timescale becomes so short that some clouds can not survive to collapse and form stars.

preprint2015arXiv

Spatial and Spectral Modeling of the Gamma-ray Distribution in the Large Magellanic Cloud

We perform spatial and spectral analyses of the Large Magellanic Cloud (LMC) gamma-ray emission collected over 66 months by the Fermi Gamma-ray Space Telescope. In our spatial analysis, we model the LMC cosmic-ray distribution and gamma-ray production using observed maps of the LMC interstellar medium, star formation history, interstellar radiation field, and synchrotron emission. We use bootstrapping of the data to quantify the robustness of spatial model performance. We model the LMC gamma-ray spectrum using fitting functions derived from the physics of $π^0$ decay, bremsstrahlung, and inverse Compton scattering. We find the integrated gamma-ray flux of the LMC from 200 MeV to 20 GeV to be $1.38 \pm 0.02 \times 10^{-7}$ ph cm$^{-2}$ s$^{-1}$, of which we attribute about 10% to inverse Compton scattering and 40% to bremsstrahlung. From our work, we conclude that the spectral index of the LMC cosmic-ray proton population is 2.4$\pm$0.2, and we find that cosmic-ray energy loss through gamma-ray production is concentrated within a few hundred pc of acceleration sites. Assuming cosmic-ray energy equipartition with magnetic fields, we estimate LMC cosmic rays encounter an average magnetic field strength ~3 $μ$G.

preprint2014arXiv

Dust and Gas in the Magellanic Clouds from the HERITAGE Herschel Key Project. I. Dust Properties and Insights into the Origin of the Submm Excess Emission

The dust properties in the Large and Small Magellanic Clouds are studied using the HERITAGE Herschel Key Project photometric data in five bands from 100 to 500 micron. Three simple models of dust emission were fit to the observations: a single temperature blackbody modified by a power- law emissivity (SMBB), a single temperature blackbody modified by a broken power-law emissivity (BEMBB), and two blackbodies with different temperatures, both modified by the same power-law emissivity (TTMBB). Using these models we investigate the origin of the submm excess; defined as the submillimeter (submm) emission above that expected from SMBB models fit to observations < 200 micron. We find that the BEMBB model produces the lowest fit residuals with pixel-averaged 500 micron submm excesses of 27% and 43% for the LMC and SMC, respectively. Adopting gas masses from previous works, the gas-to-dust ratios calculated from our the fitting results shows that the TTMBB fits require significantly more dust than are available even if all the metals present in the interstellar medium (ISM) were condensed into dust. This indicates that the submm excess is more likely to be due to emissivity variations than a second population of colder dust. We derive integrated dust masses of (7.3 +/- 1.7) x 10^5 and (8.3 +/- 2.1) times 10^4 M(sun) for the LMC and SMC, respectively. We find significant correlations between the submm excess and other dust properties; further work is needed to determine the relative contributions of fitting noise and ISM physics to the correlations.

preprint2014arXiv

Dust and Gas in the Magellanic Clouds from the HERITAGE Herschel Key Project. II. Gas-to-Dust Ratio Variations across ISM Phases

The spatial variations of the gas-to-dust ratio (GDR) provide constraints on the chemical evolution and lifecycle of dust in galaxies. We examine the relation between dust and gas at 10-50 pc resolution in the Large and Small Magellanic Clouds (LMC and SMC) based on Herschel far-infrared (FIR), H I 21 cm, CO, and Halpha observations. In the diffuse atomic ISM, we derive the gas-to-dust ratio as the slope of the dust-gas relation and find gas-to-dust ratios of 380+250-130 in the LMC, and 1200+1600-420 in the SMC, not including helium. The atomic-to-molecular transition is located at dust surface densities of 0.05 Mo pc-2 in the LMC and 0.03 Mo pc-2 in the SMC, corresponding to AV ~ 0.4 and 0.2, respectively. We investigate the range of CO-to-H2 conversion factor to best account for all the molecular gas in the beam of the observations, and find upper limits on XCO to be 6x1020 cm-2 K-1 km-1 s in the LMC (Z=0.5Zo) at 15 pc resolution, and 4x 1021 cm-2 K-1 km-1 s in the SMC (Z=0.2Zo) at 45 pc resolution. In the LMC, the slope of the dust-gas relation in the dense ISM is lower than in the diffuse ISM by a factor ~2, even after accounting for the effects of CO-dark H2 in the translucent envelopes of molecular clouds. Coagulation of dust grains and the subsequent dust emissivity increase in molecular clouds, and/or accretion of gas-phase metals onto dust grains, and the subsequent dust abundance (dust-to-gas ratio) increase in molecular clouds could explain the observations. In the SMC, variations in the dust-gas slope caused by coagulation or accretion are degenerate with the effects of CO-dark H2. Within the expected 5--20 times Galactic XCO range, the dust-gas slope can be either constant or decrease by a factor of several across ISM phases. Further modeling and observations are required to break the degeneracy between dust grain coagulation, accretion, and CO-dark H2.

preprint2014arXiv

From Gas to Stars in Energetic Environments: Dense Gas Clumps in the 30 Doradus Region Within the Large Magellanic Cloud

We present parsec scale interferometric maps of HCN(1-0) and HCO$^{+}$(1-0) emission from dense gas in the star-forming region 30 Doradus, obtained using the Australia Telescope Compact Array. This extreme star-forming region, located in the Large Magellanic Cloud (LMC), is characterized by a very intense ultraviolet ionizing radiation field and sub-solar metallicity, both of which are expected to impact molecular cloud structure. We detect 13 bright, dense clumps within the 30 Doradus-10 giant molecular cloud. Some of the clumps are aligned along a filamentary structure with a characteristic spacing that is consistent with formation via the varicose fluid instability. Our analysis shows that the filament is gravitationally unstable and collapsing to form stars. There is a good correlation between HCO$^{+}$ emission in the filament and signatures of recent star formation activity including H$_{2}$O masers and young stellar objects (YSOs). YSOs seem to continue along the same direction of the filament toward the massive compact star cluster R136 in the southwest. We present detailed comparisons of clump properities (masses, linewidths, sizes) in 30Dor-10 to those in other star forming regions of the LMC (N159, N113, N105, N44). Our analysis shows that the 30Dor-10 clumps have similar mass but wider linewidths and similar HCN/HCO$^{+}$ (1-0) line ratios as clumps detected in other LMC star-forming regions. Our results suggest that the dense molecular gas clumps in the interior of 30Dor-10 are well-shielded against the intense ionizing field that is present in the 30Doradus region.

preprint2014arXiv

Molecular hydrogen emission in the interstellar medium of the Large Magellanic Cloud

We present the detection and analysis of molecular hydrogen emission toward ten interstellar regions in the Large Magellanic Cloud. We examined low-resolution infrared spectral maps of twelve regions obtained with the Spitzer infrared spectrograph (IRS). The pure rotational 0--0 transitions of H$_2$ at 28.2 and 17.1${\,\rm μm}$ are detected in the IRS spectra for ten regions. The higher level transitions are mostly upper limit measurements except for three regions, where a 3$σ$ detection threshold is achieved for lines at 12.2 and 8.6${\,\rm μm}$. The excitation diagrams of the detected H$_2$ transitions are used to determine the warm H$_2$ gas column density and temperature. The single-temperature fits through the lower transition lines give temperatures in the range $86-137\,{\rm K}$. The bulk of the excited H$_2$ gas is found at these temperatures and contributes $\sim$5-17% to the total gas mass. We find a tight correlation of the H$_2$ surface brightness with polycyclic aromatic hydrocarbon and total infrared emission, which is a clear indication of photo-electric heating in photodissociation regions. We find the excitation of H$_2$ by this process is equally efficient in both atomic and molecular dominated regions. We also present the correlation of the warm H$_2$ physical conditions with dust properties. The warm H$_2$ mass fraction and excitation temperature show positive correlations with the average starlight intensity, again supporting H$_2$ excitation in photodissociation regions.

preprint2014arXiv

The PdBI Arcsecond Whirlpool Survey (PAWS): Environmental Dependence of Giant Molecular Cloud Properties in M51

Using data from the PdBI Arcsecond Whirlpool Survey (PAWS), we have generated the largest extragalactic Giant Molecular Cloud (GMC) catalog to date, containing 1,507 individual objects. GMCs in the inner M51 disk account for only 54% of the total 12CO(1-0) luminosity of the survey, but on average they exhibit physical properties similar to Galactic GMCs. We do not find a strong correlation between the GMC size and velocity dispersion, and a simple virial analysis suggests that 30% of GMCs in M51 are unbound. We have analyzed the GMC properties within seven dynamically-motivated galactic environments, finding that GMCs in the spiral arms and in the central region are brighter and have higher velocity dispersions than inter-arm clouds. Globally, the GMC mass distribution does not follow a simple power law shape. Instead, we find that the shape of the mass distribution varies with galactic environment: the distribution is steeper in inter-arm region than in the spiral arms, and exhibits a sharp truncation at high masses for the nuclear bar region. We propose that the observed environmental variations in the GMC properties and mass distributions are a consequence of the combined action of large-scale dynamical processes and feedback from high mass star formation. We describe some challenges of using existing GMC identification techniques for decomposing the 12CO(1-0) emission in molecule-rich environments, such as M51's inner disk.

preprint2014arXiv

The PdBI Arcsecond Whirlpool Survey (PAWS): Multi-phase cold gas kinematic of M51

The kinematic complexity and the favorable position of M51 on the sky make this galaxy an ideal target to test different theories of spiral arm dynamics. Taking advantage of the new high resolution PdBI Arcsecond Whirlpool Survey (PAWS) data, we undertake a detailed kinematic study of M51 to characterize and quantify the origin and nature of the non-circular motions. Using a tilted-ring analysis supported by several other archival datasets we update the estimation of M51's position angle (PA=(173 +/- 3) deg) and inclination (i=(22 +/- 5) deg). Harmonic decomposition of the high resolution (40 pc) CO velocity field shows the first kinematic evidence of an m=3 wave in the inner disk of M51 with a corotation at R(CR,m=3)=1.1 +/- 0.1 kpc and a pattern speed of Omega_p(m=3) = 140 km/(s kpc). This mode seems to be excited by the nuclear bar, while the beat frequencies generated by the coupling between the m=3 mode and the main spiral structure confirm its density-wave nature. We observe also a signature of an m=1 mode that is likely responsible for the lopsidedness of M51 at small and large radii. We provide a simple method to estimate the radial variation of the amplitude of the spiral perturbation (Vsp) attributed to the different modes. The main spiral arm structure has <Vsp>=50-70 km/s, while the streaming velocity associated with the m=1 and m=3 modes is, in general, 2 times lower. Our joint analysis of HI and CO velocity fields at low and high spatial resolution reveals that the atomic and molecular gas phases respond differently to the spiral perturbation due to their different vertical distribution and emission morphology.

preprint2013arXiv

A Comparative Study of Giant Molecular Clouds in M51, M33 and the Large Magellanic Cloud

We compare the properties of giant molecular clouds (GMCs) in M51 identified by the Plateau de Bure Interferometer Whirlpool Arcsecond Survey (PAWS) with GMCs identified in wide-field, high resolution surveys of CO emission in M33 and the Large Magellanic Cloud (LMC). We find that GMCs in M51 are larger, brighter and have higher velocity dispersions relative to their size than equivalent structures in M33 and the LMC. These differences imply that there are genuine variations in the average mass surface density of the different GMC populations. To explain this, we propose that the pressure in the interstellar medium surrounding the GMCs plays a role in regulating their density and velocity dispersion. We find no evidence for a correlation between size and linewidth in any of M51, M33 or the LMC when the CO emission is decomposed into GMCs, although moderately robust correlations are apparent when regions of contiguous CO emission (with no size limitation) are used. Our work demonstrates that observational bias remains an important obstacle to the identification and study of extragalactic GMC populations using CO emission, especially in molecule-rich galactic environments.

preprint2013arXiv

ALMA Resolves 30 Doradus: Sub-parsec Molecular Cloud Structure Near the Closest Super-Star Cluster

We present ALMA observations of 30 Doradus -- the highest resolution view of molecular gas in an extragalactic star formation region to date (~0.4pc x 0.6pc). The 30Dor-10 cloud north of R136 was mapped in 12CO 2-1, 13CO 2-1, C18O 2-1, 1.3mm continuum, the H30alpha recombination line, and two H2CO 3-2 transitions. Most 12CO emission is associated with small filaments and clumps (<1pc, ~1000 Msun at the current resolution). Some clumps are associated with protostars, including "pillars of creation" photoablated by intense radiation from R136. Emission from molecular clouds is often analyzed by decomposition into approximately beam-sized clumps. Such clumps in 30 Doradus follow similar trends in size, linewidth, and surface density to Milky Way clumps. The 30 Doradus clumps have somewhat larger linewidths for a given size than predicted by Larson's scaling relation, consistent with pressure confinement. They extend to higher surface density at a given size and linewidth compared to clouds studied at 10pc resolution. These trends are also true of clumps in Galactic infrared-dark clouds; higher resolution observations of both environments are required. Consistency of clump masses calculated from dust continuum, CO, and the virial theorem reveals that the CO abundance in 30 Doradus clumps is not significantly different from the LMC mean, but the dust abundance may be reduced by ~2. There are no strong trends in clump properties with distance from R136; dense clumps are not strongly affected by the external radiation field, but there is a modest trend towards lower dense clump filling fraction deeper in the cloud.

preprint2013arXiv

Clumping and the Interpretation of kpc-Scale Maps of the Interstellar Medium: Smooth HI and Clumpy, Variable H2 Surface Density

Many recent models consider the structure of individual interstellar medium (ISM) clouds as a way to explain observations of large parts of galaxies. To compare such models to observations, one must understand how to translate between surface densities observed averaging over large (~kpc) scales and surface densities on the scale of individual clouds (~pc scale), which are treated by models. We define a "clumping factor" that captures this translation as the ratio of the mass-weighted surface density, which is often the quantity of physical interest, to the area-weighted surface density, which is observed. We use high spatial resolution (sub-kpc) maps of CO and HI emission from nearby galaxies to measure the clumping factor of both atomic and molecular gas. The molecular and atomic ISM exhibit dramatically different degrees of clumping. As a result, the ratio H2/HI measured at ~kpc resolution cannot be trivially interpreted as a cloud-scale ratio of surface densities. HI emission appears very smooth, with a clumping factor of only ~1.3. Based on the scarce and heterogeneous high resolution data available, CO emission is far more clumped with a widely variable clumping factor, median ~7 for our heterogeneous data. Our measurements do not provide evidence for a universal mass-weighted surface density of molecular gas, but also cannot conclusively rule out such a scenario. We suggest that a more sophisticated treatment of molecular ISM structure, one informed by high spatial resolution CO maps, is needed to link cloud-scale models to kpc-scale observations of galaxies.

preprint2013arXiv

Gas Kinematics on GMC scales in M51 with PAWS: cloud stabilization through dynamical pressure

We use the high spatial and spectral resolution of the PAWS CO(1-0) survey of the inner 9 kpc of the iconic spiral galaxy M51 to examine the effect of gas streaming motions on the star-forming properties of individual GMCs. We compare our view of gas flows in M51 -- which arise due to departures from axi-symmetry in the gravitational potential (i.e. the nuclear bar and spiral arms) -- with the global pattern of star formation as traced by Halpha and 24μm emission. We find that the dynamical environment of GMCs strongly affects their ability to form stars, in the sense that GMCs situated in regions with large streaming motions can be stabilized, while similarly massive GMCs in regions without streaming go on to efficiently form stars. We argue that this is the result of reduced surface pressure felt by clouds embedded in an ambient medium undergoing large streaming motions, which prevents collapse. Indeed, the variation in gas depletion time expected based on the observed streaming motions throughout the disk of M51 quantitatively agrees with the variation in observed gas depletion time scale. The example of M51 shows that streaming motions, triggered by gravitational instabilities in the form of bars and spiral arms, can alter the star formation law; this can explain the variation in gas depletion time among galaxies with different masses and morphologies. In particular, we can explain the long gas depletion times in spiral galaxies compared to dwarf galaxies and starbursts. We suggest that adding a dynamical pressure term to the canonical free-fall time produces a single star formation law that can be applied to all star-forming regions and galaxies, across cosmic time.

preprint2013arXiv

Probability Distribution Functions OF 12CO(J = 1-0) Brightness and Integrated Intensity in M51: The PAWS View

We analyse the distribution of CO brightness temperature and integrated intensity in M51 at ~40 pc resolution using new CO data from the Plateau de Bure Arcsecond Whirlpool Survey (PAWS). We present probability distribution functions (PDFs) of the CO emission within the PAWS field, which covers the inner 11 x 7 kpc of M51. We find variations in the shape of CO PDFs within different M51 environments, and between M51 and M33 and the Large Magellanic Cloud (LMC). Globally, the PDFs for the inner disk of M51 can be represented by narrow lognormal functions that cover 1 to 2 orders of magnitude in CO brightness and integrated intensity. The PDFs for M33 and the LMC are narrower and peak at lower CO intensities. However, the CO PDFs for different dynamical environments within the PAWS field depart from the shape of the global distribution. The PDFs for the interarm region are approximately lognormal, but in the spiral arms and central region of M51, they exhibit diverse shapes with a significant excess of bright CO emission. The observed environmental dependence of the shape of the CO PDFs is qualitatively consistent with changes that would be expected if molecular gas in the spiral arms has a larger range of average densities, gas temperatures and velocity fluctuations, though further work is required to disentangle the importance of large-scale dynamical effects versus star formation feedback in regulating these properties. We show that the shape of the CO PDFs for different M51 environments is only weakly related to global properties of the CO emission, but is strongly correlated with some properties of the local giant molecular cloud (GMC) and young stellar cluster populations. For galaxies with strong spiral structure such as M51, our results indicate that galactic-scale dynamical processes play a significant role in the formation and evolution of GMCs and stellar clusters.(abridged)

preprint2013arXiv

The PdBI Arcsecond Whirlpool Survey (PAWS). I. A Cloud-Scale/Multi-Wavelength View of the Interstellar Medium in a Grand-Design Spiral Galaxy

The PdBI (Plateau de Bure Interferometer) Arcsecond Whirlpool Survey (PAWS) has mapped the molecular gas in the central ~9kpc of M51 in its 12CO(1-0) line emission at cloud-scale resolution of ~40pc using both IRAM telescopes. We utilize this dataset to quantitatively characterize the relation of molecular gas (or CO emission) to other tracers of the interstellar medium (ISM), star formation and stellar populations of varying ages. Using 2-dimensional maps, a polar cross-correlation technique and pixel-by-pixel diagrams, we find: (a) that (as expected) the distribution of the molecular gas can be linked to different components of the gravitational potential, (b) evidence for a physical link between CO line emission and radio continuum that seems not to be caused by massive stars, but rather depend on the gas density, (c) a close spatial relation between the PAH and molecular gas emission, but no predictive power of PAH emission for the molecular gas mass,(d) that the I-H color map is an excellent predictor of the distribution (and to a lesser degree the brightness) of CO emission, and (e) that the impact of massive (UV-intense) young star-forming regions on the bulk of the molecular gas in central ~9kpc can not be significant due to a complex spatial relation between molecular gas and star-forming regions that ranges from co-spatial to spatially offset to absent. The last point, in particular, highlights the importance of galactic environment -- and thus the underlying gravitational potential -- for the distribution of molecular gas and star formation.

preprint2012arXiv

The Resolved Radio--FIR Correlation in Nearby Galaxies with Herschel and Spitzer

We investigate the correlation between the far-infrared (FIR) and radio continuum emission from NGC6946 on spatial scales between 0.9 and 17 kpc. We use the Herschel PACS (70, 100, 160$μ$m) and SPIRE (250$μ$m) data from the KINGFISH project. Separating the free-free and synchrotron components of the radio continuum emission, we find that FIR is better correlated with the free-free than the synchrotron emission. Compared to a similar study in M33 and M31, we find that the scale dependence of the synchrotron--FIR correlation in NGC6946 is more similar to M31 than M33. The scale dependence of the synchrotron--FIR correlation can be explained by the turbulent-to-ordered magnetic field ratio or, equivalently, the diffusion length of the cosmic ray electrons in these galaxies.

preprint2011arXiv

The Magellanic Mopra Assessment (MAGMA). I. The Molecular Cloud Population of the Large Magellanic Cloud

We present the properties of an extensive sample of molecular clouds in the Large Magellanic Cloud (LMC) mapped at 11 pc resolution in the CO(1-0) line. We identify clouds as regions of connected CO emission, and find that the distributions of cloud sizes, fluxes and masses are sensitive to the choice of decomposition parameters. In all cases, however, the luminosity function of CO clouds is steeper than dN/dL \propto L^{-2}, suggesting that a substantial fraction of mass is in low-mass clouds. A correlation between size and linewidth, while apparent for the largest emission structures, breaks down when those structures are decomposed into smaller structures. We argue that the correlation between virial mass and CO luminosity is the result of comparing two covariant quantities, with the correlation appearing tighter on larger scales where a size-linewidth relation holds. The virial parameter (the ratio of a cloud's kinetic to self-gravitational energy) shows a wide range of values and exhibits no clear trends with the CO luminosity or the likelihood of hosting young stellar object (YSO) candidates, casting further doubt on the assumption of virialization for molecular clouds in the LMC. Higher CO luminosity increases the likelihood of a cloud harboring a YSO candidate, and more luminous YSOs are more likely to be coincident with detectable CO emission, confirming the close link between giant molecular clouds and massive star formation.

preprint2010arXiv

Dust/gas correlations from Herschel Observations

Spitzer and IRAS observations of the LMC suggest an excess of FIR emission with respect to the gas surface density traced by 12CO and HI 21 cm emission lines. This "FIR excess" is noticeable near molecular clouds in the LMC, and has usually been interpreted as the presence of a self-shielded H2 component not traced by CO molecular clouds' envelopes. Based on Herschel observations, we examine the correlation between gas and dust at higher resolution than previously achieved. We consider three additional causes for the FIR excess: X factor, FIR dust emissivity, and gas-to-dust ratio variations between the diffuse and dense phases of the ISM. We examine the structure of NT80 and NT71, two molecular clouds detected in the NANTEN 12CO survey of the LMC. Dust surface density maps are derived from the HERITAGE data. The gas phase is traced by MAGMA 12CO and ATCA HI 21 cm observations of the LMC. The dust emissivity, gas-to-dust ratio, and X factor required to match the dust and gas surface densities are derived, and their correlations with the dust surface density are examined. The dust surface density is spatially correlated with the atomic and molecular gas phases. The dust temperature is consistently lower in the dense phase of the ISM than in the diffuse phase. We confirm variations in the ratio of FIR emission to gas surface density derived from HI and CO observations. There is an excess of FIR emission, spatially correlated with regions of intermediate HI and dust surface densities (Av = 1-2), and little or no CO. While there is no significant trend in the dust emissivity or gas-to-dust ratio with dust surface density, the X factor is enhanced at Av = 1-2. We conclude that H2 envelopes not traced by CO and X factor variations close to the CO boundary are more likely to cause these deviations between FIR emission and gas surface density than gas-to-dust ratio or emissivity variations.

preprint2009arXiv

The Influence of Far-Ultraviolet Radiation on the Properties of Molecular Clouds in the 30 Dor Region of the Large Magellanic Cloud

We present a complete 12CO J = 1-0 map of the prominent molecular ridge in the Large Magellanic Cloud (LMC) obtained with the 22-m ATNF Mopra Telescope. The region stretches southward by ~2 degrees, (or 1.7kpc) from 30 Doradus, the most vigorous star-forming region in the Local Group. The location of this molecular ridge is unique insofar as it allows us to study the properties of molecular gas as a function of the ambient radiation field in a low-metallicity environment. We find that the physical properties of CO-emitting clumps within the molecular ridge do not vary with the strength of the far-ultraviolet (FUV) radiation field. Since the peak CO brightness of the clumps shows no correlation with the radiation field strength, the observed constant value for CO-to-H_2 conversion factor along the ridge seems to require an increase in the kinetic temperature of the molecular gas that is offset by a decrease in the angular filling factor of the CO emission. We find that the difference between the CO-to-H_2 conversion factor in the molecular ridge and the outer Milky Way is smaller than has been reported by previous studies of the CO emission: applying the same cloud identification and analysis methods to our CO observations of the LMC molecular ridge and CO data from the outer Galaxy survey by Dame et al. (2001), we find that the average CO-to-H_2 conversion factor in the molecular ridge is X_CO = (3.9+/-2.5)10^20 cm^-2 (K km s^-1)^-1, approximately twice the value that we determine for the outer Galaxy clouds. The mass spectrum and the scaling relations between the properties of the CO clumps in the molecular ridge are similar, but not identical, to those that have been established for Galactic molecular clouds.

preprint2009arXiv

Warm and Dense Molecular Gas in the N159 Region: 12CO J=4-3 and 13CO J=3-2 Observations with NANTEN2 and ASTE

New 12CO J=4-3 and 13CO J=3-2 observations of the N159 region in the Large Magellanic Cloud have been made. The 12CO J=4-3 distribution is separated into three clumps. These new measurements toward the three clumps are used in coupled calculations of molecular rotational excitation and line radiation transfer, along with other transitions of the 12CO as well as the isotope transitions of 13CO. The temperatures and densities are determined to be ~70-80K and ~3x10^3 cm-3 in N159W and N159E and ~30K and ~1.6x10^3 cm-3 in N159S. These results are compared with the star formation activity. The N159E clump is associated with embedded cluster(s) as observed at 24 micron and the derived high temperature is explained as due to the heating by these sources. The N159E clump is likely responsible for a dark lane in a large HII region by the dust extinction. The N159W clump is associated with embedded clusters mainly toward the eastern edge of the clump only. These clusters show offsets of 20"-40" from the 12CO J=4-3 peak and are probably responsible for heating indicated by the derived high temperature. The N159W clump exhibits no sign of star formation toward the 12CO J=4-3 peak position and its western region. We suggest that the N159W peak represents a pre-star-cluster core of ~105M_sol which deserves further detailed studies. Note that recent star formation took place between N159W and N159E as indicated by several star clusters and HII regions, while the natal molecular gas toward the stars have already been dissipated by the ionization and stellar winds of the OB stars. The N159S clump shows little sign of star formation as is consistent with the lower temperature and somewhat lower density. The N159S clump is also a candidate for future star formation.

preprint2008arXiv

The Molecular Ridge Close to 30 Doradus in the Large Magellanic Cloud

With the ATNF Mopra telescope we are performing a survey in the 12CO(1-0) line to map the molecular gas in the Large Magellanic Cloud (LMC). For some regions we also obtained interferometric maps of the high density gas tracers HCO+ and HCN with the Australia Telescope Compact Array (ATCA). Here we discuss the properties of the elongated molecular complex that stretches about 2 kpc southward from 30 Doradus. Our data suggests that the complex, which we refer to as the ``molecular ridge,'' is not a coherent feature but consists of many smaller clumps that share the same formation history. Likely molecular cloud formation triggers are shocks and shearing forces that are present in the surrounding south-eastern HI overdensity region, a region influenced by strong ram pressure and tidal forces. The molecular ridge is at the western edge of the the overdensity region where a bifurcated velocity structure transitions into a single disk velocity component. We find that the 12CO(1-0) and HI emission peaks in the molecular ridge are typically near each other but never coincide. A likely explanation is the conversion of warmer, low-opacity HI to colder, high-opacity HI from which H2 subsequently forms. On smaller scales, we find that very dense molecular gas, as traced by interferometric HCO+ and HCN maps, is associated with star formation along shocked filaments and with rims of expanding shell-like structures, both created by feedback from massive stars.