The EDGE-CALIFA survey provides spatially resolved optical integral-field unit and CO spectroscopy for 125 galaxies selected from the Calar Alto Legacy Integral Field Area Survey (CALIFA) Data Release 3 sample. The Extragalactic Database for Galaxy Evolution (EDGE) presents the spatially resolved products of the survey as pixel tables that reduce the oversampling in the original images and facilitate comparison of pixels from different images. By joining these pixel tables to lower-dimensional tables that provide radial profiles, integrated spectra, or global properties, it is possible to investigate the dependence of local conditions on large-scale properties. The database is freely accessible and has been utilized in several publications. We illustrate the use of this database and highlight the effects of CO upper limits on the inferred slopes of the local scaling relations between the stellar mass, star formation rate (SFR), and H _2 surface densities. We find that the correlation between H _2 and SFR surface density is the tightest among the three relations.

A strong correlation between atmospheric pressure and molecular gas mass is found in central cluster galaxies and early-type galaxies. This trend and a similar trend with atmospheric gas density would naturally arise if the molecular clouds condensed from hot atmospheres. Limits on the ratio of molecular to atomic hydrogen in these systems exceed unity. The data are consistent with ambient pressure being a significant factor in the rapid conversion of atomic hydrogen into molecules as found in normal spiral galaxies.

We present 13CO(1-0) observations for the EDGE-CALIFA survey, which is a mapping survey of 126 nearby galaxies at a typical spatial resolution of 1.5 kpc. Using detected 12CO(1-0) emission as a prior, we detect 13CO(1-0) in 41 galaxies via integrated line flux over the entire galaxy, and in 30 galaxies via integrated line intensity in resolved synthesized beams. Incorporating our CO observations and optical IFU spectroscopy, we perform a systematic comparison between the line ratio R12/13 and the properties of the stars and ionized gas. Higher R12/13 values are found in interacting galaxies than in non-interacting galaxies. The global R12/13 slightly increases with infrared color F60/F100, but appears insensitive to other host galaxy properties such as morphology, stellar mass, or galaxy size. We also present annulus-averaged R12/13 profiles for our sample up to a galactocentric radius of 0.4r25 ( 6 kpc), taking into account the 13CO(1-0) non-detections by spectral stacking. The radial profiles of R12/13 are quite flat across our sample. Within galactocentric distances of 0.2r25, azimuthally-averaged R12/13 increases with star formation rate. However, the Spearman rank correlation tests show the azimuthally-averaged R12/13 does not strongly correlate with any other gas or stellar properties in general, especially beyond 0.2r25 from the galaxy centers. Our findings suggest that in the complex environments in galaxy disks, R12/13 is not a sensitive tracer for ISM properties. Dynamical disturbances, like galaxy interactions or the presence of a bar, also have an overall impact on R12/13, which further complicate the interpretations of R12/13 variations.

We present high spatial resolution (12pc) Atacama Large Millimeter/sub-millimeter Array CO(J=3-2) observations of the nearby lenticular galaxy NGC4429. We identify 217 giant molecular clouds within the 450pc radius molecular gas disc. The clouds generally have smaller sizes and masses but higher surface densities and observed linewidths than those of Milky Way disc clouds. An unusually steep size - line width relation and large cloud internal velocity gradients (0.05 - 0.91 km s^-1 pc^-1) and observed Virial parameters (alpha_obs,vir = 4.0) are found, that appear due to internal rotation driven by the background galactic gravitational potential. Removing this rotation, an internal Virial equilibrium appears to be established between the self-gravitational (Usg) and turbulent kinetic (Eturb) energies of each cloud, i.e. alpha_sg,vir=Usg/Eturb = 1.3. However, to properly account for both self and external gravity (shear and tidal forces), we formulate a modified Virial theorem and define an effective Virial parameter alpha_eff,vir = alpha_sg,vir + Usg/Eext (and associated effective velocity dispersion). The NGC4429 clouds then appear to be in a critical state in which the self-gravitational energy and the contribution of external gravity to the cloud's energy budget (Eext) are approximately equal, i.e. Eext/Usg 1. As such, alpha_eff,vir = 2.2 and most clouds are not virialised but remain marginally gravitationally bound. We show this is consistent with the clouds having sizes similar to their tidal radii and being generally radially elongated. External gravity is thus as important as self-gravity to regulate the clouds of NGC4429.

