91̽��

Abstract:

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(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 Very Large Array 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, 91̽��ed 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.

Abstract:

In this paper we examine the factors that shape the distribution of molecular gas surface densities on the 150 pc scale across 67 morphologically diverse star-forming galaxies in the PHANGS-ALMA CO (2-1) survey. Dividing each galaxy into radial bins, we measure molecular gas surface density contrasts, defined here as the ratio between a fixed high percentile of the CO distribution and a fixed reference level in each bin. This reference level captures the level of the faint CO floor that extends between bright filamentary features, while the intensity level of the higher percentile probes the structures visually associated with bright, dense interstellar medium features like spiral arms, bars, and filaments. We compare these contrasts to matched percentile-based measurements of the 3.6 μm emission measured using Spitzer/IRAC imaging, which trace the underlying stellar mass density. We find that the logarithms of CO contrasts on 150 pc scales are 3-4 times larger than, and positively correlated with, the logarithms of 3.6 μm contrasts probing smooth nonaxisymmetric stellar bar and spiral structures. The correlation appears steeper than linear, consistent with the compression of gas as it flows supersonically in response to large-scale stellar structures, even in the presence of weak or flocculent spiral arms. Stellar dynamical features appear to play an important role in setting the cloud-scale gas density in our galaxies, with gas self-gravity perhaps playing a weaker role in setting the 150 pc scale distribution of gas densities.

Abstract:

We present a comparative study of four physical dust models and two single-temperature modified blackbody models by fitting them to the resolved WISE, Spitzer, and Herschel photometry of M101 (NGC 5457). Using identical data and a grid-based fitting technique, we compare the resulting dust and radiation field properties derived from the models. We find that the dust mass yielded by the different models can vary by up to a factor of 3 (factor of 1.4 between physical models only), although the fits have similar quality. Despite differences in their definition of the carriers of the mid-IR aromatic features, all physical models show the same spatial variations for the abundance of that grain population. Using the well-determined metallicity gradient in M101 and resolved gas maps, we calculate an approximate upper limit on the dust mass as a function of radius. All physical dust models are found to exceed this maximum estimate over some range of galactocentric radii. We show that renormalizing the models to match the same Milky Way high-latitude cirrus spectrum and abundance constraints can reduce the dust mass differences between models and bring the total dust mass below the maximum estimate at all radii.

Abstract:

We apply the Tremaine-Weinberg method to 19 nearby galaxies using stellar mass surface densities and velocities derived from the PHANGS-MUSE survey, to calculate (primarily bar) pattern speeds (ΩP). After quality checks, we find that around half (10) of these stellar-mass-based measurements are reliable. For those galaxies, we find good agreement between our results and previously published pattern speeds, and we use rotation curves to calculate major resonance locations (corotation radii and Lindblad resonances). We also compare these stellar-mass-derived pattern speeds with Hα (from MUSE) and CO(J = 2 - 1) emission from the PHANGS-ALMA survey. We find that in the case of these clumpy interstellar medium (ISM) tracers, this method erroneously gives a signal that is simply the angular frequency at a representative radius set by the distribution of these clumps (Ωclump), and that this Ωclump is significantly different from ΩP (∼20% in the case of Hα, and ∼50% in the case of CO). Thus, we conclude that it is inadvisable to use "pattern speeds"derived from ISM kinematics. Finally, we compare our derived pattern speeds and corotation radii, along with bar properties, to the global parameters of these galaxies. Consistent with previous studies, we find that galaxies with a later Hubble type have a larger ratio of corotation radius to bar length, more molecular-gas-rich galaxies have higher ΩP, and more bulge-dominated galaxies have lower ΩP. Unlike earlier works, however, there are no clear trends between the bar strength and ΩP, nor between the total stellar mass surface density and the pattern speed.

Abstract:

Supermassive black hole (SMBH) masses can be measured by resolving the dynamical influences of the SMBHs on tracers of the central potentials. Modern long-baseline interferometers have enabled the use of molecular gas as such a tracer. We present here Atacama Large Millimeter/submillimeter Array observations of the elliptical galaxy NGC 7052 at 0′′.11 (⁠37pc) resolution in the 12CO(2-1) line and 1.3 mm continuum emission. This resolution is sufficient to resolve the region in which the potential is dominated by the SMBH. We forward model these observations, using a multi-Gaussian expansion of a Hubble Space Telescope F814W image and a spatially constant mass-to-light ratio to model the stellar mass distribution. We infer an SMBH mass of 2.5±0.3×109M⊙ and a stellar I-band mass-to-light ratio of 4.6±0.2M⊙/L⊙,I (3σ confidence intervals). This SMBH mass is significantly larger than that derived using ionized gas kinematics, which however appears significantly more kinematically disturbed than the molecular gas. We also show that a central molecular gas deficit is likely to be the result of tidal disruption of molecular gas clouds due to the strong gradient in the central gravitational potential.