Abstract
Bars may induce morphological features, such as rings, through their resonances. Previous studies suggested that the presence of "dark gaps," or regions of a galaxy where the difference between the surface brightness along the bar major axis and that along the bar minor axis is maximal, can be attributed to the location of bar corotation. Here, using GALAKOS, a high-resolution N-body simulation of a barred galaxy, we test this photometric method's ability to identify the bar corotation resonance. Contrary to previous work, our results indicate that "dark gaps" are a clear sign of the location of the 4:1 ultraharmonic resonance instead of bar corotation. Measurements of the bar corotation can indirectly be inferred using kinematic information, e.g., by measuring the shape of the rotation curve. We demonstrate our concept on a sample of 578 face-on barred galaxies with both imaging and integral field observations and find that the sample likely consists primarily of fast bars.
Export citation and abstract BibTeX RIS
1. Introduction
Observations reveal that 50% to over 70% of nearby disk galaxies host a bar (e.g., Eskridge et al. 2000; Aguerri et al. 2009; Nair & Abraham 2010; Masters et al. 2011), including our own Milky Way (Binney et al. 1991; Weiland et al. 1994; Hammersley et al. 2001; Benjamin et al. 2005). Bars are long-lived features with significant effects on the stellar and gas distribution and kinematics throughout galaxies (Athanassoula et al. 2005; Barazza et al. 2008; Sheth et al. 2008). To inform a complete picture of galaxy evolution, we must understand the formation of bars and their impact on a galaxy's dynamics and evolutionary track. Bars form both as a result of interactions (Noguchi 1987; Elmegreen et al. 1991; Romano-Díaz et al. 2008; Martinez-Valpuesta et al. 2016) and spontaneously from gravitational instabilities (Toomre 1964; Polyachenko 2013). The strong nonaxisymmetric nature of bars allows for the transport of angular momentum, energy, and mass across large radial scales within a galaxy (e.g., Debattista & Sellwood 2000; Athanassoula 2003).
Additionally, bars are thought to rotate as a solid body with an angular velocity called the pattern speed (Ωp ) and can thus drive resonances with the rotating matter in the disk beyond the physical extent of the bar itself. The pattern speed of bars is one of the fundamental parameters necessary to understand the dynamics of disk galaxies. Corotation is defined as the radius at which the galaxy's circular angular rotation is equal to the bar pattern speed. The bar pattern speed also often determines the locations of the 2:1 Lindblad resonance (Lindblad 1941) and the 4:1 ultraharmonic resonance (UHR), which can both induce the formation of structures, such as rings (e.g., Buta 1986).
These dynamical parameters and bar resonance locations typically require stellar kinematics measurements to derive, with Tremaine & Weinberg (1984) outlining the most accurate and only model-independent method. Some attempts to tie morphological features, such as rings, to the locations of resonances have been made (e.g., Buta 1986, 2017). The locations of such features could then be used to infer resonance locations and pattern speeds but are model dependent. In particular, Buta (2017) proposed that the location of "dark gaps," or regions of a galaxy along the bar minor axis that are dimmer in surface brightness than the equivalent radius along the bar major axis, corresponds to the location of bar corotation.
Buta (2017) used multiband images of ringed, barred galaxies to identify these "dark gaps" and compared the results with an N-body simulation from Schwarz (1984). "Dark gaps" were attributed to a deficiency of material toward L4,5, the often stable Lagrangian points arising from the gravitational potential of a bar that may become unstable in the presence of a strong bar (Binney & Tremaine 1987). Buta (2017) argues that identifying gaps in the galaxy images along the bar minor axis should correspond to the L4,5 Lagrange points, where unstable "banana" orbits of the stars result in a decrease in the number density of stars surrounding these points. Kim et al. (2016) analyzed the presence of dark gaps with other bar properties and discovered that the light deficit produced by dark gaps was larger in longer and stronger bars. They interpreted these findings as evidence for dark gaps produced by the redistribution of mass due to bar evolution processes.
Here we use state-of-the-art Milky Way–like N-body simulation results to test this proposed method directly on a more general case of a barred galaxy. Using the simulations, we construct synthetic, face-on maps of the stellar surface density and identify bar resonance locations with these gaps. We show that the location of "dark gaps" in this nonringed galaxy does not correspond to the L4,5 Lagrangian points but instead occurs close to the location of the 4:1 ultraharmonic resonance.
In this work, an alternative to the Buta (2017) corotation interpretation of dark gaps is suggested based on results of high-resolution N-body simulations (GALAKOS; D'Onghia & Aguerri 2020). This paper is structured as follows. Section 2 briefly describes the numerical simulation and observational data used to calibrate and apply this method. Section 3 describes our new method and tests its accuracy across different physical resolutions and epochs of the simulation. In Section 4, we apply this method to a sample of 578 nearby galaxies to analyze the bar dynamics in this population of barred galaxies. Lastly, we discuss the broad implications of our results in comparison to other methods in Section 5, and we summarize in Section 6.
2. Data
We use results of an N-body simulation of a Milky Way–like galaxy described in D'Onghia & Aguerri (2020) and data associated with the Sloan Digital Sky Survey (SDSS; York et al. 2000). The details of these data are described below. Throughout this work, we assume a standard cosmology of WMAP9 (Hinshaw et al. 2013) implemented via astropy 15 (Robitaille et al. 2013).
2.1. N-body Simulations
GALAKOS is a ∼90 million N-body simulation run with GADGET3. The numerical experiment follows the motion of disk stars in a Milky Way–like galaxy. The stellar disk has structural parameters configured to correspond with the current Milky Way (Bland-Hawthorn & Gerhard 2016). The simulation develops self-consistent spiral patterns and a stellar bar that arises spontaneously from the stars themselves (D'Onghia et al. 2013). After 2.5 Gyr of evolution, the simulated galaxy develops prominent spiral structure and a bar with 4.5 kpc length but does not develop a ring.
The GADGET3 cosmological code allows for fast computation of the long-range gravitational field on a particle mesh, with short-range forces calculated on a tree-based hierarchical multipole expansion. A spline kernel with scale length, hs
, softens pairwise particle interactions so that interactions beyond that scale are strictly Newtonian. This is similar to Plummer softening with a scale length = hs
/2.8. In this simulation, hs
= 40, 28, and 80 pc for the dark matter halo, stellar disk, and bulge, respectively. Full details of this simulation and the code used can be found in D'Onghia & Aguerri (2020, and their Appendix). Our interpretation of the bar and spiral structure properties is based on the density wave theory and is supported by our numerical simulations (see D'Onghia et al. 2013; D'Onghia & Aguerri 2020).
2.2. SDSS-IV MaNGA
The SDSS-IV (Blanton et al. 2017) Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey observes nearly 10,000 galaxies with integral field spectroscopy (Bundy et al. 2015). MaNGA observations employ the BOSS spectrograph on the 2.5 m telescope at Apache Point Observatory (Gunn et al. 2006), achieving a spectral resolution of R ∼ 2000 for 3600 Å < λ < 10300 Å and a typical (S/N) level of 36 (20) at a fiducial fiber magnitude of 21 (22) in the red (blue) spectrograph arm (Bundy et al. 2015). Optical fibers subtend 2'' on the sky (Smee et al. 2013) and are bundled into integral field units (IFUs) with sizes in the range 12''–32'' in diameter (19–127 fibers; Drory et al. 2015). Sky subtraction and flux calibration are accomplished using simultaneous observations of the sky and standard stars (Yan et al. 2016). The median point-spread function of the resulting data cubes is 25 and roughly corresponds to kiloparsec physical scales at the targeted redshift range (0.01 < z < 0.15). Observations are dithered and mapped onto 0
5 spectroscopic pixels (or spaxels).
The MaNGA sample is selected to have a flat distribution of i-band absolute magnitude and uniform radial coverage. This parent sample is composed of three main components, a primary sample where 80% of galaxies are covered out to 1.5 Re , a secondary sample where 80% of galaxies are covered out to 2.5 Re , and a color-enhanced supplement to improve coverage of poorly sampled regions of the near-UV − i versus Mi color–magnitude plane (Wake et al. 2017). All MaNGA data in this work are reduced using v2_5_3 of the MaNGA Data Reduction Pipeline (Law et al. 2016) and employ the Data Analysis Pipeline (Belfiore et al. 2019; Westfall et al. 2019) from the internal eighth MaNGA product launch (MPL-8), containing 6507 galaxies.
2.2.1. Barred Galaxy Sample—Galaxy Zoo:3D
Galaxy Zoo:3D (Masters et al. 2021) is a citizen science project in which masks were drawn onto galaxy images to outline morphological features such as spiral arms, bars, galaxy centers, and foreground stars. After at least 15 citizen scientists have classified a galaxy, the combined masks allow for quick identification of the stellar bar and estimates of its length and position angle. These masks effectively separate MaNGA spaxels dominated by bar light from the rest of the galaxy. Spaxels are considered to be in the bar if so flagged by 40% of participants. This threshold is stricter than those used in previous work (Fraser-McKelvie et al. 2019; Krishnarao et al. 2020b) but provides more accurate bar length measurements (see Krishnarao et al. 2020b, their Appendix A). Details of this citizen science project can be found on their website. 16
Our barred galaxy sample is composed of 578 face-on, barred galaxies identified in Galaxy Zoo:3D, with minor-to-major-axis ratios b/a > 0.6 as estimated in the NASA-Sloan Atlas (NSA; Blanton et al. 2011).
The bar's orientation is determined by finding a minimum area bounding box around bar spaxels in the deprojected galaxy plane. Each galaxy is also recentered based on this bounding box, and the bar length is determined using the length of this box. Additionally, classifications from Galaxy Zoo (Lintott et al. 2011) are used to match the orientation of spiral arms in these galaxies such that the counterclockwise direction is the direction of rotation of the spiral arm pattern, assuming that all galaxies have trailing spiral arms. Figure 1 shows a polar plot of all spaxels within the bar of the entire sample, with the bar axis oriented horizontally.
Figure 1. 2D Gaussian kernel density estimate of all bar masks in polar coordinates, normalized by bar radius.
Download figure:
Standard image High-resolution image2.2.2. PCA Resolved Galaxy Parameters
We use i-band stellar mass-to-light ratio maps for MPL-8 galaxies derived using a principal component analysis (PCA) method of stellar continuum fitting (Pace et al. 2019a, 2019b). These maps are found to be robust across a wide range of signal-to-noise ratios (2 ≤ S/N ≤ 30) and the full range of realistic stellar metallicities and foreground dust attenuations (τ ≲ 4). These mass-to-light maps allow maps of the resolved stellar mass surface density to be created after deprojecting the MaNGA spaxels using the minor-to-major-axis ratio from the NSA (Blanton et al. 2011) and have been previously used in Krishnarao et al. (2020b) and Schaefer et al. (2019).
3. Identifying Bar Resonances
Buta (2017) attributed "dark gaps" along the bar minor axis to the location of bar corotation. These "dark gaps" are defined as where the difference in the surface brightness along the bar major axis and that along the minor axis is maximal. In what follows, we use an N-body simulation to test the validity of this proposed method in identifying the bar corotation, assuming that the GALAKOS model is a realistic simulation of a barred spiral having dark gaps.
3.1. Resonances in Simulated Barred Galaxies
In the simulated galaxy, we measure the bar pattern speed, Ωp , and locations of resonances using a spectrogram as a function of radius and frequencies for the m = 2 Fourier harmonic in the stellar disk (see D'Onghia & Aguerri 2020, their Figure 11). The spectrogram is constructed using snapshots sampled at 5 Myr intervals, with the bar pattern speed corresponding to the greatest power frequency. Corotation corresponds to the radius at which the bar pattern speed is equal to the circular angular frequency Ω = Vc /R. Similarly, the Lindblad and ultraharmonic resonances occur at the radii at which the bar pattern speed is equal to Ω ± κ/2 and Ω ± κ/4, respectively, where κ is the epicyclic frequency. Figure 2 displays the spectrogram for m = 2 applied to GALAKOS that shows the presence of a long bar with a pattern speed of 39.4 ± 0.4 km s−1 kpc−1 as measured in the time interval between 2.1 < t/Gyr < 2.6.
Figure 2. m = 2 Fourier harmonic spectrogram from the GALAKOS simulation as a function of radius and frequencies between 2.1 < t/Gyr < 2.6, showing the bar pattern speed Ωp = 39.4 ± 0.4 km s−1 kpc−1. The solid line shows Ω = Vc /R, used to identify corotation, and the dashed and dotted lines show similar tracks to identify the Lindblad and ultraharmonic resonances. In this interval, RCR = 6.40 ± 0.08 kpc and the bar patter speed, Ωp , is marked in red along the y-axis label.
Download figure:
Standard image High-resolution imageThe radius of maximal power is identified using a quadratic fitting approach implemented using the Python package bettermoments 17 (Teague & Foreman-Mackey 2018) after smoothing the data using a Savitsky–Golay filter (Savitzky & Golay 1964) with a width of 3 radial bins. This allows for accurate identification, with uncertainties, of the bar pattern speed (Ωp ), radius of corotation (RCR), radii of inner and outer ultraharmonic resonance (RiUHR, RoUHR), and radii of inner and outer Lindblad resonance (RILR, ROLR). We perform this analysis on 23 time intervals of 0.25 Gyr between 1.1 < t/Gyr < 4.1.
3.2. Simulation Stellar Surface Density
To search for the presence of dark gaps toward the L4,5 Lagrange points and verify their correspondence with corotation, we compute the stellar surface density in the simulation, averaged across the same time intervals used to measure the pattern speeds and resonances above. We degrade the simulation resolution to simulate the effects of dithered MaNGA IFU observations by first taking three mock observations within dithered "beams" that sample a specified physical scale. Then, the resulting mock observations are interpolated to finer "spaxels" that are 1/4 the scale of the "beams," similar to the 2'' fibers and 05 spaxels of MaNGA. The stellar surface density snapshots are converted to bar-centered polar coordinates, with the bar aligned along 0°–180°. The bar position angle is determined using a bounding box fit around spaxels with a stellar surface density greater than 10−0.5 pc−2. Next, the mock observations are binned across a radial range of 15 kpc, with 76 overlapping radial bins of width 0.5 kpc. Similarly, the data are binned azimuthally with a width of 2
5, assuming rotational symmetry. We normalize these maps for each snapshot by first considering the mean value across all azimuth in a radial bin and finding the deviation from this mean value for each azimuthal bin relative to that bin's standard deviation. The snapshot maps are then averaged together using 10,000 bootstrap samples, with a resulting azimuthal variation map displayed in Figure 3 for the same time interval as in Figure 2.
Figure 3. Polar plot with the same orientation as in Figure 1 showing increases (red) and decreases (blue) of the stellar surface density between 2.1 < t/Gyr < 2.6 in the GALAKOS simulation. The differences from the mean values are scaled by their standard deviation, σ. A decrease in the stellar surface density is seen in the simulation along the bar minor axis, similar to the dark gaps used in Buta (2017).
Download figure:
Standard image High-resolution imageFigure 4 shows the relative increase and decrease of the stellar surface density within 30° of the bar major and minor axes, respectively, again for the same time interval. The difference between the major and minor axes' radial stellar mass profiles and the peak difference location is shown in this figure's right panel. Buta (2017) attributed this effect to the L4,5 and L1,2 Lagrange points, and thus approximately at the radius of corotation. However, this radius does not correspond to corotation in the simulation but instead is the inner ultraharmonic bar resonance location. In this example, corotation is located at a larger galactocentric radius, near the innermost radius where the differences between the major and minor axes are zero, which we call RCross,CR. Note that this metric is not always accurate and tends to underestimate the correct corotation radius slightly.
Figure 4. Left: stellar surface density deviation as a function of radius between 2.1 < t/Gyr < 2.6 and within 30° of the bar major (blue) and minor (pink) axis from the GALAKOS simulation. Right: the difference between the major- and minor-axis surface density deviations as a function of radius. The radii with a maximal difference and a null difference are marked with an X and O, respectively. In both panels, the solid, dotted, and dashed black vertical lines mark the locations of corotation, the Lindblad resonances, and the ultraharmonic resonances, respectively.
Download figure:
Standard image High-resolution imageThese synthetic observations of a simulated barred disk galaxy provide an interesting alternative to the dark-gap/corotation interpretation described by Buta (2017). The residual gap in the GALAKOS model is found to coincide with the model's inner ultraharmonic, not its corotation. Athanassoula et al. (1982) showed that the ratio of the radii of corotation and the ultraharmonic resonance can be written as
![Equation (1)](https://content.cld.iop.org/journals/0004-637X/929/2/112/revision1/apjac5d55eqn1.gif)
where and δ describes the galaxy rotation curve where V ∼ r−δ+1. δ = 1 corresponds to a flat rotation curve, and Athanassoula et al. (1982) notes that 0.7 ≤ δ ≤ 1.0 is representative of most barred galaxies. We therefore can estimate the radius of corotation based on the ultraharmonic resonance, which we call RRatio,CR, as
![Equation (2)](https://content.cld.iop.org/journals/0004-637X/929/2/112/revision1/apjac5d55eqn2.gif)
Alternatively, corotation may also correspond with the innermost radius where the major-to-minor-axis differences are zero (RCross,CR; labeled as "First Zero Intersection" in Figure 4), though we later show that this method is not always accurate.
We then test the ability of our dark-gap-based method to identify the UHR throughout the dynamical evolution of our simulated galaxy. Figure 5 shows the estimated and true radii of the UHR across 23 time intervals of 0.25 Gyr that overlap in the range of 1.1 < t/Gyr < 4.1 and for different physical resolutions of the stellar density images. Over this time interval, the bar evolves and decreases its pattern speed through angular momentum transport. Our method tends to work best when considered at the approximate physical resolution of MaNGA spaxels, where, at higher resolution, finer density variations add significant noise and local extrema. Figure 6 shows the same methodology but when applied for corotation as estimated using Equation (2). The uncertainties associated with RRatio,CR are large here because we assume a large range in the δ parameter that describes the rotation curve shape.
Figure 5. The estimated vs. true UHR radius across 23 time intervals of 09.25 Gyr in the range of 1.1 < t/Gyr < 4.1 in the GALAKOS simulation, with points color-coded by the bar pattern speed at each snapshot. Each panel corresponds to a different physical spaxel resolution (labeled in the upper left corner of each panel) of the stellar density images used in the measurements, with red dashed outlines highlighting resolution scales corresponding to MaNGA galaxies. The red dotted line shows a one-to-one relation.
Download figure:
Standard image High-resolution imageFigure 6. Our estimated vs. true corotation radius assuming a range of rotation curve shapes across 23 time intervals of 09.25 Gyr in the range of 1.1 < t/Gyr < 4.1 in the GALAKOS simulation, with points color-coded by the bar pattern speed at each snapshot. Each panel corresponds to a different physical resolution (labeled in the upper left corner of each panel) of the stellar density images used in the measurements, with red dashed outlines highlighting resolution scales corresponding to MaNGA galaxies. The red dotted line shows a one-to-one relation.
Download figure:
Standard image High-resolution imageFigure 7 shows the same again, but for an estimate of corotation using RCross,CR, the innermost radius where the major-to-minor-axis differences are zero. This method does not adequately predict the corotation radius's valid location across different evolutionary phases of the galaxy. Instead, the location of this cross point is strongly affected by the bar pattern speed.
Figure 7. The estimated vs. true corotation radius using the innermost radius where the major-to-minor-axis differences are zero across 23 time intervals of 09.25 Gyr in the range of 1.1 < t/Gyr < 4.1 in the GALAKOS simulation, with points color-coded by the bar pattern speed at each snapshot. Each panel corresponds to a different physical resolution (labeled in the upper left corner of each panel) of the stellar density images used in the measurements, with red dashed outlines highlighting resolution scales corresponding to MaNGA galaxies. The red dotted line shows a one-to-one relation.
Download figure:
Standard image High-resolution imageFigure 8 shows all three of these metrics in terms of a fractional error (True/Estimated) across different time intervals and physical resolution scales. Within the range of resolutions corresponding with MaNGA galaxies, we can describe the accuracy and precision for identifying the UHR as . Similarly, by assuming a rotation curve shape, we find
. Although RCross,CR is not an ideal estimator, it returns an averaged accuracy of
. We adopt these simulation-derived ratios as a correction factor, with uncertainties, to adjust the measured locations and improve the accuracy of the predicted resonance locations. With this method validated on simulation data and with an understanding of the corrections needed to be made and uncertainties involved, we can apply it to a sample of observed galaxies, both ringed and nonringed, to better understand its validity.
Figure 8. Left: true/estimated radii of the ultraharmonic resonance (RUHR; top panel) and corotation as measured using a ratio (RRatio,CR; middle panel) and using the first zero crossing point (RCross,CR; bottom panel) at different physical resolution scales in 23 intervals of 0.25 Gyr between 1.1 < t/Gyr < 4.1. The pattern speed of the bar colors the points in the corresponding simulation snapshot. The approximate physical scale of the MaNGA fibers spans the resolution scales in bold. A precise and accurate estimate would be narrowly distributed and centered on one (dashed red lines). Right: histograms of the true/estimated radii for the resolutions comparable to MaNGA galaxies. The median ratios and uncertainties estimated using the 16th and 84th percentiles of the distributions are displayed for each metric. These will be used as a correction factor.
Download figure:
Standard image High-resolution image4. Observational Results
The 578 barred galaxies with bar masks and resolved galaxy parameters provide a unique statistical view into the impact of a bar on galaxies' distribution of matter. To consider our sample as a whole and understand the significance of rings in our barred galaxies, we split our sample into two subsamples of 211 ringed and 367 nonringed galaxies using a 50% unbiased vote fraction on the Galaxy Zoo question "t08_odd_feature_a19_ring" from Hart et al. (2016). We then project all ringed and nonringed galaxies onto a bar-oriented cylindrical coordinate system, with the radius in units of the bar radius (Rbar) and azimuth angle oriented such that the bar major axis spans from 0° to 180°. We use visual classifications from Galaxy Zoo (Lintott et al. 2011) to orient all spiral arms in a counterclockwise spiral pattern. In this bar reference frame, we assume a 180° rotational symmetry so that spaxels at an azimuth of 30° and 210° are equivalent. To ensure an equal sampling of different radii in galaxies, we restrict each galaxy to only extend out to the maximum bar-oriented radius found in all azimuthal directions. Additionally, we remove spaxels with a g-band S/N < 2 to only consider spaxels that lie on the galaxy disk. We can then search for azimuthal and radial surface density variations in the context of bar-driven resonances.
4.1. Surface Density Variations
Maps of the stellar mass surface density from Pace et al. (2019a, 2019b) provide a resolved measure of the distribution of stars in MaNGA galaxies. The maps are binned in a similar radial and azimuthal schema as the N-body simulations: radial bins are computed relative to the bar radius, overlap with a width of 0.2R/Rbar, and have a step size of 0.04R/Rbar, out to a maximum radius R = 3Rbar. We normalize these maps for each galaxy by first considering the mean value across all azimuth in a radial bin and finding the deviation from this mean value for each azimuthal bin in terms of the standard deviation computed across all azimuth. We also do the same using SDSS g-band imaging to test the ability to perform this analysis using only imaging data. This process produces maps for all 578 galaxies, showing azimuthal variations across different radial bins, scaled by each radial bin's standard deviation. We combine the resulting maps using 10,000 bootstrap resamples to create Figure 9. The stellar mass surface density and g-band magnitude show similar behavior in the simulations and Figure 3.
Figure 9. Same as in Figure 3, but for the 578 barred MaNGA galaxies using the PCA stellar mass surface density of Pace et al. (2019a, 2019b) (top half; 0°–180°) and SDSS G-band imaging (bottom half; 180°–360°), separated into 211 ringed (left) and 367 nonringed (right) galaxies. The radial axis is normalized by the bar radius.
Download figure:
Standard image High-resolution imageFigure 10 follows the same method used for Figure 4, but with the observed stellar mass surface density and g-band magnitudes of MaNGA galaxies. The peak difference is identified in the same manner, and the correction factors are applied to adjust for the expected offsets based on the simulation (see Figure 8). We find the average location of the inner ultraharmonic resonance to be , with estimated corotation radii with both methods of
and
. Bars with RCR/Rbar ≥ 1.4 are considered "slow," while bars with RCR/Rbar < 1.4 are considered "fast" (Debattista & Sellwood 2000). Bars with RCR/Rbar < 1 are often referred to as "ultrafast," but it is unclear whether these bars are a real feature or a result of imperfect estimates of either RCR or Rbar (e.g., Buta 2017). Our results suggest that slow bars are the most likely scenario for this sample of galaxies as an ensemble, since RCR/Rbar > 1.4. The sample of galaxies considered in Buta (2017) are also primarily slow, with RCR/Rbar ≥ 1.58. However, the relatively large uncertainties in this population-level diagnostic do not categorically exclude the fast rotator alternative.
Figure 10. Same as in Figure 4, but for the 578 barred MaNGA galaxies, separated into 211 ringed (top panel) and 367 nonringed (bottom panel) galaxies using the PCA stellar mass surface density of Pace et al. (2019a, 2019b) (solid lines) and SDSS g-band imaging (dotted lines). The inner ultraharmonic resonance (RUHR) is approximated at the location of maximal difference between the major and minor axes, and corotation is approximated as RRatio,CR/RUHR = 1.8 ± 0.3 (solid shading) and as the first zero intersection (RCross,CR; hatched shading). All estimates have also been adjusted to include the correction factors found in the simulations (see Figure 8).
Download figure:
Standard image High-resolution imageWe also perform the same analysis on the 578 individual galaxies considered here, with 185 galaxies returning reasonable estimates of RUHR and RCR using the PCA stellar mass surface density as shown in Figure 11. Four example galaxy images, with their inferred ultraharmonic resonance location marked with shaded rings, are shown in Figure 12. While many galaxies seem to be classified as slow, the large uncertainties with the corotation radius estimates leave no galaxies with a 2σ estimate of RCR/Rbar that is slow. Combined with the fact that the peak of the RCR/Rbar distribution for these 185 galaxies is within the fast threshold, it seems likely that most bars in MaNGA galaxies in our sample are, in fact, fast rotators. This is in general agreement with an independent study using a smaller sample of MaNGA galaxies with the Tremaine & Weinberg (1984) method, which found (Garma-Oehmichen et al. 2020).
Figure 11. Estimated radius of the inner ultraharmonic resonance (RUHR) and radius of corotation (RCR) for individual ringed (left) and nonringed (right) MaNGA galaxies. Data are limited to galaxies where measurement errors are under 0.2Rbar for identifying the peak difference. Error bars are shown for one-fifth of the data points. The dotted horizontal lines mark the boundaries between ultrafast, fast, and slow bars. Marginalized Gaussian kernel density estimates are shown along both axes, with the dashed lines including the associated errors. A total of 74 ringed and 111 nonringed galaxies are plotted, with 5 ultrafast, 12 fast, and 57 slow bars in the ringed sample and 22 ultrafast, 30 fast, and 59 slow bars in the nonringed sample as diagnosed using RCross,CR. When considering the errors, the number of galaxies with ultrafast and slow bars drops significantly.
Download figure:
Standard image High-resolution imageFigure 12. Example SDSS gri images with the inferred location of the inner ultraharmonic resonance and corotation shown in the right-hand image pair in blue and red, respectively, for ringed (left) and nonringed (right) galaxies. The image axes are scaled based on the measured bar radius. The MANGA-ID of each galaxy is displayed above the left-hand image pair. The purple hexagons mark the footprint of the MaNGA survey. These images are exactly what are used in the citizen science classifications of Masters et al. (2021).
Download figure:
Standard image High-resolution image5. Discussion
Rings are often suggested to be linked to resonances with the pattern speed of bars. However, measurements of the pattern speed and location of resonances are difficult to observe in large samples of galaxies, especially for face-on galaxies. Zhang & Buta (2007) described a phase shift method to locate corotation radii in face-on barred and spiral galaxies, but this method assumes that spirals and bars are quasi-steady modes in galactic disks. Later, Buta (2017) described a novel method for detecting resonances, attributing "dark gaps" in galaxies to the location of corotation. However, if we accept that GALAKOS is a reliable model of a barred galaxy, then it provides an interesting alternative interpretation of the gaps. These same "dark gaps" may indicate the location of the 4:1 ultraharmonic resonance with the bar, in both ringed and nonringed galaxies. This finding now opens up a new avenue to explore the effects of bar resonances on galaxies' structure using large samples of galaxy images.
Kim et al. (2016) analyzed the relation of the dark gaps with the bar properties for a sample of barred galaxies. They found that larger and stronger bars produced stronger dark gaps, indicating that the growth process of the bar is linked with the formation or the dark gaps. Similar to our cases, their dark gaps are located at radius smaller than the bar radius (see their Figure 5). In addition, they run a galaxy numerical model to understand the formation of the dark gaps. In this simulation the dark gaps formed were also located at a radius smaller than the bar radius. This is similar to what we obtain in our GALAKOS simulation. All these facts imply that if the dark gaps would be located at the corotation radius (Buta 2017), all bars would be ultrafast rotators. This would be in contradiction with the classical view of bars supported by the orbits of the x1-family (e.g., Contopoulos & Papayannopoulos 1980; Athanassoula 1992).
The most accurate method for measuring the bar pattern speed, which is used to derive the corotation radius, is using the kinematic and model-independent method proposed by Tremaine & Weinberg (1984). The method is, however, affected by multiple error sources, notably the determination of geometric parameters and modeling of the rotation curve (Garma-Oehmichen et al. 2020). Additionally, it is often best suited to work with early-type galaxies. Guo et al. (2019) used this method to measure the bar pattern speed and corotation radii of a sample of 53 MaNGA galaxies. Of these, 10 galaxies are also in our barred sample and have reliable estimates from our new method without considering any kinematic information. Figure 13 shows the difference between our estimated corotation radius and the measurements from Guo et al. (2019) using their kinematic and photometric methods. While large discrepancies are present in some individual galaxies, the uncertainties combined with the kinematic method cannot rule out an overall agreement. The median value of the ratio between our estimate using RRatio,CR and the Guo et al. (2019) mass-weighted kinematic method is , showing that our method tends to estimate the corotation radius to be farther out than using the Tremaine & Weinberg (1984) method. Figure 13 also shows differences between our bar length estimates from Galaxy Zoo:3D and the estimates derived in Guo et al. (2019). In the future, applying the Tremaine & Weinberg (1984) method to the full sample of MaNGA galaxies will allow for more statistically robust validation of our method. The kinematic approach to identifying the corotation radius is likely still the most accurate method, but our new method explored in this work may be preferable for more face-on galaxies, whose kinematics are less reliable in projection and whose morphological features can be better characterized.
Figure 13. The ratio of the estimated corotation radii and the estimates from the mass-weighted (cyan circle) and light-weighted (pink cross) measurements from Guo et al. (2019) for galaxies overlapping in both samples. The top and middle panels are for the kinematic- and photometric-based estimates from Guo et al. (2019), respectively, and shading encompasses both measurements' errors. Bottom panel: the ratio between our bar length estimate (Galaxy Zoo:3D) and the estimates in Guo et al. (2019; units of arcsec). The "Galaxy Index" is an arbitrary index describing each overlapping galaxy in our sample and that of Guo et al. (2019).
Download figure:
Standard image High-resolution imageThough the method used to measure bar length can impact estimates of resonance location, it is difficult to gauge which method is best. We note that the Galaxy Zoo:3D bar length estimates are generally consistent with estimates using other methods, including in our overlap with the work of Guo et al. (2019) with a median of Rbar,GZ:3D/Rbar,Guo = 0.9 ± 0.1. Guo et al. (2019) used a combination of three methods to estimate bar lengths, including the ellipticity radial profile, the ellipticity position angle profile, and the Fourier decomposition. In Krishnarao et al. (2020, their Appendix), the Galaxy Zoo:3D bar lengths were tested more extensively with a Fourier-based method of Kraljic et al. (2012), showing consistent results.
In addition, we can explore the correlation between different galaxy morphological types and their bar speeds. Using the most likely morphological types as diagnosed in Galaxy Zoo 2 (Hart et al. 2016), we compare differences in the corotation radii. Figure 14 shows the distribution of measured resonance radii separated by their morphological types. Several galaxies have a likely morphology from Galaxy Zoo 2 as SA, or not containing a bar. These galaxies likely have weaker bars than those with a likely morphology of SB, which would have bars more clearly visible in SDSS imaging. Here we use this difference as an approximate proxy for bar strength. With our current sample, we cannot robustly diagnose any correlation with morphology or bar strength and bar speeds. Understanding the impact of the bar strength on the bar speed could be studied in the future with more detailed analysis, but it is beyond the scope of this paper.
Figure 14. Same as in Figure 11, but separated by morphological types from Galaxy Zoo 2 (Hart et al. 2016), with SB vs. SA (left panel; proxy for strong vs. weak bar) and spirals vs. lenticulars (right panel).
Download figure:
Standard image High-resolution imageIn our analysis of individual galaxies, only 32% of our full sample of 578 galaxies showed a clear "dark-gap" signature in the PCA stellar mass surface density maps. While this is relatively low, it may still be possible to improve the success rate using deeper imaging data and to vary the amount of smoothing. In particular, it may be possible to identify "dark gaps" using the deep multiband imaging from the DESI legacy imaging survey (Dey et al. 2019). This approach also does not involve the more complicated effects of spatial covariances resulting from dithered IFU observations.
In this paper, we only consider the stellar mass and light concerning bar-driven resonances because the GALAKOS simulation lacks gas. It has been shown that gas can significantly impact the stellar dynamics of bars, which our test simulations do not account for (e.g., Athanassoula 2003; Combes 2008). However, cosmological simulations with the same level of detail as GALAKOS that include gas are not yet available but in the future will provide an ideal way to calibrate this gap method. In addition, including simulations of ringed galaxies and other types will help determine the universality of the gap method on all barred galaxies. Initial examination of the MaNGA emission-line data suggests that the gas distribution may also closely track the azimuthal variations of the stellar density shown in Figure 9 based on a similar map using the Hα equivalent width. However, the precise location of the "dark gaps" is slightly shifted from the stellar matter. Fraser-McKelvie et al. (2020) investigated the correlation between the Hα morphology of bars and their star formation rates, finding correspondence between different gas morphologies and the amount of quenching in bars. To fully understand the nature of the gas with resonances, we must, in the future, similarly calibrate our metrics using hydrodynamic simulations.
6. Conclusions
This paper has proposed a new method to measure the radius of both corotation and the ultraharmonic resonance in galaxies without measurements of their internal kinematics. Our method is very similar to the "dark-gap" method of Buta (2017), but with a significant change in the interpretation, prompted by tests on the N-body simulation, GALAKOS (D'Onghia & Aguerri 2020). We summarize our main findings below:
- 1.In the GALAKOS simulation, dark gaps are found in the stellar mass surface density but appear not to be related to the L4,5 Lagrangian points as was proposed by Buta (2017). Instead, the model favors the gaps to be linked to the inner 4:1 ultraharmonic resonance.
- 2.If the rotation of the galaxy hosting the "dark gap" is known, then the corotation resonance can be precisely identified using Equation (1).
- 3.Alternatively, assuming a generally flat rotation curve, the radius of corotation can be predicted as RRatio,CR/RUHR = 1.8 ± 0.3 but with decreased precision.
- 4.This method works best for galaxies with imaging, or IFU spaxels, at ∼0.2–0.4 kpc scales, similar to MaNGA galaxies, so that the strong nonaxisymmetric variations from the bar in the innermost regions do not dominate.
- 5.Applying this method to a sample of 578 barred MaNGA galaxies reveals that, on average, the MaNGA sample represents a population of bars that are slow rotators, but it cannot be ruled out that the bars are fast rotators given the uncertainties, especially on the determination of the bar radius.
- 6.
- 7.
- 8.None of the individual bars can be confidently determined to be slow rotators when considering their 2σ uncertainties.
In particular, our method does not require kinematic information to determine the fundamental parameters describing the dynamics of bars. It can be used on large samples of barred galaxies from imaging surveys, such as the Legacy Survey (Dey et al. 2019). This method allows for the nature of rings and other structures in barred galaxies to be reexamined in terms of bar-driven resonances, paving the way for better diagnostics of galaxies' internal dynamical structures. With the continued use of citizen science projects like Galaxy Zoo to build larger training samples, it may soon be possible to identify bar resonances with automated machine-vision techniques or neural networks.
D.K. and Z.J.P. acknowledge support from the NSF CAREER award AST-1554877. D.K. is supported by an NSF Astronomy and Astrophysics Postdoctoral Fellowship under award AST-2102490. J.A.L.A. is supported by the Spanish MINECO grant AYA2017-83204-P.
J.G.F.-T. gratefully acknowledges the grant support provided by Proyecto Fondecyt Iniciación No. 11220340, and also from ANID Concurso de Fomento a la Vinculación Internacional para Instituciones de Investigación Regionales (Modalidad corta duración) Proyecto No. FOVI210020, and from the Joint Committee ESO-Government of Chile 2021 (ORP 023/2021).
Funding for the Sloan Digital Sky Survey IV has been provided by the Alfred P. Sloan Foundation, the U.S. Department of Energy Office of Science, and the Participating Institutions. SDSS acknowledges support and resources from the Center for High-Performance Computing at the University of Utah. The SDSS website is www.sdss.org.
SDSS is managed by the Astrophysical Research Consortium for the Participating Institutions of the SDSS Collaboration including the Brazilian Participation Group, the Carnegie Institution for Science, Carnegie Mellon University, the Chilean Participation Group, the French Participation Group, Harvard-Smithsonian Center for Astrophysics, Instituto de Astrofísica de Canarias, The Johns Hopkins University, Kavli Institute for the Physics and Mathematics of the Universe (IPMU)/University of Tokyo, the Korean Participation Group, Lawrence Berkeley National Laboratory, Leibniz Institut für Astrophysik Potsdam (AIP), Max-Planck-Institut für Astronomie (MPIA Heidelberg), Max-Planck-Institut für Astrophysik (MPA Garching), Max-Planck-Institut für Extraterrestrische Physik (MPE), National Astronomical Observatories of China, New Mexico State University, New York University, University of Notre Dame, Observatório Nacional/MCTI, The Ohio State University, Pennsylvania State University, Shanghai Astronomical Observatory, United Kingdom Participation Group, Universidad Nacional Autónoma de México, University of Arizona, University of Colorado Boulder, University of Oxford, University of Portsmouth, University of Utah, University of Virginia, University of Washington, University of Wisconsin, Vanderbilt University, and Yale University.
Facility: Sloan. -
Software: astropy (Robitaille et al. 2013; Astropy Collaboration et al. 2018), matplotlib (Hunter 2007), seaborn (Waskom et al. 2020), scipy (Virtanen et al. 2020), sdss-marvin (Cherinka et al. 2019), bettermoments (Teague & Foreman-Mackey 2018).
Footnotes
- 15
- 16
- 17