We measure the star formation rate (SFR) per unit gas mass and the star formation efficiency (SFEgas for total gas, SFEmol for the molecular gas) in 81 nearby galaxies selected from the EDGE-CALIFA survey, using 12CO (J = 1–0) and optical IFU data. For this analysis we stack CO spectra coherently by using the velocities of Hα detections to detect fainter CO emission out to galactocentric radii r gal ∼ 1.2r 25 (∼3R e) and include the effects of metallicity and high surface densities in the CO-to-H2 conversion. We determine the scale lengths for the molecular and stellar components, finding a close to 1:1 relation between them. This result indicates that CO emission and star formation activity are closely related. We examine the radial dependence of SFEgas on physical parameters such as galactocentric radius, stellar surface density Σ⋆, dynamical equilibrium pressure P DE, orbital timescale τ orb, and the Toomre Q stability parameter (including star and gas Q star+gas). We observe a generally smooth, continuous exponential decline in the SFEgas with r gal. The SFEgas dependence on most of the physical quantities appears to be well described by a power law. Our results also show a flattening in the SFEgas–τ orb relation at log[τorb]∼7.9–8.1 and a morphological dependence of the SFEgas per orbital time, which may reflect star formation quenching due to the presence of a bulge component. We do not find a clear correlation between SFEgas and Q star+gas.

Understanding how galaxies cease to form stars represents an outstanding challenge for galaxy evolution theories. This process of “star formation quenching” has been related to various causes, including active galactic nuclei activity, the influence of large-scale dynamics, and the environment in which galaxies live. In this paper, we present the first results from a follow-up of CALIFA survey galaxies with observations of molecular gas obtained with the APEX telescope. Together with the EDGE-CARMA observations, we collected 12CO observations that cover approximately one effective radius in 472 CALIFA galaxies. We observe that the deficit of galaxy star formation with respect to the star formation main sequence (SFMS) increases with the absence of molecular gas and with a reduced efficiency of conversion of molecular gas into stars, which is in line with the results of other integrated studies. However, by dividing the sample into galaxies dominated by star formation and galaxies quenched in their centres (as indicated by the average value of the Hα equivalent width), we find that this deficit increases sharply once a certain level of gas consumption is reached, indicating that different mechanisms drive separation from the SFMS in star-forming and quenched galaxies. Our results indicate that differences in the amount of molecular gas at a fixed stellar mass are the primary drivers for the dispersion in the SFMS, and the most likely explanation for the start of star formation quenching. However, once a galaxy is quenched, changes in star formation efficiency drive how much a retired galaxy differs in its star formation rate from star-forming ones of similar masses. In other words, once a paucity of molecular gas has significantly reduced star formation, changes in the star formation efficiency are what drives a galaxy deeper into the red cloud, hence retiring it.

