Search

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

The double red clump (RC) observed in color-magnitude diagrams of the Milky Way bulge is at the heart of the current debate on the structure and formation origin of the bulge. This feature can be explained by the difference between the two RCs either in distance ("X-shaped scenario") or in chemical composition ("multiple-population scenario"). Here we report our high-resolution spectroscopy for the RC and red giant branch stars in a high-latitude field (b ~ -8.5$\deg$) of the bulge. We find a difference in [Fe/H] between the stars in the bright and faint RC regimes, in the sense that the bright stars are enhanced in [Fe/H] with respect to the faint stars by 0.149 $\pm$ 0.036 dex. The stars on the bright RC are also enhanced in [Na/Fe] but appear to be depleted in [Al/Fe] and [O/Fe], although more observations are required to confirm the significance of these differences. Interestingly, these chemical patterns are similar to those observed among multiple stellar populations in the metal-rich bulge globular cluster Terzan 5. In addition, we find a number of Na-rich stars, which would corroborate the presence of multiple populations in the bulge. Our results support an origin of the double RC from dissolved globular clusters that harbor multiple stellar populations. Thus, our study suggests that a substantial fraction of the outer bulge stars would have originated from the assembly of such stellar systems in the early phase of the Milky Way formation., Comment: 17 pages, 15 figures, accepted for publication in ApJ

