Galaxy Zoo: Kinematics of strongly and weakly barred galaxies

We study the bar pattern speeds and corotation radii of 225 barred galaxies, using IFU data from MaNGA and the Tremaine-Weinberg method. Our sample, which is divided between strongly and weakly barred galaxies identiﬁed via Galaxy Zoo, is the largest that this method has been applied to. We ﬁnd lower pattern speeds for strongly barred galaxies than for weakly barred galaxies. As simulations show that the pattern speed decreases as the bar exchanges angular momentum with its host, these results suggest that strong bars are more evolved than weak bars. Interestingly, the corotation radius is not diﬀerent between weakly and strongly barred galaxies, despite being proportional to bar length. We also ﬁnd that the corotation radius is signiﬁcantly diﬀerent between quenching and star forming galaxies. Additionally, we ﬁnd that strongly barred galaxies have signiﬁcantly lower values for R , the ratio between the corotation radius and the bar radius, than weakly barred galaxies, despite a big overlap in both distributions. This ratio classiﬁes bars into ultrafast bars ( R < 1.0; 11% of our sample), fast bars (1.0 < R < 1.4; 27%) and slow bars ( R > 1.4; 62%). Simulations show that R is correlated with the bar formation mechanism, so our results suggest that strong bars are more likely to be formed by diﬀerent mechanisms than weak bars. Finally, we ﬁnd a lower fraction of ultrafast bars than most other studies, which decreases the recently claimed tension with Λ CDM. However, the median value of R is still lower than what is predicted by simulations.


INTRODUCTION
Bars are a relatively common structure in galaxies, with about 30%-60% of nearby galaxies hosting a bar, depending on the redshift and wavelength range of the study (Marinova & Jogee 2007;Menéndez-Delmestre et al. 2007; Barazza et al. 2008;Sheth et al. 2008;Nair & Abraham 2010b;Masters et al. 2011).Bars can drive angular momentum outwards and funnel gas to the centre of the galaxy (Athanassoula 1992b;Davoust & Contini 2004; Rodriguez-★ E-mail: tobias.geron@physics.ox.ac.uk (TG) Fernandez & Combes 2008;Athanassoula et al. 2013;Villa-Vargas et al. 2010;Fragkoudi et al. 2016;Vera et al. 2016;Spinoso et al. 2017;George et al. 2019;Seo et al. 2019).This will result in significant "secular" evolution of the host, caused directly by its bar (Kormendy & Kennicutt 2004;Cheung et al. 2013;Sellwood 2014;Díaz-García et al. 2016b;Kruk et al. 2018).Moreover, multiple studies have found that bars appear more often in massive, red and gas-poor galaxies (Masters et al. 2012;Cervantes Sodi 2017;Vera et al. 2016;Fraser-McKelvie et al. 2020).For example, Kruk et al. (2018) found that bars are redder than disks and that the disks of barred galaxies are redder than the disks in unbarred galaxies.These results suggest that bars might be linked to the quenching of their host.This can be either the result of triggering a starburst in the centre of the galaxy, after extensive inflows of gas (Alonso-Herrero & Knapen 2001;Sheth et al. 2005;Jogee et al. 2005;Hunt et al. 2008;Carles et al. 2016), or by making the gas too dynamically hot for star formation (Zurita et al. 2004;Haywood et al. 2016;Khoperskov et al. 2018;Athanassoula 1992b;Reynaud & Downes 1998;Sheth et al. 2000).In any case, a bar is a very common and important structure in a galaxy, so understanding bars is fundamental to understanding galaxy evolution.
Bars are historically classified into weak or strong.Since de Vaucouleurs (1959,1963), three subclasses are recognised: unbarred (SA), strongly barred (SB) and weakly barred (SAB).Weakly barred galaxies were thought to be an intermediate class between unbarred and strongly barred, having lengths and contrast in between SA and SB bars.Bars classified as weak were usually small and faint, whereas bars classified as strong were long and obvious de Vaucouleurs (1959Vaucouleurs ( , 1963)).Morphological arguments are still used to determine bar type.Nair & Abraham (2010a) produced a catalogue of detailed visual morphological classifications and they distinguished between weak and strong by looking at whether the bar dominated the light distribution.The bar strength can also be estimated using the maximum ellipticity and boxiness of the isophotes (Athanassoula 1992a;Laurikainen & Salo 2002;Erwin 2004;Gadotti 2011).One can also look at the surface brightness profiles of bars.Previous work has shown that stronger bars have flat profiles, while weaker bars have exponential profiles (Elmegreen & Elmegreen 1985;Elmegreen et al. 1996;Kim et al. 2015;Kruk et al. 2018).Clearly, there are many ways to characterise bars and their strength.However, the community has yet to reach a consensus on how to best define weak and strong bars and on which detection method is superior.
This problem was addressed more recently by the Galaxy Zoo (GZ) team, who combined the efforts of citizen scientists and machine learning to provide morphological classifications of galaxies (Lintott et al. 2008;Walmsley et al. 2022).These morphological classifications included a distinction between weak and strong bars based on visual morphology.Volunteers are shown examples of weak and strong bars prior to classification.The strong bars are typically large and obvious structures, whereas the weak bars can be smaller and fainter.Géron et al. (2021) used the morphological classifications from GZ based on images from the Dark Energy Camera Legacy Survey (DECaLS, Dey et al. 2019) to study weak and strong bars, and found that around 28% of all disk galaxies have a weak bar, while 16% had a strong bar.They also found that, when correcting for bar length, any difference they observed between weak and strong bars disappeared.Thus, they suggested that weak and strong bars are not fundamentally different physical phenomena.Instead, they proposed the existence of a continuum of bar types, which varies from 'weakest' to 'strongest'.Most research on bars has traditionally been focussed on stronger bars, as they are more obvious and clearer structures.However, weaker bars are still very common structures in galaxies and need to be included in more studies to obtain a more complete picture.
The bar pattern speed (Ω bar ), or the rotational frequency of the bar, is one of the most important dynamical parameters that describe a bar.It is intrinsically linked to the evolution of the bar and its host.It is typically found that, as the bar exchanges angular momentum with its host, the bar grows and the pattern speed decreases (Debattista & Sellwood 2000;Athanassoula 2003;Martinez-Valpuesta et al. 2006;Okamoto et al. 2015).
If the bar pattern speed and galaxy kinematics are known, one can calculate the corotation radius (R CR ), which is the radius at which the angular speed of the stars in the disc is equal to the pattern speed of the bar.Additionally, one can also calculate the dimensionless corotation radius-to-bar radius ratio, R. A large value for R implies that the point of corotation is far outside the bar region.This ratio is typically used to separate bars into 'fast' (1.0 < R < 1.4) and 'slow' (R > 1.4) bars (Debattista & Sellwood 2000;Rautiainen et al. 2008;Aguerri et al. 2015).There is a known correlation between the formation of the bar and R. Bars that are triggered by tidal interactions tend to be in the slow regime for a longer time and have higher values for R than bars formed by global bar instabilities (Sellwood 1981;Miwa & Noguchi 1998;Martinez-Valpuesta et al. 2016, 2017).
There is also a known tension between simulations and observations on the distribution of the ratio R. Cosmological simulations predict that bars slow down significantly due to dynamical friction with their dark matter halo, which results in high values for R.However, observations typically find lower values of R, which has been highlighted as a challenge for the ΛCDM cosmology used in these cosmological simulations (Algorry et al. 2017;Peschken & Łokas 2019;Roshan et al. 2021b).
It is suggested that bars cannot extend beyond their corotation radius (Contopoulos 1980(Contopoulos , 1981;;Athanassoula 1992b).This implies that bars with R < 1 should not exist.However, these so-called 'ultrafast' bars have been repeatedly observed (Buta & Zhang 2009;Aguerri et al. 2015;Cuomo et al. 2019;Guo et al. 2019;Garma-Oehmichen et al. 2020;Krishnarao et al. 2022).This discrepancy between observation and theory remains an open question, although some suggest that the cause for this problem is rooted in incorrect estimates of the bar radius (Cuomo et al. 2021;Roshan et al. 2021a).
It is becoming clear that the pattern speed and the parameters derived from it (such as corotation radius and R) are important.However, it is also quite challenging to correctly estimate the bar pattern speed.Nevertheless, various methods exist to measure this dynamical parameter.For example, one can match the observed surface gas distribution or gas velocity field with simulations where Ω bar is a free parameter (Sanders & Tubbs 1980;Hunter et al. 1988;Lindblad & Kristen 1996;Weiner et al. 2001;Rautiainen et al. 2008;Treuthardt et al. 2008).Alternatively, one can subtract a rotation model from the gas velocity field and look at the morphology of the residuals to estimate the pattern speed (Sempere et al. 1995;Font et al. 2011Font et al. , 2017)).Other morphological features are helpful to determine the bar pattern speed, such as rings (Buta 1986;Rautiainen & Salo 2000;Muñoz-Tuñón et al. 2004;Pérez et al. 2012), the shape and offset of dust lanes (Athanassoula 1992b;Sánchez-Menguiano et al. 2015) and the morphology of spiral arms (Puerari & Dottori 1997;Aguerri et al. 1998;Buta & Zhang 2009;Sierra et al. 2015).
However, all these methods require some sort of modelling.The only reliable direct and model-independent method to determine the bar pattern speed is the Tremaine-Weinberg (TW) method (Tremaine & Weinberg 1984).It has been used extensively in the past to study bar pattern speeds (Aguerri et al. 2015;Cuomo et al. 2019;Guo et al. 2019;Garma-Oehmichen et al. 2020).The TW method uses surface brightness and line-of-sight (LOS) velocity data to estimate the pattern speed.
In this paper, we use the TW method on integral-field spectroscopy data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (Bundy et al. 2015) to estimate bar pattern speeds, corotation radii and the dimensionless ratio R for a sample of 225 galaxies.This is the largest sample to date measured with the TW method and includes both weakly and strongly barred galaxies, identified using GZ.
The structure of the paper is as follows: in Section 2, we explain the Tremaine-Weinberg method in detail.The data and sample selection is explained in Section 3. Section 4 shows our results, which are discussed in Section 5. Finally, our conclusions are summarised in Section 6.Where necessary, we assumed a standard flat cosmological model with H 0 = 70 km s −1 Mpc −1 , Ω m = 0.3 and Ω Λ = 0.7.