We present an empirical relation between the cold gas surface density (Sigma(gas)) and the optical extinction (A(V)) in a sample of 103 galaxies from the Extragalactic Database for Galaxy Evolution (EDGE) survey. This survey provides CARMA interferometric CO observations for 126 galaxies included in the Calar Alto Legacy Integral Field Area (CALIFA) survey. The matched, spatially resolved nature of these data sets allows us to derive the Sigma(gas)-A(V) relation on global, radial, and kpc (spaxel) scales. We determine A(V) from the Balmer decrement (H alpha/H beta). We find that the best fit for this relation is Sigma(gas) (M-circle dot pc(-2)) similar to 26 x A(V) (mag), and that it does not depend on the spatial scale used for the fit. However, the scatter in the fits increases as we probe smaller spatial scales, reflecting the complex relative spatial distributions of stars, gas, and dust. We investigate the Sigma(gas)/A(V) ratio on radial and spaxel scales as a function of EW(H alpha). We find that at larger values of EW(H alpha) (i.e. actively star-forming regions) this ratio tends to converge to twice the value expected for a foreground dust screen geometry (similar to 30 M-circle dot pc(-2) mag(-1)). On radial scales, we do not find a significant relation between the Sigma(gas)/A(V) ratio and the ionized gas metallicity. We contrast our estimates of Sigma(gas) using A(V) with compilations in the literature of the gas fraction on global and radial scales as well as with well-known scaling relations such as the radial star formation law and the Sigma(gas)-Sigma(*) relation. These tests show that optical extinction is a reliable proxy for estimating Sigma(gas) in the absence of direct sub/millimeter observations of the cold gas.

We present a multilinear analysis to determine the significant predictors of star formation in galaxies using the combined EDGE-CALIFA sample of galaxies. We analyze 1845 kpc-scale lines of sight across 39 galaxies with molecular line emission measurements from EDGE combined with optical IFU data drawn from CALIFA. We use the Least Absolute Shrinkage and Selection Operator (LASSO) to identify significant factors in predicting star formation rates. We find that the local star formation rate surface density is increased by higher molecular gas surface densities and stellar surface densities. In contrast, we see lower star formation rates in systems with older stellar populations, higher gas- and stellar-phase metallicities and larger galaxy masses. We also find a significant increase in star formation rate with galactocentric radius normalized by the disk scale length, which suggests additional parameters regulating star formation rate not explored in this study.

Dark-matter-only simulations predict that dark matter halos have steep, cuspy inner density profiles, while observations of dwarf galaxies find a range of inner slopes that are often much shallower. There is debate whether this discrepancy can be explained by baryonic feedback or if it may require modified dark matter models. In Paper I of this series, we obtained high-resolution integral field Hα observations for 26 dwarf galaxies with M* = 108.1−109.7 . We derived rotation curves from our observations, which we use here to construct mass models. We model the total mass distribution as the sum of a generalized Navarro–Frenk–White (NFW) dark matter halo and the stellar and gaseous components. Our analysis of the slope of the dark matter density profile focuses on the inner 300–800 pc, chosen based on the resolution of our data and the region resolved by modern hydrodynamical simulations. The inner slope measured using ionized and molecular gas tracers is consistent, and it is additionally robust to the choice of stellar mass-to-light ratio. We find a range of dark matter profiles, including both cored and cuspy slopes, with an average of , shallower than the NFW profile, but steeper than those typically observed for lower-mass galaxies with M* ∼ 107.5 . Simulations that reproduce the observed slopes in those lower-mass galaxies also produce slopes that are too shallow for galaxies in our mass range. We therefore conclude that supernova feedback models do not yet provide a fully satisfactory explanation for the observed trend in dark matter slopes.

Dark-matter-only simulations predict that dark matter halos have cusp-like inner density profiles, while observations of low-mass galaxies have found a range of inner slopes that are typically much shallower. It is still not well established whether this discrepancy can be explained by baryonic feedback or if it may require modified dark matter models. To better understand the diversity of dark matter profiles in dwarf galaxies, we undertook a survey of 26 low-mass galaxies (logM(*)/M-circle dot = 8.4-9.8, v(max) = 50-140 km s(-1)) within 30 Mpc using the Palomar Cosmic Web Imager, which is among the largest integral field spectroscopic surveys of its type. In this paper, we derive H alpha velocity fields for the full sample with a typical spatial resolution of similar to 160 pc. We extract rotation curves and verify their robustness to several choices in the analysis. We present a method for improving the velocity precision obtained from image slicing spectrographs using narrowband H alpha images. For 11 galaxies, we compare the H alpha velocity fields to CO kinematics measured using CARMA, finding the maps to be in good agreement. The standard deviation of the difference is typically similar to 7 km s(-1), comparable to the level of turbulence in the interstellar medium, showing that the two tracers have substantially the same bulk kinematics. In a companion paper, we will use the rotation curves produced here to construct mass models of the galaxies and determine their dark matter density profiles.