It is a major open question which physical processes stop the accretion of gas onto giant molecular clouds (GMCs) and limit the efficiency at which gas is converted into stars within these GMCs. While feedback from supernova explosions has been the popular feedback mechanism included in simulations of galaxy formation and evolution, `early' feedback mechanisms such as stellar winds, photoionisation and radiation pressure are expected to play an important role in dispersing the gas after the onset of star formation. These feedback processes typically take place on small scales ($\sim 10-100$ pc) and their effects have therefore been difficult to constrain in environments other than the Milky Way. We apply a novel statistical method to $\sim 1$" resolution maps of CO and Ha emission across a sample of nine nearby disc galaxies, in order to measure the time over which GMCs are dispersed by feedback from young, high-mass stars, as a function of the galactic environment. We find that GMCs are typically dispersed within $\sim$ 3 Myr after the emergence of unembedded high-mass stars, showing no significant trend with galactocentric radius. Comparison with analytical predictions demonstrates that, independently of the environment, early feedback mechanisms (particularly photoionisation and stellar winds) play a crucial role in dispersing GMCs and limiting their star formation efficiency in nearby galaxies. Finally, we show that the efficiency at which the energy injected by these early feedback mechanisms couples with the parent GMC is relatively low (a few tens of per cent), such that the vast majority of momentum and energy emitted by the young stellar populations escapes the parent GMC., Comment: 17 pages, 8 figures, 3 tables; submitted to MNRAS (October 23, 2020)

The spatial distribution of metals reflects, and can be used to constrain, the processes of chemical enrichment and mixing. Using PHANGS-MUSE optical integral field spectroscopy, we measure the gas phase oxygen abundances (metallicities) across 7,138 HII regions in a sample of eight nearby disc galaxies. In Paper I (Kreckel et al. 2019) we measure and report linear radial gradients in the metallicities of each galaxy, and qualitatively searched for azimuthal abundance variations. Here, we examine the two-dimensional variation in abundances once the radial gradient is subtracted, Delta(O/H), in order to quantify the homogeneity of the metal distribution and to measure the mixing scale over which HII region metallicities are correlated. We observe low (0.03--0.05 dex) scatter in Delta(O/H) globally in all galaxies, with significantly lower (0.02--0.03 dex) scatter on small (<600 pc) spatial scales. This is consistent with the measurement uncertainties, and implies the two-dimensional metallicity distribution is highly correlated on scales of <600 pc. We compute the two point correlation function for metals in the disc in order to quantify the scale lengths associated with the observed homogeneity. This mixing scale is observed to correlate better with the local gas velocity dispersion (of both cold and ionized gas) than with the star formation rate. Selecting only HII regions with enhanced abundances relative to a linear radial gradient, we do not observe increased homogeneity on small scales. This suggests that the observed homogeneity is driven by the mixing introducing material from large scales rather than by pollution from recent and on-going star formation., Comment: 17 pages, 14 figures. Accepted for publication in MNRAS

Using the PHANGS-ALMA CO (2-1) survey, we characterize molecular gas properties on ${\sim}$100 pc scales across 102,778 independent sightlines in 70 nearby galaxies. This yields the best synthetic view of molecular gas properties on cloud scales across the local star-forming galaxy population obtained to date. Consistent with previous studies, we observe a wide range of molecular gas surface densities (3.4 dex), velocity dispersions (1.7 dex), and turbulent pressures (6.5 dex) across the galaxies in our sample. Under simplifying assumptions about sub-resolution gas structure, the inferred virial parameters suggest that the kinetic energy of the molecular gas typically exceeds its self-gravitational binding energy at ${\sim}$100 pc scales by a modest factor (1.3 on average). We find that the cloud-scale surface density, velocity dispersion, and turbulent pressure (1) increase towards the inner parts of galaxies, (2) are exceptionally high in the centers of barred galaxies (where the gas also appears less gravitationally bound), and (3) are moderately higher in spiral arms than in inter-arm regions. The galaxy-wide averages of these gas properties also correlate with the integrated stellar mass, star formation rate, and offset from the star-forming main sequence of the host galaxies. These correlations persist even when we exclude regions with extraordinary gas properties in galaxy centers, which contribute significantly to the inter-galaxy variations. Our results provide key empirical constraints on the physical link between molecular cloud populations and their galactic environment., Comment: 9 pages + appendices, ApJL in press. Data tables available at https://www.canfar.net/storage/list/phangs/RELEASES/Sun_etal_2020b

Recent stellar evolution models for globular clusters (GCs) in multiple population paradigm suggest that horizontal-branch (HB) morphology and mean period of type ab RR Lyrae variables are mostly determined by He & CNO abundances and relative ages for subpopulations. These parameters are also provided by chemical evolution models constructed to reproduce the Na-O anti-correlation. Therefore, a consistency check is possible between the synthetic HB and chemical evolution models. Furthermore, by combining them, a better constraint might be attained for star formation history and chemical abundances of subpopulations in GCs. We find, from such efforts made for four GCs, M4, M5, M15, and M80, that consistent results can be obtained from these two independent studies. In our unified model, He and Na abundances gradually increase over the generation, and therefore, the various extensions observed in both HB morphology and Na-O chemical pattern depend on the presence of later generation stars after the second generation. It is schematically shown that this observed diversity, however, would not be naturally explained by the models requiring dilution. Further spectroscopic observations are required, for metal-poor GCs in particular, to obtain a more detailed constraint from this approach., Comment: Accepted for publication in ApJ

Recent investigations of the double red clump in the color-magnitude diagram of the Milky Way bulge cast serious doubts on the structure and formation origin of the outer bulge. Unlike previous interpretation based on an X-shaped bulge, stellar evolution models and CN-band observations have suggested that this feature is another manifestation of the multiple stellar population phenomenon observed in globular clusters (GCs). This new scenario requires a significant fraction of the outer bulge stars with chemical patterns uniquely observed in GCs. Here we show from homogeneous high-quality spectroscopic data that the red giant branch stars in the outer bulge (> 5.5 degrees from the Galactic center) are clearly divided into two groups according to Na abundance in the [Na/Fe] - [Fe/H] plane. The Na-rich stars are also enhanced in Al, while the differences in O and Mg are not observed between the two Na groups. The population ratio and the Na and Al differences between the two groups are also comparable with those observed in metal-rich GCs. The only plausible explanation for these chemical patterns and characteristics appears to be that the outer bulge was mostly assembled from disrupted proto-GCs in the early history of the Milky Way., Comment: Accepted for publication in ApJL

In order to investigate the origin of multiple stellar populations found in globular clusters (GCs) in the halo and bulge of the Milky Way, we have constructed chemical evolution models for their putative low-mass progenitors. In light of recent theoretical developments, we assume that supernova blast waves undergo blowout without expelling the pre-enriched ambient gas, while relatively slow winds of massive stars, together with the winds and ejecta from low to high mass asymptotic-giant-branch stars, are all locally retained in these less massive systems. Interestingly, we find that the observed Na-O anti-correlations in metal-poor GCs can be reproduced when multiple episodes of starburst and enrichment are allowed to continue in these subsystems. A specific form of star formation history with decreasing time intervals between the successive stellar generations, however, is required to obtain this result, which is in good agreement with the parameters obtained from synthetic horizontal-branch models. The "mass budget problem" is also much alleviated by our models without ad-hoc assumptions on star formation efficiency, initial mass function, and the preferential loss of first-generation stars. We also applied these models to investigate the origin of super-He-rich red clump stars in the metal-rich bulge suggested by Lee et al. (2015). We find that chemical enrichment by the winds of massive stars can naturally reproduce the required strong He enhancement in metal-rich subsystems. Our results further underscore that gas expulsion or retention is a key factor in understanding the multiple populations in GCs., Comment: Accepted for publication in ApJ

The two red clumps (RCs) observed in the color-magnitude diagram of the Milky Way bulge is widely accepted as evidence for an X-shaped structure originated from the bar instability. A drastically different interpretation has been suggested, however, based on the He-enhanced multiple stellar population phenomenon as is observed in globular clusters (GCs). Because these two scenarios imply very different pictures on the formation of the bulge and elliptical galaxies, understanding the origin of the double RC is of crucial importance. Here we report our discovery that the stars in the two RCs show a significant (> 5.3 {\sigma}) difference in CN-band strength, in stark contrast to that expected in the X-shaped bulge scenario. The difference in CN abundance and the population ratio between the two RCs are comparable to those observed in GCs between the first- and later generation stars. Since CN-strong stars trace a population with enhanced N, Na, and He abundances originated in GCs, this is direct evidence that the double RC is due to the multiple population phenomenon, and that a significant population of stars in the Milky Way bulge were assembled from disrupted proto-GCs. Our result also calls for the major revision of the 3D structure of the Milky Way bulge given that the current view is based on the previous interpretation of the double RC phenomenon., Comment: Accepted for publication in ApJL