Theory
The Tremaine-Weinberg (TW) method is a model-independent method to determine the pattern speed of a galaxy (Tremaine & Weinberg 1984).The main assumptions of the TW method are that there is a well-defined pattern speed and that the tracer used (i.e.stars or gas) satisfies the continuity equation.To illustrate the different steps of the TW method, we will apply the TW method to one strongly barred galaxy and one weakly barred galaxy, shown in Figure 1.
Take a Cartesian coordinate system (,  ) in the sky plane with the origin in the centre of the galaxy and the -axis aligned with the line of nodes (LON), which is defined as the intersection of the sky plane and the disc plane, so it is effectively the major axis of the galaxy.Then, the Tremaine-Weinberg method can be formulated as: where Ω b is the bar pattern speed,  is the inclination of the galaxy,  LOS is the line of sight velocity, Σ is the surface brightness of the galaxy and ℎ( ) is a weight function.A delta function like ℎ( ) =  ( −  0 ) is typically used here, so that the integration happens in pseudo-slits across the IFU parallel to the LON.Multiple integrations are usually done with different offset distances  0 to ensure reliable measurement of the pattern speed (Tremaine & Weinberg 1984).In this case, Equation 1 can be simplified to: where  is called the photometric integral and  the kinematic integral.They are defined as: is effectively the luminosity-weighted mean position and  is the luminosity-weighted mean line of sight velocity.These photometric and kinematic integrals are calculated for the multiple different pseudo-slits across the IFU.These pseudo-slits are visualised on top of the MaNGA stellar flux and velocity maps in Figure 2. Every pseudo-slit is carefully placed next to each other so that they do not overlap.Every slit has a width of 0.5 arcsec, which is the same width that was used in Guo et al. (2019).Using different slit widths does not have a significant impact on the final measurement (Guo et al. 2019;Zou et al. 2019).To make optimal use of the data, we make the pseudo-slits as long as the data allows.However, each slit should be centred on the disc minor axis.This implies that a target can have multiple slits with sightly different slit lengths.However, this variation is minimal and this approach is similar to what is done in Garma-Oehmichen et al. (2020, 2022).We place as many slits as we can fit within the bar, but impose a minimum of three slits.The maximum amount of slits placed on one galaxy was 48, and the median is 10.
The limits of the integration technically go from −∞ to +∞, as shown in Equation 3.However, this which is not possible with real data.To make sure that the pseudo-slits are long enough, we test the convergence of every pseudo-slit, as suggested by Zou et al. (2019) and Zou et al. (in prep.).This is done by increasing the length of each pseudo-slit by 1 pixel until its maximum length is reached.A slit has converged if the median of the absolute value of the change in Ω b sin () in the last 5 slit lengths tested was less than 1 km s −1 arcsec −1 .Any pseudo-slit that did not meet this threshold was discarded.
We do not calculate the ratio of  and  directly.Instead, we plot  against  for the different pseudo-slits and the slope of the best-fit line going through these points will then be equal to Ω b sin ().This is done to help avoid centering errors and account for incorrect estimates of the systemic velocity (Guo et al. 2019).An example of such a plot can be found in Figure 3.
It is important to note that for an axisymmetric disc, the weighted mean position and velocity integrals will equal zero.This means that any non-zero values will be due to additional structures such as the bar (if it is not aligned or perpendicular to the LON).
The code used in this work to calculate the bar pattern speed is publicly available here 1 .

Concerns and limitations
As alluded to in Section 2.1, the TW method requires that the chosen tracer satisfies the continuity equation.Multiple studies have used gas as the tracer and have been successful in determining the pattern speed using the TW method (Zimmer et al. 2004;Hernandez et al. 2005;Emsellem et al. 2006;Fathi et al. 2009;Gabbasov et al. 2009).
Many studies have also successfully determined pattern speeds by using stars as the tracer (Merrifield & Kuĳken 1995;Debattista et al. 2002;Aguerri et al. 2003;Corsini et al. 2007), although they usually limited their sample to early-type barred galaxies.This is because there were concerns that dust obscuration and star formation  in late-type galaxies could cause the surface brightness to not trace the mass distribution properly.However, other papers show that it is possible to use the TW method on late-type galaxies, despite initial concerns (Gerssen et al. 2003;Gerssen & Debattista 2007;Treuthardt et al. 2007;Aguerri et al. 2015;Cuomo et al. 2019;Guo et al. 2019;Garma-Oehmichen et al. 2020).More recently, Williams et al. (2021) applied the TW method to stellar and gaseous tracers (using both CO and H) and found significantly different results.They attributed this inconsistency to the clumpy nature of the gaseous tracers they used, which resulted in incorrect pattern speed measurements.Thus, in this work we decided to use stars as our tracer.
As Garma-Oehmichen et al. (2020) show, centering issues are not negligible and it is crucial that the LON goes through the centre of the galaxy.In this work, we find the centre by smoothing the stellar flux data with a Gaussian filter and finding the brightest pixel in the smoothened data.
The TW method is also only applicable to galaxies with regular kinematics and on galaxies with intermediate inclinations (20 • <  < 70 • ) (Tremaine & Weinberg 1984;Aguerri et al. 2015;Cuomo et al. 2019;Garma-Oehmichen et al. 2020).Edge-on galaxies do not have enough spatial data, while the stellar velocity is not well constrained in face-on galaxies.Additionally, detecting bars in edgeon galaxies is very difficult and unreliable.It is known that the TW method is very sensitive to incorrect estimates of the PA of the galaxy (Debattista 2003;Zou et al. 2019;Garma-Oehmichen et al. 2020).Thus, a correct and reliable estimate of the position angle is crucial.We try to account for this sensitivity by performing a Monte Carlo (MC) simulation over the uncertainty of the PA (see Section 3.3 for more details).For the TW method to work, it is also important that the bar is not aligned with the major or minor axis of the galaxy, as otherwise the integrals will cancel out.We also need to be able to place a sufficient amount of slits, otherwise the straight line in the  over  plot that is used to determine Ω b is not well constrained.As the slits have to be placed on top of the bar, the TW method is not ideal for the shortest of bars, where only very few slits can be placed.This will disproportionally affect weak bars, which should be kept in mind.The specific thresholds we impose are detailed in Section 3.5.