We utilize the multi-wavelength data of M33 to study the origin of turbulence in its interstellar medium. We find that the H i turbulent energy surface density inside 8 kpc is ∼1–3 × 10 46 erg pc −2 , and has no strong dependence on galactocentric radius because of the lack of variation in H i surface density and H i velocity dispersion. Then, we consider the energies injected by supernovae (SNe), the magneto-rotational instability (MRI), and the gravity-driven turbulence from accreted materials as the sources of turbulent energy. For a constant dissipation time of turbulence, the SNe energy can maintain turbulence inside ∼4 kpc radius (equivalent to ∼0.5 R 25 ), while the MRI energy is always smaller than the turbulent energy within 8 kpc radius. However, when we let the dissipation time to be equal to the crossing time of turbulence across the H i scale height, the SNe energy is enough to maintain turbulence out to 7 kpc radius, and the sum of SNe and MRI energies is able to maintain turbulence out to 8 kpc radius. Due to lack of constraint in the mass accretion rate through the disk of M33, we cannot rule out the accretion driven turbulence as a possible source of energy. Furthermore, by resolving individual giant molecular clouds in M33, we also show that the SNe energy can maintain turbulence within individual molecular clouds with ∼1% of coupling efficiency. This result strengthens the proposition that stellar feedback is an important source of energy to maintain turbulence in nearby galaxies.

We investigate the prevalence, properties, and kinematics of extraplanar diffuse ionized gas (eDIG) in a sample of 25 edge-on galaxies selected from the CALIFA survey. We measure ionized gas scale heights from H alpha and find that 90% have measurable scale heights with a median of 0.8(-0.4)(+0.7) kpc. From the Ha kinematics, we find that 60% of galaxies show a decrease in the rotation velocity as a function of height above the midplane. This lag is characteristic of eDIG, and we measure a median lag of 21 km s(-1) kpc(-1), which is comparable to lags measured in the literature. We also investigate variations in the lag with radius. H I lags have been reported to systematically decrease with galactocentric radius. We find both increasing and decreasing ionized gas lags with radius, as well as a large number of galaxies consistent with no radial lag variation, and investigate these results in the context of internal and external origins for the lagging ionized gas. We confirm that the [S II]/H alpha and [N II]/H alpha line ratios increase with height above the midplane as is characteristic of eDIG. The ionization of the eDIG is dominated by star-forming complexes (leaky H II regions). We conclude that the lagging ionized gas is turbulent ejected gas likely resulting from star formation activity in the disk as opposed to gas in the stellar thick disk or bulge. This is further evidence for the eDIG being a product of stellar feedback and for the pervasiveness of this WIM-like phase in many local star-forming galaxies.

We employ the earlier published proper motions of the newly discovered Antlia 2 dwarf galaxy derived from Gaia data to calculate its orbital distribution in the cosmologically recent past. Using these observationally motivated orbits, we calculate the effect of the Antlia 2 dwarf galaxy on the outer H i disk of the Milky Way, using both test particle and smoothed particle hydrodynamics simulations. We find that orbits with low pericenters, ?10 kpc, produce disturbances that match the observed outer H i disk perturbations. We have independently recalculated the proper motion of the Antlia 2 dwarf from Gaia data and found a proper motion of (?(?)cos?, ?(?))=(?0.068, 0.032)(0.023, ?0.031) mas yr(?1), which agrees with results from Torrealba et al. within the errors, but gives lower mean pericenters, e.g., ?15 kpc for our fiducial model of the Milky Way. We also show that the Sagittarius dwarf galaxy interaction does not match the observed perturbations in the outer gas disk. Thus, Antlia 2 may be the driver of the observed large perturbations in the outer gas disk of the Galaxy. The current location of the Antlia 2 dwarf galaxy closely matches that predicted by an earlier dynamical analysis of the dwarf galaxy that drove ripples in the outer Galaxy, and, in particular, its orbit is nearly coplanar to the Galactic disk. If the Antlia 2 dwarf galaxy is responsible for the perturbations in the outer Galactic disk, it would have a specific range of proper motions that we predict here; this can be tested soon with Gaia DR-3 and Gaia DR-4 data.

