Issue |
A&A
Volume 666, October 2022
|
|
---|---|---|
Article Number | A175 | |
Number of page(s) | 23 | |
Section | Extragalactic astronomy | |
DOI | https://doi.org/10.1051/0004-6361/202244306 | |
Published online | 24 October 2022 |
The EDGE-CALIFA survey: The role of spiral arms and bars in driving central molecular gas concentrations
1
Max-Planck-Institut für Radioastronomie, Auf dem Hügel 69, 53121 Bonn, Germany
e-mail: syu@mpifr-bonn.mpg.de
2
Department of Astronomy, University of Maryland, College Park, MD 20742, USA
3
Department of Astronomy, University of Illinois, Urbana, IL 62801, USA
4
Steward Observatory, University of Arizona, Tucson, AZ 85721, USA
5
Instituto de Astronomía, Universidad Nacional Autónoma de México, A.P. 70-264, 04510 México, D.F., Mexico
6
Kavli Institute for Astronomy and Astrophysics, Peking University, Beijing 100871, PR China
7
Departamento de Astronomía, Universidad de Chile, Casilla 36-D, Santiago, Chile
Received:
20
June
2022
Accepted:
29
August
2022
Shocks and torques produced by non-axisymmetric structures such as spiral arms and bars may transport gas to galaxy central regions. We test this hypothesis by studying the dependence of the concentration of CO luminosity (CCO) and molecular gas (Cmol) and the star formation rate (CSFR) in the central ∼2 kpc on the strength of non-axisymmetric disk structure using a sample of 57 disk galaxies selected from the EDGE-CALIFA survey. The Cmol is calculated using a CO-to-H2 conversion factor that decreases with higher metallicity and higher stellar surface density. We find that Cmol is systematically 0.22 dex lower than CCO. We confirm that high Cmol and strong non-axisymmetric disk structure are more common in barred galaxies than in unbarred galaxies. However, we find that spiral arms also increase Cmol. We show that there is a good correlation between Cmol and the strength of non-axisymmetric structure (which can be due to a bar, spiral arms, or both). This suggests that the stronger the bars and spirals, the more efficient the galaxy is at transporting cold gas to its center. Despite the small subsample size, the Cmol of the four Seyferts are not significantly reduced compared to inactive galaxies of similar disk structure, implying that the active galactic nucleus feedback in Seyferts may not notably affect the molecular gas distribution in the central ∼2 kpc. We find that CSFR tightly correlates with Cmol in both unbarred and barred galaxies. Likewise, elevated CSFR is found in galaxies with strong disk structure. Our results suggest that the disk structure, either spirals or bars, can transport gas to the central regions, with higher inflow rates corresponding to stronger structure, and consequently boost central star formation. Both spirals and bars play, therefore, an essential role in the secular evolution of disk galaxies.
Key words: galaxies: evolution / galaxies: nuclei / galaxies: spiral / galaxies: star formation / galaxies: starburst / ISM: molecules
© S.-Y. Yu et al. 2022
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This article is published in open access under the Subscribe-to-Open model.
Open Access funding provided by Max Planck Society.
1. Introduction
Molecular gas is a key material for star formation. Observations of molecular gas provide an important tool for exploring how structures or processes influence galaxy evolution. Molecular gas is typically traced through CO observations (e.g., Sakamoto et al. 1999b; Sheth et al. 2005; Kennicutt et al. 2007; Bigiel et al. 2008; Leroy et al. 2009, 2021; Wilson et al. 2012; Bolatto et al. 2017; Sorai et al. 2019). Spatially resolved CO data furnish essential information about the global gaseous properties of galaxies, such as the radial distribution of molecular gas.
Spiral arms and bars are the most common features in disk galaxies. In the local universe, nearly 60% of disk galaxies have a bar (e.g., Aguerri et al. 2009; Díaz-García et al. 2016). Models and simulations show that the non-axisymmetric bar gravitational potential drives gas flow toward the galaxy central region along the bar dust lane (e.g., Athanassoula 1992a,b; Regan et al. 1999; Sheth et al. 2000, 2002; Regan & Teuben 2004; Kim et al. 2012; Combes et al. 2014; Fragkoudi et al. 2016; Tress et al. 2020). In a galaxy with a well-defined inner Lindblad resonance (ILR), the gas transported by the bar accumulates at the ILR at a distance of a few kiloparsecs from the center (Athanassoula 1984; Jenkins & Binney 1994); later, the response of gas to forcing by the bar can give rise to a nuclear spiral or bar pattern in the gas, which can transport the gas to the proximity of the central black hole (Englmaier & Shlosman 2000, 2004). Observations have confirmed bar-driven gas transport (e.g., Mundell & Shone 1999; Combes 2003; Zurita et al. 2004; Fathi et al. 2006; Jogee 2006; Haan et al. 2009; Querejeta et al. 2016). The rates of gas inflow range from 0.01 to 50 M⊙ yr−1 (Regan et al. 1997; Haan et al. 2009; Querejeta et al. 2016). If the gas accumulation in the center is faster than the gas consumption, the inflow of gas will result in a centrally concentrated distribution of gas. Based on 20 galaxies from the Nobeyama Radio Observatory-Owens Valley Radio Observatory (NRO-OVRO) survey (Sakamoto et al. 1999b), Sakamoto et al. (1999a) found that the molecular gas concentrations in the central 1 kpc are systematically higher in barred spirals than in unbarred spirals. Consistently, with a larger and more diverse sample of 44 galaxies, the Berkeley-Illinois-Maryland Association (BIMA) CO Survey of Nearby Galaxies (SONG) confirmed that high central gas concentrations are more common in barred galaxies (Sheth et al. 2005). Further, the central gas concentrations are higher in galaxies hosting a stronger bar (Kuno et al. 2007), suggesting that stronger bars drive more gas inflow (Regan & Teuben 2004; Hopkins & Quataert 2011; Kim et al. 2012). Komugi et al. (2008) showed that early-type galaxies harbor larger central concentrations, a trend they attribute to the effect of bulges.
The gas inflow driven by the bar leads to an increase in gas density at the center, which enhances the central star formation activity (Ho et al. 1997; Martinet & Friedli 1997; Aguerri 1999; Sheth et al. 2005; Regan et al. 2006; Schinnerer et al. 2006; Wang et al. 2012; Combes et al. 2014; Zhou et al. 2015; Lin et al. 2017; Chown et al. 2019; Lin et al. 2020; Wang et al. 2020; Tress et al. 2020; Sormani et al. 2020). Stronger bars tend to have higher levels of central enhancements of star formation rates (SFRs; Zhou et al. 2015; Lin et al. 2017). Nevertheless, some strongly barred galaxies could have normal or suppressed central star formation (Martinet & Friedli 1997; Wang et al. 2012, 2020; Consolandi et al. 2017; Díaz-García et al. 2020). In these cases, the bar may have transported gas to the center, where the gas has been consumed by star formation.
Since high central gas concentrations are more common in barred galaxies, the gas inflow driven by spiral arms should typically be weaker compared to the bar effect. A signature of spiral-driven gas inflow has been detected. Regan et al. (2006) showed that there are two unbarred galaxies among the known six disk galaxies that have central excess in the 8 μm and CO emission above the inward extrapolation of an exponential disk; hence, their molecular gas content is centrally concentrated. In Sheth et al. (2005), a few unbarred spirals show relatively high gas concentrations, although none of them reach the very high values seen in some barred galaxies. Recently, Yu et al. (2022) analyzed central star formation in 2779 nearby unbarred star-forming (SF) disk galaxies; they find higher central SFRs in galaxies with strong spiral arms, implying that strong spiral arms may transport gas to the center. The strength of spiral arms, analogous to the role of bar strength, may be a key factor in studying the impact of spiral arms. This point of view is supported by models and simulations, which have shown that the non-axisymmetric spiral potential provides an efficient mechanism for transferring angular momentum, causing a radial inflow of gas. Firstly, the quasi-static density waves (Lin & Shu 1964) predict that spiral potential generated by old stars induces large-scale galactic shocks on a gas cloud as the gas cloud crosses the arm (Roberts 1969). The large-scale galactic shocks dissipate angular momentum, causing the gas cloud in orbital motions to move radially inward inside the corotation resonance (Kalnajs 1972; Roberts 1972; Lubow et al. 1986; Kim & Kim 2014; Baba et al. 2016; Kim et al. 2020). Secondly, gravitational torque of the non-axisymmetric spiral potential generated by old stars drives gas inflow (Kim & Kim 2014; Baba et al. 2016; Kim et al. 2020). Thirdly, gravitational torque of the gaseous component has an additional minor contribution to gas inflow (Kim & Kim 2014). Kim & Kim (2014) find that the rate of gas inflow to the central region, caused by a combination of the above three processes driven by spiral arms, is higher in galaxy models with stronger and more slowly rotating arms.
The dependence of the gas inflow rate on the spiral arms strength suggests that the stronger the spiral arms, the higher the molecular gas concentration. If this is true, and given the known correlation between bar strength and gas concentration (Kuno et al. 2007), the impact of spirals and bars may be uniformly described by the strength of the non-axisymmetric disk structure. The arms and bars jointly influence the radial distribution of molecular gas. In barred galaxies, the arms first transport gas to a radial extent within the bar and the bar successively drives the gas toward the center.
A constant CO-to-H2 conversion factor (αCO) has usually been adopted to convert CO integrated intensities to molecular gas mass surface densities and then derive molecular gas concentrations (Sakamoto et al. 1999a; Sheth et al. 2005; Kuno et al. 2007). However, studies have shown that the αCO depends on the physical properties of the environment where the gas clouds are embedded (e.g., Narayanan et al. 2012; Bolatto et al. 2013). Hydrodynamical simulations suggest that the αCO could vary by orders of magnitude in different environments (Feldmann et al. 2012; Gong et al. 2020). The αCO tends to be lower if molecular clouds have decreased density and/or increased temperature (e.g., Bolatto et al. 2013), have increased velocity dispersion (e.g., Watanabe et al. 2011), or have increased metallicity (e.g., Israel 1997; Wolfire et al. 2010). Galaxy centers tend to have lower αCO, perhaps due to higher temperatures, higher metallicity, and/or dynamical effects in the centers (e.g., Strong et al. 2004; Sandstrom et al. 2013; Sánchez et al. 2014; Israel 2020; Teng et al. 2022). A central αCO depression reduces the contribution of molecular gas in the central region when deriving molecular gas concentrations, and the true molecular gas concentrations are therefore lower than those derived in the literature. The degree of central depression in αCO varies from galaxy to galaxy, exerting different levels of influence on the molecular gas concentrations. It is not known how a variable αCO quantitatively influences the molecular gas concentrations and the related results.
In this work we aim to understand the influence of variable αCO on molecular gas concentrations and study the impact of disk (spiral+bar) structure on the radial distribution of molecular gas and star formation. The paper is structured as follows. Section 2 describes the data. Section 3 presents the methods for measuring the central concentrations of CO luminosity, molecular gas and the star formation rate, and for calculating the strength of non-axisymmetric disk structure. Sections 4 and 5 show the results and discussions, respectively. A summary of the main conclusions is given in Sect. 6.
2. Observational material
2.1. Sample and data
The Extragalactic Database for Galaxy Evolution survey (EDGE; Bolatto et al. 2017) observes 12CO J = 1–0 and 13CO J = 1–0 in 126 nearby galaxies using a combination of D and E configurations (D+E) of the Combined Array for Research in Millimeter-wave Astronomy (CARMA) interferometer. It has a typical synthesized beam of ∼45, corresponding to ∼1.5 kpc for the EDGE sample. EDGE galaxies are selected from the Calar Alto Legacy Integral Field Area (CALIFA) surveys (Sánchez et al. 2012), and these galaxies form a representative sample of the SF galaxies in CALIFA, although most of them were selected for their high mid-infrared flux and convenience of scheduling observations (Bolatto et al. 2017). The average RMS noise in EDGE data cubes gives ∼50 mK, corresponding to a typical 3σ sensitivity of molecular gas mass surface density: ∼11 M⊙ pc−2 (Bolatto et al. 2017). We used the 12CO J = 1–0 moment 0 maps from D+E data cubes to quantify the central concentrations of CO luminosity and molecular gas probed in the present paper. The significant noise in the maps was rejected by applying a blanking mask through the smooth masking approach (Bolatto et al. 2017). We used smoothed masks instead of dilated masks, as the smoothed-mask moment maps provide a more complete accounting of the total CO flux (Bolatto et al. 2017). The maps were resampled so that each resolution element has ∼4 pixels (Utomo et al. 2017).
The integral field unit (IFU) data from the third CALIFA data release was used to characterize properties of metallicity, stars, and star formation. We used maps of emission line of Hα, Hβ, [O I] λ6300, N II] λ6584, [O III] λ5007, and [S II] λλ6717, 31, and maps of stellar surface density (Σ⋆), obtained from the PIPE3D pipeline (Sánchez et al. 2016). The oxygen abundance (O/H) derived based on O3N2 ratio is acquired from Marino et al. (2013). The maps were re-gridded and smoothed to match the pixel scale and resolution of the EDGE data (Utomo et al. 2017). Pixels with signal-to-noise ratio (S/N)< 2 and those occupied by foreground stars were set to blank.
From the parent EDGE sample, we excluded galaxies that were not suitable for studying disk structures based on visual inspection. First, we excluded elliptical galaxies, where a large-scale disk is not present. Second, we rejected galaxies whose optical image is severely contaminated by bright foreground stars or nearby galaxies. Galaxies presenting tidal tails, which are thus interacting galaxies, were rejected. This process was to ensure the robustness of quantifying the non-axisymmetric disk structure, that is, spiral arms and bars. Finally, nearly edge-on galaxies were removed. This selection process was done by visual inspection, rather than using a cut in inclination angle derived from their axis ratio, because the axis ratio would be underestimated for early-type edge-on galaxies with spherical stellar haloes. Disk structures are invisible and unrecoverable in edge-on galaxies. Furthermore, three galaxies without CO detection in their central regions (R < 0.2Re, where R is the galactocentric radius in the face-on viewing angle and Re is the optical half-light radius) were removed because their gas concentrations are not available (Sect. 3.2). Our final sample consists of 57 galaxies (Table A.1).
Figure 1 compares Hubble types, distances, and global molecular gas masses (Bolatto et al. 2017) of our derived sample and the EDGE parent sample. Our sample favors late-type (Sb–Sd) galaxies, showing a similar Hubble-type distribution to late-type galaxies in the EDGE parent sample. We include one S0 galaxy, NGC 5784. S0s are disk galaxies without spiral arms by definition, but some S0s could have faint arms (Kormendy & Bender 2012; Yu & Ho 2020). NGC 5784 has molecular gas mass of 109.4 M⊙ (Bolatto et al. 2017). Its gas content was not removed, contrary to the scenario that progenitors of S0s should lose gas and then structurally subside (van den Bergh 1976). The vanished arm structure in NGC 5784 therefore represents weak effect of spiral arms on radial distribution of gas. This galaxy is included in the sample. The histograms of distance and total molecular gas mass are close to those of the EDGE parent sample, indicating that there is no significant bias in these two parameters for our sample. All galaxies have available 2D multicomponent photometric decomposition from Méndez-Abreu et al. (2017). We use the bulge half-light radius (Rbul) measured at r-band from Méndez-Abreu et al. (2017) when determining disk-dominated regions in Sect. 3. Other parameters are derived in Sect. 3. We use the bar identification from Walcher et al. (2014). 34 (60%) objects are barred galaxies (SAB or SB), and 23 (40%) objects are unbarred galaxies (SB).
![]() |
Fig. 1. Histogram of Hubble types, distance, and total molecular gas mass (Mmol) of galaxies in the EDGE parent sample and our derived sample. |
2.2. Nuclear activity
The proposed relationship between central concentrations of molecular gas and strengths of non-axisymmetric disk structure may provide a tool to investigate effects of active galactic nucleus (AGN) feedback on the radial distribution of molecular gas. Measurement of the SFR based on Hα emission suffers from contamination by emission from the AGN, and AGN hosts should therefore not be involved when probing galaxy central star formation. For the above two reasons, we classify the dominant energy source for the galaxies in our sample.
We use the Baldwin-Phillips-Terlevich (BPT) diagram (Baldwin et al. 1981) to search for AGN candidates. The BPT diagram makes use of the emission line ratios [O III]/Hβ, [N II]/Hα, [S II]/Hα, and [O I]/Hα. We classify galaxies as pure SF galaxies if the central (R < 0.2Re) pixels locate below the pure star formation line of Kauffmann et al. (2003) on the [N II]/Hα versus [O III]/Hβ diagram. The reason for using 0.2Re as inner boundary is to be consistent with the definition of central surface density of molecular gas mass and SFR (see Sect. 3). Composite line ratios are caused by a combination of star formation and AGN activity (Kewley et al. 2006; Yuan et al. 2010). We classify galaxies as composite galaxies if the central pixels lie between the pure star formation line of Kauffmann et al. (2003) and the extreme starburst classification line of Kewley et al. (2001).
The remaining galaxies are further classified as candidate Seyferts, low-ionization nuclear emission line regions (LINERs), and ambiguous galaxies, using the Seyfert-LINER demarcation lines on the [S II]/Hα versus [O III]/Hβ and [O I]/Hα versus [O III]/Hβ diagram (Kewley et al. 2006). Candidate Seyferts and LINERs are those galaxies whose central pixels lie above and below the two demarcation lines, respectively. This procedure results in 6 candidate Seyferts. No LINERs are found, perhaps due to the EDGE sample selection requiring the galaxies to be mid-infrared bright. If galaxies do not meet the above criteria, we classify them as ambiguous galaxies. The BPT diagram has been widely adopted to identify AGNs (e.g., Kewley et al. 2006; Ho et al. 2008; Koss et al. 2010). However, diffuse ionization associated with old-stellar populations (Binette et al. 1994; Singh et al. 2013; Gomes et al. 2016) may also populate the regions for AGNs in BPT diagram. Lacerda et al. (2018) find that regions where the ionization is dominated by hot low-mass evolved stars have Hα equivalent widths (EWHα) < 3Å. It is therefore possible to distinguish between presence of an AGN from old-stellar ionization by introducing a cut in the EWHα (Lacerda et al. 2018, 2020; Levy et al. 2019; Kalinova et al. 2021). Following the strategy in Lacerda et al. (2020), we require that AGNs have central EWHα larger than 3 Å. We then classify the candidate AGNs with central EWHα > 3 Å as real AGNs. Candidates with central EWHα < 3 Å could have ionization contaminated by old-stellar populations and are classified as ambiguous galaxies.
The above selection results to 30 pure SF galaxies, 15 composite galaxies, 4 Seyferts (NGC 2410, NGC 2639, NGC 6394, UGC 3973), and 8 ambiguous galaxies. Our AGN classification is consistent with the classification in Lacerda et al. (2020). Specifically, UGC 3973 is a type 1 AGN, and NGC 2410, NGC 2639, and NGC 6394 are type 2 AGNs (Lacerda et al. 2020).
3. Methods
In this section we first describe how to measure the strength of non-axisymmetric disk structure, which traces average strength of spiral arms for unbarred galaxies and average strength of spirals and bars for barred galaxies. Later, we present the procedure to quantify central concentrations of CO luminosity, molecular gas, and SFR.
3.1. Strength of non-axisymmetric disk structures
We analyze the disk (spiral+bar) structure using the r-band images from the Data Release 12 (Alam et al. 2015) of the Sloan Digital Sky Survey (SDSS; York et al. 2000). To reduce the image data, we use SEP1 (Bertin & Arnouts 1996; Barbary 2016) to automatically generate a preliminary mask of foreground stars, and then manually mask out the stars, mainly inside the galaxy, that SEP has missed. After applying the mask to reject contamination, we average the fluxes over the region where the flux profile flattened to calculate background and then subtract this value from the image.
For each galaxy we run IRAF task ellipse with an exponential step of 0.05 to obtain isophotal parameters such as isophotal ellipticity (ϵ) and position angle (PA). We estimate the intrinsic galaxy light affected by the masked region using the task bmodel, which builds a smooth representation of the galaxy light based on the isophotes, and then use these values to fill in the masked region to produce the star-cleaned images. The top-left and bottom-left panel of Fig. 2 illustrates the star-cleaned images of NGC 447 and UGC 9476, respectively. The galaxy disk-averaged ellipticity ϵdisk and position angle PAdisk are obtained by averaging the profiles of ϵ and PA over the outer disk or by minimizing the real part of the m = 2 mode of 2D Fourier spectra at a radial wavenumber of zero (for details, see Grosbøl et al. 2004 and Yu et al. 2018). The measured ϵdisk and PAdisk are listed in Table A.1. Our measured ϵdisk and PAdisk are consistent with those from HyperLEDA (Makarov et al. 2014) with a standard deviation of 0.07 and 10°, respectively. The 0.07 and 10° are taken as the measurement uncertainty of ϵdisk and PAdisk, respectively. We then apply elliptical apertures with the ϵdisk and PAdisk to derive R50% (denoted as Re hereafter) and R90%, which contains 20% and 90% of total flux of the galaxy.
![]() |
Fig. 2. Illustration of how to measure the strength of non-axisymmetric disk structures. The left two panels show the star-cleaned r-band image of NGC 447 and UGC 9476. The dashed blue ellipse in the image of NGC 447 denotes its bar (Appendix B). The right panels show radial profiles of Atot. The galactocentric radius at which the bulge terminates and the disk starts to dominate (1.5 × Rbul/(1 − ϵdisk), where ϵdisk is the disk-averaged ellipticity) is indicated by the dotted line. The bar radius is indicated by the dashed line for NGC 447. The R90%, enclosing the majority of structure, is indicated by the solid line. We compute sdisk as the mean Atot over the disk-dominated region between 1.5 × Rbul/(1 − ϵdisk) and R90%, and define the strength of non-axisymmetric disk structure as 2 + log sdisk. The disk-dominated region for barred galaxies contains spirals and bars, but that for unbarred galaxies contains only spirals. |
In order to quantify the non-axisymmetric disk (spiral+bar) structure, we first define a disk-dominated region for each galaxy. A galaxy is bulgeless if a bulge component is not included in the 2D multicomponent photometric decomposition in Méndez-Abreu et al. (2017). The inner boundary of the disk-dominated region is set to for bulgeless galaxies. The
is the typical SDSS r-band point spread function (PSF) full width at half medium. For galaxies hosting a bulge, the inner boundary is set to 1.5 × Rbul/(1 − ϵdisk). The factor of (1 − ϵdisk) is to take into account the geometric difference between the round bulge or PSF and the inclined disk. The factor of 1.5 is to ensure that the majority of bulge flux is excluded. The outer boundary of the disk-dominated region is set to R90%, which by definition encloses the majority of the disk (Yu et al. 2018).
The strength is the most fundamental property of spirals and bars, and therefore of the disk structure. We adopt a homogeneous method to quantify the disk (spiral+bar) structure based on 1D Fourier decomposition of azimuthal intensity distribution along each isophote (Elmegreen et al. 1989, 2011; Laurikainen et al. 2004; Rix & Zaritsky 1995; Grosbøl et al. 2004; Durbala et al. 2009; Baba 2015; Kendall et al. 2011, 2015; Yu et al. 2018; Yu & Ho 2020). We run the ellipse task with ellipticity fixed to ϵdisk, position angle fixed to PAdisk, and with a linear step of 14. The resulting azimuthal intensity distribution, I(R, θ), along an isophote at radius R is decomposed through
where I0 is the azimuthally averaged intensity and represents the axisymmetric disk, and Im is the Fourier amplitude, which measures the structure amplitude. The relative Fourier amplitude is defined as
The m = 2 is often used to quantify the strength of bars or spirals (e.g., Grosbøl et al. 2004; Elmegreen et al. 2007; Durbala et al. 2009). However, since the shape of I(R, θ) is not an exact cosine, a bar or spiral would contribute to higher order modes (Rix & Zaritsky 1995). Following the strategy in Yu et al. (2018), we define relative amplitude of structure as the quadratic sum of the relative Fourier amplitude of m = 2, 3, and 4 modes:
The top-right and bottom-right panels of Fig. 2 illustrate the Atot(R) for NGC 447 and UGC 9476, respectively. The disk-dominated regions of NGC 447 and UGC 9476 are marked by the region between the vertical dotted (the end of bars or bulges) and the vertical solid (the end disks) lines in Fig. 2. The bar of NGC 447, quantified using isophotal analysis described in Appendix B, is illustrated by the dashed ellipse in the top-left panel and the bar radius is shown as the vertical dashed line in the top-right panel. The disk-dominated region for barred galaxies contains spirals and bars, but that for unbarred galaxies contains only spirals. The average relative amplitude of the disk structure, sdisk, is defined as the mean Atot over the disk dominated region:
where Rd belongs to the disk-dominated region. sdisk = 0.49 for NGC 447 and sdisk = 0.29 for UGC 9476. The strength of non-axisymmetric disk structure is then defined as 2 + log sdisk. The logarithmic format is used as the spiral or bar effect may be highly nonlinear (Yu et al. 2021). The value of 2 is added to make the quantity greater than zero. Uncertainty of the structure strength is estimated by considering the standard error of the mean value, error of ϵdisk, and error of PAdisk. The derived 2 + log sdisk and uncertainties are listed in Table A.1. In strongly barred galaxies, the disk structure strength is driven by bars, because the strong bar is long and dominates the disk-dominated region. However, in weakly barred galaxies, the disk structure strength is driven by spirals, because the weak bar is short (Elmegreen et al. 2007) and the spiral arms rule the disk-dominated region. Despite the different radial extent occupied by bars and spirals, the disk may facilitate the formation of arms and bars to a similar level of strength (Salo et al. 2010; Díaz-García et al. 2019, but see Athanassoula et al. 2009). The disk structure strength (2 + log sdisk) traces average spiral strength for unbarred galaxies, and average spiral and bar strength for barred galaxies.
The bar strength has been quantified in the literature by its ellipticity or length (Martin 1995; Martinet & Friedli 1997; Aguerri 1999; Marinova & Jogee 2007; Barazza et al. 2008, 2009; Gadotti 2009; Li et al. 2011). It is instructive to compare bar strengths defined in our framework with these classic indicators. We calculate average relative Fourier amplitude over the region occupied by the bar: sbar and then compute the bar strength: 2 + log sbar. We find a good relation where the value of 2 + log sbar gets higher with increasing bar ellipticity (Fig. B.2) and bar length (Fig. B.3). A more sophisticated measure of bar strength is to estimate the gravitational torque (Block et al. 2004; Laurikainen et al. 2004; Buta et al. 2005). Nevertheless, it has been shown that bar torques are very well correlated with bar ellipticity (Laurikainen & Salo 2002; Block et al. 2004). We therefore argue that our approach for quantifying the bars is validated.
3.2. Central concentrations of CO luminosity and molecular gas
We use the EDGE 12CO J = 1–0 integrated intensity to derive the central concentrations of CO luminosity and molecular gas. Pixels without CO detections or with a S/N of less than 1 are blanked. Using S/N of 3 as the threshold yields consistent gas concentrations (Pearson correlation coefficient ρ > 0.9). Sakamoto et al. (1999a) and Sheth et al. (2005) defined the molecular gas concentration as the ratio of the central 1 kpc molecular gas surface density to the density average within the diameter at 25 mag arcsec−2 (D25; R25 = D25/2). Kuno et al. (2007) suggested using a fraction of the galaxy diameter as the inner aperture, instead of a fixed radius, to avoid a possible dependence on the galaxy size. Larger galaxies tend to have longer bars (Gadotti 2009), which may have a more prominent effect on the gas distribution (Kuno et al. 2007). Chown et al. (2019) used a different definition, the ratio of optical half-light radius to CO flux half-light radius. This definition measures gas concentrations in central ∼8 kpc for our sample. However, an aperture of this size is too large and would cancel out the difference in central gas concentration between barred and unbarred galaxies (Kuno et al. 2007), making this definition undesirable.
The molecular gas concentrations calculated by these authors are based on a constant αCO, and therefore have bias and a certain amount of uncertainty if αCO varies from location to location within the galaxy. To be precise, they measured central concentrations of CO luminosity. Using the projection parameters, ϵdisk and PAdisk, to take into account the inclination effect and, similar to the strategy in Kuno et al. (2007), we define the central concentration of CO luminosity (CCO) as the ratio of the average central CO intensity density (ICO(R ≤ 0.2Re)) to the disk-averaged value (ICO(R ≤ Re)):
The galaxy size is characterized by the Re. The inner radius of 0.2Re is chosen to ensure that the aperture for each galaxy is larger than the beam size of 45. The average 0.2Re of galaxies in our sample corresponds to 1.1 kpc (2.2 kpc in diameter), consistent with the inner radius used in Kuno et al. (2007). Equation (5) therefore measures CO luminosity concentration in the central ∼2 kpc. The CCO should be slightly smaller than the concentration computed using 1 kpc as the inner aperture diameter, especially for highly concentrated CO distributions. We refrain from doing the same calculation as in Sakamoto et al. (1999a) and Sheth et al. (2005), because the central 1 kpc is smaller than the beam size for most of the galaxies in our sample.
The integrated intensity was converted to molecular surface density using a variable, αCO. We considered two functional forms of αCO (i.e., αCO(R) and αCO(Z, Σ⋆)) to probe their effect on the calculation of molecular gas concentrations. Sandstrom et al. (2013) solved for αCO by assuming that dust and gas are well mixed and the dust-to-gas ratio is approximately constant within a given region of ∼1 kpc. They find that the αCO in the center is on average a factor of 2 lower than the rest of the galaxy. We adopted the empirical normalized average radial profile of αCO from Sandstrom et al. (2013) and multiplied the profile by the Milky Way αCO, MW = 4.4 M⊙ pc−2 (K km s−1)−1 to obtain the un-normalized profile αCO(R) (Fig. 3). We fit the data using a broken function:
![]() |
Fig. 3. αCO plotted against galactocentric radius normalized by R25. Top panel:αCO(R) obtained by multiplying the empirical normalized average profile of αCO adopted from Sandstrom et al. (2013) by Milky Way αCO, MW = 4.4 M⊙ pc−2 (K km s−1)−1; the errors are also from Sandstrom et al. (2013) and the black curve is the best-fit function. Bottom panel:αCO(Z, Σ⋆), marked in gray, derived from Eq. (8), which considers the decrease in αCO driven by higher metallicity (Z) and stellar surface density (Σ⋆); the black curve denotes the mean αCO for a given radius. |
The best-fit parameters are: k = 1.2, b = 0.3, and Rc = 0.3. The center of αCO(R) gives 2.1 M⊙ pc−2 (K km s−1)−1. The αCO(R) increases with radius, flattens out at r = 0.3R25, and then is fixed at 4.6 M⊙ pc−2 (K km s−1)−1 (black curve in the top panel of Fig. 3). Compared with the outer disk, the central αCO(R) is a factor of 2.2 (0.34 dex) lower. The same best-fit function is used to convert the CO intensity to molecular surface density (Σmol) for each galaxy, and then we calculate the molecular gas concentration (Cmol(αCO[R])) as the ratio of central molecular surface density (Σmol(αCO[R];R ≤ 0.2Re)) to the disk-averaged surface density (Σmol(αCO[R];R ≤ Re)):
Cmol(αCO[R]) is only for understanding the effect of the same central αCO depression, not for exploring the impact of disk structure on central concentrations of molecular gas.
The degree of central depression in αCO varies from galaxy to galaxy (Sandstrom et al. 2013), imposing different level of influence on molecular gas concentrations. Considering two primary dependences of αCO, Bolatto et al. (2013) suggested a prescription of αCO (their Eq. (31)), which involves metallicity (Z) to take into account the CO-faint molecular gas and involves total surface density (Σtotal) to take into account the effects of temperature and velocity dispersion. We replace the Σtotal with Σ⋆, because the stellar component is dominant in our sample. We assume that giant molecular clouds have a characteristic mean surface density of 100 M⊙ pc−2, as did Bolatto et al. (2013) to derive prediction. We therefore have
with γ = 0.5 for Σ⋆ > 100 M⊙ pc−2 and γ = 0 otherwise. Z⊙ is the solar metallicity. Since the O/H measurements are not available at some locations in the galaxy, we fit a straight line to the 12 + log(O/H) versus R, and use the best-fit line to regenerate a 2D map of 12 + log(O/H). The solar oxygen abundance of 8.7 is adopted. The metallicity (Z) is given by: log(Z/Z⊙)=12 + log(O/H)−8.7 (Marino et al. 2013). To avoid the masked region, we do azimuthal averaging of Σ⋆ at each radius to obtain the radial profile of Σ⋆, then regenerate a 2D map of Σ⋆. These two maps are entered into Eq. (8) to calculate αCO(Z, Σ⋆). The derived αCO(Z, Σ⋆) for each galaxy are marked by gray curves in Fig. 3, and the average values for a given normalized radius are marked by thick black curves. The derived αCO(Z, Σ⋆) are consistent with the observations in Sandstrom et al. (2013) within uncertainty. The average αCO(Z, Σ⋆) yields 1.3 M⊙ pc−2 (K km s−1)−1 in the center (R = 0), and 5.8 M⊙ pc−2 (K km s−1)−1 in the outer part (0.8 < R/R25 < 1). Compared with the outer disk, the central αCO(Z, Σ⋆) is therefore lower by a factor of 4.5 (0.65 dex). We convert CO intensity to molecular gas surface density using the derived αCO(Z, Σ⋆) for each galaxy, and define the molecular gas concentration:
The derived Cmol(αCO[Z, Σ⋆]) is used to explore the impact of disk structure on the molecular gas concentration. The calculation of molecular gas concentration (e.g., Cmol(αCO[R]) or Cmol(αCO[Z, Σ⋆])) does not depend on the normalization of αCO, because any multiplicative constant term in the numerator and denominator cancels out. Instead, the shape of the αCO profile influence the molecular gas concentrations.
To estimate the uncertainty, we perform a Monte Carlo test on the CCO, Cmol(αCO[R]), and Cmol(αCO[Z, Σ⋆]), by adding random noise to CO integrated intensities, ϵdisk, and PAdisk, according to their statistical errors. We repeat the calculations with the randomly perturbed data values 1000 times and find the standard deviation of the results.
The purpose of deriving Cmol(αCO[R]) with αCO(R) fixed the same for all the galaxy is to first understand the effect of the same central αCO(R) depression, which helps reveal the effect of αCO(Z, Σ⋆), which changes from galaxy to galaxy. Figure 4 compares CCO with Cmol(αCO[R]) and Cmol(αCO[Z, Σ⋆]). The CCO correlates strongly with Cmol(αCO[R]) with a Pearson correlation coefficient ρ = 0.99, and Cmol(αCO[R]) are systematically 0.19 dex lower. Similarly, in addition to being strongly correlated (ρ = 0.96), Cmol(αCO[Z, Σ⋆]) are systematically 0.22 dex lower than CCO. The difference between CCO and Cmol(αCO[Z, Σ⋆]) is larger than between CCO and Cmol(αCO[R]), because the central depression is deeper in αCO(Z, Σ⋆) than in αCO(R). The scatter in the CCO-Cmol(αCO[R]) relation is 0.05 dex, which is lower than the scatter of 0.08 dex in the CCO-Cmol(αCO[Z, Σ⋆]) relation. This is because the same αCO(R) is used, whereas αCO(Z, Σ⋆) varies from galaxy to galaxy. By comparing CCO-Cmol(αCO[R]) and CCO-Cmol(αCO[Z, Σ⋆]) relation, we can be sure that the lower value in Cmol(αCO[Z, Σ⋆]) compared to CCO is caused by the central depression in αCO(Z, Σ⋆).
![]() |
Fig. 4. Correlations between CO luminosity concentration (CCO) and molecular gas mass concentration (Cmol). Top panel:Cmol(αCO[R]) derived based on αCO(R), which is a function of radius and is set the same for all the galaxies (Eq. (6)). Bottom panel:Cmol(αCO[Z, Σ⋆]) computed based on αCO[Z, Σ⋆], which is a function of Z and Σ⋆ and varies from galaxy to galaxy (Eq. (8)). The dashed lines mark the 1 : 1 relations. |
The caveat is that there is still a fair dispersion around the average prescription αCO(Z, Σ⋆), representing the variation in local parameters such as temperature and surface density of giant molecular clouds, although αCO(Z, Σ⋆) has considered the main drivers of αCO variations (see the detailed discussion in Bolatto et al. 2013). We use the molecular mass map derived using αCO(Z, Σ⋆) and probe the connection between Cmol(αCO[Z, Σ⋆]) and disk structure. We adopt
in the remainder of this paper. We have verified that the resulting CCO and Cmol are almost identical to those computed using dilated-mask moment maps (see Bolatto et al. 2017 for the derivation of dilated-mask moment maps), so using the dilated-mask moment maps does not affect our results. The derived CCO, Cmol and their uncertainties are listed in Table A.1.
3.3. Central concentrations of star formation rate
We use the extinction-corrected SFR derived from Hα flux using a Salpeter initial mass function and the nebular extinction based on the Balmer decrement (Bolatto et al. 2017). SFR in pixels with EWHα less than 6 Å are blanked, because they are primarily ionized by evolved stars (Sánchez et al. 2014). The derived SFR traces ongoing star formation averaged over the past ∼10 Myr (Kennicutt & Evans 2012). The SFR surface densities (ΣSFR) are then computed by dividing the SFR by the physical area of each pixel.
Measuring SFRs in active galaxies is challenging because tracers of young stars are to some extent contaminated with radiation from the AGN. We removed the four Seyferts and eight ambiguous galaxies, where the dominant power source may be an AGN, to avoid contamination in quantifying the SFR concentration. This leads to 45 galaxies. Cid Fernandes et al. (2011) showed that a value of EWHα ≤ 3 Å corresponds to retired regions in galaxies, which have ceased star formations and are now ionized by hot low-mass evolved stars. The EWHα ≥ 6 Å are primarily dominated by recent star formation, while those with 3 Å < EWHα < 6 Å are mixed (Sánchez et al. 2014; Lacerda et al. 2018). We require that the average central EWHα(R ≤ 0.2Re) ≥ 6 Å to ensure robust measurement of SFR. This leads to the final sample of 38 galaxies for measuring the SFR concentration.
We define the concentration of SFR (CSFR) as the ratio of the nuclear SFR surface density (ΣSFR(R ≤ 0.2Re)) to the disk-averaged SFR surface density (ΣSFR(R ≤ Re)):
If the central SFR is enhanced with respect to global SFR, the ΣSFR(R ≤ 0.2Re) increases for a given ΣSFR(R ≤ Re), and thus CSFR is elevated. Uncertainties of CSFR are estimated from the errors of SFR, ϵdisk, and PAdisk through a Monte Carlo test. The derived CSFR and uncertainties are listed in Table A.1.
4. Results
In this section we first demonstrate the ability of CCO and Cmol to measure central concentrations, then explore distributions of CCO and Cmol in barred and unbarred galaxies, and finally study the dependence of CCO and Cmol on the strengths of disk (spiral+bar) structure. We use the Pearson correlation coefficients to analyze the relationships. The Pearson correlation coefficient measures tightness of a linear relation between two set of data. The Spearman correlation coefficient measures relation tightness, regardless of whether the relation is linear or not. Still, there is not significant difference between the two correlation coefficients of the relations presented in this paper. We use the Pearson correlation coefficient (ρ) in this work due to the perspicuity of linearity.
4.1. Profiles of the molecular gas surface density
To demonstrate the ability of CCO and Cmol to measure central concentrations, we extract profiles of CO intensity (ICO(R)) and molecular gas surface mass (Σmol(R); αCO(Z, Σ⋆) is used) using a series of elliptical rings with measured ϵdisk and PAdisk. and
are respectively the 33th and 66th percentile of the CCO in the sample.
and
are respectively the 33th and 66th percentile of the Cmol in the sample. We separate the data into four groups: high-CCO galaxies (
), low-CCO galaxies (
), high-Cmol galaxies (
), and low-Cmol galaxies (
).
To visualize their shape, we normalize the ICO(R) and Σmol(R) to their values at R = Re, that is, ICO(R)/ICO(Re) and Σmol(R)/Σmol(Re), and plot them as a function of R/Re in Fig. 5 (top two panels). The average ICO(R)/ICO(Re) profiles of high-CCO galaxies and average Σmol(R)/Σmol(Re) profiles of high-Cmol galaxies are approximately exponential, but they slightly bends upward at R ≈ 0.6 Re. The average profiles for low-CCO and low-Cmol galaxies flatten below Re and become almost flat at the center. The difference between low and high concentrations is the most pronounced at the center, consistent with the definition of CCO and Cmol, which compare the value averaged within 0.2Re to the value averaged within Re. In contrast, there is no significant difference in the profiles beyond Re between galaxies with low and high concentrations. We point out that the CCO and Cmol reflect the global shape rather than a central precipitous rise within 0.2Re. At the center, the Σmol(0)/Σmol(Re) is lower than ICO(0)/ICO(Re), because the αCO(Z, Σ⋆) is on average 0.65 dex lower in the center than in the outer disk. This is consistent with Fig. 4, where the Cmol is systematically lower than CCO.
The definitions of CCO and Cmol include a normalization term of galaxy size and therefore remove the dependence on the size. Nevertheless, we also visualize the profiles normalized by the size of the CO disk. The R20%,CO (R20%,mol) and R80%,CO (R80%,mol), which contains 20% and 80% of total CO flux (total molecular mass) are derived. We, respectively, normalize the ICO(R) and Σmol(R) to their values at R = R80%,CO and R = R80%,mol, that is, ICO(R)/ICO(R80%,CO) and Σmol(R)/Σmol(R80%,mol) (bottom two panels in Fig. 5).
![]() |
Fig. 5. Average normalized profiles of CO integrated intensity (ICO(R)/ICO(Re) and ICO(R)/ICO(R80%,CO)) and molecular mass surface density (Σmol(R)/Σmol(Re) and Σmol(R)/Σmol(R80%,mol)). Blue profiles are for high-CCO or high-Cmol galaxies, while red profiles are for low-CCO or low-Cmol galaxies. The shaded regions mark the standard deviation in y-axis value for a given normalized radius. |
Compared to those normalized to the values at R = Re, the difference in ICO(R)/ICO(R80%,CO) between high- and low-CCO galaxies or the difference in Σmol(R)/Σmol(R80%,mol) between high- and low-Cmol galaxies becomes smaller, especially at the center. However, other main features as in the profiles normalized using Re are still present. Inspired by the definition of galaxy light concentration (Conselice 2003), we additional test two kinds of concentrations, without consideration on the stellar components, which are log(R80%,CO/R20%,CO) and log(R80%,mol/R20%,mol). We find that the CCO are consistent with log(R80%,CO/R20%,CO) with ρ = 0.83 and p < 0.01, and that Cmol are consistent with log(R80%,mol/R20%,mol) with ρ = 0.73 and p < 0.01. We therefore confirm that the CCO and Cmol truly reflect the central concentrations of CO intensity and molecular gas mass.
4.2. Molecular gas concentrations in barred and unbarred galaxies
Figure 6 presents histograms of CCO (top) and Cmol (bottom) for unbarred and barred galaxies. Two key features are revealed. The CCO and Cmol in unbarred galaxies are mainly moderate to low. The unbarred galaxies with the highest concentration is NGC 7819, which has CCO = 0.85 and Cmol = 0.69, meaning that the average central CO intensity and molecular mass surface density are, respectively, a factor of 7 and 5 higher than the disk-averaged values. The concentration in NGC 7819 is comparable to the highest gas concentration of barred galaxies. In contrast, there is a wide range of CCO and Cmol in barred galaxies. They can be as low as those in unbarred galaxies with the lowest concentrations, or they can reach very high values, far exceeding the concentrations in most of the unbarred galaxies. The barred galaxy with the highest concentration is NGC 447, which has CCO = 1.01 and Cmol = 0.89, meaning that the average central CO intensity and molecular mass surface density are a factor of 10 and 8 higher than the disk-averaged values, respectively. The histogram of Cmol is similar to that of CCO, albeit shifting toward lower values.
![]() |
Fig. 6. Histograms of the CO luminosity concentrations (CCO) and molecular gas concentrations (Cmol). The results for barred and unbarred galaxies are marked in black and gray, respectively. The p value of the Kolmogorov-Smirnov test (pKS) is presented. |
We perform a Kolmogorov-Smirnov (KS) test for the CCO in barred and unbarred galaxies, and a test for the Cmol. We obtain the same pKS value of 0.18 for the two tests. The pKS is small because the barred galaxies tend to have higher CCO or Cmol, or, in other words, high CCO or Cmol are more common in barred galaxies. But the pKS is not small enough to reject the null hypothesis that the CCO or Cmol in barred and unbarred galaxies are drawn from the same distribution. This is because there are a few unbarred spiral with relatively high concentrations (see also Figs. 1 and 2 in Sheth et al. 2005).
Assuming a constant αCO, Sheth et al. (2005) showed that molecular gas concentrations in the central 1 kpc is systematically higher in barred spirals than unbarred spirals (Sakamoto et al. 1999a; Kuno et al. 2007; see also Komugi et al. 2008). We derived Cmol using the αCO(Z, Σ⋆) that varies from galaxy to galaxy and confirmed the previous findings in the literature. The αCO(Z, Σ⋆) therefore does not significantly affect the main conclusions in Sakamoto et al. (1999a), Sheth et al. (2005), and Kuno et al. (2007). By comparing the molecular gas concentrations in barred and unbarred galaxies, it has been concluded that the bar is an efficient mechanism for driving gas to the galaxy center (Sakamoto et al. 1999a; Sheth et al. 2005; Kuno et al. 2007). On the other hand, this implies that the effect of spiral arms in driving gas inflow is in general moderate to weak.
Instructively, bars tend to be stronger than spiral arms by using Fourier amplitude, arm/inter-arm contrast, or gravitational torque as a measure of their strength (Buta et al. 2005; Durbala et al. 2009; Bittner et al. 2017; Yu & Ho 2020). We confirm this behavior in Fig. 7, where histogram of the strengths of disk structure is shown. The disk structure in unbarred galaxies equals the spiral arms, whereas the disk structure of barred galaxy includes both spiral arms and bars. The disk structure in barred galaxies are systematically stronger than in unbarred galaxies. A KS test yields pKS < 0.01, suggesting that we can reject the null hypothesis that the two sets of 2 + log sdisk are drawn from the same parent distribution. In strongly barred galaxies, the disk structure strength is driven by bars, as the strong bar is long and dominates the disk-dominated region. In weakly barred galaxies, the disk structure strength is actually driven by spirals, as the weak bar is short (Elmegreen et al. 2007) and the spirals dominate the disk-dominated region. Stronger bars tend to be associated with stronger spiral arms, either because the bars drive spirals (e.g., Yuan & Kuo 1997; Athanassoula et al. 2009) or because the disk facilitates the formation of arms and bars to a similar extent (Salo et al. 2010; Díaz-García et al. 2019). As in Figs. 6 and 7, barred spirals tend to have more centrally concentrated molecular gas distribution and, meanwhile, have stronger disk structures. The similarities between Figs. 6 and 7 imply a connection between gas concentration and disk (spiral+bar) structure strength.
![]() |
Fig. 7. Histogram of the strength of non-axisymmetric disk (spiral+bar) structures (2 + log sdisk). The 2 + log sdisk traces average spiral strength for unbarred galaxies, but traces average spiral and bar strength for barred galaxies. The results for barred and unbarred galaxies are marked in black and gray, respectively. The p value of the Kolmogorov-Smirnov test (pKS) is presented. |
4.3. Molecular gas concentrations and structure strength
We study the dependence of CCO and Cmol on the strength of non-axisymmetric disk structure using the parent sample of 57 galaxies to understand the impact of disk (spiral+bar) structure. The CCO and Cmol increase as the disk structure strengths, regardless of whether the galaxies are barred or unbarred (left two panels in Fig. 8). The Pearson correlation coefficients of CCO-(2 + log sdisk) relation yields ρ = 0.79 (p < 0.01) and ρ = 0.53 (p < 0.01) for unbarred and barred galaxies, respectively. Those of Cmol-(2 + log sdisk) relation yields ρ = 0.78 (p < 0.01) and ρ = 0.51 (p < 0.01) for unbarred and barred galaxies, respectively. Figures 6 and 7 are projections of Fig. 8. In unbarred galaxies, we reveal a novel trend that galaxies with stronger disk (spiral) structure tend to have more centrally concentrated CO or molecular gas distribution. In barred galaxies, our reported relationship is consistent with the previous study by Kuno et al. (2007). By assuming a constant αCO and using gravitational torque as a measure of bar strength calculated from Laurikainen & Salo (2002), Kuno et al. (2007) presented a fairly tight positive correlation between molecular gas concentration in the central ∼2 kpc and bar strength.
![]() |
Fig. 8. Dependence of CO luminosity concentrations (CCO) and molecular gas concentrations (Cmol) on the strengths of disk structure (2 + log sdisk). The left two panels show the results for the parent sample, whereas the right two panels show the results for galaxies with i ≤ 65°. The formulas indicate linear fits considering both barred and unbarred galaxies. They are derived using PCA and marked with solid lines. The Pearson correlation coefficients, ρ, for unbarred and barred galaxies are given at the top of each panel. |
The relations for unbarred and barred galaxies are similar (Fig. 8), although unbarred galaxies tend to have lower gas concentration and weaker structure compared to barred galaxies (Figs. 6 and 7). It may imply that the relations of gas concentration against structure strength for unbarred and barred galaxies are drawn from a unified correlation. We derive best-fit straight lines for the correlation considering both unbarred and barred spirals using principal component analysis (PCA; Pearson 1901). Unlike the least squares fit, which only considers the deviation of y-axis values from a model, the PCA fit simultaneously takes into account the x- and y-axis values. The PCA fit is preferred in this work because there should be some degree of intrinsic dispersion in the x-axis values of the Cmol- or CCO-(2 + log sdisk) relations. PCA transforms the data set, for example, CCO versus 2 + log sdisk, into a set of coefficients for two new orthonormal bases (which replace the old bases: unit CCO and 2 + log sdisk along the x and y direction, respectively.). The two new orthonormal bases are the eigenvectors of the data covariant matrix and are ordered according to their eigenvalues to successively maximize variance of the data. The slope of the straight line is determined by the direction of the eigenvector with the largest eigenvalue, along which the data show the greatest variance. The straight line is the best-fit linear relation between the two parameters. The best-fit line for the CCO-(2 + log sdisk) relation yields
The uncertainties are obtained by repeating 1000 times PCA using samples of the same size, regenerated using bootstrapping with replacement, and computing standard deviations of the results. The scatter in CCO for a given 2 + log sdisk is 0.24 dex. The best-fit line for the Cmol-(2 + log sdisk) relation yields
The scatter in Cmol for a given 2 + log sdisk is 0.29 dex. The slopes of the CCO-(2 + log sdisk) and Cmol-(2 + log sdisk) relations are consistent with each other within uncertainty. The intercept of the Cmol-(2 + log sdisk) relation is lower then the intercept of the CCO-(2 + log sdisk) relation with statistical significance, because Cmol is systematically lower than CCO (Fig. 4).
Although we have used an elliptical aperture, with measured ϵdisk and PAdisk, to account for the inclination effect, we are not able to completely eliminate this effect if the inclination angle (i) is too high. Pixels along the minor axis of the galaxy observed CO emission from a larger range of radius than those along the major axis. This effect becomes more significant as i increases, making CCO or Cmol in galaxies with high i more uncertain. The inclination effect has been mitigated because edge-on galaxies have been removed during the sample selection (Sect. 2). The measurement of 2 + log sdisk is less affected as the resolution in SDSS image is better. 12 galaxies in our sample are highly inclined (65° < i < 75°). In Fig. 8 (right two panels), we plot the relationships for a subsample of 45 galaxies with i ≤ 65°. The new CCO- and Cmol-(2 + log sdisk) relations for barred and unbarred galaxies remains almost unchanged.
The best-fit line derived for the CCO-(2 + log sdisk) relation considering both the unbarred and barred galaxies with i ≤ 65° yields
The scatter in CCO for a given structure strength is 0.26 dex. The same for the Cmol-(2 + log sdisk) relation yields
The scatter in Cmol for a given structure strength is 0.31 dex. The best-fit lines (Eqs. (14) and (15)) are consistent with the results considering all galaxies (Eqs. (12) and (13)) within uncertainty. These results show that although the inclination affects the calculated gas concentrations in fact, this effect is not significant in our sample.
4.4. Uncertainty from the bar identification
Our results for barred and unbarred galaxies are obtained using bar identification based on visual inspection from Walcher et al. (2014). Accordingly, 34 (60%) galaxies are barred galaxies and 23 (40%) galaxies are unbarred galaxies. The bar fraction is consistent with previous studies that nearly 60% of disk galaxies in the local universe have a bar (e.g., Sheth et al. 2008; Aguerri et al. 2009; Díaz-García et al. 2016; Erwin 2018). Still, there is uncertainty in identifying bars and it may influence our results. To understand the uncertainty and make our results more robust, we use two more independent bar identifications. The first independent bar identification is based on isophotal analysis (Sect. B). Isophotes are widely adopted to identify and quantify bars (e.g., Athanassoula & Misiriotis 2002; Laine et al. 2002; Erwin & Sparke 2003; Menéndez-Delmestre et al. 2007; Aguerri et al. 2009; Li et al. 2011). As a result, there are 25 (44%) barred galaxies and 32 (56%) unbarred galaxies. This derived bar fraction is underestimated, likely because the isophote-based approach is less sensitive to weak bars. The short and weak bars observed in SDSS images may not have an apparent effect on the isophotes due to limited image resolution (Erwin 2018). Another independent result is from Méndez-Abreu et al. (2017), who performed 2D multicomponent decomposition for the CALIFA galaxies. A galaxy is considered to be a barred galaxy if a bar is included in the decomposition. The inclusion of a bar model was determined by iteratively fitting different composite models and finding the one that best matches the observed structural properties (Méndez-Abreu et al. 2017). There are 32 (56%) barred galaxies and 25 (44%) unbarred galaxies.
We calculate ρ for the CCO- and Cmol-(2 + log sdisk) correlations in barred and unbarred galaxies using bar identification based on isophotal analysis and the 2D decomposition (rows [1] and [2], respectively, in Table 1). If the two new bar identifications are used, the ρ for CCO- and Cmol-(2 + log sdisk) relations in unbarred galaxies ranges from 0.61 to 0.66 (p < 0.01). Likewise, the ρ for the relations in barred galaxies ranges from 0.45 to 0.56 (p < 0.01). Overall, the change in ρ is small and each correlation remains statistical significant (p < 0.05). The three independent bar identifications yield consistent results. The uncertainty from the bar identification therefore does not adversely affect our results.
Ordinary correlation coefficients.
4.5. Partial correlation coefficients
In Fig. 8 we show that galaxies with stronger disk (spiral+bar) structure tend to have more centrally concentrated distribution of CO intensity and molecular gas mass. However, there may be a network of interdependences between gas concentrations, disk structure, and bulge properties. Komugi et al. (2008) found that early-type galaxies show higher CO luminosity concentration, defined as a ratio of average CO intensity within central 1 kpc and 3 kpc, and these authors attributed the trend to the larger bulges seen in earlier Hubble types. However, their definition of concentration used a fixed aperture and had a dependence on galaxy size, as early-type galaxies are larger (e.g., Muñoz-Mateos et al. 2015). By using galaxy size as aperture, Kuno et al. (2007) found no link between gas concentrations and Hubble types. Consistent with Kuno et al. (2007), we find that the correlation between CCO and Hubble types (ρ = 0.04 and p = 0.77; see Fig. 9) and the correlation between Cmol and Hubble types (ρ = 0.11 and p = 0.43) are not apparent.
![]() |
Fig. 9. Comparison of CCO and Cmol with Hubble type. Black points show the data. The Hubble types are from Walcher et al. (2014), and the corresponding numerical Hubble stages are used to calculate the Pearson correlation coefficients between the two parameters. |
We note that the Hubble classification system is not only based on the visual inspection of relative bulge size, but also the character and apparent resolution of spiral arms (Hubble 1926; Sandage et al. 1975). In addition to larger bulges, early type galaxies tend to have a stronger bar (Gadotti 2011), and more tightly wound (Yu & Ho 2019) and weaker (Yu & Ho 2020) spiral arms. The Hubble type may not reveal any unique galaxy properties. We therefore do not apply the analysis of partial correlation coefficients described below to the Hubble types.
One possible concern is that the central bulge may adversely affect the reported relationships, or drive them in some way. Both spiral arms and bars may be associated with the bulges. Bars would lose angular momentum to the disk, bulge, and halo during the secular evolution (Tremaine & Weinberg 1984; Sellwood 1980; Debattista & Sellwood 1998; Athanassoula 2003), and thus slow down, which makes the star obits in the bar more elongated. The bar become longer and, meanwhile, stronger as the stars are more concentrated in the bar (Elmegreen et al. 2007). A more prominent bulge indicates a more evolved system and thus a more evolved bar, resulting in an increase in bar strengths (Laurikainen et al. 2007; Zhou et al. 2015; Gadotti 2011). A massive classical bulge weakens spiral arms (Bittner et al. 2017; Yu & Ho 2020), as massive bulges would stabilize the disk and therefore reduce the dynamically active disk mass that responds to the spiral perturbation (Bertin et al. 1989). In addition to the connection to classical bulges, bars and spirals contribute to the growth of central disky pseudo bulges (Kormendy & Kennicutt 2004). We use the bulge-to-total light ratio (B/T) from Méndez-Abreu et al. (2017) as a proxy to gauge the degree of bulge prominence. To understand the possible influence from the bulge, we compute partial Pearson correlation coefficients (ρ′) by removing common dependence on B/T using the python package pingouin (Vallat 2018) (rows [1] and [2] in Table 2). In both barred and unbarred galaxies, the ρ′ of CCO- and Cmol-(2 + log sdisk) relations remain almost unchanged compared to the original ordinary Pearson correlation coefficients ρ. It suggests that our results are not induced by the central bulges.
Partial correlation coefficients.
5. Discussion
Non-axisymmetric disk structures are essential to understand the galaxy secular evolution (Kormendy & Kennicutt 2004). As the two most generic structures, spirals and bars may play an important role in driving gas into the central regions. Here the gas provides the fuel for the central star formation, which leads to the growth of pseudo bulges (e.g., Kormendy & Kennicutt 2004). While the effects of bars has been extensively explored, the effects of spirals are less understood.
5.1. Bar-driven gas inflow
Theoretical models and simulations show that a stellar bar imposes a non-axisymmetric gravitational potential on the disk to drive the gas clouds flow down the dust lanes along the leading edge of the bar (Athanassoula 1992a,b; Regan et al. 1999; Sheth et al. 2000, 2002; Regan & Teuben 2004; Combes et al. 2014; Fragkoudi et al. 2016; Tress et al. 2020; Sormani et al. 2020). In a galaxy where there is a well-defined ILR introduced by the bar, the bar-driven inward flowing gas accumulates at the ILR at a distance of a few kiloparsecs from the center (Athanassoula 1984; Jenkins & Binney 1994); later the response of gas to forcing by the bar can give rise to a nuclear spiral or bar pattern in the gas, which can transport the gas to the proximity of the central black hole (Englmaier & Shlosman 2000, 2004). Observations have confirmed the bar-driven gas transport (Mundell & Shone 1999; Combes 2003; Zurita et al. 2004; Fathi et al. 2006; Jogee 2006; Haan et al. 2009; Querejeta et al. 2016). The rates of bar-driven gas inflow range from 0.01 to 50 M⊙ yr−1 (Regan et al. 1997; Haan et al. 2009; Querejeta et al. 2016), leading to an accumulation of gas in the center (Sheth et al. 2005). The molecular gas concentrations in the central 1 kpc are systematically higher in barred spirals than in unbarred spirals (Sakamoto et al. 1999a; Sheth et al. 2005; Kuno et al. 2007). By comparing the molecular gas concentrations in barred and unbarred galaxies, it has been concluded that the bar is an efficient tool for driving gas to the galaxy center. Furthermore, this bar effect is more efficient in stronger bars (Regan & Teuben 2004; Hopkins & Quataert 2011; Kim et al. 2012). Consistently, galaxies hosting a stronger bar are found to have a more centrally concentrated gas distribution (Kuno et al. 2007).
These previous observational studies calculated molecular gas concentrations by assuming a constant αCO and they therefore actually measured CO luminosity concentrations. However, αCO depends on the conditions of gas clouds, such as, density, temperature, metallicity, and velocity dispersion (e.g., Narayanan et al. 2012; Bolatto et al. 2013). Likely as a result of these effects, galaxy centers tend to have lower αCO. The Milky Way has a αCO near its center that is 3–10 times lower than in the solar neighborhood (e.g., Sodroski et al. 1995; Dahmen et al. 1998; Strong et al. 2004). Observations have also found substantially lower αCO in nearby galaxy centers than the Milky Way disk value (e.g., Sandstrom et al. 2013; Israel 2020). The αCO depression in the center reduces the molecular gas concentration compared to that derived using a constant αCO. The degree of central depression in αCO varies from galaxy to galaxy, imposing different level of influence on molecular gas concentrations. To understand the effect of a variable αCO on molecular gas concentrations, we follow the suggestion in Bolatto et al. (2013) to derive the αCO(Z, Σ⋆) according to the local metallicity and surface density (Sect. 3). The resulting αCO(Z, Σ⋆) is on average a factor of 4.5 lower at the center than the outer disk. We use the αCO(Z, Σ⋆) to derive molecular gas concentrations (Cmol). We find that the CO luminosity concentrations (CCO) tightly correlate with Cmol (Fig. 4), but Cmol are systematically 0.22 dex lower. By accounting for the primary dependences of αCO on galaxy local properties, we confirm that high Cmol are more common in barred than unbarred galaxies and that barred galaxies with a stronger bar tend to have higher Cmol (Fig. 8). A variable αCO therefore does not affect the conclusions drawn in Sakamoto et al. (1999a), Sheth et al. (2005) and Kuno et al. (2007).
Nevertheless, the KS test shows that Cmol in barred galaxies are not statistically significantly different from Cmol in unbarred galaxies. This is because there are some unbarred galaxies with relatively high Cmol. To thoroughly explain why the high Cmol are more common in barred than in unbarred spirals, the spiral effect should be taken into account.
5.2. Spiral-driven gas inflow
Gas clouds are subjected to large-scale shocks upon passage across spiral arms that can trigger gravitational collapse and accelerate the formation of new stars (Roberts 1969; Kalnajs 1972; Roberts 1972), which successfully explains the observed bright SF regions along the arms. Hydrodynamical simulations of isolated disk with a prior input spiral density wave confirm the spiral shock on gas (Kim & Kim 2014; Kim et al. 2014, 2020; Baba et al. 2016; Pettitt et al. 2020). The influence of spiral arms on gas clouds is implicit in the Silk-Elmegreen star formation law (Elmegreen 1997; Silk 1997), which proposes that multiple passages of gas clouds through the spiral gravitational potential favor the growth of gas clouds via cloud-cloud collision. In agreement with the compression of gas by shocks, Colombo et al. (2014) show that, in M51, the giant molecular clouds along the arms are brighter than those between the arms, and have a higher number density of more massive clouds, suggesting the arms of M51 promoting the cloud formation. Recently, Meidt et al. (2021) show that the logarithm of molecular arm-interarm contrasts, traced by CO emission observation, are higher than that of old stars, traced by 3.6 μm observation, in a correlation steeper than linear even in the presence of weak or flocculent spiral arms. Stronger spirals are expected to trigger stronger shocks (Roberts 1969; Kim & Kim 2014). Consistently, Seigar & James (2002), Kendall et al. (2015), and Yu et al. (2021) have found that galaxies with strong spiral arms are found to have more intense global SFR or specific SFR than galaxies with weak arms.
Studies on spatially resolved star formation efficiency (SFE) present mixed results. The SFE across the density wave arms of M51 is reduced due to the streaming motion that introduces an additional stabilization effect on gas clouds (Meidt et al. 2013). The SFE enhancement along the arms of NGC 6946 is relatively well defined (Rebolledo et al. 2012), while the enhancement is absence in NGC 628 (Kreckel et al. 2016; see also Foyle et al. 2010). Recently, Querejeta et al. (2021) show that spiral arms accumulate gas and star formation, without systematically increasing the SFE. Unlike the case of M51 (Colombo et al. 2014), the giant molecular cloud (GMC) properties change slightly but not significantly between arms and inter-arms in noninteracting galaxies (Rosolowsky et al. 2021; Colombo et al. 2022). Although the spiral arms of our Galaxy seem have little or no effect on the enhancement of the SFE (Eden et al. 2015; Ragan et al. 2016; Urquhart et al. 2021), evidence that supports the arms influencing the dense gas distribution or kinematics has been reported (Sakai et al. 2015; Urquhart et al. 2021). The complexity of star formation in the arms and interarm regions may be attributed to the turbulent and streaming motions induced when the shock increases the gas density, which prevents gas cloud collapse and reduces the SFE (Meidt et al. 2013; Leroy et al. 2017; Yajima et al. 2019; Querejeta et al. 2019). Spiral shocks are also essential for explaining the azimuthal offset between arm ridges of different components or for explaining the smaller winding angles of younger stellar arms (Grosbol & Patsis 1998; Gittins & Clarke 2004; Egusa et al. 2009; Martínez-García et al. 2014; Yu & Ho 2018).
The funneling of the gas toward the center is expected, as the shock causes the gas to lose angular momentum (Kalnajs 1972; Roberts 1972; Lubow et al. 1986). Lubow et al. (1986) included the back reaction of gas on stellar disk to study spiral shocks and found that parameters representing the solar neighborhood produced a gas accretion rate of ∼(0.2–0.4) M⊙ yr−1, which is consistent with mass drift derived from chemical modeling in the Milky Way (Lacey & Fall 1985) and consistent with those in nearby galaxies derived based on gravitational torques (García-Burillo et al. 2009; Haan et al. 2009). The mass inflow is actually caused by a combination of three processes: angular momentum loss at spiral shocks (∼50% contribution), stellar spiral gravitational torque (∼40% contribution), and gaseous spiral gravitational torque (∼10% contribution; Kim & Kim 2014). By examining the nonlinear gas responses to an imposed stellar spiral potential with various properties in thin disks, Kim & Kim (2014) show that the arms drive mass inflows at rates of ∼0.05–3.0 M⊙ yr−1, with larger values corresponding to stronger and more slowly rotating arms. Although study of individual galaxies has shown that a few unbarred spiral galaxies have relatively high gas concentration (Sheth et al. 2005; Regan et al. 2006), the effect of spiral arms did not receive enough attention. By using a statistical large sample selected from the EDGE, we show that unbarred galaxies with stronger spiral arms tend to have more centrally concentrated molecular gas distribution, providing significant evidence for the spiral-driven transport of molecular gas to the central 2 kpc of galaxies. We verify that it is not caused by the method adopted to identify the bars or induced by central bulges.
5.3. Gas inflow driven by both spirals and bars
By quantifying the strength of the disk (spiral+bar) structure, we find that the molecular gas concentration increases as the structure strengthens, regardless of whether they are unbarred or barred (Fig. 8). Despite the higher gas concentrations and stronger structure found in barred galaxies, there is no significant difference between the correlations for barred and unbarred spirals, implying that the gas concentrations may be uniformly driven by shocks and torques produced by non-axisymmetric disk (spiral+bar) structures. The scatter in molecular gas concentration for given structure strength is ∼0.3 dex. The scatter arises from several aspects. There are degeneracies between spiral strength and spiral pattern speed and degeneracies between bar strength and bar orbit morphologies. For a given spiral strength, spiral arms with a slower pattern speed trigger a stronger shock and a higher gas infall rate (Kim & Kim 2014), and a higher gas concentration is expected. For a given bar strength, a bar could have different orbit morphologies, leading to various gas inflow rates and hence a variation in gas concentrations (Regan & Teuben 2004). The inflow of gas to the center may trigger a starburst at the center, which in turn consumes gas in a short timescale. This occurs periodically and is regulated by a balance between the inflow rate and central accumulation of the gas (Krumholz & Kruijssen 2015), perhaps leading to periodic changes in the molecular gas concentrations. The AGN-driven outflow, if present, may remove central gas content to reduce gas concentrations and cause some degree of variation in the CCO- or Cmol-(2 + log sdisk) relation. The possible AGN feedback is discussed in Sect. 5.4.
We reinterpret the previous finding that barred galaxies tend to have higher molecular gas concentrations than unbarred galaxies (Sakamoto et al. 1999a; Sheth et al. 2005; Kuno et al. 2007; Komugi et al. 2008). By assuming that the spiral- and bar-driven instabilities are the key drivers of gas inflow, the finding should be traced to the fact that barred spirals generally have stronger disk non-axisymmetric structure than unbarred spirals (Fig. 7) and therefore have more gas transported to the center (Fig. 8), leading to higher molecular gas concentrations.
5.4. Gas concentrations in AGN hosts
AGNs may transfer radiative or mechanical energy to the surrounding gas and blow cold gas out of the galaxy (e.g., Silk & Rees 1998; Croton et al. 2006; Zinger et al. 2020). However, studies of low redshift AGNs have shown that, compared to inactive galaxies of similar mass and morphology, the AGN hosts have similar amounts of atomic (Ho et al. 2008; Fabello et al. 2011; Ellison et al. 2019), molecular (Maiolino et al. 1997; Saintonge et al. 2017; Shangguan et al. 2020; Yesuf & Ho 2020; Koss et al. 2021), and total gas content (Vito et al. 2014; Shangguan et al. 2018; Shangguan & Ho 2019). These results suggest that low redshift AGN feedback is inefficient at removing gas, or, alternatively, AGN feedback effects may be limited to the galaxy central regions. Outflows within a few kiloparsecs from the center have been identified in AGN hosts (Karouzos et al. 2016; Kang & Woo 2018; Fluetsch et al. 2019), although the contribution of AGN outflow to feedback is unclear (Woo et al. 2017). Ellison et al. (2021) show that the gas fractions of central AGN regions are lower than those in SF regions and they suggest that the expelling of gas by AGN feedback may happen on central sub-galactic kiloparsec scales instead of global scale. If this is true, AGN hosts should have notably lower gas concentrations compared to inactive galaxies of similar mass or morphology.
We test the effect of AGNs on the radial distribution of molecular gas using the CCO- and Cmol-(2 + log sdisk) relations. In Fig. 10 we plot the CCO- and Cmol-(2 + log sdisk) relations separately for 30 pure SF, 15 composite, and 4 Seyfert galaxies. The solid lines in each panel all mark the best-fit linear function of the relation for pure SF galaxies. The CCO- and Cmol-(2 + log sdisk) relations for pure SF galaxies are consistent within uncertainty with those constructed using the full sample in Fig. 8. The composite galaxies have similar ranges of CCO and Cmol, and do not have significantly lower CCO or Cmol compared to that in SF galaxies. They lie around the SF best-fit functions.
![]() |
Fig. 10. Dependence of CCO and Cmol on 2 + log sdisk for pure SF, composite, and Seyfert galaxies. The solid line marks the best-fit linear function of the relation constructed using the sample of pure SF galaxies. |
Among our 4 Seyferts, there is one type 1 AGN (UGC 3973) and three type 2 AGNs (NGC 2410, NGC 2639, NGC 6394), as identified by Lacerda et al. (2020). NGC 2410 and NGC 2639 are strong AGNs (Kalinova et al. 2021). There are one Seyferts lie above the SF best-fit function, and three lie below it. Ellison et al. (2021) show that NGC 2639 has significantly lower local gas fraction in the AGN region than in SF regions. Consistent with their results, this galaxy (bottom-left corner on the diagrams for Seyferts) has the lowest gas concentrations: CCO = −0.22 and Cmol = −0.52. However, it also has almost the weakest structure strength: 2 + log sdisk = 0.83. The reduction of the central gas content in this galaxy may be caused by its AGN feedback (Ellison et al. 2021). However, because NGC 2639 lacks a bar and its spiral arms almost vanish, the inflow of gas driven by spiral arms and bars ceases. The low molecular gas concentration observed in NGC 2639 can therefore be explained by the lack of disk instabilities that transport gas inward. Taking the best-fit function for the pure SF galaxies as a reference, the CCO or Cmol in Seyferts tend to be lower but not significant.
We calculate the difference (ΔCCO) between the measured CCO and the best-fit function constructed using pure SF galaxies at the given measured 2 + log sdisk, and the difference for Cmol (ΔCmol) in the same way. Figure 11 presents the histogram of ΔCCO and ΔCmol for pure SF, composite, and Seyfert galaxies. The KS test between ΔCCO in pure SF and composite galaxies yields a pKS = 0.94, and the same for ΔCmol yields a pKS = 0.63. This suggests that we cannot reject the null hypothesis that the ΔCCO in pure SF and composite galaxies come from the same parent distribution. Likewise, the KS test between ΔCCO in pure SF and Seyfert galaxies yields a pKS = 0.34, and the same for ΔCmol yields a pKS = 0.12. This suggests that we cannot reject the null hypothesis that the ΔCCO in pure SF and Seyfert galaxies come from the same parent distribution. It is important to note that the size of the Seyferts subsample is small, which makes our results less reliable. Future analyses of more Seyferts may strengthen our findings.
![]() |
Fig. 11. Histogram of ΔCCO and ΔCmol for pure SF, composite, and Seyfert galaxies. ΔCCO (ΔCmol) is the difference between the measured CCO (Cmol) and, at the measured strength, the value of the best-fit CCO-(2 + log sdisk) (Cmol-(2 + log sdisk)) relation constructed using pure SF galaxies as shown in Fig. 10. |
We therefore do not reveal a significant drop in gas concentrations of Seyferts compared to non-AGN galaxies of similar strength in our sample. If the shocks and torques produced by non-axisymmetric structures such as spiral arms and bars are key factors to regulate the molecular gas distribution, our results suggest that the AGN feedback in Seyferts do not have a notable effect on the radial distribution of molecular gas on central ∼2 kpc scales. Still, the AGN feedback of Seyferts in our sample may not be powerful enough or the feedback may work on smaller scale, study of which requires higher resolution observation.
5.5. Molecular gas concentration and star formation rate concentration
Star formation on global and sub-galactic scales couples directly with the molecular gas (e.g., Kennicutt 1989; Bigiel et al. 2008; Leroy et al. 2008, 2013; Schruba et al. 2011; Bolatto et al. 2017). A centrally concentrated distribution of SFR is therefore expected for a galaxy with centrally concentrated molecular gas. We focus on the subsample of 38 non-AGN galaxies with robust SFR measurements (Sect. 3.3). In this subsample, AGNs, ambiguous galaxies, and galaxies with central EWHα < 6 Å are excluded. By defining CCO, Cmol, and CSFR over the same radial extent, we find that both CCO and Cmol tightly correlate with CSFR, regardless of whether the galaxies are barred or unbarred (ρ = 0.84–0.91), as illustrated in Fig. 12. The variable αCO(Z, Σ⋆) does not influence the relationship strength between molecular gas concentrations and SFR concentrations.
![]() |
Fig. 12. Dependence of CSFR on CCO and Cmol. AGNs, LINERs and galaxies with a quenched center have been excluded. The Pearson correlation coefficients, ρ, for unbarred, barred, and all galaxies are given at the top of each panel. The formulas indicate the linear fits, derived using PCA and marked with solid lines. |
Chown et al. (2019) find a loose correlation between gas concentrations and central upturn of EWHα profile, a measure of centrally enhanced SFR, in barred galaxies (ρ = 0.64), and cannot find a similar one in unbarred galaxies (ρ = 0.09). The failure to find a tight correlation in Chown et al. (2019) is due to the following two reasons. First, they adopted different radial scales to calculate the two quantities. Their gas concentration was defined over a mean radial extent R50%,CO ≈ 4 kpc, while they measured the upturn of EWHα at R = 0. Second, determining by eye the turning point of the EWHα profile to calculate its central upturn introduces uncertainty. For instance, the EWHα profile of NGC 7819 has at least two turn-up and one turnover, and the turn-up at R ≈ is chosen by Chown et al. (2019). However, this chosen turning is the weakest one among others. We use the change point analysis described in Watkins et al. (2019) to estimate the amplitude of the turning. Instead, if the strongest turning at R ≈ 16″ is used, NGC 7819 has a strong EWHα central upturn ∼1 rather than 0.3 reported by Chown et al. 2019. The EWHα central upturn of ∼1 is consistent with our result that NGC 7819 has a highly centrally concentrated SFR distribution.
By performing PCA, we find the best-fit straight lines for the CCO- and Cmol-CSFR relations (Fig. 12) and the two relations have a small scatter of 0.04 dex. Together with the CCO- and Cmol-(2 + log sdisk) relations (Fig. 8), it is expected to see a connection between the central concentration of SFR and the strength of non-axisymmetric disk structure.
5.6. Central concentrations of SFR driven by spiral arms and bars
We probe the dependence of CSFR on the disk structure strength in Fig. 13. First, similar to CCO and Cmol shown in Fig. 8, CSFR for unbarred galaxies are mainly moderate to weak, whereas those for barred galaxies span a wide range and can reach particularly high values. Highly centrally concentrated SFRs are more common in barred than in unbarred galaxies. Second, the central concentration of SFR increases as the disk (spiral+bar) structure strengthens in both barred and unbarred galaxies. The relation for unbarred and barred galaxies has a correlation coefficient ρ = 0.71 (p < 0.01) and ρ = 0.70 (p < 0.01), respectively. We perform PCA to obtain the best-fit line:
![]() |
Fig. 13. Dependence of central concentration of the SFR (CSFR) on the strength of non-axisymmetric disk structure (2 + log sdisk). AGNs, LINERs and galaxies with a quenched center have been excluded. The Pearson correlation coefficients for unbarred and barred galaxies are given at the top of each panel. The formulas indicate the linear fits, derived using PCA and marked with solid lines. |
Uncertainties are obtained through bootstrapping with replacement. The scatter in CSFR for a given structure strength is 0.18 dex. Bar identification based on isophotal analysis (Appendix B) or 2D decomposition from Méndez-Abreu et al. (2017) yields consistent results, although the correlation coefficients vary slightly (row [3] in Table 1). Morphological quenching proposes that the bulges stabilizes the central gaseous disk to prevent the formation of new stars (Martig et al. 2009). We test the effect of morphological quenching using partial correlation coefficients. The residual dependence of CSFR on 2 + log sdisk after removing the effect of B/T yields ρ′ = 0.62 (p < 0.01) in unbarred galaxies and ρ′ = 0.70 (p < 0.01) in barred galaxies (row [3] in Table 2). The reported correlations are not introduced by bulges or no obvious signature of morphological quenching is detected, perhaps because galaxies with a quenched center have been excluded by requiring central EWHα ≥ 6 Å (Sect. 3.3). After removing the effect of CCO or Cmol, no residual dependence of CSFR on 2 + log sdisk is found (rows [4] and [5] in Table 2). This suggests that the disk instabilities drive gas inflow to the centers, and the increased gas surface density raises central SFR surface density without significantly changing the central SFE. It is consistent with the result that the SFE derived based on molecular gas does not significantly change with radius (Leroy et al. 2008; Muraoka et al. 2019).
Our reported relationships are in agreement with previous findings. The pioneering works by Martinet & Friedli (1997) and Aguerri (1999) found a tight correlation between bar strength and global SF activity in noninteracting barred galaxies. Further, Zhou et al. (2015) and Lin et al. (2017) showed that star formation activity in the center is enhanced in galaxies with stronger bars (see also James et al. 2009). Strong bars may have suppressed central star formation (Wang et al. 2012, 2020; Consolandi et al. 2017; Díaz-García et al. 2020), especially in those massive gas-poor galaxies (Díaz-García et al. 2020) with particularly weak spiral arms (Wang et al. 2020). However, these gas-poor galaxies with star formation that have been quenched are not of interest to the present work because the effects of gas inflow driven by disk instabilities have been removed. Consistent with our results, Yu et al. (2022) used a ratio of central fiber specific SFR acquired from the MPA-JHU catalog (Brinchmann et al. 2004) to global specific SFR derived from Salim et al. (2018) as a measure of central enhancement of star formation, and found enhanced central SFR in galaxies with strong spiral arms.
We show that barred and unbarred galaxies with strong non-axisymmetric disk (spiral+bar) structures tend to have highly concentrated distribution of both molecular gas and SFR, providing evidence for bar-driven and, especially, spiral-driven gas transport to the central region. The accumulation of gas in the center boosts central SFR. It is consistent with the concept of secular evolution, which describes the slow rearrangement of energy and mass resulting from interactions facilitated by non-axisymmetric galaxy structures (Combes & Sanders 1981; Kormendy 1982; Pfenniger & Norman 1990; Sellwood & Wilkinson 1993; Kormendy & Kennicutt 2004). Yu et al. (2022) also shed light on this perspective. They found that the stronger the arm structure, the more enhanced central SFR, the higher the fraction of galaxies hosting a pseudo bulge. The gas inflow triggered by disk instabilities drives secular growth of the central pseudo bulges, and the spirals and bars therefore play an essential role in the galaxy secular evolution.
6. Summary and conclusions
Models and simulations have shown that the instabilities induced by spirals and bars are efficient mechanisms for transporting gas to the central regions, with the inflow rates depending on the perturbation strength (e.g., Athanassoula 1992a,b; Hopkins & Quataert 2011; Kim & Kim 2014). Observations supporting bar-driven gas inflow have shown that barred galaxies tend to have higher molecular gas concentrations than unbarred galaxies (Sakamoto et al. 1999a; Sheth et al. 2005; Kuno et al. 2007; Komugi et al. 2008) and that stronger bars are associated with higher molecular gas concentrations (Kuno et al. 2007). These previous studies used a constant CO-to-H2 conversion factor (αCO) to derive molecular gas concentrations. However, the αCO depends on the physical properties of the environment where the gas clouds are embedded (e.g., Narayanan et al. 2012; Bolatto et al. 2013). It is not clear how a parameter-dependent αCO influences the calculation of molecular gas concentrations and the conclusions drawn in previous studies. Although spiral-driven gas inflow is predicted (Kalnajs 1972; Roberts 1972; Lubow et al. 1986; Kim & Kim 2014), less is known about the spiral effect on the radial distribution of molecular gas.
We used a sample of 57 disk galaxies selected from the EDGE-CALIFA survey to investigate the connection between molecular gas concentrations and non-axisymmetric disk (spiral+bar) structures. The structure strength is defined as the logarithm of the average relative amplitude of spiral arms and bars, and it therefore traces average spiral strength for unbarred galaxies but average spiral and bar strength for barred galaxies. Molecular gas mass surface density is derived using a αCO that decreases with higher metallicity and higher stellar surface density, following the prescription in Bolatto et al. (2013).The αCO at the center is on average a factor of 4.5 lower than in the outer disk. The central concentrations of CO luminosity (CCO) and molecular gas (Cmol) and the SFR (CSFR) are defined, respectively, as the logarithmic ratio of the average central (< 0.2Re) CO intensity, the molecular mass surface density, and the SFR surface density to their disk-averaged (< Re) surface densities. The 0.2Re on average corresponds to ∼1.1 kpc in our sample, and the CCO, Cmol, and CSFR therefore measure concentrations in the central ∼2 kpc. Our results are independent of methods used to identify bars and are not caused by bulges. The main findings are as follows.
-
By construction, Cmol and CCO are tightly correlated, with Cmol 0.22 dex lower on average due to the lower central region value of αCO. The tight correlation between Cmol and CCO implies that several conclusions obtained using a constant αCO in previous studies still stand, although they overestimate central molecular gas concentrations. We therefore confirm the previous finding that high Cmol values are more common in barred galaxies. However, a few unbarred galaxies with strong spiral arms can also have high Cmol.
-
There is a good correlation between Cmol and the strength of non-axisymmetric structure, which can be due to a bar, spiral arms, or both. It supports the idea that central gas concentrations are due to bar-driven and spiral-driven gas transport.
-
Despite the small subsample size, the Cmol of four Seyferts are not significantly reduced compared to inactive galaxies of similar disk structure. This suggests that the AGN feedback in Seyferts may not notably affect the molecular gas distribution in the central ∼2 kpc.
-
Cmol tightly correlates with CSFR. Galaxies with stronger disk (spiral+bar) structure tend to have higher CSFR, regardless of whether they are barred or unbarred. After removing the effect of Cmol, no residual dependence of CSFR on structure strength is found, suggesting that spiral arms and bars transport gas to the centers but have no additional significant effect on further raising the central SFE.
Our results provide significant evidence for bar-driven and, especially, spiral-driven gas transport to the central regions of galaxies. The accumulation of gas in the center increases the central SFR, which likely facilitates the buildup of central disky pseudo bulges. Both spiral arms and bars therefore play an essential role in the secular evolution of disk galaxies.
Acknowledgments
We are grateful to the anonymous referee for helpful feedback. S.Y.Y. is indebted to Karl Menten for his great support during the pandemic. S.Y.Y. acknowledges support by the Alexander von Humboldt Foundation. R.C.L. acknowledges support provided by a National Science Foundation (NSF) Astronomy and Astrophysics Postdoctoral Fellowship under award AST-2102625. L.C.H. was supported by the National Science Foundation of China (11721303, 11991052, 12011540375) and the China Manned Space Project (CMS-CSST-2021-A04 and CMS-CSST-2021-A06). We benefited from discussions with Adam K. Leroy. Support for CARMA construction was derived from the Gordon and Betty Moore Foundation, the Eileen and Kenneth Norris Foundation, the Caltech Associates, the states of California, Illinois, and Maryland, and the NSF. Funding for CARMA development and operations were supported by NSF and the CARMA partner universities.
References
- Aguerri, J. A. L. 1999, A&A, 351, 43 [NASA ADS] [Google Scholar]
- Aguerri, J. A. L., Méndez-Abreu, J., & Corsini, E. M. 2009, A&A, 495, 491 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Alam, S., Albareti, F. D., Allende Prieto, C., et al. 2015, ApJS, 219, 12 [Google Scholar]
- Athanassoula, E. 1984, Phys. Rep., 114, 319 [Google Scholar]
- Athanassoula, E. 1992a, MNRAS, 259, 328 [NASA ADS] [CrossRef] [Google Scholar]
- Athanassoula, E. 1992b, MNRAS, 259, 345 [Google Scholar]
- Athanassoula, E. 2003, MNRAS, 341, 1179 [Google Scholar]
- Athanassoula, E., & Misiriotis, A. 2002, MNRAS, 330, 35 [Google Scholar]
- Athanassoula, E., Romero-Gómez, M., Bosma, A., & Masdemont, J. J. 2009, MNRAS, 400, 1706 [NASA ADS] [CrossRef] [Google Scholar]
- Baba, J. 2015, MNRAS, 454, 2954 [NASA ADS] [CrossRef] [Google Scholar]
- Baba, J., Morokuma-Matsui, K., Miyamoto, Y., Egusa, F., & Kuno, N. 2016, MNRAS, 460, 2472 [NASA ADS] [CrossRef] [Google Scholar]
- Baldwin, J. A., Phillips, M. M., & Terlevich, R. 1981, PASP, 93, 5 [Google Scholar]
- Barazza, F. D., Jogee, S., & Marinova, I. 2008, ApJ, 675, 1194 [Google Scholar]
- Barazza, F. D., Jablonka, P., Desai, V., et al. 2009, A&A, 497, 713 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Barbary, K. 2016, J. Open Source Software, 1, 58 [NASA ADS] [CrossRef] [Google Scholar]
- Bertin, E., & Arnouts, S. 1996, A&AS, 117, 393 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bertin, G., Lin, C. C., Lowe, S. A., & Thurstans, R. P. 1989, ApJ, 338, 78 [CrossRef] [Google Scholar]
- Bigiel, F., Leroy, A., Walter, F., et al. 2008, ApJ, 136, 2846 [Google Scholar]
- Binette, L., Magris, C. G., Stasińska, G., & Bruzual, A. G. 1994, A&A, 292, 13 [NASA ADS] [Google Scholar]
- Bittner, A., Gadotti, D. A., Elmegreen, B. G., et al. 2017, MNRAS, 471, 1070 [NASA ADS] [CrossRef] [Google Scholar]
- Block, D. L., Buta, R., Knapen, J. H., et al. 2004, AJ, 128, 183 [NASA ADS] [CrossRef] [Google Scholar]
- Bolatto, A. D., Wolfire, M., & Leroy, A. K. 2013, ARA&A, 51, 207 [CrossRef] [Google Scholar]
- Bolatto, A. D., Wong, T., Utomo, D., et al. 2017, ApJ, 846, 159 [Google Scholar]
- Brinchmann, J., Charlot, S., White, S. D. M., et al. 2004, MNRAS, 351, 1151 [Google Scholar]
- Buta, R., Vasylyev, S., Salo, H., & Laurikainen, E. 2005, AJ, 130, 506 [NASA ADS] [CrossRef] [Google Scholar]
- Chown, R., Li, C., Athanassoula, E., et al. 2019, MNRAS, 484, 5192 [NASA ADS] [CrossRef] [Google Scholar]
- Cid Fernandes, R., Stasińska, G., Mateus, A., & Vale Asari, N. 2011, MNRAS, 413, 1687 [Google Scholar]
- Colombo, D., Hughes, A., Schinnerer, E., et al. 2014, ApJ, 784, 3 [NASA ADS] [CrossRef] [Google Scholar]
- Colombo, D., Duarte-Cabral, A., Pettitt, A. R., et al. 2022, A&A, 658, A54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Combes, F. 2003, in Active Galactic Nuclei: From Central Engine to Host Galaxy, eds. S. Collin, F. Combes, & I. Shlosman, ASP Conf. Ser., 290, 411 [NASA ADS] [Google Scholar]
- Combes, F., & Sanders, R. H. 1981, A&A, 96, 164 [NASA ADS] [Google Scholar]
- Combes, F., García-Burillo, S., Casasola, V., et al. 2014, A&A, 565, A97 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Conselice, C. J. 2003, ApJS, 147, 1 [NASA ADS] [CrossRef] [Google Scholar]
- Consolandi, G., Dotti, M., Boselli, A., Gavazzi, G., & Gargiulo, F. 2017, A&A, 598, A114 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Croton, D. J., Springel, V., White, S. D. M., et al. 2006, MNRAS, 365, 11 [Google Scholar]
- Dahmen, G., Huttemeister, S., Wilson, T. L., & Mauersberger, R. 1998, A&A, 331, 959 [Google Scholar]
- Debattista, V. P., & Sellwood, J. A. 1998, ApJ, 493, L5 [Google Scholar]
- Díaz-García, S., Salo, H., Laurikainen, E., & Herrera-Endoqui, M. 2016, A&A, 587, A160 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Díaz-García, S., Salo, H., Knapen, J. H., & Herrera-Endoqui, M. 2019, A&A, 631, A94 [CrossRef] [EDP Sciences] [Google Scholar]
- Díaz-García, S., Moyano, F. D., Comerón, S., et al. 2020, A&A, 644, A38 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Durbala, A., Buta, R., Sulentic, J. W., & Verdes-Montenegro, L. 2009, MNRAS, 397, 1756 [NASA ADS] [CrossRef] [Google Scholar]
- Eden, D. J., Moore, T. J. T., Urquhart, J. S., et al. 2015, MNRAS, 452, 289 [NASA ADS] [CrossRef] [Google Scholar]
- Egusa, F., Kohno, K., Sofue, Y., Nakanishi, H., & Komugi, S. 2009, ApJ, 697, 1870 [NASA ADS] [CrossRef] [Google Scholar]
- Ellison, S. L., Brown, T., Catinella, B., & Cortese, L. 2019, MNRAS, 482, 5694 [NASA ADS] [CrossRef] [Google Scholar]
- Ellison, S. L., Wong, T., Sánchez, S. F., et al. 2021, MNRAS, 505, L46 [NASA ADS] [CrossRef] [Google Scholar]
- Elmegreen, B. G. 1997, Rev. Mex. Astron. Astrofis. Conf. Ser., 6, 165 [NASA ADS] [Google Scholar]
- Elmegreen, B. G., Elmegreen, D. M., & Seiden, P. E. 1989, ApJ, 343, 602 [CrossRef] [Google Scholar]
- Elmegreen, B. G., Elmegreen, D. M., Knapen, J. H., et al. 2007, ApJ, 670, L97 [NASA ADS] [CrossRef] [Google Scholar]
- Elmegreen, D. M., Elmegreen, B. G., Yau, A., et al. 2011, ApJ, 737, 32 [NASA ADS] [CrossRef] [Google Scholar]
- Englmaier, P., & Shlosman, I. 2000, ApJ, 528, 677 [NASA ADS] [CrossRef] [Google Scholar]
- Englmaier, P., & Shlosman, I. 2004, ApJ, 617, L115 [NASA ADS] [CrossRef] [Google Scholar]
- Erwin, P. 2018, MNRAS, 474, 5372 [Google Scholar]
- Erwin, P., & Sparke, L. S. 2003, ApJS, 146, 299 [NASA ADS] [CrossRef] [Google Scholar]
- Fabello, S., Kauffmann, G., Catinella, B., et al. 2011, MNRAS, 416, 1739 [NASA ADS] [CrossRef] [Google Scholar]
- Fathi, K., Storchi-Bergmann, T., Riffel, R. A., et al. 2006, ApJ, 641, L25 [NASA ADS] [CrossRef] [Google Scholar]
- Feldmann, R., Gnedin, N. Y., & Kravtsov, A. V. 2012, ApJ, 758, 127 [NASA ADS] [CrossRef] [Google Scholar]
- Fluetsch, A., Maiolino, R., Carniani, S., et al. 2019, MNRAS, 483, 4586 [NASA ADS] [Google Scholar]
- Foyle, K., Rix, H. W., Walter, F., & Leroy, A. K. 2010, ApJ, 725, 534 [NASA ADS] [CrossRef] [Google Scholar]
- Fragkoudi, F., Athanassoula, E., & Bosma, A. 2016, MNRAS, 462, L41 [NASA ADS] [CrossRef] [Google Scholar]
- Gadotti, D. A. 2009, MNRAS, 393, 1531 [Google Scholar]
- Gadotti, D. A. 2011, MNRAS, 415, 3308 [Google Scholar]
- García-Burillo, S., Fernández-García, S., Combes, F., et al. 2009, A&A, 496, 85 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gittins, D. M., & Clarke, C. J. 2004, MNRAS, 349, 909 [CrossRef] [Google Scholar]
- Gomes, J. M., Papaderos, P., Vílchez, J. M., et al. 2016, A&A, 585, A92 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gong, M., Ostriker, E. C., Kim, C.-G., & Kim, J.-G. 2020, ApJ, 903, 142 [CrossRef] [Google Scholar]
- Grosbøl, P., Patsis, P. A., & Pompei, E. 2004, A&A, 423, 849 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Grosbol, P. J., & Patsis, P. A. 1998, A&A, 336, 840 [Google Scholar]
- Haan, S., Schinnerer, E., Emsellem, E., et al. 2009, ApJ, 692, 1623 [NASA ADS] [CrossRef] [Google Scholar]
- Ho, L. C., Filippenko, A. V., & Sargent, W. L. W. 1997, ApJ, 487, 591 [NASA ADS] [CrossRef] [Google Scholar]
- Ho, L. C., Darling, J., & Greene, J. E. 2008, ApJ, 681, 128 [NASA ADS] [CrossRef] [Google Scholar]
- Hopkins, P. F., & Quataert, E. 2011, MNRAS, 415, 1027 [NASA ADS] [CrossRef] [Google Scholar]
- Hubble, E. P. 1926, ApJ, 64, 321 [Google Scholar]
- Israel, F. P. 1997, A&A, 328, 471 [NASA ADS] [Google Scholar]
- Israel, F. P. 2020, A&A, 635, A131 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- James, P. A., Bretherton, C. F., & Knapen, J. H. 2009, A&A, 501, 207 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Jenkins, A., & Binney, J. 1994, MNRAS, 270, 703 [Google Scholar]
- Jogee, S. 2006, in Physics of Active Galactic Nuclei at all Scales, ed. D. Alloin, 693, 143 [NASA ADS] [CrossRef] [Google Scholar]
- Kalinova, V., Colombo, D., Sánchez, S. F., et al. 2021, A&A, 648, A64 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kalnajs, A. J. 1972, Astrophys. Lett., 11, 41 [NASA ADS] [Google Scholar]
- Kang, D., & Woo, J.-H. 2018, ApJ, 864, 124 [NASA ADS] [CrossRef] [Google Scholar]
- Karouzos, M., Woo, J.-H., & Bae, H.-J. 2016, ApJ, 819, 148 [NASA ADS] [CrossRef] [Google Scholar]
- Kauffmann, G., Heckman, T. M., Tremonti, C., et al. 2003, MNRAS, 346, 1055 [Google Scholar]
- Kendall, S., Kennicutt, R. C., & Clarke, C. 2011, MNRAS, 414, 538 [NASA ADS] [CrossRef] [Google Scholar]
- Kendall, S., Clarke, C., & Kennicutt, R. C. 2015, MNRAS, 446, 4155 [NASA ADS] [CrossRef] [Google Scholar]
- Kennicutt, R. C. Jr. 1989, ApJ, 344, 685 [CrossRef] [Google Scholar]
- Kennicutt, R. C., & Evans, N. J. 2012, ARA&A, 50, 531 [NASA ADS] [CrossRef] [Google Scholar]
- Kennicutt, R. C. Jr., Calzetti, D., Walter, F., et al. 2007, ApJ, 671, 333 [NASA ADS] [CrossRef] [Google Scholar]
- Kewley, L. J., Dopita, M. A., Sutherland, R. S., Heisler, C. A., & Trevena, J. 2001, ApJ, 556, 121 [Google Scholar]
- Kewley, L. J., Groves, B., Kauffmann, G., & Heckman, T. 2006, MNRAS, 372, 961 [Google Scholar]
- Kim, Y., & Kim, W.-T. 2014, MNRAS, 440, 208 [NASA ADS] [CrossRef] [Google Scholar]
- Kim, W.-T., Seo, W.-Y., & Kim, Y. 2012, ApJ, 758, 14 [NASA ADS] [CrossRef] [Google Scholar]
- Kim, W.-T., Kim, Y., & Kim, J.-G. 2014, ApJ, 789, 68 [NASA ADS] [CrossRef] [Google Scholar]
- Kim, W.-T., Kim, C.-G., & Ostriker, E. C. 2020, ApJ, 898, 35 [Google Scholar]
- Komugi, S., Sofue, Y., Kohno, K., et al. 2008, ApJS, 178, 225 [NASA ADS] [CrossRef] [Google Scholar]
- Kormendy, J. 1982, ApJ, 257, 75 [NASA ADS] [CrossRef] [Google Scholar]
- Kormendy, J., & Bender, R. 2012, ApJS, 198, 2 [Google Scholar]
- Kormendy, J., & Kennicutt, R. C. Jr. 2004, ARA&A, 42, 603 [NASA ADS] [CrossRef] [Google Scholar]
- Koss, M., Mushotzky, R., Veilleux, S., & Winter, L. 2010, ApJ, 716, L125 [Google Scholar]
- Koss, M. J., Strittmatter, B., Lamperti, I., et al. 2021, ApJS, 252, 29 [CrossRef] [Google Scholar]
- Kreckel, K., Blanc, G. A., Schinnerer, E., et al. 2016, ApJ, 827, 103 [NASA ADS] [CrossRef] [Google Scholar]
- Krumholz, M. R., & Kruijssen, J. M. D. 2015, MNRAS, 453, 739 [Google Scholar]
- Kuno, N., Sato, N., Nakanishi, H., et al. 2007, PASJ, 59, 117 [NASA ADS] [Google Scholar]
- Lacerda, E. A. D., Cid Fernandes, R., Couto, G. S., et al. 2018, MNRAS, 474, 3727 [Google Scholar]
- Lacerda, E. A. D., Sánchez, S. F., Cid Fernandes, R., et al. 2020, MNRAS, 492, 3073 [Google Scholar]
- Lacey, C. G., & Fall, S. M. 1985, ApJ, 290, 154 [NASA ADS] [CrossRef] [Google Scholar]
- Laine, S., Shlosman, I., Knapen, J. H., & Peletier, R. F. 2002, ApJ, 567, 97 [Google Scholar]
- Laurikainen, E., & Salo, H. 2002, MNRAS, 337, 1118 [Google Scholar]
- Laurikainen, E., Salo, H., Buta, R., & Vasylyev, S. 2004, MNRAS, 355, 1251 [Google Scholar]
- Laurikainen, E., Salo, H., Buta, R., & Knapen, J. H. 2007, MNRAS, 381, 401 [Google Scholar]
- Leroy, A. K., Walter, F., Brinks, E., et al. 2008, AJ, 136, 2782 [Google Scholar]
- Leroy, A. K., Walter, F., Bigiel, F., et al. 2009, AJ, 137, 4670 [Google Scholar]
- Leroy, A. K., Walter, F., Sandstrom, K., et al. 2013, AJ, 146, 19 [Google Scholar]
- Leroy, A. K., Schinnerer, E., Hughes, A., et al. 2017, ApJ, 846, 71 [Google Scholar]
- Leroy, A. K., Schinnerer, E., Hughes, A., et al. 2021, ApJS, 257, 43 [NASA ADS] [CrossRef] [Google Scholar]
- Levy, R. C., Bolatto, A. D., Sánchez, S. F., et al. 2019, ApJ, 882, 84 [NASA ADS] [CrossRef] [Google Scholar]
- Li, Z.-Y., Ho, L. C., Barth, A. J., & Peng, C. Y. 2011, ApJS, 197, 22 [Google Scholar]
- Lin, C. C., & Shu, F. H. 1964, ApJ, 140, 646 [Google Scholar]
- Lin, L., Li, C., He, Y., Xiao, T., & Wang, E. 2017, ApJ, 838, 105 [NASA ADS] [CrossRef] [Google Scholar]
- Lin, L., Li, C., Du, C., et al. 2020, MNRAS, 499, 1406 [NASA ADS] [CrossRef] [Google Scholar]
- Lubow, S. H., Balbus, S. A., & Cowie, L. L. 1986, ApJ, 309, 496 [NASA ADS] [CrossRef] [Google Scholar]
- Maiolino, R., Ruiz, M., Rieke, G. H., & Papadopoulos, P. 1997, ApJ, 485, 552 [NASA ADS] [CrossRef] [Google Scholar]
- Makarov, D., Prugniel, P., Terekhova, N., Courtois, H., & Vauglin, I. 2014, A&A, 570, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Marino, R. A., Rosales-Ortega, F. F., Sánchez, S. F., et al. 2013, A&A, 559, A114 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Marinova, I., & Jogee, S. 2007, ApJ, 659, 1176 [Google Scholar]
- Martig, M., Bournaud, F., Teyssier, R., & Dekel, A. 2009, ApJ, 707, 250 [Google Scholar]
- Martin, P. 1995, AJ, 109, 2428 [NASA ADS] [CrossRef] [Google Scholar]
- Martinet, L., & Friedli, D. 1997, A&A, 323, 363 [NASA ADS] [Google Scholar]
- Martínez-García, E. E., Puerari, I., Rosales-Ortega, F. F., et al. 2014, ApJ, 793, L19 [CrossRef] [Google Scholar]
- Meidt, S. E., Schinnerer, E., García-Burillo, S., et al. 2013, ApJ, 779, 45 [NASA ADS] [CrossRef] [Google Scholar]
- Meidt, S. E., Leroy, A. K., Querejeta, M., et al. 2021, ApJ, 913, 113 [NASA ADS] [CrossRef] [Google Scholar]
- Méndez-Abreu, J., Ruiz-Lara, T., Sánchez-Menguiano, L., et al. 2017, A&A, 598, A32 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Menéndez-Delmestre, K., Sheth, K., Schinnerer, E., Jarrett, T. H., & Scoville, N. Z. 2007, ApJ, 657, 790 [Google Scholar]
- Muñoz-Mateos, J. C., Sheth, K., Regan, M., et al. 2015, ApJS, 219, 3 [Google Scholar]
- Mundell, C. G., & Shone, D. L. 1999, MNRAS, 304, 475 [NASA ADS] [CrossRef] [Google Scholar]
- Muraoka, K., Sorai, K., Miyamoto, Y., et al. 2019, PASJ, 71, S15 [NASA ADS] [CrossRef] [Google Scholar]
- Narayanan, D., Krumholz, M. R., Ostriker, E. C., & Hernquist, L. 2012, MNRAS, 421, 3127 [NASA ADS] [CrossRef] [Google Scholar]
- Pearson, K. 1901, London. Edinburgh Dublin Philos. Mag. J. Sci., 2, 559 [Google Scholar]
- Pettitt, A. R., Dobbs, C. L., Baba, J., et al. 2020, MNRAS, 498, 1159 [Google Scholar]
- Pfenniger, D., & Norman, C. 1990, ApJ, 363, 391 [Google Scholar]
- Querejeta, M., Meidt, S. E., Schinnerer, E., et al. 2016, A&A, 588, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Querejeta, M., Schinnerer, E., Schruba, A., et al. 2019, A&A, 625, A19 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Querejeta, M., Schinnerer, E., Meidt, S., et al. 2021, A&A, 656, A133 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ragan, S. E., Moore, T. J. T., Eden, D. J., et al. 2016, MNRAS, 462, 3123 [NASA ADS] [CrossRef] [Google Scholar]
- Rebolledo, D., Wong, T., Leroy, A., Koda, J., & Donovan Meyer, J. 2012, ApJ, 757, 155 [CrossRef] [Google Scholar]
- Regan, M. W., & Teuben, P. J. 2004, ApJ, 600, 595 [NASA ADS] [CrossRef] [Google Scholar]
- Regan, M. W., Vogel, S. N., & Teuben, P. J. 1997, ApJ, 482, L143 [NASA ADS] [CrossRef] [Google Scholar]
- Regan, M. W., Sheth, K., & Vogel, S. N. 1999, ApJ, 526, 97 [NASA ADS] [CrossRef] [Google Scholar]
- Regan, M. W., Thornley, M. D., Vogel, S. N., et al. 2006, ApJ, 652, 1112 [NASA ADS] [CrossRef] [Google Scholar]
- Rix, H.-W., & Zaritsky, D. 1995, ApJ, 447, 82 [Google Scholar]
- Roberts, W. W. 1969, ApJ, 158, 123 [NASA ADS] [CrossRef] [Google Scholar]
- Roberts, W. W. Jr., & Shu, F. H. 1972, Astrophys. Lett., 12, 49 [NASA ADS] [Google Scholar]
- Rosolowsky, E., Hughes, A., Leroy, A. K., et al. 2021, MNRAS, 502, 1218 [NASA ADS] [CrossRef] [Google Scholar]
- Saintonge, A., Catinella, B., Tacconi, L. J., et al. 2017, ApJS, 233, 22 [Google Scholar]
- Sakai, N., Nakanishi, H., Matsuo, M., et al. 2015, PASJ, 67, 69 [NASA ADS] [CrossRef] [Google Scholar]
- Sakamoto, K., Okumura, S. K., Ishizuki, S., & Scoville, N. Z. 1999a, ApJS, 124, 403 [NASA ADS] [CrossRef] [Google Scholar]
- Sakamoto, K., Okumura, S. K., Ishizuki, S., & Scoville, N. Z. 1999b, ApJ, 525, 691 [NASA ADS] [CrossRef] [Google Scholar]
- Salim, S., Boquien, M., & Lee, J. C. 2018, ApJ, 859, 11 [Google Scholar]
- Salo, H., Laurikainen, E., Buta, R., & Knapen, J. H. 2010, ApJ, 715, L56 [NASA ADS] [CrossRef] [Google Scholar]
- Sánchez, S. F., Kennicutt, R. C., Gil de Paz, A., et al. 2012, A&A, 538, A8 [Google Scholar]
- Sánchez, S. F., Rosales-Ortega, F. F., Iglesias-Páramo, J., et al. 2014, A&A, 563, A49 [CrossRef] [EDP Sciences] [Google Scholar]
- Sánchez, S. F., García-Benito, R., Zibetti, S., et al. 2016, A&A, 594, A36 [Google Scholar]
- Sandage, A. 1975, in Galaxies and the Universe, eds. A. Sandage, M. Sandage, & J. Kristian, 1 [Google Scholar]
- Sandstrom, K. M., Leroy, A. K., Walter, F., et al. 2013, ApJ, 777, 5 [Google Scholar]
- Schinnerer, E., Böker, T., Emsellem, E., & Lisenfeld, U. 2006, ApJ, 649, 181 [NASA ADS] [CrossRef] [Google Scholar]
- Schruba, A., Leroy, A. K., Walter, F., et al. 2011, AJ, 142, 37 [NASA ADS] [CrossRef] [Google Scholar]
- Seigar, M. S., & James, P. A. 2002, MNRAS, 337, 1113 [NASA ADS] [CrossRef] [Google Scholar]
- Sellwood, J. A. 1980, A&A, 89, 296 [NASA ADS] [Google Scholar]
- Sellwood, J. A., & Wilkinson, A. 1993, Rep. Prog. Phys., 56, 173 [NASA ADS] [CrossRef] [Google Scholar]
- Shangguan, J., & Ho, L. C. 2019, ApJ, 873, 90 [NASA ADS] [CrossRef] [Google Scholar]
- Shangguan, J., Ho, L. C., & Xie, Y. 2018, ApJ, 854, 158 [NASA ADS] [CrossRef] [Google Scholar]
- Shangguan, J., Ho, L. C., Bauer, F. E., Wang, R., & Treister, E. 2020, ApJ, 899, 112 [NASA ADS] [CrossRef] [Google Scholar]
- Sheth, K., Regan, M. W., Vogel, S. N., & Teuben, P. J. 2000, ApJ, 532, 221 [NASA ADS] [CrossRef] [Google Scholar]
- Sheth, K., Vogel, S. N., Regan, M. W., et al. 2002, AJ, 124, 2581 [NASA ADS] [CrossRef] [Google Scholar]
- Sheth, K., Vogel, S. N., Regan, M. W., Thornley, M. D., & Teuben, P. J. 2005, ApJ, 632, 217 [NASA ADS] [CrossRef] [Google Scholar]
- Sheth, K., Elmegreen, D. M., Elmegreen, B. G., et al. 2008, ApJ, 675, 1141 [Google Scholar]
- Silk, J. 1997, ApJ, 481, 703 [CrossRef] [Google Scholar]
- Silk, J., & Rees, M. J. 1998, A&A, 331, L1 [NASA ADS] [Google Scholar]
- Singh, R., van de Ven, G., Jahnke, K., et al. 2013, A&A, 558, A43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Sodroski, T. J., Odegard, N., Dwek, E., et al. 1995, ApJ, 452, 262 [Google Scholar]
- Sorai, K., Kuno, N., Muraoka, K., et al. 2019, PASJ, 71, S14 [Google Scholar]
- Sormani, M. C., Tress, R. G., Glover, S. C. O., et al. 2020, MNRAS, 497, 5024 [Google Scholar]
- Strong, A. W., Moskalenko, I. V., Reimer, O., Digel, S., & Diehl, R. 2004, A&A, 422, L47 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Teng, Y.-H., Sandstrom, K. M., Sun, J., et al. 2022, ApJ, 925, 72 [NASA ADS] [CrossRef] [Google Scholar]
- Tremaine, S., & Weinberg, M. D. 1984, MNRAS, 209, 729 [Google Scholar]
- Tress, R. G., Sormani, M. C., Glover, S. C. O., et al. 2020, MNRAS, 499, 4455 [NASA ADS] [CrossRef] [Google Scholar]
- Urquhart, J. S., Figura, C., Cross, J. R., et al. 2021, MNRAS, 500, 3050 [Google Scholar]
- Utomo, D., Bolatto, A. D., Wong, T., et al. 2017, ApJ, 849, 26 [NASA ADS] [CrossRef] [Google Scholar]
- Vallat, R. 2018, J. Open Source Software, 3, 1026 [NASA ADS] [CrossRef] [Google Scholar]
- van den Bergh, S. 1976, ApJ, 206, 883 [NASA ADS] [CrossRef] [Google Scholar]
- Vito, F., Maiolino, R., Santini, P., et al. 2014, MNRAS, 441, 1059 [NASA ADS] [CrossRef] [Google Scholar]
- Walcher, C. J., Wisotzki, L., Bekeraité, S., et al. 2014, A&A, 569, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Wang, J., Kauffmann, G., Overzier, R., et al. 2012, MNRAS, 423, 3486 [NASA ADS] [CrossRef] [Google Scholar]
- Wang, J., Athanassoula, E., Yu, S.-Y., et al. 2020, ApJ, 893, 19 [NASA ADS] [CrossRef] [Google Scholar]
- Watanabe, Y., Sorai, K., Kuno, N., & Habe, A. 2011, MNRAS, 411, 1409 [NASA ADS] [CrossRef] [Google Scholar]
- Watkins, A. E., Laine, J., Comerón, S., Janz, J., & Salo, H. 2019, A&A, 625, A36 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Wilson, C. D., Warren, B. E., Israel, F. P., et al. 2012, MNRAS, 424, 3050 [Google Scholar]
- Wolfire, M. G., Hollenbach, D., & McKee, C. F. 2010, ApJ, 716, 1191 [Google Scholar]
- Woo, J.-H., Son, D., & Bae, H.-J. 2017, ApJ, 839, 120 [Google Scholar]
- Yajima, Y., Sorai, K., Kuno, N., et al. 2019, PASJ, 71, S13 [NASA ADS] [CrossRef] [Google Scholar]
- Yesuf, H. M., & Ho, L. C. 2020, ApJ, 900, 107 [NASA ADS] [CrossRef] [Google Scholar]
- York, D. G., Adelman, J., Anderson, J. E. Jr., et al. 2000, AJ, 120, 1579 [Google Scholar]
- Yu, S.-Y., & Ho, L. C. 2018, ApJ, 869, 29 [NASA ADS] [CrossRef] [Google Scholar]
- Yu, S.-Y., & Ho, L. C. 2019, ApJ, 871, 194 [NASA ADS] [CrossRef] [Google Scholar]
- Yu, S.-Y., & Ho, L. C. 2020, ApJ, 900, 150 [NASA ADS] [CrossRef] [Google Scholar]
- Yu, S.-Y., Ho, L. C., Barth, A. J., & Li, Z.-Y. 2018, ApJ, 862, 13 [NASA ADS] [CrossRef] [Google Scholar]
- Yu, S.-Y., Ho, L. C., & Wang, J. 2021, ApJ, 917, 88 [NASA ADS] [CrossRef] [Google Scholar]
- Yu, S.-Y., Xu, D., Ho, L. C., Wang, J., & Kao, W.-B. 2022, A&A, 661, A98 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Yuan, C., & Kuo, C.-L. 1997, ApJ, 486, 750 [NASA ADS] [CrossRef] [Google Scholar]
- Yuan, T. T., Kewley, L. J., & Sanders, D. B. 2010, ApJ, 709, 884 [NASA ADS] [CrossRef] [Google Scholar]
- Zhou, Z.-M., Cao, C., & Wu, H. 2015, AJ, 149, 1 [Google Scholar]
- Zinger, E., Pillepich, A., Nelson, D., et al. 2020, MNRAS, 499, 768 [CrossRef] [Google Scholar]
- Zurita, A., Relaño, M., Beckman, J. E., & Knapen, J. H. 2004, A&A, 413, 73 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
Appendix A: Table of derived quantities
Table A.1 presents the quantities derived in Sect. 3.
Table of derived quantities
Appendix B: Bar identification based on isophotal analysis
We perform an additional bar identification based on isophotes. It has been found that the ϵ profile tends to increase with semimajor axis (SMA) and the PA profile tend to keep constant within the region dominated by a bar, beyond which the ϵ would drop and PA would twist (Athanassoula & Misiriotis 2002; Laine et al. 2002; Erwin & Sparke 2003; Menéndez-Delmestre et al. 2007; Aguerri et al. 2009; Li et al. 2011). We thus search for drops in the ϵ profile and twists in the PA profile as bar candidates. We inspect the images to ensure that dust lanes, star formation, rings, or spiral arms are not responsible for these features in the ϵ and PA profiles. A candidate bar is identified as a true bar only if there is a bar-like structure in the r-band image consistent with its ϵ, PA, and SMA. We give preference to the bar candidate with maximum ϵ. Figure B.1 illustrates the bar identification for NGC 447. The ellipse in the left panel denotes the bar with measured ϵbar, PAbar, and SMAbar; the arrow in the middle and right panel denotes the ϵ drop and PA twist caused by the bar. Then thePAbar and SMAbar are used to calculate the intrinsic bar radius:
where ΔPA = PAdisk − PAbar. The vertical dashed line in the top-right panel of Fig. 2 denotes the Rbar for NGC 447. The isophotal analysis yields a bar fraction of only 42% and it may have missed some weak bars. The isophote-based bar identification is used to understand the uncertainty in our results introduced by using different methods to identify bars (Sect. 4.4).
The ellipticity and relative length are measures of bar strength (Martin 1995; Martinet & Friedli 1997; Aguerri 1999; Marinova & Jogee 2007; Barazza et al. 2008, 2009; Gadotti 2009; Li et al. 2011). It is instructive to compare bar strengths defined in our framework with the measured ϵbar and Rbar/R25. We calculate average relative Fourier amplitude over the region occupied by the bar sbar and then compute the bar strength: 2 + log sbar. Since the ϵbar are projected values, we focus on galaxies with inclination angle less than 65° to avoid projection effects. As shown in Figs. B.2 and B.3, the 2 + log sbar are well correlated with the ϵbar and Rbar/R25. The value of 2 + log sbar gets higher with increasing bar ellipticity (ρ = 0.67 and p < 0.01) and bar length (ρ = 0.87 and p < 0.01).
![]() |
Fig. B.1. Illustration of how to identify a bar using profiles of ϵ and PA. The ellipse in the left panel denotes the bar with measured ϵbar, PAbar and SMAbar. The arrows in the bottom-right and top-right panels denotes the ϵ drop and the PA twist caused by the bar, respectively. |
![]() |
Fig. B.2. Relation between the bar strength (2 + log sbar) and the bar ellipticity (ϵbar). The Pearson correlation coefficient, ρ, and the p value are shown at the top. |
![]() |
Fig. B.3. Relation between the bar strength (2 + log sbar) and the relative length of the bar (Rbar/R25). The Pearson correlation coefficient, ρ, and the p value are shown at the top. |
All Tables
All Figures
![]() |
Fig. 1. Histogram of Hubble types, distance, and total molecular gas mass (Mmol) of galaxies in the EDGE parent sample and our derived sample. |
In the text |
![]() |
Fig. 2. Illustration of how to measure the strength of non-axisymmetric disk structures. The left two panels show the star-cleaned r-band image of NGC 447 and UGC 9476. The dashed blue ellipse in the image of NGC 447 denotes its bar (Appendix B). The right panels show radial profiles of Atot. The galactocentric radius at which the bulge terminates and the disk starts to dominate (1.5 × Rbul/(1 − ϵdisk), where ϵdisk is the disk-averaged ellipticity) is indicated by the dotted line. The bar radius is indicated by the dashed line for NGC 447. The R90%, enclosing the majority of structure, is indicated by the solid line. We compute sdisk as the mean Atot over the disk-dominated region between 1.5 × Rbul/(1 − ϵdisk) and R90%, and define the strength of non-axisymmetric disk structure as 2 + log sdisk. The disk-dominated region for barred galaxies contains spirals and bars, but that for unbarred galaxies contains only spirals. |
In the text |
![]() |
Fig. 3. αCO plotted against galactocentric radius normalized by R25. Top panel:αCO(R) obtained by multiplying the empirical normalized average profile of αCO adopted from Sandstrom et al. (2013) by Milky Way αCO, MW = 4.4 M⊙ pc−2 (K km s−1)−1; the errors are also from Sandstrom et al. (2013) and the black curve is the best-fit function. Bottom panel:αCO(Z, Σ⋆), marked in gray, derived from Eq. (8), which considers the decrease in αCO driven by higher metallicity (Z) and stellar surface density (Σ⋆); the black curve denotes the mean αCO for a given radius. |
In the text |
![]() |
Fig. 4. Correlations between CO luminosity concentration (CCO) and molecular gas mass concentration (Cmol). Top panel:Cmol(αCO[R]) derived based on αCO(R), which is a function of radius and is set the same for all the galaxies (Eq. (6)). Bottom panel:Cmol(αCO[Z, Σ⋆]) computed based on αCO[Z, Σ⋆], which is a function of Z and Σ⋆ and varies from galaxy to galaxy (Eq. (8)). The dashed lines mark the 1 : 1 relations. |
In the text |
![]() |
Fig. 5. Average normalized profiles of CO integrated intensity (ICO(R)/ICO(Re) and ICO(R)/ICO(R80%,CO)) and molecular mass surface density (Σmol(R)/Σmol(Re) and Σmol(R)/Σmol(R80%,mol)). Blue profiles are for high-CCO or high-Cmol galaxies, while red profiles are for low-CCO or low-Cmol galaxies. The shaded regions mark the standard deviation in y-axis value for a given normalized radius. |
In the text |
![]() |
Fig. 6. Histograms of the CO luminosity concentrations (CCO) and molecular gas concentrations (Cmol). The results for barred and unbarred galaxies are marked in black and gray, respectively. The p value of the Kolmogorov-Smirnov test (pKS) is presented. |
In the text |
![]() |
Fig. 7. Histogram of the strength of non-axisymmetric disk (spiral+bar) structures (2 + log sdisk). The 2 + log sdisk traces average spiral strength for unbarred galaxies, but traces average spiral and bar strength for barred galaxies. The results for barred and unbarred galaxies are marked in black and gray, respectively. The p value of the Kolmogorov-Smirnov test (pKS) is presented. |
In the text |
![]() |
Fig. 8. Dependence of CO luminosity concentrations (CCO) and molecular gas concentrations (Cmol) on the strengths of disk structure (2 + log sdisk). The left two panels show the results for the parent sample, whereas the right two panels show the results for galaxies with i ≤ 65°. The formulas indicate linear fits considering both barred and unbarred galaxies. They are derived using PCA and marked with solid lines. The Pearson correlation coefficients, ρ, for unbarred and barred galaxies are given at the top of each panel. |
In the text |
![]() |
Fig. 9. Comparison of CCO and Cmol with Hubble type. Black points show the data. The Hubble types are from Walcher et al. (2014), and the corresponding numerical Hubble stages are used to calculate the Pearson correlation coefficients between the two parameters. |
In the text |
![]() |
Fig. 10. Dependence of CCO and Cmol on 2 + log sdisk for pure SF, composite, and Seyfert galaxies. The solid line marks the best-fit linear function of the relation constructed using the sample of pure SF galaxies. |
In the text |
![]() |
Fig. 11. Histogram of ΔCCO and ΔCmol for pure SF, composite, and Seyfert galaxies. ΔCCO (ΔCmol) is the difference between the measured CCO (Cmol) and, at the measured strength, the value of the best-fit CCO-(2 + log sdisk) (Cmol-(2 + log sdisk)) relation constructed using pure SF galaxies as shown in Fig. 10. |
In the text |
![]() |
Fig. 12. Dependence of CSFR on CCO and Cmol. AGNs, LINERs and galaxies with a quenched center have been excluded. The Pearson correlation coefficients, ρ, for unbarred, barred, and all galaxies are given at the top of each panel. The formulas indicate the linear fits, derived using PCA and marked with solid lines. |
In the text |
![]() |
Fig. 13. Dependence of central concentration of the SFR (CSFR) on the strength of non-axisymmetric disk structure (2 + log sdisk). AGNs, LINERs and galaxies with a quenched center have been excluded. The Pearson correlation coefficients for unbarred and barred galaxies are given at the top of each panel. The formulas indicate the linear fits, derived using PCA and marked with solid lines. |
In the text |
![]() |
Fig. B.1. Illustration of how to identify a bar using profiles of ϵ and PA. The ellipse in the left panel denotes the bar with measured ϵbar, PAbar and SMAbar. The arrows in the bottom-right and top-right panels denotes the ϵ drop and the PA twist caused by the bar, respectively. |
In the text |
![]() |
Fig. B.2. Relation between the bar strength (2 + log sbar) and the bar ellipticity (ϵbar). The Pearson correlation coefficient, ρ, and the p value are shown at the top. |
In the text |
![]() |
Fig. B.3. Relation between the bar strength (2 + log sbar) and the relative length of the bar (Rbar/R25). The Pearson correlation coefficient, ρ, and the p value are shown at the top. |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.