Calculation of corotation radius and R
After the bar pattern speed is obtained using the TW method, one can calculate the corotation radius ( CR ).This is where the centrifugal and gravitational forces balance each other in the rest frame of the bar, which means that the stars in the disc will have the same angular velocity as the bar pattern speed at the corotation radius (Cuomo et al. 2019;Guo et al. 2019).Various papers calculate this by doing  CR =  c /Ω b , where  c is the circular velocity in the flat part of the rotation curve (Aguerri et al. 2015;Cuomo et al. 2019;Guo et al. 2019).However, this assumes that the corotation radius lies in the region where the rotation curve has flattened.This is not necessarily the case and can lead to incorrect estimates of  CR and R.
Instead, we will use the rotation curve of the galaxy to calculate the corotation radius.The rotation curve can be obtained by using the MaNGA stellar velocity data (see Section 3.4).The bar pattern speed is multiplied by a radius range, which effectively indicates how fast the tracer moves at any radius for that particular pattern speed.The radius at which this curve intersects with the galaxy rotation curve, is the corotation radius.An example of this can be found in Figure 4.
Garma-Oehmichen et al. ( 2020) used a similar approach to ours and compared their results to results obtained by using the  CR =  c /Ω b method.On average, they found a relative difference of ∼15%, indicating that the simplified approach introduces a significant bias.
The corotation radius can be used to calculate the dimensionless parameter R, defined as R =  CR / bar , where  bar is the deprojected bar radius.We obtain estimates for the uncertainty on Ω b ,  CR and R by performing a Monte Carlo simulation using the errors on the input variables and assuming Gaussianity (see Section 3.3 for more details).The posterior distributions of the final pattern speed, corotation radius and R for our example galaxies are shown in Figure 5.

