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, $\Sigma_\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 $\Sigma_\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.

We analyze the 1D spatial power spectra of dust surface density and mid to far-infrared emission at $24-500\,\mu$m in the LMC, SMC, M31, and M33. By forward-modelling the point-spread-function (PSF) on the power spectrum, we find that nearly all power spectra have a single power-law and point source component. A broken power-law model is only favoured for the LMC 24 $\mu$m MIPS power spectrum and is due to intense dust heating in 30 Doradus. We also test for local power spectrum variations by splitting the LMC and SMC maps into $820$ pc boxes. We find significant variations in the power-law index with no strong evidence for breaks. The lack of a ubiquitous break suggests that the spatial power spectrum does not constrain the disc scale height. This contradicts claims of a break where the turbulent motion changes from 3D to 2D. The power spectrum indices in the LMC, SMC, and M31 are similar ($2.0-2.5$). M33 has a flatter power spectrum ($1.3$), similar to more distant spiral galaxies with a centrally-concentrated H$_{2}$ distribution. We compare the power spectra of HI, CO, and dust in M31 and M33 and find that HI power spectra are consistently flatter than CO power spectra. These results cast doubt on the idea that the spatial power spectrum traces large scale turbulent motion in nearby galaxies. Instead, we find that the spatial power spectrum is influenced by (1) the PSF on scales below $\sim3$ times the FWHM, (2) bright compact regions (30 Doradus), and (3) the global morphology of the tracer (an exponential CO disc).

We present an atlas of ultraviolet and infrared images of ~15,750 local (d < 50 Mpc) galaxies, as observed by NASA's WISE and GALEX missions. These maps have matched resolution (FWHM 7.5'' and 15''), matched astrometry, and a common procedure for background removal. We demonstrate that they agree well with resolved intensity measurements and integrated photometry from previous surveys. This atlas represents the first part of a program (the z=0 Multi-wavelength Galaxy Synthesis) to create a large, uniform database of resolved measurements of gas and dust in nearby galaxies. The images and associated catalogs are publicly available at the NASA/IPAC Infrared Science Archive. This atlas allows us estimate local and integrated star formation rates (SFRs) and stellar masses (M$_\star$) across the local galaxy population in a uniform way. In the appendix, we use the population synthesis fits of Salim et al. (2016, 2018) to calibrate integrated M$_\star$ and SFR estimators based on GALEX and WISE. Because they leverage an SDSS-base training set of >100,000 galaxies, these calibrations have high precision and allow us to rigorously compare local galaxies to Sloan Digital Sky Survey results. We provide these SFR and M$_\star$ estimates for all galaxies in our sample and show that our results yield a "main sequence" of star forming galaxies comparable to previous work. We also show the distribution of intensities from resolved galaxies in NUV-to-WISE1 vs. WISE1-to-WISE3 space, which captures much of the key physics accessed by these bands.

We present maps of the dust properties in the Small and Large Magellanic Clouds (SMC, LMC) from fitting Spitzer and Herschel observations with the \citet{DL07} dust model. We derive the abundance of the small carbonaceous grain (or polycyclic aromatic hydrocarbon; PAH) component. The global PAH fraction (q_pah, the fraction of the dust mass in the form of PAHs) is smaller in the SMC (1.0$^{+0.3}_{-0.3}$%) than in the LMC (3.3$^{+1.4}_{-1.3}$%). We measure the PAH fraction in different gas phases (H II regions, ionized gas outside of H II regions, molecular gas, and diffuse neutral gas). H II regions appear as distinctive holes in the spatial distribution of the PAH fraction. In both galaxies, the PAH fraction in the diffuse neutral medium is higher than in the ionized gas, but similar to the molecular gas. Even at equal radiation field intensity, the PAH fraction is lower in the ionized gas than in the diffuse neutral gas. We investigate the PAH life-cycle as a function of metallicity between the two galaxies. The PAH fraction in the diffuse neutral medium of the LMC is similar to that of the Milky Way ($\sim4.6$%), while it is significantly lower in the SMC. Plausible explanations for the higher PAH fraction in the diffuse neutral medium of the LMC compared to the SMC include: a more effective PAH production by fragmentation of large grains at higher metallicity, and/or the growth of PAHs in molecular gas.

We utilize archival far-infrared maps from the Herschel Space Observatory in four Local Group galaxies (Small and Large Magellanic Clouds, M31, and M33). We model their Spectral Energy Distribution (SED) from 100 to 500 $\mu$m using a single-temperature modified blackbody emission with a fixed emissivity index of $\beta = 1.8$. From the best-fit model, we derive the dust temperature, $T_{\rm d}$, and the dust mass surface density, $\Sigma_{\rm d}$, at 13 parsec resolution for SMC and LMC, and at 167 parsec resolution for all targets. This measurement allows us to build the distribution of dust mass and luminosity as functions of dust temperature and mass surface density. We compare those distribution functions among galaxies and between regions in a galaxy. We find that LMC has the highest mass-weighted average $T_{\rm d}$, while M31 and M33 have the lowest mass-weighted average $T_{\rm d}$. Within a galaxy, star forming regions have higher $T_{\rm d}$ and $\Sigma_{\rm d}$ relative to the overall distribution function, due to more intense heating by young stars and higher gas mass surface density. When we degrade the resolutions to mimic distant galaxies, the mass-weighted mean temperature gets warmer as the resolution gets coarser, meaning the temperature derived from unresolved observation is systematically higher than that in highly resolved observation. As an implication, the total dust mass is lower (underestimated) in coarser resolutions. This resolution-dependent effect is more prominent in clumpy star-forming galaxies (SMC, LMC, and M33), and less prominent in more quiescent massive spiral (M31).

The dust-to-metals ratio describes the fraction of the heavy elements contained in dust grains, and its variation provides key insights into the life cycle of dust. We measure the dust-to-metals ratio in M101, a nearby galaxy with a radial metallicity (Z) gradient spanning $\sim$1 dex. We fit the dust spectral energy distribution from 100 to 500 $\mu m$ with five variants of the modified blackbody dust emission model in which we vary the temperature distribution and how emissivity depends on wavelength. Among them, the model with a single temperature blackbody modified by a broken power-law emissivity gives the statistically best fit and physically most plausible results. Using these results, we show that the dust-to-gas ratio is proportional to $\rm Z^{1.7}$. This implies that the dust-to-metals ratio is not constant in M101, but decreases as a function of radius, equivalent to a lower fraction of metals trapped in dust at low metallicity (large radius). The dust-to-metals ratio in M101 remains at or above what would be predicted by the minimum depletion level of metals observed in the Milky Way. Our current knowledge of metallicity-dependent CO-to-H$_2$ conversion factor suggests that variations in the conversion factor cannot be responsible for the dust-to-metals ratio trends we observe. This change of dust-to-metals ratio is significantly correlated with molecular hydrogen fraction, which suggests that the accretion of gas phase metals onto existing dust grains could be a mechanism contributing to a variable dust-to-metals ratio.