We present a kpc-scale analysis of the relationship between the molecular depletion time (tau(mol)(dep)) and the orbital time (tau(orb)) across the field of 39 face-on local galaxies, selected from the EDGE-CALIFA sample. We find that, on average, 5 per cent of the available molecular gas is converted into stars per orbital time, or tau(mol)(dep) similar to 20 tau(orb). The resolved relation shows a scatter of similar to 0.5 dex. The scatter is ascribable to galaxies of different morphologies that follow different tau(mol)(dep)-tau(orb) relations which decrease in steepness from early-to late types. The morphologies appear to be linked with the star formation rate surface density, the molecular depletion time, and the orbital time, but they do not correlate with the molecular gas content of the galaxies in our sample. We speculate that in our molecular gas rich, early-type galaxies, the morphological quenching (in particular the disc stabilization via shear), rather than the absence of molecular gas, is the main factor responsible for their current inefficient star formation.

Deriving circular velocities of galaxies from stellar kinematics can provide an estimate of their total dynamical mass, provided a contribution from the velocity dispersion of the stars is taken into account. Molecular gas (e.g. CO), on the other hand, is a dynamically cold tracer and hence acts as an independent circular velocity estimate without needing such a correction. In this paper, we test the underlying assumptions of three commonly used dynamical models, deriving circular velocities from stellar kinematics of 54 galaxies (S0-Sd) that have observations of both stellar kinematics from the Calar Alto Legacy Integral Field Area (CALIFA) survey, and CO kinematics from the Extragalactic Database for Galaxy Evolution (EDGE) survey. We test the asymmetric drift correction (ADC) method, as well as Jeans, and Schwarzschild models. The three methods each reproduce the CO circular velocity at 1R(e) to within 10 per cent. All three methods show larger scatter (up to 20 per cent) in the inner regions (R < 0.4R(e)) that may be due to an increasingly spherical mass distribution (which is not captured by the thin disc assumption in ADC), or non-constant stellar M/L ratios (for both the JAM and Schwarzschild models). This homogeneous analysis of stellar and gaseous kinematics validates that all three models can recover M-dyn at 1R(e) to better than 20 per cent, but users should be mindful of scatter in the inner regions where some assumptions may break down.

We present a comparative study of molecular and ionized gas kinematics in nearby galaxies. These results are based on observations from the EDGE survey, which measured spatially resolved 12CO(J = 1–0) in 126 nearby galaxies. Every galaxy in EDGE has corresponding resolved ionized gas measurements from CALIFA. Using a sub-sample of 17 rotation-dominated, star-forming galaxies where precise molecular gas rotation curves could be extracted, we derive CO and Hα rotation curves using the same geometric parameters out to ≳1 Re. We find that ∼75% of our sample galaxies have smaller ionized gas rotation velocities than the molecular gas in the outer part of the rotation curve. In no case is the molecular gas rotation velocity measurably lower than that of the ionized gas. We suggest that the lower ionized gas rotation velocity can be attributed to a significant contribution from extraplanar diffuse ionized gas in a thick, turbulence-supported disk. Using observations of the Hγ transition, also available from CALIFA, we measure ionized gas velocity dispersions and find that these galaxies have sufficiently large velocity dispersions to support a thick ionized gas disk. Kinematic simulations show that a thick disk with a vertical rotation velocity gradient can reproduce the observed differences between the CO and Hα rotation velocities. Observed line ratios tracing diffuse ionized gas are elevated compared to typical values in the midplane of the Milky Way. In galaxies affected by this phenomenon, dynamical masses measured using ionized gas rotation curves will be systematically underestimated.