MaNGA survey
We need resolved stellar velocity and stellar flux data in order to implement the TW method, which we obtain from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (Bundy et al. 2015).MaNGA is part of the Sloan Digital Sky Survey IV (SDSS-IV) collaboration (Blanton et al. 2017).More specifically, we used data from the seventeenth data release of SDSS (Abdurro'uf et al. 2022).MaNGA used the Baryon Oscillation Spectroscopic Survey (BOSS) Spectrograph, which has a resolution of R ∼ 2000 and a wavelength coverage of 3600 -10,000 Å (Smee et al. 2013), on the 2.5m Sloan Telescope at Apache Point Observatory (Gunn et al. 2006).Every integral field unit (IFU) consists of 19-127 optical fibers, stacked hexagonally (Drory et al. 2015).Most galaxies are covered out to 1.5 effective radii (R  ), while a third are covered out to 2.5 R  .We make use of the maps that are binned to S/N ∼ 10 using the Voronoi binning algorithm (Westfall et al. 2019).
For more information on the observing strategy, survey design, data reduction process, sample selection and the data analysis pipeline, please refer to Law et al. (2015); Yan et al. (2016); Law et al. (2016); Wake et al. (2017); Belfiore et al. (2019); Westfall et al. (2019).All the stellar masses and SFRs in this paper come from the Pipe3D value added catalog (Sánchez et al. 2016a,b).The SFRs in Pipe3D are estimated from the H flux and is dust and aperture corrected.
For more details, please refer to (Sánchez et al. 2016b).Finally, this paper made extensive use of the Marvin software in order to access MaNGA data (Cherinka et al. 2019).

Galaxy Zoo and the Legacy Survey
We have used the Galaxy Zoo (GZ) project to obtain morphological classifications and find weak and strong bars.Here, citizen scientists classify galaxies according to a decision tree (Lintott et al. 2008(Lintott et al. , 2011)).We made use of the latest iteration of Galaxy Zoo, namely Galaxy Zoo DESI (GZ DESI, Walmsley et al., in prep.).GZ DESI sources images from the DESI Legacy Imaging Surveys2 (Dey et al. 2019), which consists of three individual projects: the Dark Energy Camera Legacy Survey (DECaLS), the Beĳing-Arizona Sky Survey (BASS) and the Mayall z-band Legacy Survey (MzLS), which covers ∼ 14, 000 deg 2 of sky.As shown by Géron et al. (2021), the DESI Legacy Imaging Surveys are sufficiently deep so that weak bars are visible and can be identified by the volunteers (the median 5 point source depth of DECaLS is r = 23.6,Dey et al. 2019).GZ DESI uses classifications from citizen scientists to train machine classifications based on the Bayesian convolutional neural networks described in Walmsley et al. (2022), which we rely on for our morphology measurements.The decision tree of GZ DESI, up to the bar question, is shown in Figure 6.Note that volunteers will only reach the bar question after they identified the target as being a disk that is not edge-on.

Inclination, bar length and position angles
In order to perform the TW method, we need multiple additional parameters.We need the inclination of the galaxy, the kinematic position angle of the galaxy, the length of the bar and the position angle of the bar.

Position angle
Various papers using the TW method often use photometric approaches to obtain the position angle of the galaxy, such as fitting ellipses in the sky plane (Aguerri et al. 2015;Cuomo et al. 2019;Guo et al. 2019).The outermost isophotes are then used to estimate the position angle.However, multiple issues are associated with this method.First of all, there is no systematic way to determine which and how many outer isophotes to use to determine the position angle.Secondly, the presence of a bar, especially a strong bar, will influence these estimates.Additionally, spiral arms, rings, companion galaxies and foreground stars will all also affect these measurements significantly.However, MaNGA allows us to use kinematic position angles rather than photometric ones.We obtained the global kinematic position angle and its uncertainty by using the Python package PaFit on the stellar velocity maps from MaNGA.Small-scale disturbances in the velocity field are removed by smoothing the stellar velocity maps using a 10x10 pixel sliding window that calculates the median at every position, before the kinematic position angle is calculated.This package is based on the method detailed in Appendix C of Krajnović et al. (2006) 3 .It constructs a bi-anti-symmetric map based on the original input.The position angle that minimises the difference between the original velocity map and the symmetrised map is considered to be the best-fit global kinematic position angle.The error on the best-fit position angle is defined as the range of angles for which the difference in  2 with the best-fit angle is less than 9 (Krajnović et al. 2006), which corresponds to a 3 confidence limit.

Galaxy Zoo DECaLS
Figure 6.The decision tree of GZ DESI up to the bar question.It is worth noting that volunteers will only reach the bar question after they said the target is a disk galaxy that is not viewed edge-on.The full decision tree is shown in Walmsley et al. (2022).
Multiple other papers have successfully used this code to study galaxy kinematics before (Cappellari et al. 2007;Krajnović et al. 2011).The kinematic position angle does not suffer from the issues that plague the photometric position angle, which is why we used the kinematic one in this work.However, the kinematic position angle is not infallible.The bar can twist the inner parts of the disc velocity field, which can affect the measurement of the kinematic position angle.Additionally, the coverage of the IFU, inclination of the galaxy and the difference in position angle between the bar and disc will have some influence as well.These effects are described in more detail in Appendix A3 of Guo et al. (2019).

Inclination
Estimating the inclination of a barred galaxy, especially a strongly barred one, is not straight forward.As a strong bar is a very obvious component, it will make the galaxy appear more inclined.Therefore, we carefully measured the inclinations of our barred galaxies ourselves, using the elliptical isophote analysis technique described by Jedrzejewski (1987) using the Python package photutils4 on the r-band images from the Legacy Survey.We typically averaged the ellipticity of the outermost 5% of fitted isophotes, which usually corresponded to 5 isophotes.However, to guarantee the bar does not affect our measurement, we excluded any isophotes that are within the bar region.This meant that we used less than 5 isophotes for some targets that had long bars.The ellipticity profiles and r-band images of all our targets were inspected individually to make sure that the final value was correct.However, spiral arms, foreground stars and rings will all bias this measurement to some degree.To estimate the error on this ellipticity measurement, we correctly combine the errors associated with the isophotes used to calculate the ellipticity.

Bar length and bar position angle
There are multiple ways to determine the length of the bar.One possibility involves ellipse fitting again (Laine et al. 2002;Erwin 2005;Marinova & Jogee 2007;Aguerri et al. 2009), but as bars are associated with spiral arms, rings and ansae, this method is prone to inconsistencies.Other methods include Fourier decomposition (Aguerri et al. 2000) and using explainable artificial intelligence and saliency mapping techniques (Bhambra et al. 2022).Additionally, the Galaxy Zoo:3D project provides bar masks for galaxies in MaNGA, based on SDSS images, which can be used to estimate bar length (Masters et al. 2021).
A more straight-forward approach is to manually measure bar lengths.Manual bar length measurements have been successfully used in various studies (Erwin 2019;Géron et al. 2021).Additionally, Hoyle et al. (2011) have found that manual bar length measurements between different volunteers agree within 10% of each other and they show that manual measurements can be unbiased and robust against systematic effects.Additionally, Díaz-García et al. (2016a) have found that their manual bar lengths agree with bar lengths determined by various other automated techniques.
Thus, manual bar length measurements are used in this work.The bar lengths were measured manually by one of the authors (TG) on grz-images obtained from the Legacy Survey.A number of measures were put in place to make sure these measurements were done as consistently and correctly as possible.For example, the order of the measurements was completely randomised, so that it was was not known whether the bar that was being measured was classified as a strong or weak bar by GZ.The measurements themselves were done in DS9 (Joye & Mandel 2003) with a measurement tool that automatically records the distance measured.Additionally, every bar was measured twice and the final bar length distribution is modelled by a Gaussian centered around the average of the two measurements and with an uncertainty equal to half the difference between the two measurements.Finally, all measurements were inspected again afterwards to make sure no mistakes were made.These bar lengths were successfully used before in Géron et al. (2021), where they were compared to another bar length catalog (Hoyle et al. 2011).The bar lengths are deprojected using the method described by Gadotti et al. (2007): where  is the inclination of the galaxy,  is the difference between the position angle of the bar and of the galaxy and R b,obs and R b,deproj are the observed and deprojected bar lengths, respectively.The position angles of the bar were also obtained from these manual measurements and are similarly modelled by a Gaussian centered around the average of the measurements and with an uncertainty equal to half the difference between the two measurements.
The uncertainties on the inclination, disk PA, bar length and bar PA are used to estimate the uncertainty on the bar pattern speed, corotation radius and R.This is achieved by assuming Gaussianity over these input parameters and performing a Monte Carlo simulation with 1,000 iterations.
An overview of all the input parameters for 50 randomly selected targets is given in Table 1.The full table can be found online here 5 .

Rotation curve
As mentioned in Section 2.3, we need the rotation curve of the galaxy in order to obtain the corotation radius.The rotation curve can be determined from the stellar velocity IFU data from MaNGA.We look at the spaxels in a 5 arcsec aperture along the position angle 5 https://doi.org/10.5281/zenodo.7567945 of the galaxy.The true stellar velocity in every spaxel is calculated from the observed stellar velocity by doing: where V obs and V rot are the observed and true velocity in that spaxel,  is the inclination of the galaxy and  is the azimuthal angle measured relative to the position angle of the galaxy.The distance to the centre of the galaxy is deprojected using equation 4. The corrected velocities and deprojected distances are used to fit a two parameter arctan function, described in Courteau (1997): where V sys is the systemic velocity, V c is the asymptotic velocity, r 0 is the spatial centre of the galaxy and r  is the transition radius.The rotation curve flattens at r  and goes towards V c in this model.For our purposes, V sys and r 0 are assumed to equal zero.

Sample selection
We use the machine classifications from GZ DESI (Walmsley et al. in prep.).GZ works based on a decision tree structure.As you can see in Figure 6, this means that the question "Is there a bar feature through the centre of the galaxy?" is only answered when the galaxy is a not edge-on disk galaxy.To guarantee reliable bar classifications, we must apply additional thresholds on the fraction of people that would have been asked the bar question (N bar6 ) and the fraction of people that would have voted for a certain answer (e.g.,  strong bar ), as predicted by the automated classifications.We choose to apply  features/disk ≥ 0.27,  not edge-on ≥ 0.68 and  bar ≥ 0.5.For more information on these thresholds, please refer to Géron et al. (2021) and Walmsley et al. (2022).These thresholds resulted in a sample of 3,106 galaxies that consists of relatively face-on disk galaxies with reliable bar classifications.The same classifications are also used to assign a bar type (no bar, weak bar or strong bar) to every galaxy.The galaxy had no bar if  strong bar + weak bar < 0.5.If this was not the case and if  weak bar ≥  strong bar , then the galaxy had a weak bar.Otherwise, it had a strong bar.This classification scheme was used before in Géron et al. (2021) and is shown in Table 2.The galaxies that are identified as unbarred were removed from our sample, which reduced the sample size to 1,679 barred galaxies.
In order to avoid selection effects, we work with a volumelimited sample by imposing additional thresholds on the redshift (0.01 <  < 0.05) and absolute r-band magnitude ( r < −18.96, values obtained from the NASA-Sloan Atlas), which removed 519 galaxies from our sample.Limitations of the TW method (see Section 2.2) also impose a few additional thresholds on our sample selection.The TW method is also not developed for galaxies with irregular kinematics, as one of the main assumptions of the TW method is the existence of a well-defined pattern speed.The stellar velocity field of every galaxy was inspected by eye and 475 irregular galaxies were removed.Additionally, the bar cannot align with the disc major or minor axis.Thus, galaxies where the PA of the bar was within 10 • of the major or minor axis of the galaxy were Table 1.The plate-ifu number, right ascension, declination, inclination, position angle of the disk, position angle of the bar, the (projected) bar radius, redshift and the bar type for 50 randomly selected galaxies.The full table can be found online here.removed from our sample, which affected 193 galaxies.The TW method only works on galaxies with intermediate inclination, so we limit our sample to galaxies with inclinations between 20 • and 70 • , which removed a further 32 galaxies.

Plate
As our methodology requires to reliably perform a linear fit in the  against  plots, we require each galaxy to have at least three pseudo-slits.Similarly, to ensure the robustness of the linear fit, we use normalized root mean squared error (NRMSE) to estimate the fit quality.We only included targets that had a median NRMSE of all the MC iterations lower than 0.2.Finally, we fit the rotation curve of the galaxy with a two-parameter arctan function (see Section 3.4).However, some galaxies are not described correctly by this function, especially galaxies with highly irregular kinematics (which have mostly already been removed by this point).Thus, a threshold of median NRMSE < 0.2 is imposed on this fit as well.These threshold values for NRMSE were chosen after careful visual inspection of their fits.Additionally, targets where more than 10% of the MC iterations were unable to provide a value for pattern speed (e.g., due to not being able to place enough pseudo-slits), were excluded as well.Applying these last thresholds result in a final sample that contains 225 galaxies, with 122 strongly barred and 103 weakly barred galaxies.
Reliable bar pattern speed estimates were obtained for all these targets.However, it was found that for a small subset of these targets, especially for those with low pattern speeds, estimating the corotation radius and R is difficult.This was because the corotation radius was so high, that it fell far outside the MaNGA field of view.It was judged that extrapolating the velocity curves too much results in unreliable estimates for the corotation radius.Therefore, we excluded any targets where we had to extrapolate by more than a factor of two.This affected 15 of our 225 galaxies.

Bar pattern speeds, corotation radii and R
The final bar pattern speeds of all our weakly and strongly barred galaxies are shown in Figure 7.As mentioned in Section 3.3.2,inclination measurements are prone to biases, so to be cautious we show both Ω b sin () (top row) and Ω b (bottom row).Additionally, the pattern speed is measured in observational units (km s −1 arcsec −1 , left column), which are then converted to physical units (km s −1 kpc −1 , right column).
Despite the distributions for the weakly and strongly barred samples overlapping considerably, an Anderson-Darling test reveals that they are still significantly different, with strongly barred galaxies having lower average pattern speeds than weak barred galaxies.The p-values for the Ω b sin () distributions are 0.003 (which corresponds to 3.0) and <0.001 (>3.3) for the plots using observational units and physical units, respectively.The p-values for the Ω b plots are 0.013 (2.5) and 0.002 (3.1), respectively.
We can conclude that strongly barred galaxies have significantly lower bar pattern speeds than weakly barred galaxies.The median, together with the 16 th and 84 th percentiles, is Ω b = 23.36 +9.25 −8.1 km s −1 kpc −1 for strongly barred galaxies and Ω b = 25.91 +10.42   −7.26   km s −1 kpc −1 for weakly barred galaxies.
The final corotation radii for our target galaxies are shown in Figure 8.The median values are R CR = 8.33 +4.57−3.31 kpc and R CR = 7.19 +3.82  −2.96 kpc for the strongly barred and weakly barred sample, respectively.The results from a two-sample Anderson-Darling test show that the distributions of the corotation radii in physical units are not significantly different between weak and strong bars (p-value = 0.012; 2.5), as the significance is below 3.
With a p-value of 0.001, which corresponds to 3.3, we conclude that strong bars have significantly lower values for R than weak bars, as shown in Figure 9.However, please note that the distributions still overlap significantly.The median value for strongly barred galaxies is R = 1.53 +0.74  −0.53 and R = 1.88 +1.08 −0.75 for weakly barred galaxies.
As mentioned above, R is used to divide bars into ultrafast (R < 1.0), fast (1.0 < R < 1.4) and slow (R > 1.4).Most bars in our sample seem to be slow bars (62% of our sample), as shown in Table 4.This fraction is higher than what most other studies find.Conversely, we find less ultrafast and fast bars (11% and 27%, respectively) than most other studies.
The final values for pattern speeds, corotation radii and R of 50 randomly selected galaxies is shown in Table 3.The full table can be found online here7 .

Relationship between the parameters
The bar pattern speed and corotation radius should be inversely proportional to each other.A higher pattern speed will result in a steeper gradient of the straight line shown in Figure 4, resulting in it intersecting with the rotation curve at a shorter distance, which produces a lower corotation radius (see Section 2.3 for more details).This is shown explicitly in Figure 10, where the inversely proportional relationship becomes very clear.
We can also see that galaxies with the highest values of R tend to have lower values for the bar pattern speed.This also makes sense, as low pattern speeds will result in larger corotation radii, which increases R. Conversely, galaxies with lower R tend to have lower values for the corotation radius.
Figure 11 shows the corotation radius plotted against the bar radius.Lines of equal values of R are found diagonally over this Figure, as R is defined as the ratio between the corotation radius and the bar radius.Thus, we can divide this figure into three regions, one with all the ultrafast bars, one with all the fast bars and one with all the slow bars.The bar pattern speed is shown with the colour gradient.The galaxies with the fastest pattern speeds mostly have low values for the bar radius and corotation radius.The galaxies with the lowest pattern speeds tend to have higher values for the corotation radius, as well as higher values for R.
The bar radius is plotted against the bar pattern speed, corotation radius and R in Figure 12.We find that the bar pattern speed decreases as bar radius increases.Though the Spearman correlation index is quite small (R = -0.31)due to the high amounts of scatter, its significance is high (4.7).A more careful look reveals that all the largest bars have lower values for the bar pattern speeds.The sample is divided into strongly barred (orange) and weakly barred (blue).The MC is done over observational units (which include arcsec), which are afterwards converted to kpc.The left column shows the results for the observational units, while the right column shows the results for the physical units.The vertical dashed lines show the median values for every histogram.The full lines are kernel density estimates of these histograms, using a Gaussian kernel.The p-value of a two-sample Anderson-Darling test is shown inside each subplot, with the null hypothesis being that the two samples are drawn from the same population.We see that, on average, strongly barred galaxies have significantly lower bar pattern speeds, despite there being significant overlap between the two populations.The p-value of the comparisons in physical units is <0.001-0.002(which corresponds to 3.1-3.3).
Figure 8.The final median values for R CR for every galaxy, after doing a Monte Carlo simulation of 1,000 iterations.The sample is divided into strongly barred (orange) and weakly barred (blue).The MC is done over observational units (which is in arcsec), which are afterwards converted to kpc.The left column shows the results for the observational units, while the right column shows the results for the physical units.The vertical dashed lines show the median values for every histogram.The full lines are kernel density estimates of these histograms, using a Gaussian kernel.The p-value of a two-sample Anderson-Darling test is shown inside each subplot, with the null hypothesis being that the two samples are drawn from the same population.As the p-value of the comparison in physical units is 0.012 (which corresponds to 2.5), we conclude that we see no significant difference between weak and strong bars in terms of their corotation radii.Conversely, all bars with higher values for their pattern speed are relatively short.
The corotation radius increases with bar radius (R = 0.5; 7.75).This is because larger bars need to have a larger corotation radius, as a bar can only grow up to its corotation radius.R is observed to decrease with bar radius (R = -0.46;6.97).However, R is very sensitive to correct bar length estimates, so this trend could merely be a reflection of that.Interestingly, the median trends for the weakly and strongly barred subsamples are very similar to each other.
Many properties of galaxies vary with stellar mass (Brinchmann & Ellis 2000;Brinchmann et al. 2004;Noeske et al. 2007;Lara-López et al. 2010) and bars, especially stronger bars, are known to appear more often in massive galaxies (Masters et al. 2012;Cervantes Sodi 2017;Géron et al. 2021).In Figure 13, we plot the stellar mass against the bar pattern speed, corotation radius and R to see if any of these parameters are correlated with stellar mass as well.
We see that the pattern speed and R do not correlate with stellar mass (R = -0.05;0.75 and R = -0.12;1.73, respectively).This shows that the differences we observed in pattern speed in Figures 7 and 9 are not due to differences in stellar mass of our targets.Interestingly, the corotation radius does increase with stellar mass (R = 0.48; 7.44).This is because more massive galaxies tend to host stronger and longer bars, which tend to have larger corotation radii, as shown in the middle panel of Figure 12.

Quenching
Our previous results reveal a complicated interplay between bar pattern speed, corotation radii, R, bar length and bar type.It is also known that strong bars are more often found in red sequence galaxies (Masters et al. 2012;Vera et al. 2016;Géron et al. 2021).This suggest a potential link between these dynamical parameters and quenching.Galaxies can be classified as star forming or quiescent based on their location on the SFR-stellar mass plane.We can use the star formation main sequence (SFMS) defined in Belfiore et al. ( and assume that all galaxies that are 1 (= 0.39 dex) below this line are undergoing quenching and everything else is star forming (Belfiore et al. 2018).57% of our barred galaxies are quenching, whereas 43% are star forming.The bar pattern speeds, corotation radii and values for R for all the star forming and quenching galaxy subsamples are shown in Figure 14.An Anderson-Darling test between the subsamples shows that the pattern speed and R are not significantly different (both are <3).However, with a p-value of <0.001, which corresponds to >3.3, the corotation radii are significantly different.Thus, our results suggest that quenching galaxies tend to have significantly higher corotation radii than star forming galaxies.
Table 3.The bar pattern speeds, corotation radii and R for 50 randomly selected bars.The full table can be found online here.

Plate-ifu
11014-12705 9.40 +0.99All strongly barred galaxies are coloured orange, while all weakly barred galaxies are coloured blue.The median trend for the weakly and strongly barred galaxies is shown with the blue and orange full lines, respectively.Additionally, the general median trend of all barred galaxies is shown in the dashed black line.The Spearman correlation coefficient, R, and its significance, , are shown in every subplot.The median error on the x and y axis is shown in the top-right corner.We see that the pattern speed and R decrease with bar length, while the corotation radius increases.[-] R = -0.12,= 1.73 Figure 13.The bar pattern speed (left panel), corotation radius (middle panel) and R (right panel) against the stellar mass.All strongly barred galaxies are coloured orange, while all weakly barred galaxies are coloured blue.The median trend for the weakly and strongly barred galaxies is shown with the blue and orange full lines, respectively.Additionally, the general median trend of all barred galaxies is shown in the dashed black line.The Spearman correlation coefficient, R, and its significance, , are shown in every subplot.The median error on the x and y axis is shown in the top-right corner.We see that the pattern speed and R do not show a significant trend with stellar mass, while the corotation radius seems to increase with stellar mass.

Comparison with other work
Various other studies have also tried to measure bar pattern speeds, corotation radii and R. Rautiainen et al. (2008) determine pattern speeds, corotation radii and R for a sample of 38 galaxies with data from the Ohio State University Bright Spiral Galaxy Survey (Eskridge et al. 2002).Aguerri et al. (2015) used the Calar Alto Legacy Integral Field Area (CALIFA, Sánchez et al. 2012) survey on 15 galaxies and found that all of their bars were consistent with being fast.Font et al. (2017) combined Spitzer images of 68 barred galaxies with previously determined corotation radii to estimate values for R. Cuomo et al. (2019) looked at 16 weakly barred galaxies using data from CALIFA.Guo et al. (2019) used MaNGA data to obtain estimates for pattern speeds, corotation radii and R for a total of 53 barred galaxies.Garma-Oehmichen et al. ( 2020) combined data from MaNGA and CALIFA to study a sample of 18 galaxies.Finally, Garma-Oehmichen et al. ( 2022) used MaNGA to study 97 barred galaxies.
We compare our results with these studies in Figure 15.Our distribution of the bar pattern speed, corotation radius and R falls well within the range that is usually observed and we see no obvious deviations.Our bar pattern speed distribution agrees especially well with the other studies that have larger sample sizes (n > 50).An interesting trend is observed when looking at the various distributions of R. The median value of R seems to be moving upwards as the sample sizes increase.This could be attributed to larger samples typically being more representative of large variety of galaxy and bar types, which could have an effect on the observed distribution of R. The sample is divided into quenching galaxies (red) and star forming galaxies (blue).The vertical dashed lines show the median values for every histogram, while the full lines are kernel density estimates of these histograms, using a Gaussian kernel.The p-value of a two-sample Anderson-Darling test is shown in the top-right corner of every subplot.The null hypothesis is that the two samples in each subplot are drawn from the same population.We can see that the quenching and star forming subsamples are not significantly different in terms of pattern speed and R, but are in terms of corotation radius.
Please refer to Appendix A for a similar comparison to various other studies, but distinguishing between weakly and strongly barred galaxies as well.

Are strong bars older than weak bars?
We have found in Figure 7 that strongly barred galaxies have significantly lower bar pattern speeds than weakly barred galaxies, especially in terms of Ω b sin () (p-value < 0.001; >3.3).However, there is still a large overlap between the two samples.It is worth noting that this difference is not due to differences in stellar mass, as shown in Figure 13.Additionally, we found that pattern speed is negatively correlated to bar length in Figure 12.This is in agreement with Cuomo et al. (2020), who used the CALIFA and MaNGA surveys and also found that stronger bars have lower bar pattern speeds.Font et al. (2017) also found that the largest bars have the lowest pattern speeds, while the bars with the largest pattern speeds are all very small.Using CALIFA, MaNGA and Pan-STARRS DR1 (PS1), Lee et al. (2022) also found that the bar pattern speed is negatively correlated to bar length and strength.
Simulations suggest that the bar grows in size and the pattern speed slows down as the bar exchanges angular momentum with its host (Debattista & Sellwood 2000;Athanassoula 2003;Martinez-Valpuesta et al. 2006;Okamoto et al. 2015).Thus, our results suggest that strong bars are older and more evolved structures than weakly barred galaxies.Alternatively, stronger bars could be simply more efficient at redistributing angular momentum, as this depends on various parameters, such as the velocity dispersion and mass distribution of the emitting and absorbing components (Athanassoula 2003).

How are strong and weak bars triggered?
We found that strongly barred galaxies have statistically significantly lower values for R (p-value = 0.001; 3.3), which was defined as R =  CR / bar , than weakly barred galaxies, as shown in Figure 9.However, it is important to note that there is still a big overlap in both distributions.Additionally, we find a relationship between R and bar length in Figure 12, which shows that long bars have lower values for R.However, R is very dependent on the bar length and inclination estimates (Cuomo et al. 2021;Roshan et al. 2021a), so this trend could be primarily caused by differences in bar length between weak and strong bars.A more detailed study using different metrics to measure bar length could help clarify this issue.
We know from simulations that R depends on the formation of the bar.Bars triggered by tidal interactions tend to have higher values for R than bars formed by global bar instabilities.Additionally, tidally induced bars stay in the slow regime for a longer time (Sellwood 1981;Miwa & Noguchi 1998;Martinez-Valpuesta et al. 2016, 2017).This seems to suggest that strong bars are more likely to be triggered by bar instabilities, whereas weak bars are more likely to be formed by tidal interactions.This statement could possibly be tested observationally by looking at the environment of a large mass-matched sample of strongly and weakly barred galaxies.Interestingly, this is not in agreement with Cuomo et al. (2020) and Guo et al. (2019), who found no relationship between bar strength and R, although this could be caused by the lower sample sizes used in these studies.

Why and where do we see ultrafast bars?
Ultrafast bars, which are bars that have R < 1, should not exist according to our current theoretical understanding.This is because bars are thought to not be able to extend beyond the corotation radius of the galaxy (Contopoulos 1980(Contopoulos , 1981;;Athanassoula 1992b).Nevertheless, they have been found observationally.Multiple studies find that 26 -67% of their bars have R < 1, while 7-40% have a 1 upper limit that has R < 1 (Aguerri et al. 2015;Cuomo et al. 2019;Guo et al. 2019;Garma-Oehmichen et al. 2020).In our sample, ∼11% of galaxies have R < 1, while ∼2% have a 1 upper limit that has R < 1.This is significantly lower than what others studies typically find.A more detailed breakdown can be found in Table 4.
Given that R is defined as R =  CR / bar , a low value for R can  arise either because  CR is underestimated or  bar is overestimated.Both seem to happen simultaneously: even though most ultrafast bars have a lower corotation radius (see Figure 10), we see in Figure 11 that ultrafast bars still have a relatively broad range of bar and corotation radii.Bars are often associated with spiral arms and rings, so it is often not straight-forward to measure the bar length correctly (Hilmi et al. 2020;Cuomo et al. 2021;Roshan et al. 2021a).An underestimation of the corotation radius can be a consequence of either an overestimation of the pattern speed, or because the rotation curve is not fitted properly.Additionally, if the inclination of the galaxy is not measured properly, the line of sight velocities in the rotation curve will be corrected incorrectly, which will affect the corotation radius as well.Interestingly, we see slightly more ultrafast bars among strong bars than weak bars (13.9% and 7.8%, respectively).

Slow bars
As noted in Table 4, we find that ∼62% of galaxies have a slow bar.Our fraction of slow bars is higher than what other studies have typically found.This can be due to multiple factors, such as our larger and more representative sample, which includes weak and strong bars from a volume-limited sample, from a wide range of magnitudes.Additionally, a correct measurement of the bar length is crucial to a correct estimate of R (Cuomo et al. 2021).Differ-ent authors use different methods of measuring bar length, which will change the final distribution of R (for more details on our bar length measurements, see Section 3.3).However, perhaps most importantly, Guo et al. (2019) have shown that estimates of the pattern speed will be systematically lower when using kinematic position angles, compared to photometric position angles.This is because the method to calculate the kinematic position angle works by minimising asymmetry in the velocity field, which will reduce the values for the kinematic integrals.This will lower the estimates for the patterns speed, which will, in turn, increase the estimates for R and produce more slow bars.We used kinematic position angles in this work, which could partially explain the higher fraction of slow bars.Finally, if we take the errors into account, we find that ∼35% of our targets have a 1 lower limit that is greater than 1.4.Thus, we can confidently exclude the fast regime for only these targets, which is more consistent with other studies.

Strong and weak: part of a continuum
We have found that strong bars tend to have lower bar pattern speeds.However, as Figure 7 shows, their distributions overlap significantly.There is no clear threshold in pattern speed, corotation radius or R that separates weak and strong bars.A closer look at Figure 12, where we plot these parameters against bar radius, reveals that the differences in pattern speed and R are driven by the smallest and largest bars.At intermediate bar radii, the distributions of the two populations overlap.Additionally, the median trend of the weakly and strongly barred subsamples are almost identical in these figures.Figure 8 shows that weakly and strongly barred galaxies do not have significantly different corotation radii.As bars are able to grow up until their corotation radius (Contopoulos 1980(Contopoulos , 1981;;Athanassoula 1992b), this result suggests that weak bars still have the possibility to grow up to the same length as strong bars.Either they have not had to time to do so yet or something else is preventing them.
These results are consistent with the idea of a bar continuum, proposed by Géron et al. (2021), who found that any distinction in fibre SFR, gas mass and depletion timescale between weak and strong bars disappeared when correcting for bar length.They suggested that weak and strong bars are not fundamentally distinct physical phenomena.Instead, bar types are continuous, and vary from 'weakest' to 'strongest'.Our measurements of the dynamical parameters of weak and strong bars support this conclusion as well.
As mentioned in Section 2.2, we require that at least three pseudo-slits can be placed on the bar.This means that this method does not work for the shortest of bars, which are typically weak bars.Thus, it should be kept in mind when interpreting our results that the shortest and weakest bars are not included in our analysis.This manifests itself as well in the errors associated with the measured pattern speeds.The highest errors are associated with bars with the fewest pseudo-slits.This means that weak bars typically have higher errors for the bar pattern speed than strong bars: 6.83 km s −1 kpc −1 and 4.89 km s −1 kpc −1 , respectively.

Dark matter haloes and tension with ΛCDM
The fraction of ultrafast and fast bars found in this work is less than what some other studies tend to find, but combined they still comprise 38% of our sample.Observing a high fraction of ultrafast and fast bars has been raised as a tension for the ΛCDM cosmological paradigm, as cosmological simulations predict that bars should slow down significantly (Algorry et al. 2017;Peschken & Łokas 2019;Fragkoudi et al. 2021;Roshan et al. 2021a,b;Frankel et al. 2022).This slowdown of the bar and increase of R is typically attributed to the dynamical friction applied to the bar by the DM halo (Debattista & Sellwood 1998, 2000;Fragkoudi et al. 2021).
As we find fewer ultrafast and fast bars than other studies, this tension is decreased somewhat, but the median value of R in our sample (R = 1.66 +1.05 −0.62 ) is still significantly lower than that predicted from simulations, whose average values at z∼0 are typically R > 2.5 (Algorry et al. 2017;Peschken & Łokas 2019;Roshan et al. 2021b).
Other studies have tried to relieve the tension as well.Frankel et al. (2022) has recently shown that simulations obtain higher values of R than observations, mostly because simulations predict shorter bars, rather than slower bars.Additionally, Fragkoudi et al. (2021) actually do find fast bars in their cosmological simulations in baryon-dominated discs and claim that the DM fraction is too high in other simulations.A lower DM fraction or lower central DM density will lower the dynamical friction, and thus, allow fast bars to exist.Finally, Beane et al. (2022) have shown that the gas phase of the disk can help to stabilise the bar pattern speed and prevent it from slowing down.
As mentioned above, R is significantly lower for strong bars than for weak bars.This suggests that the DM fractions of strong and weak bars are different as well.Studying the relationship between the DM halo and bars will help us understand the evolution of bars in general.This can be done, for example, with Jeans Anisotropic Modelling (JAM) of these galaxies (Cappellari 2008).This would provide estimates for the DM fraction and allow us to study the intrinsic connection between the DM halo and the dynamical parameters of bars in greater detail.

Effect on quenching
It is also known that strong bars are more often found in red sequence galaxies (Masters et al. 2012;Vera et al. 2016;Cervantes Sodi 2017;Fraser-McKelvie et al. 2020).In addition, Géron et al. (2021) showed that stronger bars have the ability to facilitate quenching, whereas weaker bars do not.Figure 14 shows that the bar pattern speed and R are not significantly different between star forming and quenching galaxies.These results suggest that how fast bars rotate, both in terms of pattern speed and R, has no significant or measurable impact on quenching.This seems odd, but one way to look at this is by looking at the timescales involved.A bar with a pattern speed of ∼25 km s −1 kpc −1 will make a full rotation once every ∼250 Myr.As (secular) quenching usually happens on ∼Gyr timescales (Smethurst et al. 2015), this means that a bar will usually have made multiple full rotations before the galaxy is quenched.
Interestingly, the corotation radius is significantly higher for quenching galaxies.It is known that a bar can grow up to its corotation radius (Contopoulos 1980(Contopoulos , 1981;;Athanassoula 1992b), so this result relates back to longer and stronger bars having more of an effect on quenching, which agrees with the findings of Géron et al. (2021).

CONCLUSION
We have used the TW method on MaNGA IFU data for a sample of 225 barred galaxies, which is the largest sample this method has been applied to so far.The TW method produces bar pattern speeds, which is used to calculate corotation radii and R, the ratio between the corotation radius and bar radius.We have used Galaxy Zoo morphological classifications to distinguish between strongly and weakly barred galaxies.This allows us to study the bar pattern speed, corotation radius and R for a statistically significant sample of strongly and weakly barred galaxies, and compare them with each other.We have found that: • Though there is significant overlap, we find that the bar pattern speeds between weakly and strongly barred galaxies are significantly different, as the p-value of the comparisons in physical units is <0.001-0.002(which corresponds to 3.1-3.3).The median bar pattern speed of strongly barred galaxies (Ω b = 23.36 +9.25  −8.1 km s −1 kpc −1 ) is lower than that of weakly barred galaxies (Ω b = 25.91 +10.42  −7.26 km s −1 kpc −1 ).Additionally, we find that the bar pattern speed is inversely proportional to bar length.We also show that this difference is not due to differences in stellar mass in our targets.Simulations suggest that the pattern speed goes down as the bar evolves and exchanges angular momentum, so our results suggest that strong bars are older and more evolved structures than weak bars.
• We could not find evidence that the corotation radius between weakly (R CR = 7.19 +3.82  −2.96 kpc) and strongly barred galaxies (R CR = 8.33 +4.57−3.31 kpc) is significantly different, as the p-value of the comparison in physical units is 0.012 (which corresponds to 2.5).As bars can grow up until their corotation radius, this result suggests that weak bars still have the possibility of becoming as long as strong bars.
• Despite the significant overlap in the distributions, we find that R is statistically significantly lower for strong bars than for weak bars (p-value = 0.001; 3.3).The median value for strong bars is 1.53 +0.74  −0.53 , while it is 1.88 +1.08 −0.75 for weak bars.Additionally, we find that R is inversely proportional to bar length and that these differences are not caused by differences in stellar mass.As R is related to the formation of bars, this suggests that weak bars are more likely to be formed by tidal interactions, whereas strong bars are more likely to be triggered by global bar instabilities.
• We do not see a distinct cutoff or threshold in pattern speed, corotation radius or R that separates strong and weak bars from each other.In fact, the overlap is still quite significant.This is consistent with Géron et al. (2021), who stated that strong and weak bars are not distinct physical phenomena, but rather lie on a continuum of bar types that vary from 'weakest' to 'strongest'.
• 11% of our sample host ultrafast bars, 27% host fast bars and 62% of all bars are slow.We have a slightly higher fraction of slow bars and a lower fraction of ultrafast bars than most other studies.This can be attributed to various factors, such as the bigger and more representative sample used in this study.
• However, only ∼2% of our targets have a 1 upper limit that has R < 1.Similarly, we can only confidently exclude the (ultra)fast regime for ∼35% of our galaxies (i.e. they have a 1 lower limit that has R > 1.4) • The lower fraction of ultrafast bars among our sample decreases the recent tension with ΛCDM.However, the median value of R in our sample is still significantly lower than what is predicted from simulations.
• We do not see any significant difference between the star forming and quenching subsamples in terms of pattern speed or R.However, quenching galaxies do have significantly higher corotation radii than star forming galaxies.

Figure 1 .
Figure 1.DECaLS postage stamps (64x64 arcsec) of a strongly (left) and weakly barred galaxy (right), on which we will apply the Tremaine-Weinberg method.

Figure 2 .
Figure 2. The stellar flux (top row) and stellar velocity (bottom row) for our strongly barred (left column) and weakly barred (right column) example galaxies.The different pseudo-slits, over which the kinematic and photometric integrals are calculated, are visualised on top of the maps in white outlines.

Figure 3 .
Figure3.The kinematic integral,  , is plotted against the photometric integral,  , for all the pseudo-slits for a strongly barred galaxy (left) and a weakly barred galaxy (right).The equation of the best-fit line going through these points is shown in each plot.The slope of these lines is equal to Ω b sin ().

Figure 4 .
Figure 4. Visualisation of how the corotation radius is obtained, for a strongly barred galaxy (left) and a weakly barred galaxy (right).The blue dots are stellar velocity measurement from MaNGA in a 5 arcsec aperture along the major axis of the galaxy.The black line is the best-fit rotation curve, please refer to Section 3.4 for more details on how the rotation curve is calculated.The black dashed line is obtained by multiplying Ω b with a radius range.The distance where this line and the rotation curve intersect (indicated by the orange cross), defines the corotation radius.The dotted vertical line is the deprojected bar radius.

Figure 5 .
Figure 5.The posterior distributions of the pattern speed (left column), corotation radius (middle column) and R (right column) for a strongly barred galaxy (top row) and weakly barred galaxy (bottom row), obtained from performing Monte Carlo (MC) simulations of a 1,000 iterations in order to characterise the uncertainty on each measurement.The median value is indicated in every plot by a black vertical line, while the 16 th and 84 th percentile are shown by the dashed vertical lines.These values are also printed in each subplot.

Figure 7 .
Figure 7.The final median values for Ω b sin () (top row) and Ω b (bottom row) for every galaxy, after doing a Monte Carlo simulation of 1,000 iterations.The sample is divided into strongly barred (orange) and weakly barred (blue).The MC is done over observational units (which include arcsec), which are afterwards converted to kpc.The left column shows the results for the observational units, while the right column shows the results for the physical units.The vertical dashed lines show the median values for every histogram.The full lines are kernel density estimates of these histograms, using a Gaussian kernel.The p-value of a two-sample Anderson-Darling test is shown inside each subplot, with the null hypothesis being that the two samples are drawn from the same population.We see that, on average, strongly barred galaxies have significantly lower bar pattern speeds, despite there being significant overlap between the two populations.The p-value of the comparisons in physical units is <0.001-0.002(which corresponds to 3.1-3.3).

Figure 9 .Figure 10 .
Figure 9.The final median values for R for every galaxy, after doing a Monte Carlo simulation of 1,000 iterations.The sample is divided into strongly barred (orange) and weakly barred (blue).The vertical dashed lines show the median values for every histogram.The full lines are kernel density estimates of these histograms, using a Gaussian kernel.The p-value of a twosample Anderson-Darling test is shown inside each subplot, with the null hypothesis being that the two samples are drawn from the same population.We see that strong bars have significantly lower values of R than weak bars, despite the big overlap (p-value = 0.001; 3.3).

Figure 11 .
Figure 11.The corotation radius is plotted against the deprojected bar length.As R =  CR / bar , this figure is divided into three regions: the region with slow bars (R > 1.4), the region with fast bars (1 < R < 1.4) and the region with ultrafast bars (R < 1).The colour indicates the bar pattern speed.The median error on the x and y axis is shown in the top-left corner.To aid visualisation, the colours used to indicate the bar pattern speed were capped at the 16 th and 84 th percentile.

Figure 12 .
Figure12.The bar pattern speed (left panel), corotation radius (middle panel) and R (right panel) against the bar radius.All strongly barred galaxies are coloured orange, while all weakly barred galaxies are coloured blue.The median trend for the weakly and strongly barred galaxies is shown with the blue and orange full lines, respectively.Additionally, the general median trend of all barred galaxies is shown in the dashed black line.The Spearman correlation coefficient, R, and its significance, , are shown in every subplot.The median error on the x and y axis is shown in the top-right corner.We see that the pattern speed and R decrease with bar length, while the corotation radius increases.

Figure 14 .
Figure 14.The final median values for bar pattern speed (left panel), corotation radius (middle panel) and R (right panel) after doing a Monte Carlo simulation of 1,000 iterations.The sample is divided into quenching galaxies (red) and star forming galaxies (blue).The vertical dashed lines show the median values for every histogram, while the full lines are kernel density estimates of these histograms, using a Gaussian kernel.The p-value of a two-sample Anderson-Darling test is shown in the top-right corner of every subplot.The null hypothesis is that the two samples in each subplot are drawn from the same population.We can see that the quenching and star forming subsamples are not significantly different in terms of pattern speed and R, but are in terms of corotation radius.

Figure 15 .
Figure 15.A comparison of estimates of the bar pattern speeds (left), corotation radii (middle) and R (right) found in various works.All the histograms are normalised and offset from each other vertically to facilitate comparison.The median, 25 th and 75 th percentile for every distribution are indicated by the short full and dashed lines.The studies are ordered by sample size, which is also shown in the left panel.The values from Rautiainen et al. (2008) and Aguerri et al.(2015)  were converted from the observational units cited in their papers to physical units using redshifts obtained from the NASA/IPAC Extragalactic Database (NED).The values obtained fromGuo et al. (2019) were converted using their own cited redshifts.Similarly, the values fromFont et al. (2017) were converted using the distances cited.The values for the bar pattern speeds fromRautiainen et al. (2008) are not publicly available, hence the empty histogram.

Figure A1 .
Figure A1.A comparison of estimates of the bar pattern speeds (left), corotation radii (middle) and R (right) found in various works.The data is split between weakly barred galaxies (blue) and strongly barred galaxies (orange).Font et al. (2017) have both weak and strong bars, whereas the other studies focus only on either weak or strong bars.All the histograms are normalised and offset from each other vertically to facilitate comparison.The median, 25 th and 75 th percentile for every distribution are indicated by the short full and dashed lines.The studies are ordered by sample size, which is also shown in the left panel.The values from Aguerri et al. (2015) were converted from the observational units cited in the papers to physical units using redshifts obtained from the NASA/IPAC Extragalactic Database (NED).The values obtained fromGuo et al. (2019) were converted using their own cited redshifts.Similarly, the values fromFont et al. (2017) were converted using the distances cited in the paper.

Table 2 .
Géron et al. (2021)of the galaxy are used to determine its bar type (no bar, weak bar or strong bar), according to the following scheme.This method of classification is identical to the one inGéron et al. (2021).