Free Access
Issue
A&A
Volume 631, November 2019
Article Number A94
Number of page(s) 23
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/201936000
Published online 28 October 2019

© ESO 2019

1. Introduction

The most conspicuous stellar structures of spiral galaxies are the arms that sweep out from the centre of the disc in unbarred galaxies, or from near the ends of the stellar bar. Spiral arms are sites rich in gas with intense star formation, H II regions, and dust (e.g. Schweizer 1976; Elmegreen et al. 2011). This makes them more prominent in blue bands of the spectrum, but their backbone is typically composed of old stars (e.g. Elmegreen 1981; Elmegreen & Elmegreen 1989; Knapen & Beckman 1996; Eskridge et al. 2002). Roughly two thirds of all massive galaxies are spirals (e.g. Lintott et al. 2011; Willett et al. 2013).

There are three main types of spiral arms (e.g. Elmegreen 1990): grand-design, flocculent, and multi-armed. The fraction of each of these classes in disc galaxies is ∼18%, ∼50%, and ∼32%, respectively, based on the morphological classifications made by Elmegreen et al. (2011) and Buta et al. (2015) to the Spitzer Survey of Stellar Structure in Galaxies (S4G; Sheth et al. 2010), that comprises over 2000 nearby galaxies. Grand-design galaxies (Lin & Shu 1966) are characterised by two long and well-defined arms. Flocculent spirals (Elmegreen 1981) present short and fragmented spiral arm sections. They are more common in faint galaxies (Elmegreen & Elmegreen 1985). The fairly symmetric multi-armed spirals constitute a category between grand-design and flocculent, closer to the former, characterised by a central two-armed pattern that develops long ramifications in the outer parts of the optical disc. In this work we analyse the shapes of these three types of spirals in the S4G survey.

How spiral arms are formed in galactic discs remains a matter of debate. One of the most widely accepted explanations of the formation and maintenance of spirals is the density wave theory (Lin & Shu 1964). The arms are considered as density waves that rotate with a radius-independent angular speed through the stars and gas in a shearing, differentially rotating galactic disc. In this scenario, the stars in the inner parts of the disc rotate faster than the spiral and can overtake it, whereas in the outer parts they lag with respect to the spiral pattern.

Spirals might also be tidally triggered in interactions with companion galaxies (Kormendy & Norman 1979), as it is the case of the Whirlpool Galaxy (or M 51) (see e.g. Salo & Laurikainen 2000a; Watkins et al. 2015, and references therein), the first identified spiral nebula in history (Rosse 1850) using the Leviathan of Parsonstown telescope (1.8 m aperture).

In the swing amplification mechanism (Goldreich & Lynden-Bell 1965a,b; Toomre 1981), shear and self-gravity transform initially leading waves into trailing ones, so that the wave is amplified when the stabilising effects of random motions are negligible. A constraint on the theory of the formation of spirals is the number of arms expected by swing amplification. D’Onghia (2015) showed that the total number of spiral arms anticorrelates with the disc-to-total mass fraction at 2.2 disc scalelengths: submaximal discs are expected to be multi-armed galaxies. High resolution N-body simulations by Fujii et al. (2018) confirm that the disc-mass fraction controls the number of spirals.

Grand-design spiral arms and the symmetric inner parts of multi-armed galaxies have been interpreted as long-lived spiral density waves, while flocculent galaxies and the outer parts of multi-armed galaxies have been associated to random local gravitational instabilities in the disc, swing amplified into transient spiral arms (e.g. Elmegreen et al. 2011, and references therein). According to numerical models (e.g. Sellwood 2011, and references therein), the spiral arms in flocculent galaxies are short-lived (of the order of ∼100 Myr), whereas grand-design spirals exhibit arms that last longer (∼1 Gyr). High-resolution N-body simulations by D’Onghia et al. (2013) showed that ragged spiral structure can survive long after the original perturbing influence has been removed.

How tightly wound the spiral arms are can be quantified by measuring the pitch angle (Binney & Tremaine 1987), defined as the angle subtended by the tangent to the spiral arm, relative to the tangent to a circle at a point at a given galactocentric radius. By definition, it lies between −90° and 90°. In theoretical modelling, positive (negative) pitch angles indicate trailing (leading) pattern with respect to rotation.

The pitch angle is not necessarily constant, regardless of the class of the spiral arms (e.g. Kennicutt 1981). It has been claimed that pitch angles depend on the central mass concentration and atomic gas density (e.g. Kennicutt 1981; Block et al. 1994; Davis et al. 2015; Yu & Ho 2019), on the galactic shear rate (e.g. Seigar et al. 2006; Grand et al. 2013), or on the steepness of the rotation curves (Seigar et al. 2005, 2014): spiral arms have been found to be more open in galaxies with rising rotation curves, whereas those with falling rotation curves are generally tightly wound. Recent work by Davis et al. (2017) reports a very tight scaling relation between the mass of the central supermassive black hole and the spiral pitch angle which, if true, hints at a surprisingly intimate link between the large-scale structure of discs and the mass in the nucleus (see also Seigar et al. 2008; Berrier et al. 2013).

Spiral arms are known to contribute to the rearrangement of gas that leads to the formation of disc-like bulges (e.g. Kormendy & Kennicutt 2004). This makes them important agents for the secular evolution of disc galaxies, a process in which bars also play a fundamental role. Whether spirals are bar-driven or not has been debated over the past decades, ever since Kormendy & Norman (1979) suggested that stellar bars can lead the formation of spiral density waves. The latter has been mainly tested from the correlations between gravitational torques and density amplitudes associated to bars and the spirals (e.g. Seigar et al. 2003; Block et al. 2004; Buta et al. 2005, 2009; Durbala et al. 2009; Salo et al. 2010).

A strong coupling between bars and spirals is also expected from the manifold theory, which is a relatively new paradigm for the formation of spirals and rings (e.g. Romero-Gómez et al. 2006, 2007; Patsis 2006; Voglis et al. 2006; Athanassoula et al. 2009a). Numerical simulations show that galactic material gets confined in tubes (invariant manifolds) that extend from the two unstable Lagrangian points at the end of the bar (L1 and L2). This theory predicts a dependence of the pitch angle of spirals on the bar perturbation strength (Athanassoula et al. 2009b).

The main goals of this paper are (i) to reassess the existence of scaling relations involving the pitch angle of the arms in the S4G, (ii) to analyse the coupling between the amplitudes of bars and spiral arms in all kinds of spirals galaxies with unprecedented sample size, (iii) to test the predicted correlation between bar strength and pitch angle, and (iv) to check if there is any difference in the distribution of the pitch angles for barred and non-barred galaxies.

The paper is organised as follows: in Sect. 2 we present the sample and the data used in this work. In Sects. 3 and 4 we explain our methodology to derive the global pitch angle of the galaxies as well as the strength of spirals and bars. In Sect. 5 we study the distribution of the pitch angle across the Hubble sequence, as well as the distribution of spiral types and bar frequency in our sample versus T-type. In Sect. 6 we assess the radial variation of the pitch angle, in Sect. 7 we test whether the pitch angle and amplitude of spirals are controlled by the strength of the bar, and in Sect. 8 we analyse the dependence of the pitch angle of the spirals on the central stellar mass concentration and on the mass of the central supermassive black hole. In Sect. 9 we discuss the properties of spiral arms, their dependence on disc fundamental parameters, and their hypothetical coupling with bars. Section 10 summarises the most important results of this work.

2. Data and sample selection

2.1. The S4G survey

The S4G survey (Sheth et al. 2010) consists of 2352 nearby galaxies (distance ≲40 Mpc) observed in the 3.6 μm and 4.5 μm bands with the InfraRed Array Camera (IRAC; Fazio et al. 2004) installed on board the Spitzer Space Telescope (Werner et al. 2004). Data taken from HyperLEDA were used to define the S4G sample, which is composed of bright and large galaxies (extinction-corrected total blue magnitude mBcorr <  15.5 mag and blue light isophotal angular diameter D25 >  1′) located away from the Milky Way plane (Galactic latitude |b|> 30°).

The S4G is not complete in any quantitative sense, but it is comprised of galaxies of all Hubble types (T) with stellar masses spanning ∼5 orders of magnitude. Due to selection based on H I recessional velocities, it is deficient in gas-poor early-type galaxies: this is not a serious drawback for the current study because we are not interested in S0s, that lack spiral arms. In addition, the S4G missed around 400 late-type galaxies without 21 cm systemic velocity measurements listed in HyperLEDA: we are currently obtaining i-band photometry with ground-based telescopes for the ∼50% of those galaxies that lack high-resolution deep near-IR and optical archival imaging. The current work is based on the original S4G.

Bulge-to-total mass ratios (B/T) and disc scale-lengths (hR) of S4G galaxies are available from the 2D photometric decomposition models by Salo et al. (2015). We take the 3.6 μm isophotal radii at the surface brightness 25.5 mag arcsec−2 (R25.5) from Muñoz-Mateos et al. (2015).

2.2. Morphological classification

We use the morphological classification of the galaxies in the S4G sample made by Buta et al. (2015), that includes families (SB, SA B ̲ $ \rm SA\underline{B} $, SAB, S A ̲ B $ \rm S\underline{A}B $, SA) and the revised Hubble stages (T). We also use arm classes from Buta et al. (2015) and Elmegreen et al. (2011).

2.3. Pitch angles of spiral arms

Herrera-Endoqui et al. (2015) measured the pitch angles of spiral arm segments in S4G galaxies using 3.6 μm photometry, which are analysed in this paper. Only measurements with “ok” quality flags (=1 − 2) are used.

Herrera-Endoqui et al. (2015) identified the different spiral segments in unsharp-masked images, created by dividing the original images by their smoothed versions. In unsharp-masked images one can easily identify subtle galactic structures that appear against a bright and diffuse background. These images were displayed in different scales to make the spiral arms stand out. They visually marked points tracing the spiral segments (and thus there are uncertainties in the parametrisation associated to human errors). Then, they performed a linear fit in the disc plane using polar coordinates log (r) versus θ, where logarithmic arms appear as straight lines. Thus, several measurements of the pitch angle can appear for a single galaxy, corresponding to logarithmic segments at different radial distances. In Fig. 1 we illustrate these measurements for M 51 and NGC 3992. For further details, the reader is referred to Herrera-Endoqui et al. (2015).

thumbnail Fig. 1.

3.6 μm image of NGC 5194 (upper panels) and NGC 3992 (lower panels) in sky plane (left), with fitted spiral arm segments, and logarithmic polar plot of the same galaxies de-projected to the disc plane (right). Images and measurements are taken from Herrera-Endoqui et al. (2015).

2.4. Sample of not-highly inclined spiral galaxies

For this study, we use a sample of 391 galaxies with inclinations lower than 65° (according to Salo et al. 2015), with reliable measurements of pitch angles of individual spiral segments1 (from Herrera-Endoqui et al. 2015), and classified spiral types (from Buta et al. 2015). Of these, 269 galaxies are barred (i.e. classified as SB, SA B ̲ $ \rm SA\underline{B} $, SAB, or S A ̲ B $ \rm S\underline{A}B $), 76 are grand-design, 157 are multi-armed, and 158 host flocculent spirals.

3. Global pitch angles of spirals

In order to parameterise the winding of the spirals, so that it can be studied as a function of galaxy parameters, we need to determine the global pitch angle of the galaxies (ϕ hereafter). We will assume that the spiral segments are all trailing, and will therefore treat ϕ always as a positive quantity. Also, we do not care whether the apparent spiral opening angles in sky are clockwise or counter-clockwise (the winding direction of arms in spiral galaxies is distributed at random in the local Universe, see e.g. Hayes et al. 2017).

3.1. Mean pitch angle of individual logarithmic segments

We first calculate the mean and median of the absolute value of the pitch angle measurements of logarithmic segments (|ϕ|mean and |ϕ|median, respectively). Uncertainties in the determination of |ϕ| are estimated by computing the standard deviation of the mean of the measurements of different segments. We compare |ϕ|mean and |ϕ|median in the top panel of Fig. 2: the mean absolute difference is 3.49° ±0.17°.

thumbnail Fig. 2.

Comparison between mean and (a) median pitch angle, (b) mean of the innermost logarithmic segments, and (c) mean weighted by the arc length of the arms. The y = x straight line is shown in red, and the dotted green line shows the linear fit to the cloud of points. Error bars indicate the standard deviation of the mean calculated from the segments fitted in each galaxy.

We also calculate the mean of the absolute value of the two innermost spiral segments (or one for the 6 galaxies with a single measurement), |ϕ|inner, to trace the pitch angle in the central parts of galaxies and in the vicinity of the bar. The mean absolute difference between |ϕ|mean and |ϕ|inner is 4.70° ±0.24° (central panel of Fig. 2).

In addition, we propose a measurement of the global pitch angle of disc galaxies by weighting the average by the relative arc length of the logarithmic segments, |ϕ|weighted. The equation for a logarithmic spiral in polar coordinates (r, θ) can be written as follows (e.g. Lin & Shu 1964; Elmegreen & Elmegreen 1987):

r ( θ ) = R · e θ tan ϕ , $$ \begin{aligned} r(\theta )=R \cdot e^{\theta \, \mathrm{tan} \phi }, \end{aligned} $$(1)

where θ = 0 at r = R ∈ ℝ and ϕ refers to the pitch angle. For a given galaxy, the arc length (s) of the logarithmic spiral segment i extending between (ri, θi) and ( r i , θ i ) $ (r_{i}^\prime,\theta_{i}^\prime) $ with pitch angle ϕi can be approximated as (e.g. Zwikker 1963):

s i = θ i θ i r 2 + ( d r / d θ ) 2 · d θ | r i r i | · 1 + tan | ϕ i | 2 tan | ϕ i | , $$ \begin{aligned} s_{i}=\int _{\theta _{i}}^{\theta _{i}{\prime }} \sqrt{r^2+(\mathrm{d}r/\mathrm{d}\theta )^{2}} \cdot \mathrm{d}\theta \approx |r_{i}{\prime }-r_{i}| \cdot \dfrac{\sqrt{1+\mathrm{tan} |\phi _{i}|^{2}}}{\mathrm{tan} |\phi _{i}|}, \end{aligned} $$(2)

and then

| ϕ | weighted = 1 S i = 1 n | ϕ i | · s i , $$ \begin{aligned} |\phi |_{\rm weighted}=\dfrac{1}{S} \sum _{i=1}^{n} |\phi _{i}| \cdot s_{i}, \end{aligned} $$(3)

where n is the total number of measured spiral arm segments, and S = i = 1 n s i $ S=\sum\nolimits_{i=1}^{n}s_{i} $. This is in practice very close to weighting with the relative radial extent of the spirals segments over which measurements were performed. The mean absolute difference between |ϕ|mean and |ϕ|weighted is 2.9° ±0.14° (see bottom panel of Fig. 2).

The global pitch angles of each of the galaxies in our sample2 are listed in Table A.1. Hereafter, we will use all proxies of |ϕ| for the analysis in this paper.

3.2. Pitch angles from 2D Fourier analysis

In order to test the sensitivity of the pitch angle determination to the employed methodology, we also fit the spiral structure using Fourier transform spectral analyses of the de-projected 3.6 μm images. We apply a similar method as in Puerari & Dottori (1992) and Puerari et al. (2000), who had in turn followed the approach by Kalnajs (1975) and Considere & Athanassoula (1982) (for a more sophisticated determination, see Puerari et al. 2014).

Images are decomposed into a flux-weighted superposition of m-armed logarithmic spirals of different pitch angles (see e.g. Puerari et al. 2000; Ma 2001; Davis et al. 2012), whose amplitudes are calculated as:

A ( p , m ) = 1 D π + π R inner R outer I ( u , θ ) e i ( m θ + p u ) d u d θ , $$ \begin{aligned} A(p,m)=\frac{1}{D} \int _{-\pi }^{+\pi } \int _{R_{\rm inner}}^{R_{\rm outer}} I(u,\theta ) e^{\mathrm{-i} (m \theta +p u)} \mathrm{d}u \mathrm{d}\theta , \end{aligned} $$(4)

where u = ln (r), I(u, θ) refers to the intensity distribution in polar coordinates, the variable p ∈ ℝ defines the pitch angle of the spiral as arctan (−m/p), D = R inner R outer π + π I ( r , θ ) d r d θ $ D=\int_{R_{\mathrm{inner}}}^{R_{\mathrm{outer}}}\int_{-\pi}^{+\pi} I(r,\theta)\, \mathrm{d}r \mathrm{d}\theta $ is the normalisation factor, and Rinner and Router are the limits in the radial range for the fit.

Choosing Rinner and Router is the greatest source of human error for the fit (Davis et al. 2012). Where available, we use the minimum and maximum radii where Herrera-Endoqui et al. (2015) fit the logarithmic segments. By doing so, the central bulge is excluded from the fit and the outer radius reaches the limit of the arms in the images. Whenever clear outer arm breaks are seen (e.g. M 51), we limit the fit to the inner segments. When the galaxies are barred, the masked region has a radius ≳rbar (at times, we mask regions well beyond the bar radius to avoid short spurious spiral segments that can contaminate the fit). Occasionally, when the spiral pattern presents multiple breaks, we limit the fit to the range [rbar, 2rbar]. By visual inspection of the images, we also choose the range of pitch angles for the fit (positive or negative depending on winding direction).

We use the limits p ∈ [ − 50, 50] and step-size (0.25) advocated by Puerari et al. (2000). Finally, for m = 2, we select the value of p with the highest amplitude, which defines the final pitch angle of the spiral arms, called as |ϕ|Fourier (for NGC 5054 we used m = 3 to fit the three arms). Examples of the resulting fits for M 51 and NGC 3992 are shown in Fig. 3.

thumbnail Fig. 3.

Upper panels: logarithmic fit to the spiral arms in the de-projected 3.6 μm image of NGC 5194 (left) (only central parts considered) (18–23 μ3.6μm (AB) magnitude scale) and NGC 3992 (right) (21–25 μ3.6μm (AB) magnitude scale), using 2D Fourier transform spectral analysis (see text). Lower panels: logarithm of the m = 2 Fourier amplitude as a function of the pitch angle for the two galaxies shown above. The vertical line denotes the maximum value, that determines the global pitch angle assigned to the galaxy.

We could reliably identify the maximum of the pitch angle p associated to the spiral pattern (|ϕ|Fourier) for a subsample of 32 of the galaxies in our master sample (∼10%), of which most (25) are barred (i.e. they are not SA) according to Buta et al. (2015). In Table 1 we list the measurements of the pitch angle using the 2D Fourier method. We also indicate the radial range where the fit was done. A rough estimate of the uncertainties in |ϕ|Fourier is calculated by reducing by an arbitrary 20 % the radial range for the fit: an offset of (Router − Rinner)/10 is added (subtracted) to Rinner (from Router).

Table 1.

Pitch angle of a subsample of S4G galaxies determined based on 2D Fourier analysis.

3.3. Comparison between methods and with values in the literature

Recent work by Yu et al. (2018) characterised the pitch angle of the spiral arms in the Carnegie-Irvine Galaxy Survey (Ho et al. 2011) using Fourier methods3: there is an overlap of 49 galaxies between their survey and our sample of 391 galaxies. In addition, there are seven S4G galaxies from our sample with measurements available from Davis et al. (2017)4. Of those 56 galaxies with pitch angle measurements in the literature, we could reliably apply the Fourier method to 11 (Sect. 3.2): in Fig. 4 we show a comparison of the obtained values, showing a good agreement.

thumbnail Fig. 4.

For a subsample of 11 S4G galaxies, comparison between the pitch angles obtained in this work and in the literature (see legend) applying Fourier analysis (see text). The y = x straight line is shown in red.

In general, there is also good agreement between the average pitch angle obtained in this work – based on the measurements of individual logarithmic segments – and the values of |ϕ|Fourier from the literature (Fig. 5): we obtain a mean absolute difference of 3.9° ±1.0° and 8.4° ±1.4° for |ϕ|mean and |ϕ|weighted, respectively. The scatter and outliers in the comparison are, most likely, due to the employed methodology to estimate the pitch angle.

thumbnail Fig. 5.

Comparison between mean and weighted mean pitch angles calculated in this work, and the pitch angle obtained by Yu et al. (2018) (black) and Davis et al. (2017) (blue diamonds) via Fourier analysis (see the text). The dashed red line corresponds to the straight line y = x, while the dotted green line shows the linear fit to the data cloud.

Differences in |ϕ| for the two methods applied in this paper (Fourier analysis and average of the pitch angles of individual logarithmic segments) over 3.6 μm images for a subsample of 32 galaxies are shown in Fig. 6. We compare |ϕ|Fourier with |ϕ|mean, |ϕ|inner, and |ϕ|weighted, obtaining mean absolute differences of 6.8° ±1.0°, 10.0° ±2.2°, and 5.9° ±0.7°, respectively. Discrepancies are associated to decoupled spiral modes that are not fitted with the Fourier method.

thumbnail Fig. 6.

For a subsample of S4G galaxies, comparison between the pitch angles determined with Fourier methods (from this work, see text) and (a) mean pitch angle, (b) inner mean pitch angles, and (c) mean pitch angle weighted by the arc length of the segments. We show measurements with very reliable (=1, filled circle) and acceptable (=2, empty diamonds) flags in Herrera-Endoqui et al. (2015). As in Fig. 5, the red dashed line corresponds to the straight line y = x, while the green dotted line shows the linear fit to the cloud of points.

Obtaining a single value that characterises the winding of the spirals is not straightforward, and it is sensitive to the utilised methodology and its drawbacks (see also discussion in e.g. Elmegreen et al. 1992). Given the uncertainties inherent to the various methodologies, the analysis of the possible dependence of pitch angle on bar strength (Sect. 7.1) is done based on visual measurements of |ϕ| as well as with the Fourier transform (using values from this work and from the literature).

4. Amplitude of non-axisymmetries

4.1. Bar strength and size

We use deprojected visual measurements of bar sizes (rbar) and ellipticities (ϵ) from Herrera-Endoqui et al. (2015). We use normalised Fourier amplitudes Am = Im/I0, where I0 indicates the m = 0 surface density component, calculated by Díaz-García et al. (2016a) by applying the NIR-QB code (Salo et al. 1999; Laurikainen & Salo 2002) to 3.6 μm S4G images. They also derived radial profiles of tangential forces normalised to the mean radial force field (Combes & Sanders 1981):

Q T ( r ) = max ( | F T ( r , ϕ ) | ) | F R ( r , ϕ ) | · $$ \begin{aligned} Q_{\rm T}\,(r)=\frac{\mathrm{max}\left( |F_{\rm T}\,(r,\phi )| \right)}{\langle |F_{\rm R}\,(r,\phi )|\rangle }\cdot \end{aligned} $$(5)

The maximum of A2 and QT at the bar region are used as proxies of the bar strength and named A 2 max $ A_2^{\rm max} $ and Qb, respectively (e.g. Buta & Block 2001; Laurikainen et al. 2002, 2004).

In addition, Díaz-García et al. (2016a) obtained a first-order model of the dark matter halo rotation curve and implemented a correction on QT for the contribution of the halo to the radial forces ( Q T halo-corr $ Q_{\mathrm{T}}^{\text{halo-corr}} $), following Buta et al. (2004).

Finally, Díaz-García et al. (2016a) eliminated the spiral contribution to the bar local force by setting to zero the Fourier density amplitudes beyond the bar radius for m >  0, as done in Salo et al. (2010). We use the maximum at the bar region ( Q b bar only $ Q_{\mathrm{b}}^{\mathrm{bar-only}} $) as an intrinsic measure on the bar-only force.

4.2. Spiral strength

Here we define a method to calculate the strength of spiral arms in the same regions [ r i , r i ] $ [r_{i},r_{i}^\prime] $ where the pitch angles of logarithmic segments were measured:

Q T Segment i = max r i r r i Q T ( r ) , $$ \begin{aligned}Q_{\rm T}^{\mathrm{Segment} \,i}=\max _{r_{i} \le r \le r_{i}{\prime }} Q_{\rm T}\,(r), \end{aligned} $$(6)

Q T,halo-corr Segment i = max r i r r i Q T halo-corr (r), $$ \begin{aligned}Q_{\text{T,} \text{halo-corr}}^{\mathrm{Segment} \,i}=\max _{r_{i} \le r \le r_{i}{\prime }} Q_{\rm T}^{\text{halo-corr}}\,(r),\end{aligned} $$(7)

A 2 Segmenti = max r i r r i A 2 (r). $$ \begin{aligned}A_{2}^{\mathrm{Segment} \,i}=\max _{r_{i} \le r \le r_{i}{\prime }} A_{2}\,(r). \end{aligned} $$(8)

For each galaxy, we then calculate the mean of the amplitude of the spiral segments:

Q T Spiral = Q T Segment i i[1,n] , $$ \begin{aligned}Q_{\rm T}^\mathrm{Spiral}=\left\langle Q_{\rm T}^{\mathrm{Segment} \,i}\right\rangle _{i\in [1,{n}]}\!, \end{aligned} $$(9)

Q T,halo-corr Spiral = Q T,halo-corr Segmenti i[1,n] , $$ \begin{aligned}Q_{\text{T,halo-corr}}^\mathrm{Spiral}=\left\langle Q_{\text{T,} \text{halo-corr}}^{\mathrm{Segment} \,i}\right\rangle _{i\in [1,{n}]}\!, \end{aligned} $$(10)

A 2 Spiral = A 2 Segmenti i[1,n] , $$ \begin{aligned}A_{2}^\mathrm{Spiral}=\left\langle A_{2}^{\mathrm{Segment} \,i}\right\rangle _{i\in [1,{n}]}\!, \end{aligned} $$(11)

where n refers to the number of logarithmic segments fitted by Herrera-Endoqui et al. (2015). In Appendix B we also make use of the mean amplitude of the two innermost and two outermost spiral segments.

There is a clear correlation between all the proxies of the spiral strength (Fig. 7)5. At the radii at which these segments were measured the contribution of the dark matter halo to the overall radial force field might be significant, and thus it may be important to perform the halo correction: this explains the scatter in the comparison of Q T Spiral, halo-corr $ Q_{\mathrm{T}}^{\text{Spiral, halo-corr}} $ versus QT (lower panel of Fig. 7).

thumbnail Fig. 7.

Comparison between the proxies of the amplitude of the spirals used in this work (Eqs. (9)–(11)). We show the mean value per galaxy, and measurements over individual segments. The Spearman’s rank correlation coefficients (significances) are 0.62 (4.14 × 10−43) and 0.82 (0.0) for the upper and lower panels, respectively. With a red dashed line we show the linear fit to the red data points.

5. Spiral arms across the Hubble sequence

5.1. Frequency of spiral types versus T-type

The distribution of the different arm classes as a function of the Hubble stage in our sample is shown in Fig. 8 (upper panel). Flocculent arms are mainly found among late-type spirals (peak for Scd, i.e. T = 6), multi-armed spirals mainly have intermediate T-types (peak for Sc, i.e. T = 5), while grand-design galaxies are typically early-types (peak for Sb, i.e. T = 3). All types of spirals are found for 3 ≲ T ≲ 6.

thumbnail Fig. 8.

Upper panel: histograms of the distributions of the galaxies in our sample (showing the barred galaxies separately) and of the different arm classes (grand-design, multi-armed, and flocculent) as a function of revised Hubble stage. Small offsets (≤0.1) have been added in the x-axis for the sake of avoiding line overlapping. Lower panel: bar fraction for the different spiral types. With vertical lines we represent binomial counting errors.

In the lower panel of Fig. 8 we show the bar fraction (fbar) for the different spiral types. fbar is known to have a bimodal distribution versus T in the S4G, with a minimum (40–45%) for T ≈ 4 (Sbc galaxies) (Buta et al. 2015; Díaz-García et al. 2016a). This trend is maintained when different classes of spirals are studied separately.

The global bar fraction in our sample of spiral galaxies is 69 ± 2.3%6. It is slightly lower for multi-armed spirals (59 ± 3.9%) and gets higher for grand-design (71 ± 5.2%) and flocculent galaxies (77 ± 3.3%). Among the 133 galaxies with T >  5, which are known to be faint (M* ≲ 1010 M) (see e.g. Fig. 1 in Laine et al. 2016), the fraction of barred galaxies is very high for all kinds of spirals, namely 100% for grand-design, 96.0 ± 3.9% for multi-armed, and 90.7 ± 2.9% for flocculent spirals.

5.2. Pitch angles versus T-type

For grand-design and multi-armed spirals, the global pitch angle increases with increasing T (Fig. 9), as expected, but a large scatter is found. The pitch angle is independent of T for flocculent spirals (Fig. 10), slightly flattening out the global distribution (see also Table 2). We do not find statistical and significant differences in the distribution of the pitch angle for barred and non-barred galaxies, regardless of spiral type or method used to estimate the average |ϕ|.

thumbnail Fig. 9.

Mean pitch angle (left panel), mean of the two innermost spiral segments (middle panel), and weighted mean pitch angle (right panel) as a function of the integer value of the revised numerical Hubble stage, for all the grand-design and multi-armed spirals in our sample. The running mean and standard deviation of the mean are shown for the barred (red) and non-barred (blue) galaxies separately, whenever more than one data point lies within the bin (i.e. when T ≤ 8 for barred galaxies and T ≤ 5 for non-barred ones). We have added small random offsets (≲0.3) to the T values in the x-axis (integers) to avoid point overlapping.

thumbnail Fig. 10.

As in Fig. 9, but including flocculent spirals (green). The running mean and standard deviation of the mean are shown for the galaxies that are either barred (red) or non-barred (blue, covering only the T-type range 1 ≤ T ≤ 6 where they appear).

Table 2.

Average pitch angles (in degrees) as a function of Hubble stage.

In Fig. 11 we show the internal scatter of |ϕ| versus T for all types of spirals in our sample. We note that σ|ϕ| is calculated, for each galaxy, from the pitch angle of the segments, whenever the host has more than one reliable measurement of |ϕ| (98.5% of the galaxies in our sample). The distribution is fairly flat. The mean σ|ϕ| is 9.5° ±0.3°, regardless of the presence of a bar, with a mild decrease in the end of the Hubble sequence.

thumbnail Fig. 11.

As in Fig. 10, but for the standard deviation of |ϕ| in each of the galaxies. The lines indicate the running median per T-type for barred and non-barred galaxies, and the vertical lines correspond to the standard deviation of the mean. The horizontal orange dotted line traces the mean σ|ϕ| (≈9.5°). Mean residual differences with respect to this value for each T bin (typically ≲2°) are shown in the inner panel on the upper left corner.

6. Pitch angle as a function of radius

Pitch angles in individual spiral galaxies are not necessarily constant as a function of galactocentric radius (e.g. Kennicutt 1981) (see also Fig. 1), and spiral arms are often asymmetric (e.g. Elmegreen et al. 2011). In Fig. 12 we study the radial variation of the pitch angle for all the galaxies in our sample, showing all the individual measurements on logarithmic spiral segments (see also Tables 3 and 4). Pitch angles that were measured within the radial interval [ r i , r i ] $ [r_{i},r_{i}^\prime] $ are represented as a function of the boundaries ri and r i $ r_{i}^\prime $, normalised by the disc size (traced by R25.5 and hR).

thumbnail Fig. 12.

Upper panel: radial variation of the pitch angle based on the individual measurements of logarithmic segments by Herrera-Endoqui et al. (2015). In the x-axis we plot the boundaries of the radial interval where the pitch was measured, normalised by the 25.5 mag arcsec−2 isophotal radius. The lines correspond to the running mean (bin size of 0.1), and the error bars are the standards deviation of the mean (the values appear listed in Table 3). Barred and non-barred galaxies are plotted separately. Lower panel: same as above but using the disc scale-length for the normalization, and showing only the running mean (bin size of 1, see Table 4). The colours indicate flocculent, multi-armed, and grand-design galaxies, as indicated in the legend.

Table 3.

Average pitch angles as a function of radius.

Table 4.

As in Table 3, but defining the radial ranges with respect to the disc scale-length (hR).

On average, the distribution of pitch angles barely change with radius (Fig. 12), regardless of spiral class (lower panel). The statistical trend is the same for barred and non-barred galaxies (upper panel).

7. Spiral properties as a function of bar strength

7.1. Pitch angle versus bar strength

In Fig. 13 we probe a possible dependence of winding of the arms on the bar-induced gravitational torque, which is expected from the manifold theory (see Sect. 9.3.1). We compare |ϕ|weighted and |ϕ|inner to the tangential-to-radial forcing evaluated at the bar radius: QT (rbar) (upper and central panels). We also use measurements of |ϕ|Fourier (lower panel), from this work and from the literature, to make sure that the results are not biased by the employed methodology. Overlaid in the upper panels of Fig. 13 is the expected correlation (roughly traced from their Fig. 5) from the numerical models by Athanassoula et al. (2009b), but we note that they evaluated tangential-to-radial forces exactly at the L1 point (typically an unstable saddle point, Binney & Tremaine 1987). L1 and L2 are hard to determine accurately in our images, since the kinematic data for our sample of galaxies is scarce and their dynamical modelling is beyond the scope of this paper.

thumbnail Fig. 13.

Weighted mean pitch angle (upper panel) and mean of the innermost spiral segments (middle panel) versus tangential-to-radial forces evaluated at the bar radius. The blue diamonds and line outline the trend in the simulations by Athanassoula et al. (2009b) for QT (L1). With red triangles we show the few galaxies in our sample hosting outer rings of type R1 or R1R2. On the lower panel we show the same plot but using pitch angles obtained via 2D Fourier analysis from this work and from the literature, using the measurements from Yu et al. (2018) (black points), and from Davis et al. (2017, 2019) (green diamonds).

The average values of |ϕ| for a given QT (rbar)-bin obtained in this work are consistent with those in the simulations by Athanassoula et al. (2009b). However, we hardly find any dependence (statistical tests7 confirm a very weak correlation: ρ ≈ 0.18 − 0.2, p <  0.01 for |ϕ|weighted, |ϕ|mean, and |ϕ|inner). For a better comparison with the aforementioned models, we analyse separately the eight galaxies in our sample hosting R1 or R1R2 outer rings. Their morphology is reproduced with manifold calculations (see Fig. 2 in Athanassoula et al. 2009b) and their L1 point location is most likely close to rbar. In spite of this, we do not find a clear trend between pitch angle and bar torque8. |ϕ|Fourier is not found to correlate with bar strength either. We confirmed the lack of a strong correlation between bar strength and pitch angles when the maximum bar torque – Qb – or the halo-corrected forces evaluated at the bar end – Q T halo-corr ( r bar ) $ Q_{\mathrm{T}}^{\text{halo-corr}}\,(r_{\mathrm{bar}}) $ – are used instead of QT (rbar): we obtain correlation coefficients (significances) of 0.18 (0.006) and 0.17 (0.01), respectively, when using |ϕ|mean (and roughly the same values for |ϕ|weighted and |ϕ|inner).

7.2. Spiral strength versus bar strength

In order to shed more light on the coupling between bars and spirals, here we assess, with unprecedented sample size, the possible dependence between their strengths, which has been extensively debated in the literature to this day (see Sect. 9.3.2). The amplitude of spirals in the S4G is correlated with bar strength (Fig. 14). This is clear for all proxies of the strength of spirals ( A 2 Spiral $ A_{2}^{\mathrm{Spiral}} $, Q T Spiral $ Q_{\mathrm{T}}^{\mathrm{Spiral}} $, Q T,halo-corr Spiral $ Q_{\text{T,halo-corr}}^{\mathrm{Spiral}} $), shown in the different panels (see Spearman’s correlation coefficient and significance above the panels). In Fig. 15 we show that the correlation is clear when the bar strength is estimated from Q b bar only $ Q_{\mathrm{b}}^{\mathrm{bar-only}} $ (spiral contribution to bar forcing suppressed).

thumbnail Fig. 14.

Bar strength (from Díaz-García et al. 2016a) as a function of spiral strength, measured from m = 2 Fourier amplitudes (Eq. (11)) (left), from tangential-to-radial forces (Eq. (9)) (middle), also including the correction for the halo dilution (Eq. (10)) (right). Different types of spirals are plotted with different symbols, as indicated in the legend. The red dashed line shows the linear fit to the cloud of points (in green and black we show the fit for flocculent and grand-design+multi-armed, respectively).

thumbnail Fig. 15.

Bar-only force versus spiral strength (Eq. (9)). The colour palette and symbols are the same as in Fig. 14.

The correlation between bars and spirals holds true when different types of spirals are studied separately. Interestingly, when using QT and A2, the trend for flocculent spirals is fairly similar as for grand-design and multi-armed. We note that the effect of the correction for the halo-dilution on the force profiles is slightly larger for T ≳ 5 (Díaz-García et al. 2016a). Thus, the difference between Q T,halo-corr Spiral $ Q_{\text{T, halo-corr}}^{\mathrm{Spiral}} $ and Q T Spiral $ Q_{\mathrm{T}}^{\mathrm{Spiral}} $ is more pronounced in flocculent spirals (typically hosted by late-type galaxies, see previous Sect. 5.1) than in other types of spirals.

Likewise, we find that the strength of the spirals increases with increasing bar size (Fig. 16), regardless of whether rbar is measured in physical units or relative to the disc size. In addition, the gravitational torques associated with the spirals (with and without halo-correction) increase with increasing pitch angle (Fig. 17, see also Fig. E.1): this is somehow expected, as more tightly-wrapped spirals tend to be more axisymmetric than loosely-wound arms.

thumbnail Fig. 16.

Bar size, in physical units (left) and normalised to the disc size (central and right panels), as a function of the spiral strength, measured from the m = 2 Fourier amplitude and from the tangential-to-radial forces. The colour palette and symbols are the same as in Fig. 14.

thumbnail Fig. 17.

Mean pitch angle as a function of spiral strength, measured from tangential-to-radial forces (Eq. (6)). The colour palette and symbols are the same as in Fig. 14. The red dashed line represents the linear fit to the data cloud, while the green and black lines show the fit for flocculent and grand-design+multi-armed galaxies, respectively.

8. Pitch angle and central mass concentration

The winding of the spiral arms has been claimed to depend on the prominence of the bulge (e.g. Kennicutt 1981), the steepness of the rotation curve (e.g. Seigar et al. 2005), and on the mass of the central supermassive black holes (MBH) (Davis et al. 2017). The small scatter in the latter scaling relation would make the pitch angle suitable for the prediction of MBH even in bulge-less galaxies. Here we assess these scaling relations in the S4G survey.

8.1. Pitch angle versus MBH

We estimate MBH from the central velocity dispersion (σ*) from HyperLEDA (available for 117 galaxies in our sample) and the calibration from Gültekin et al. (2009):

log 10 ( M BH ) = 8.12 + 4.24 · log 10 ( σ / 200 ) . $$ \begin{aligned} \mathrm{log_{10}}\,(M_{\rm BH})=8.12+4.24\cdot \mathrm{log_{10}}\,(\sigma _{*}/200). \end{aligned} $$(12)

In order to have an estimate of the uncertainties in MBH, for a subsample of S4G galaxies with disc inclinations lower than 65°, we checked that the inferred black hole masses are consistent with other measurements in the literature (e.g. compilations by Cisternas et al. 2013; Graham & Scott 2013; Sahu et al. 2019; Davis et al. 2019), using direct measurements and also the MBH − σ* relation (Fig. 18), but we find scatter and outliers (e.g. NGC 5055) in the comparison that emphasise the uncertainties in our analysis.

thumbnail Fig. 18.

For a sample of S4G disc galaxies with inclinations lower than 65°, we show a comparison of direct measurements of black hole masses from the compilations by Graham & Scott (2013), Sahu et al. (2019), Davis et al. (2019) – and also from Cisternas et al. (2013) (obtained via MBH − σ* relation and using direct measurements) – with the ones estimated using LEDA’s central velocity dispersions. In red, we show the straight line y = x.

We do not find a strong dependence between pitch angle and MBH (Fig. 19): correlation coefficients (significances) of −0.24 (0.008), −0.10 (0.27), and −0.03 (0.78) for |ϕ|mean, |ϕ|inner, and |ϕ|weighted, respectively. Only a weak correlation is detected with |ϕ|mean, that becomes slightly stronger (ρ = −0.34 and p = 0.0008) when only grand-design and multi-armed galaxies are considered (for the flocculent spirals alone, no correlation is seen whatsoever). Since the correlation is weak, the scatter in MBH for a given |ϕ| (about two decades) is much too large for this relation to yield any useful predictions concerning black hole mass. This contradicts previous reports in the literature (e.g. Davis et al. 2017): this could be explained by their use of smaller samples, but we note the lack of clear correlation can be due to the uncertainty in |ϕ| or to the use of indirect estimates of MBH.

thumbnail Fig. 19.

Mass of supermassive black holes (left) (estimated from central stellar velocity dispersions) and inner slope of the stellar component of the rotation curve (right) (tracer of central stellar mass concentration) as a function of mean pitch angle. The red dashed straight line (left panel) corresponds to the fit from the tight correlation presented by Davis et al. (2017): their values are consistent with ours, but we do not find such a strong dependence of the pitch angle on MBH. Different colours correspond to different spiral types, as indicated in the legend.

8.2. Pitch angle versus dRv* (0)

Díaz-García et al. (2016a) calculated the inner gradient of the stellar component of the rotation curve (dRv* (0) hereafter). They performed a polynomial fit to the inner part of the disc+bulge component of the rotation curve (V3.6 μm) and took the linear term as an estimate of the inner slope, following Lelli et al. (2013).

Here, we use dRv* (0) as a proxy of the central stellar mass concentration9 (e.g. Erroz-Ferrer et al. 2016; Díaz-García et al. 2016a; Salo & Laurikainen 2017). We show that the pitch angle of the galaxy is not dependent on dRv* (0) (right panel of Fig. 19).

9. Discussion

Here we discuss the properties of spiral arms based on the study of their pitch angles and amplitudes in the S4G survey, as presented in previous sections. Particular attention is given to the theoretical coupling between bars and spirals, which is assessed with observations. All types of spirals are included in our analysis. We take advantage of the wealth in late-type galaxies hosting flocculent arms in the S4G.

In this work we use near-IR (3.6 μm) imaging, which is known to be an excellent tracer of the old stellar populations. The pitch angle of a galaxy can vary depending on the used wavelength (e.g. Martínez-García et al. 2014). According to Pour-Imani et al. (2016), the typical difference between the pitch angle measured at 3.6 μm and 7 μm is ∼3.75°. Recent work by Yu & Ho (2018) reported a mild decrease of the pitch from red (including 3.6 μm) to blue passbands, which is in agreement with the expectations from the density wave theory.

The classification of spiral types is not strongly dependent on wavelength: for instance, Kendall et al. (2011) showed that there is a correspondence between galaxies being grand-design in infrared and optical passbands. Elmegreen & Elmegreen (1984) had claimed that old stellar waves could not be found in patchy spirals, proposing that flocculent spirals are mainly built of star-forming regions (Seiden & Gerola 1979; Elmegreen 1981). Recent work by Elmegreen et al. (2011) studied 46 galaxies in the S4G and showed that most of the galaxies that are flocculent in optical wavelengths can also be classified as flocculent at 3.6 μm. They also showed that, in the S4G, early-type spiral galaxies tend to be multi-armed and grand-design, having underlying stellar waves, whereas their late-type counterparts host flocculent spirals. We confirm this in Sect. 5 using a factor of ∼8 bigger sample (see also Bittner et al. 2017). These trends in the Hubble sequence are probably linked to the report by Hart et al. (2016) that galaxies hosting many arms – typically multi-armed and flocculent spirals – in the Galaxy Zoo 2 (Willett et al. 2013) have bluer colours than their two-armed counterparts (i.e. the grand-design).

It is well known that the majority of galaxies present spiral arms that cannot be fitted with a single pitch angle (e.g. Savchenko & Reshetnikov 2013; Herrera-Endoqui et al. 2015). In fact, in Sect. 5.1 we showed that the internal variation of |ϕ| in galaxies is approximately 10° on average, and the differences in certain cases can be larger, even by a factor of two. Thus, by estimating global winding angles one misses information on individual segments, in the present paper and in most of the work in the literature.

In Sect. 6 we also showed that the distribution of |ϕ| does not change with increasing radius, in a statistical sense, for all types of spirals. This implies that in discs with many logarithmic spiral segments |ϕ| can either get larger or smaller with increasing radius, and thus the arms do not systematically wind up or loosen in the galaxy outskirts. We are in an eclectic position regarding the radial variation of |ϕ| discussed in the literature: Savchenko & Reshetnikov (2013) and Davis et al. (2015) showed that the spiral pitch angles increase with decreasing radius, but an opposite trend was found by Davis & Hayes (2014).

In this work we propose the average pitch angle of galaxies as a global measure of the global winding angle, including a weighting for the relative arc length of the logarithmic segments (Sect. 3), so that we can study its dependence on global properties of the host galaxies. Overall, we find a good agreement between the values of the global |ϕ| determined with different methods (Fig. 6) and from the literature (Figs. 4 and 5). However, we also show scatter and outliers in the comparison, mainly caused by the differences in the employed methodology and its sensitivity to radial variations of |ϕ|, and to asymmetrical or decoupled spiral arms.

9.1. Observational constraints on spiral formation: pitch angle and disc properties

The pitch angle of the spirals in the S4G increases with increasing T for grand-design and multi-armed galaxies, on average, unlike in the case of flocculent galaxies (Sect. 5.2). The dispersion for each T-type is quite large: this is somewhat unexpected, since the winding angle of the arms is a criterion for morphological classification in the de Vaucouleurs Hubble–Sandage system (de Vaucouleurs 1959), applied by Buta et al. (2015) to our sample. The lack of strong correlation might be due to uncertainties in the T-type determination or the asymmetry of the spiral arms. It might also be a consequence of the larger weight given to other criteria to classify galaxies (e.g. smoothness or bulge prominence). The weak, or missing, correlation between pitch angle and T-type has been previously reported in the literature, and is reviewed in this section.

Early-type galaxies are expected to host more tightly wound spirals on average (e.g. Considere & Athanassoula 1988; Puerari & Dottori 1992; García Gómez & Athanassoula 1993), but counterexamples are known. Using K-band imaging for 45 face-on spirals, Seigar & James (1998a) found that their measurements of arm pitch angle did not correlate with Hubble type (obtained from optical images). They argued that this was due to the fact that the morphology of the old stellar population bears little resemblance to the optical morphology used to classify galaxies. In this work we find that the relation between the global pitch angle and the T-type has a lot of scatter and is not very tight, even when 3.6 μm imaging was used for both the morphological classification and the pitch angle measurements.

Early work by Kennicutt (1981) using 113 galaxies showed that the pitch angle correlates only in an average sense with T-type, and there is a large spread in the values. Kennicutt (1981) also reported that the shape of the spirals is mainly determined by kinematic parameters (absolute rotation velocity). A correlation of pitch angle and the fraction of light in the disc component has been predicted by theoreticians (e.g. Lin & Shu 1964; Roberts et al. 1975). Also, Block et al. (2004) discussed that open spirals imply high disc-to-halo mass ratios. Nevertheless, we checked and confirmed that the pitch angle does not depend on the total dynamical mass or halo-to-stellar mass ratio, derived from inclination-corrected 21 cm line widths (see Figs. 11. and C.1. in Díaz-García et al. 2019).

Kennicutt (1981) also showed that the spiral arms might not be either purely logarithmic (associated to a density wave origin) or hyperbolic (linked to a tidal origin). The environment cannot explain the formation of all spirals in all disc galaxies in the local universe, since spiral structure in the near-IR has been found in isolated galaxies (Kendall et al. 2011). However, Seigar & James (1998a) argued that tidal effects from near neighbours seem to have the effect of enhancing the amplitude of spirals. Analysis of Figs. 11. and C.1. in Díaz-García et al. (2019) reveals that the pitch angle in the S4G does not depend on the tidal interaction strength between galaxies or with the projected surface density of galaxies, taken from Laine et al. (2014) (see also Hart et al. 2016). In itself this does not rule out the influence of tidal perturbations on spiral formation: Salo & Laurikainen (2000b) showed in terms of N-body simulations of M 51 that, excluding the outer bridge and tail, the transient spiral density waves excited by interaction were very similar to those the model would generate in isolation.

Observations of disc galaxies at 2.1 μm by Block et al. (1994) lead the authors to conclude that large-scale structure of spirals is mainly intrinsic (i.e. not caused by tidal interactions or by the environment). They claimed that the stellar disc temperature, the gas content, and the disc mass are the main agents determining the spiral morphology. In addition, Davis et al. (2015) reported the existence of a fundamental plane of spiral structure in disc galaxies, determined by the pitch angle, the density of H I within the disc, and the mass of the stellar bulge.

A dependence between spiral arm pitch angles and the central mass concentration is expected from the density wave theory (Lin & Shu 1964, 1966). Kennicutt (1981) found that the spiral winding angle correlates only weakly with bulge-to-disc ratio. Recent work by Yu & Ho (2019) using 79 galaxies from the CALIFA survey (Sánchez et al. 2012) shows a dependence of pitch angles on the prominence of bulges, the galaxy concentration, and the total stellar mass, with a large scatter. Masters et al. (2019) performs the analysis of the spiral arms in the Galaxy Zoo project (Lintott et al. 2011; Willett et al. 2013) and reports a lack of strong correlation between bulge prominence and spiral arm tightness: galaxies with small bulges show a wide range of winding angles, while those with prominent bulges host more tightly wound arms. The authors speculate that the spiral arms are not static density waves, but rather re-formed structures that tend to wind-up in time. A link between |ϕ| and the relative bulge mass would also imply a correlation of the former with the disc mass fraction: such a correlation was not found by Font et al. (2019), based on the analysis of a sample of 79 disc galaxies with estimates of the corotation radius and the pattern speed of bars and spirals.

In this paper we do not find a dependence of the pitch angle on central stellar mass concentration (Sect. 8, Fig. 19). We checked and confirmed that the pitch angle does not depend on bulge-to-total mass ratio and bulge mass, estimated from 2D decompositions models by Salo et al. (2015). We also reassessed the dependence of the pitch on gas fraction and gas density (Appendix C), finding a weak correlation that holds for grand-design and multi-armed galaxies, but vanishes for flocculent spirals.

A correlation between spiral arm pitch angle and galactic shear rate has been found in N-body simulations by Grand et al. (2013) and recently confirmed by Fujii et al. (2018)10. Possible observational confirmation was provided by Seigar et al. (2006), who claimed that the dark matter halo profiles could be constrained using spiral arm morphology. Recent work by Yu et al. (2018) shows that the pitch angle correlates weakly, if at all, with the shear rate. Yu et al. (2018) actually reproduced the same plot from Seigar et al. (2006) using the same sample and the same shear rate estimates, but their own measurements of pitch angle, finding no clear correlation. We note that Seigar et al. (2006) found a connection between pitch and shear only when the latter was measured at a certain physical radius (10 kpc). The connection vanished when such radii were chosen as a function of the disc scale length. A correlation between spiral pitch angle and shear would imply an intimate link between spiral winding, central mass concentration, and dark matter haloes, which we cannot confirm in the S4G. In Appendix F we checked and confirmed the lack of correlation between shear and pitch angle for a subsample of 17 galaxies with available high-resolution rotation curves.

The whole picture for the formation of spiral arms remains unclear. An interpretation that disentangles the physical mechanisms that drive the formation and evolution of spiral arms when several observational parameters come into play is far from established. We conclude that, in the S4G, the pitch angle of spirals does not obey any definite scaling relation involving disc fundamental properties such as bulge prominence, relative dark matter content, gas fraction, stellar and dynamical mass, or environmental parameters.

9.2. Relation between the mass of the central supermassive black hole and the pitch angle of the spiral arms

Supermassive black holes are present in the centres of massive galaxies (see e.g. Lynden-Bell & Rees 1971; Kormendy & Richstone 1995; Schödel et al. 2002; Gillessen et al. 2009), with masses spanning several orders of magnitude. Recent pioneering work by the Event Horizon Telescope Collaboration (2019) produced the first image of the event horizon of a supermassive black hole in the centre of M 87, caught in the act of swallowing material.

Seigar et al. (2008) found a relationship between spiral arm pitch angle and central supermassive black hole mass (MBH). Further evidence was provided by Berrier et al. (2013) and, more recently, by Davis et al. (2017), who studied 44 spiral galaxies with measurements of MBH, and reported a surprisingly tight correlation between MBH and |ϕ|. If real, this could potentially explain previous findings in the literature of the dependence of the pitch and the bulge mass (e.g. Yoshizawa & Wakamatsu 1975), given the known scaling relation between central black holes and bulges (e.g. Kormendy & Richstone 1995). The tight correlation between |ϕ| and MBH would allow for the prediction of supermassive black hole mass even in bulge-less galaxies.

In Sect. 8 we used measurements in the literature of central stellar velocity dispersions to estimate the mass of the supermassive black hole hosted by the spiral galaxies in our sample. Although these are less direct estimates of MBH, we have a factor of 2 bigger sample. For a subsample of galaxies, we confirmed the fairly good agreement between direct and indirect estimates of black hole masses. Our values are consistent with the aforementioned trend, but we find no clear link between black hole mass and pitch angle (only a weak correlation is found for |ϕ|mean, confirmed by statistical tests). Recent work by Graham et al. (2019) reassessed the trend between ϕ and central velocity dispersion of the same sample of galaxies (see open black stars and circles in their Fig. 5), and the resulting scaling relation seems to be weaker.

Although the mean pitch angle correlates weakly with supermassive black hole mass (inferred from σ*) we conclude that no tight scaling relation between |ϕ| and MBH is found in the S4G. This questions previous studies in the literature that were based on smaller samples (e.g. Seigar et al. 2008; Davis et al. 2017). However, we note that part of the scatter in the relation might be due to the uncertainty in the determination of both parameters.

9.3. Interplay between bars and spirals

The effect of bar-induced secular evolution on mixing and reshaping discs is subtle, but its importance is becoming consolidated. For instance, by analysing average stellar density profiles, Díaz-García et al. (2016b) showed that barred galaxies have larger disc scale-lengths and fainter extrapolated surface brightness (∼15% difference) than their non-barred counterparts, in agreement with Sánchez-Janssen & Gadotti (2013). This proved the effect of bars pushing out the material that lies beyond corotation, in agreement with various simulation models (e.g. Athanassoula & Misiriotis 2002; Minchev et al. 2011; Athanassoula 2012). To further support this scenario, Seidel et al. (2015) showed that the flattening of the outer stellar velocity dispersion profiles increases with increasing bar strength. The next question is whether bars play a role in driving the formation of spirals.

Kendall et al. (2011) claimed that there is no evidence that bars trigger spirals (see also Seigar & James 1998a), since there is no connection between a galaxy being barred and having spiral structure. However, we confirm that grand-design spirals are more frequently barred than non-barred (in agreement with Elmegreen & Elmegreen 1982), and we find that this is also the case for multi-armed and flocculent galaxies (Sect. 5). Among late-type galaxies, most of the spiral galaxies are barred.

Nonetheless, non-axisymmetries cannot explain the formation of spirals in all the disc galaxies in the local universe, since spiral structure in the near-IR has been found in non-barred galaxies. In particular, nearly 30% of the spirals in our sample are non-barred, regardless of the type. In these cases, it seems unlikely that the spiral pattern could have been excited by bars that dissolved afterwards, since bars are believed to be long-lasting (e.g. Athanassoula 2003; Shen & Sellwood 2004; Villa-Vargas et al. 2010; Athanassoula et al. 2013; Sánchez-Blázquez et al. 2011; Gadotti et al. 2015) unless the disc is very gas rich (e.g. Bournaud & Combes 2002).

Masters et al. (2019) reports that, at a given bulge size, strongly barred galaxies have more loosely wound arms than non-barred galaxies: this is linked by the authors to the effect of bars slowing the winding speed of the arms. In this work we do not find any difference in the distribution of the pitch angles for barred and non-barred galaxies when 1 ≤ T ≤ 5 (Sect. 5): this questions the role of bars in the formation and evolution of spirals.

9.3.1. The hypothetical link between pitch angle and bar strength: expectations from the manifold theory

The manifold theory (e.g. Romero-Gómez et al. 2006, 2007; Patsis 2006; Voglis et al. 2006; Athanassoula et al. 2009a) makes specific predictions for the shape of spirals and their coupling with bars. The models by Athanassoula et al. (2009b) show a tight correlation between the pitch angle and the tangential-to-radial forces evaluated at the unstable L1 point (see their Fig. 5). Recent observational work by Font et al. (2019) does not confirm this causality. We find a very weak correlation between gravitational torques evaluated at the bar end – QT (rbar) – and the pitch angle (Sect. 7.1).

We also checked that no clear dependence of |ϕ| on bar strength is found when other proxies of the latter are used – such as the maximum bar torque (Qb), the intrinsic ellipticity, or the m = 2 Fourier amplitude – or when the dilution of bar gravitational torques by the halo contribution to the overall radial force field was modelled. For a more direct comparison, we additionally checked a possible dependence between bar strength and the average pitch angle of the innermost spiral segments, which are closest to the unstable Lagrangian points in barred galaxies, confirming the negative result. We also analysed this trend with other Fourier-based measurements of pitch angles, both taken from the literature or from this work, to make sure that the lack of a strong correlation cannot be attributed to the employed methodology to characterise the spiral winding in observations.

We conclude that the formation of the majority of the spiral arms in the S4G cannot be explained with the manifold theory. Likewise, Díaz-García et al. (2019) did not confirm the manifold origin of rings in the S4G from the predicted correlation between the ratio of outer-to-inner ring semi-major axes and the bar strength.

9.3.2. Coupling of bars and spirals analysed from their amplitudes

The observational test of the coupling between bars and spirals has been controversial for years. Based on the comparison of equivalent angles of arms and bars (measured as in Seigar & James 1998b), Seigar et al. (2003) did not find evidence for the role of bars driving spiral structure in a sample of 17 inclined galaxies. Based on the Fourier analysis of near-IR KS-band imaging, Block et al. (2004) and Buta et al. (2009) noted that stronger bars are associated with stronger arms. Buta et al. (2005) removed the spiral contribution to the bar torque (following Buta et al. 2003) in the OSUBSGS sample (Eskridge et al. 2002) and showed a correlation between the bar and spiral torques.

On the other hand, Durbala et al. (2009) did not find a correlation between the strengths of bars and spirals in their sample of isolated galaxies. Buta et al. (2009) found the relation to be rather weak using KS-band images from the Anglo-Australian Telescope (AAT) for 23 galaxies, and concluded that some arms might be genuine bar-driven (i.e. bars and spirals are formed from the same disc instability) but other factors might explain their formation as well.

Salo et al. (2010) reassessed the bar-spiral coupling in the above-referred OSUBSGS and AAT samples and showed a correlation between the local spiral density amplitude and the local bar-only forcing – the bar-induced tangential-to-radial forces after the effect of the spirals is suppressed – out to ≈1.6 bar radii. Within this radius, the authors interpret that spirals are a long-lasting continuation of the bar mode, whereas beyond the radius the spirals might be decoupled (e.g. Sellwood & Sparke 1988) or transient (e.g. Sellwood & Kahn 1991; Laurikainen & Salo 2002).

More recently, Díaz-García et al. (2016b) showed that A2 amplitudes are larger across the discs of barred galaxies than those of their non-barred counterparts (see their Fig. 9) when 0 ≤ T <  5. Likewise, the amplitude of non-axisymmetries beyond the bar radius were more pronounced for the strongly barred galaxies than for the weakly barred ones: this is consistent with bars exciting the spiral pattern. In addition, Elmegreen et al. (2011) found that longer bars with larger m = 2 density amplitudes (more evolved) appear together with stronger arms. Bittner et al. (2017) showed that galaxies with high bar contrasts tend to have higher arm contrasts as well. Lately, Sánchez-Menguiano et al. (2017) studied the oxygen abundance of spiral arms and interarm regions and found small differences for barred galaxies, not detected in their unbarred counterparts, suggesting that bars may have a subtle effect on the chemical distribution of these galaxies.

In Sect. 7.2 we showed a clear correlation between the strength of spiral arms and bars, as well as a moderate dependency between spiral strength and bar length. In addition, the mean amplitude of the arms was found to increase with increasing pitch angle (this is also confirmed in Appendix E, for all proxies of spiral strength, using all spiral segments instead of the mean). Altogether, the analysis of the amplitudes of spirals and bars supports a coupling between these stellar structures, even though the details of this connection are far from clear.

As discussed in Elmegreen et al. (2011, and references therein), grand-design spiral arms and the symmetric inner parts of multi-armed galaxies are thought to be the result of spiral density waves, that could have been driven by bars. On the other hand, properties of flocculent galaxies and the outer parts of multi-armed galaxies can in principle be explained with random local gravitational instabilities in the stellar and gaseous component (e.g. Goldreich & Lynden-Bell 1965b; Julian & Toomre 1966), swing amplified into spiral arms (e.g. Toomre 1981). Thus, it is meaningful that the correlation between the strength of bars and spirals is fairly similar for all types of spirals. Most likely, bars are not responsible for exciting flocculent spiral structure or the outer spiral modes of multi-armed galaxies galaxies. Thus, the correlation between the amplitudes of bars and spirals may not necessarily imply causation. In fact, in Appendix B we show that the correlation holds even when the mean amplitude of the two outermost spiral segments is considered (only slightly smaller correlation coefficients than for the innermost ones).

We conclude that the correlation between bar and spiral amplitudes might simply indicate that discs that are prone to the development of bars with large amplitudes are also reactive to the formation of prominent spirals, as discussed in Salo et al. (2010). Bars might play a role in creating spirals (Sanders & Huntley 1976), but mainly in strongly barred galaxies (e.g. Sellwood & Sparke 1988).

10. Summary and conclusions

In this paper we characterise the winding angle and amplitude of the spirals arms of 391 not-highly inclined galaxies (i <  65°) from the Spitzer Survey of Stellar Structure in Galaxies (S4G; Sheth et al. 2010). We study barred and non-barred galaxies separately. We analyse all types of spirals, that is, grand-design, multi-armed, and flocculent.

We determine the pitch angle of the arms using the human-supervised measurements of spiral segments from Herrera-Endoqui et al. (2015) (Fig. 1). When multiple spiral segments appear for individual galaxies (98.5% of the cases), the global pitch angle (|ϕ|) is obtained from the mean (|ϕ|mean), median (|ϕ|median), the mean of innermost logarithmic segments (|ϕ|inner), and the mean weighted by the arc length of the segments (|ϕ|weighted) (Figs. 2, 4, 5, and 6). Values of |ϕ| are listed in Table A.1. In addition, for a subsample of 32 galaxies we could reliably determine the pitch angle using 2D Fourier analysis (|ϕ|Fourier) (Fig. 3, Table 1).

Pitch angles are analysed as a function of galaxy properties – such as morphological type (T) or the estimated mass of the central supermassive black hole – and the strength of non-axisymmetries. The latter is estimated from gravitational torques (QT) and m = 2 Fourier amplitudes (A2) from Díaz-García et al. (2016a) evaluated at the bar region (bar strength), and also at the radial ranges where the pitch angle of logarithmic segments were fitted (spiral strength) (Fig. 7).

The main results of this paper are the following:

  • The distribution of grand-design, multi-armed, and flocculent spirals peaks at T = 3 (Sb), T = 5 (Sc), and T = 6 (Scd), respectively (Fig. 8). Roughly 2/3 of the spirals in our sample are barred (i.e. classified as SB, SA B ̲ $ \rm SA\underline{B} $, SAB, or S A ̲ B $ \rm S\underline{A}B $ by Buta et al. 2015), and for T >  5 the bar fraction is larger than 90%.

  • In galaxies with more than one measured spiral segment, |ϕ| has a dispersion of σ = 9.5 ± 0.3° on average, but the differences on individual galaxies can be ≳15 − 20° (Fig. 11).

  • The distribution of the pitch angles as a function of galactocentric radius is flat (Fig. 12). Hence, statistically speaking, the pitch angles can either increase or decrease with radius.

  • On average, the pitch angle increases with increasing Hubble stage for grand-design and multi-armed spirals (Fig. 9), as expected by the definition of T, but the scatter in this relation is large. Among flocculent spirals, |ϕ| is not correlated with T (Fig. 10).

  • We confirm that the strengths of bars and spirals are correlated (Fig. 14). This either supports the role of bars driving the formation of spirals (Sanders & Huntley 1976) or indicates that the discs that are prone to the formation of stronger bars are also more reactive to the formation of spirals of larger amplitudes (see e.g. Salo et al. 2010). The latter interpretation is favoured by the observed coupling of the amplitudes of bar and arms even in flocculent spirals, which are implausible to be excited by bars, or when only the outermost segments of multi-armed galaxies are analysed.

  • We confirm that long bars (in physical size or relative to the disc size) tend to be hosted by discs with prominent spirals (Fig. 16), in agreement with Elmegreen et al. (2011).

  • The pitch angle of the spiral arms is not controlled by the bar-induced perturbation strength (Fig. 13). This questions the manifold origin of spiral arms (e.g. Athanassoula et al. 2009b).

  • The distribution of pitch angles for barred and non-barred galaxies is roughly the same when 1 ≤ T ≤ 5 (Figs. 9, 10, and 12). This questions the role of bars driving spiral density waves.

  • The mean pitch angle is not strongly correlated with the mass of the supermassive black hole (Fig. 19), estimated from central stellar velocity dispersion (Fig. 18), and does not depend on the central concentration of stellar mass or shear (Fig. F.1). This questions previous results in the literature using smaller samples (e.g. Kennicutt 1981; Davis et al. 2015) and casts doubt on the validity of predicting MBH from |ϕ| (e.g. Davis et al. 2017).


1

Only six galaxies in our sample (1.5%) have a single measurement of the pitch angle, 74 galaxies (18.9%) have two measured logarithmic segments, 86 (22.0%) have three, 105 (26.9%) have four, 58 (14.8%) have five, and 62 galaxies (15.9%) have more than five (with two cases having a maximum of nine segments).

2

Preliminary analysis of these measurements (and the properties of barred and non-barred galaxies) using unsupervised machine learning techniques (self-organising maps) has been done by Díaz-García et al. (2019).

3

We note that the mean difference between the pitch angles obtained at 3.6 μm and from I-band images by Yu et al. (2018) is very small: ≈0.1° for the 2D method and ≈0.3° for the 1D method (Yu & Ho 2018).

4

These measurements were performed using using the software Spirality (Shields et al. 2015a,b), 2DFFT software (Davis et al. 2012, 2016), and compared to computer vision software (Davis & Hayes 2014).

5

Díaz-García et al. (2016a) reported a bimodality in the comparison between Qb and A 2 max $ A_{2}^{\mathrm{max}} $, which might be due to the sensitivity of Qb to the enhanced radial forces by galactic bulges (Block et al. 2001; Laurikainen & Salo 2002; Díaz-García et al. 2016a). As a matter of fact, no bimodality is present when A 2 Spiral $ A_{2}^{\mathrm{Spiral}} $ is plotted versus Q T Spiral $ Q_{\mathrm{T}}^{\mathrm{Spiral}} $, which are measured at a radii where the bulge dilution should be much smaller.

6

We calculate binomial counting errors: Δ f bar = f bar · ( 1 f bar ) / N gals $ \Delta f_{\mathrm{bar}}=\sqrt{f_{\mathrm{bar}}\cdot (1-f_{\mathrm{bar}})/N_{\mathrm{gals}}} $, where fbar refers to the bar fraction and Ngals to the total number of galaxies.

7

The Spearman’s rank correlation assesses the existence of a monotonic relation between two variables. ρ = ( − 1)+1 implies perfect (anti)correlation. Small p-values (< 0.01) indicate significant correlation.

8

The low number of rings of type R1, R2, or R1R2 is most likely due to the bias of the S4G towards late-type galaxies. We note that we do not sample galaxies hosting rings of type R2 exclusively.

9

Díaz-García (2016) showed that dRv* (0) tightly correlates with M* (and thus with the maximum circular velocity), T, and the bulge-to-total mass ratio (see their Fig. 4.6). This was also shown by Erroz-Ferrer et al. (2016) for a subsample of 29 S4G galaxies.

10

A similar dependence between local shear rate and pitch angle has also been demonstrated in collisional N-body simulation of Saturn ring’s self-gravity wake structures, which arise from superposition of swing amplified noise in a disc with Toomre parameter Q <  2 − 3 (Salo et al. 2018): such structures agree very well with the stellar dynamical predictions of Julian & Toomre (1966) and can be viewed as analogous of flocculent spiral segments.

11

Under the assumption that the halo within the optical disc contributes approximately a constant fraction of the total halo mass (∼4%), Díaz-García et al. (2016a,c) found that the trend of the Mhalo/M* versus M* relation agreed with the prediction of ΛCDM models, fitted at z ≈ 0 based on abundance matching and halo occupation distribution methods (e.g. Moster et al. 2010).

Acknowledgments

We thank the anonymous referee for comments that improved this paper. We acknowledge financial support from the European Union’s Horizon 2020 research and innovation programme under Marie Skłodowska-Curie grant agreement No 721463 to the SUNDIAL ITN network, and from the Spanish Ministry of Economy and Competitiveness (MINECO) under grant number AYA2016-76219-P. S.D. acknowledges the financial support from the visitor and mobility program of the Finnish Centre for Astronomy with ESO (FINCA), funded by the Academy of Finland grant number 306531. J.H.K. acknowledges support from the Fundación BBVA under its 2017 programme of assistance to scientific research groups, for the project “Using machine-learning techniques to drag galaxies from the noise in deep imaging”, and from the Leverhulme Trust through the award of a Visiting Professorship at LJMU. H.S. acknowledges financial support from the Academy of Finland (grant no: 297738). M.H.E. thanks CONACyT for supporting his postdoctoral position through the grant number 252531. We thank Sébastien Comerón for excellent comments on the manuscript. We thank the S4G team for their work on the different pipelines and the data release. We thank Si Yue-Yu for providing their measurements of pitch angles for a subsample of S4G galaxies. We thank Joachim Janz, Aaron Watkins, and Elena D’Onghia for useful discussions. Facilities: Spitzer (IRAC).

References

  1. Athanassoula, E. 2003, MNRAS, 341, 1179 [NASA ADS] [CrossRef] [Google Scholar]
  2. Athanassoula, E. 2012, MNRAS, 426, L46 [NASA ADS] [Google Scholar]
  3. Athanassoula, E., & Misiriotis, A. 2002, MNRAS, 330, 35 [NASA ADS] [CrossRef] [Google Scholar]
  4. Athanassoula, E., Romero-Gómez, M., & Masdemont, J. J. 2009a, MNRAS, 394, 67 [NASA ADS] [CrossRef] [Google Scholar]
  5. Athanassoula, E., Romero-Gómez, M., Bosma, A., & Masdemont, J. J. 2009b, MNRAS, 400, 1706 [NASA ADS] [CrossRef] [Google Scholar]
  6. Athanassoula, E., Machado, R. E. G., & Rodionov, S. A. 2013, MNRAS, 429, 1949 [NASA ADS] [CrossRef] [Google Scholar]
  7. Berrier, J. C., Davis, B. L., Kennefick, D., et al. 2013, ApJ, 769, 132 [NASA ADS] [CrossRef] [Google Scholar]
  8. Binney, J., & Tremaine, S. 1987, Galactic Dynamics (Princeton, NJ: Princeton University Press) [Google Scholar]
  9. Bittner, A., Gadotti, D. A., Elmegreen, B. G., et al. 2017, MNRAS, 471, 1070 [NASA ADS] [CrossRef] [Google Scholar]
  10. Block, D. L., Bertin, G., Stockton, A., et al. 1994, A&A, 288, 365 [NASA ADS] [Google Scholar]
  11. Block, D. L., Puerari, I., Knapen, J. H., et al. 2001, A&A, 375, 761 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Block, D. L., Buta, R., Knapen, J. H., et al. 2004, AJ, 128, 183 [NASA ADS] [CrossRef] [Google Scholar]
  13. Bournaud, F., & Combes, F. 2002, A&A, 392, 83 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Buta, R., & Block, D. L. 2001, ApJ, 550, 243 [NASA ADS] [CrossRef] [Google Scholar]
  15. Buta, R., Block, D. L., & Knapen, J. H. 2003, AJ, 126, 1148 [NASA ADS] [CrossRef] [Google Scholar]
  16. Buta, R., Laurikainen, E., & Salo, H. 2004, AJ, 127, 279 [NASA ADS] [CrossRef] [Google Scholar]
  17. Buta, R., Vasylyev, S., Salo, H., & Laurikainen, E. 2005, AJ, 130, 506 [NASA ADS] [CrossRef] [Google Scholar]
  18. Buta, R. J., Knapen, J. H., Elmegreen, B. G., et al. 2009, AJ, 137, 4487 [NASA ADS] [CrossRef] [Google Scholar]
  19. Buta, R. J., Sheth, K., Athanassoula, E., et al. 2015, ApJS, 217, 32 [NASA ADS] [CrossRef] [Google Scholar]
  20. Cisternas, M., Gadotti, D. A., Knapen, J. H., et al. 2013, ApJ, 776, 50 [NASA ADS] [CrossRef] [Google Scholar]
  21. Combes, F., & Sanders, R. H. 1981, A&A, 96, 164 [NASA ADS] [Google Scholar]
  22. Considere, S., & Athanassoula, E. 1982, A&A, 111, 28 [NASA ADS] [Google Scholar]
  23. Considere, S., & Athanassoula, E. 1988, A&AS, 76, 365 [NASA ADS] [Google Scholar]
  24. Davis, D. R., & Hayes, W. B. 2014, ApJ, 790, 87 [NASA ADS] [CrossRef] [Google Scholar]
  25. Davis, B. L., Berrier, J. C., Shields, D. W., et al. 2012, ApJS, 199, 33 [NASA ADS] [CrossRef] [Google Scholar]
  26. Davis, B. L., Kennefick, D., Kennefick, J., et al. 2015, ApJ, 802, L13 [NASA ADS] [CrossRef] [Google Scholar]
  27. Davis, B. L., Berrier, J. C., & Shields, D. W. 2016, Astrophysics Source Code Library [record ascl:1608.015] [Google Scholar]
  28. Davis, B. L., Graham, A. W., & Seigar, M. S. 2017, MNRAS, 471, 2187 [NASA ADS] [CrossRef] [Google Scholar]
  29. Davis, B. L., Graham, A. W., & Cameron, E. 2019, ApJ, 873, 85 [NASA ADS] [CrossRef] [Google Scholar]
  30. de Vaucouleurs, G. 1959, Handbuch Phys., 53, 275 [Google Scholar]
  31. Díaz-García, S. 2016, Ph. D. Thesis, University of Oulu [Google Scholar]
  32. Díaz-García, S., Salo, H., Laurikainen, E., & Herrera-Endoqui, M. 2016a, A&A, 587, A160 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Díaz-García, S., Salo, H., & Laurikainen, E. 2016b, A&A, 596, A84 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Díaz-García, S., Salo, H., Laurikainen, E., & Leaman, R. 2016c, ArXiv e-prints [arXiv:1611.01844] [Google Scholar]
  35. Díaz-García, S., Díaz-Suárez, S., Knapen, J. H., & Salo, H. 2019, A&A, 625, A146 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. D’Onghia, E. 2015, ApJ, 808, L8 [NASA ADS] [CrossRef] [Google Scholar]
  37. D’Onghia, E., Vogelsberger, M., & Hernquist, L. 2013, ApJ, 766, 34 [NASA ADS] [CrossRef] [Google Scholar]
  38. Durbala, A., Buta, R., Sulentic, J. W., & Verdes-Montenegro, L. 2009, MNRAS, 397, 1756 [NASA ADS] [CrossRef] [Google Scholar]
  39. Elmegreen, D. M. 1981, ApJS, 47, 229 [NASA ADS] [CrossRef] [Google Scholar]
  40. Elmegreen, B. G. 1990, Ann. New York Acad. Sci., 596, 40 [NASA ADS] [CrossRef] [Google Scholar]
  41. Elmegreen, D. M., & Elmegreen, B. G. 1982, MNRAS, 201, 1021 [NASA ADS] [CrossRef] [Google Scholar]
  42. Elmegreen, D. M., & Elmegreen, B. G. 1984, ApJS, 54, 127 [NASA ADS] [CrossRef] [Google Scholar]
  43. Elmegreen, B. G., & Elmegreen, D. M. 1985, ApJ, 288, 438 [NASA ADS] [CrossRef] [Google Scholar]
  44. Elmegreen, D. M., & Elmegreen, B. G. 1987, ApJ, 314, 3 [NASA ADS] [CrossRef] [Google Scholar]
  45. Elmegreen, B. G., & Elmegreen, D. M. 1989, ApJ, 342, 677 [NASA ADS] [CrossRef] [Google Scholar]
  46. Elmegreen, B. G., Elmegreen, D. M., & Montenegro, L. 1992, ApJS, 79, 37 [NASA ADS] [CrossRef] [Google Scholar]
  47. Elmegreen, D. M., Elmegreen, B. G., Yau, A., et al. 2011, ApJ, 737, 32 [NASA ADS] [CrossRef] [Google Scholar]
  48. Erroz-Ferrer, S., Knapen, J. H., Leaman, R., et al. 2016, MNRAS, 458, 1199 [NASA ADS] [CrossRef] [Google Scholar]
  49. Erwin, P. 2018, MNRAS, 474, 5372 [NASA ADS] [CrossRef] [Google Scholar]
  50. Eskridge, P. B., Frogel, J. A., Pogge, R. W., et al. 2002, ApJS, 143, 73 [NASA ADS] [CrossRef] [Google Scholar]
  51. Event Horizon Telescope Collaboration (Akiyama, K., et al.) 2019, ApJ, 875, L1 [NASA ADS] [CrossRef] [Google Scholar]
  52. Fazio, G. G., Hora, J. L., Allen, L. E., et al. 2004, ApJS, 154, 10 [NASA ADS] [CrossRef] [Google Scholar]
  53. Font, J., Beckman, J. E., James, P. A., & Patsis, P. A. 2019, MNRAS, 482, 5362 [NASA ADS] [CrossRef] [Google Scholar]
  54. Fujii, M. S., Bédorf, J., Baba, J., & Portegies Zwart, S. 2018, MNRAS, 477, 1451 [NASA ADS] [CrossRef] [Google Scholar]
  55. Gadotti, D. A., Seidel, M. K., Sánchez-Blázquez, P., et al. 2015, A&A, 584, A90 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  56. García Gómez, C., & Athanassoula, E. 1993, A&AS, 100, 431 [NASA ADS] [Google Scholar]
  57. Gillessen, S., Eisenhauer, F., Trippe, S., et al. 2009, ApJ, 692, 1075 [NASA ADS] [CrossRef] [Google Scholar]
  58. Giovanelli, R., & Haynes, M. P. 1988, Extragalactic neutral hydrogen., ed. K. I. Kellermann, & G. L. Verschuur, 522–562 [Google Scholar]
  59. Goldreich, P., & Lynden-Bell, D. 1965a, MNRAS, 130, 97 [NASA ADS] [CrossRef] [Google Scholar]
  60. Goldreich, P., & Lynden-Bell, D. 1965b, MNRAS, 130, 125 [NASA ADS] [CrossRef] [Google Scholar]
  61. Graham, A. W., & Scott, N. 2013, ApJ, 764, 151 [NASA ADS] [CrossRef] [Google Scholar]
  62. Graham, A. W., Soria, R., & Davis, B. L. 2019, MNRAS, 484, 814 [NASA ADS] [Google Scholar]
  63. Grand, R. J. J., Kawata, D., & Cropper, M. 2013, A&A, 553, A77 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  64. Gültekin, K., Richstone, D. O., Gebhardt, K., et al. 2009, ApJ, 698, 198 [NASA ADS] [CrossRef] [Google Scholar]
  65. Hart, R. E., Bamford, S. P., Willett, K. W., et al. 2016, MNRAS, 461, 3663 [NASA ADS] [CrossRef] [Google Scholar]
  66. Hayes, W. B., Davis, D., & Silva, P. 2017, MNRAS, 466, 3928 [NASA ADS] [CrossRef] [Google Scholar]
  67. Herrera-Endoqui, M., Díaz-García, S., Laurikainen, E., & Salo, H. 2015, A&A, 582, A86 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  68. Ho, L. C., Li, Z.-Y., Barth, A. J., Seigar, M. S., & Peng, C. Y. 2011, ApJS, 197, 21 [NASA ADS] [CrossRef] [Google Scholar]
  69. Julian, W. H., & Toomre, A. 1966, ApJ, 146, 810 [NASA ADS] [CrossRef] [Google Scholar]
  70. Kalnajs, A. J. 1975, in La Dynamique des galaxies spirales, ed. L. Weliachew, 241, 103 [NASA ADS] [Google Scholar]
  71. Kendall, S., Kennicutt, R. C., & Clarke, C. 2011, MNRAS, 414, 538 [NASA ADS] [CrossRef] [Google Scholar]
  72. Kennicutt, Jr., R. C. 1981, AJ, 86, 1847 [NASA ADS] [CrossRef] [Google Scholar]
  73. Kennicutt, Jr., R., & Hodge, P. 1982, ApJ, 253, 101 [NASA ADS] [CrossRef] [Google Scholar]
  74. Knapen, J. H., & Beckman, J. E. 1996, MNRAS, 283, 251 [NASA ADS] [CrossRef] [Google Scholar]
  75. Kormendy, J., & Kennicutt, Jr., R. C. 2004, ARA&A, 42, 603 [NASA ADS] [CrossRef] [Google Scholar]
  76. Kormendy, J., & Norman, C. A. 1979, ApJ, 233, 539 [NASA ADS] [CrossRef] [Google Scholar]
  77. Kormendy, J., & Richstone, D. 1995, ARA&A, 33, 581 [NASA ADS] [CrossRef] [Google Scholar]
  78. Laine, J., Laurikainen, E., Salo, H., et al. 2014, MNRAS, 441, 1992 [NASA ADS] [CrossRef] [Google Scholar]
  79. Laine, J., Laurikainen, E., & Salo, H. 2016, A&A, 596, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  80. Laurikainen, E., & Salo, H. 2002, MNRAS, 337, 1118 [NASA ADS] [CrossRef] [Google Scholar]
  81. Laurikainen, E., Salo, H., & Rautiainen, P. 2002, MNRAS, 331, 880 [NASA ADS] [CrossRef] [Google Scholar]
  82. Laurikainen, E., Salo, H., & Buta, R. 2004, ApJ, 607, 103 [NASA ADS] [CrossRef] [Google Scholar]
  83. Lelli, F., Fraternali, F., & Verheijen, M. 2013, MNRAS, 433, L30 [NASA ADS] [CrossRef] [Google Scholar]
  84. Lelli, F., McGaugh, S. S., & Schombert, J. M. 2016, AJ, 152, 157 [Google Scholar]
  85. Lin, C. C., & Shu, F. H. 1964, ApJ, 140, 646 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  86. Lin, C. C., & Shu, F. H. 1966, Proc. Nat. Acad. Sci., 55, 229 [Google Scholar]
  87. Lintott, C., Schawinski, K., Bamford, S., et al. 2011, MNRAS, 410, 166 [NASA ADS] [CrossRef] [Google Scholar]
  88. Lynden-Bell, D., & Rees, M. J. 1971, MNRAS, 152, 461 [NASA ADS] [CrossRef] [Google Scholar]
  89. Ma, J. 2001, Chin. J. Astron. Astrophys., 1, 395 [NASA ADS] [CrossRef] [Google Scholar]
  90. Martínez-García, E. E., Puerari, I., Rosales-Ortega, F. F., et al. 2014, ApJ, 793, L19 [NASA ADS] [CrossRef] [Google Scholar]
  91. Masters, K. L., Lintott, C. J., Hart, R. E., et al. 2019, MNRAS, 487, 1808 [NASA ADS] [CrossRef] [Google Scholar]
  92. Minchev, I., Famaey, B., Combes, F., et al. 2011, A&A, 527, A147 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  93. Moster, B. P., Somerville, R. S., Maulbetsch, C., et al. 2010, ApJ, 710, 903 [NASA ADS] [CrossRef] [Google Scholar]
  94. Muñoz-Mateos, J. C., Sheth, K., Regan, M., et al. 2015, ApJS, 219, 3 [NASA ADS] [CrossRef] [Google Scholar]
  95. Patsis, P. A. 2006, MNRAS, 369, L56 [NASA ADS] [CrossRef] [Google Scholar]
  96. Pour-Imani, H., Kennefick, D., Kennefick, J., et al. 2016, ApJ, 827, L2 [NASA ADS] [CrossRef] [Google Scholar]
  97. Puerari, I., Block, D. L., Elmegreen, B. G., Frogel, J. A., & Eskridge, P. B. 2000, A&A, 359, 932 [NASA ADS] [Google Scholar]
  98. Puerari, I., & Dottori, H. A. 1992, A&AS, 93, 469 [NASA ADS] [Google Scholar]
  99. Puerari, I., Elmegreen, B. G., & Block, D. L. 2014, AJ, 148, 133 [NASA ADS] [CrossRef] [Google Scholar]
  100. Roberts, Jr., W. W., Roberts, M. S., & Shu, F. H. 1975, ApJ, 196, 381 [NASA ADS] [CrossRef] [Google Scholar]
  101. Romero-Gómez, M., Masdemont, J. J., Athanassoula, E., & García-Gómez, C. 2006, A&A, 453, 39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  102. Romero-Gómez, M., Athanassoula, E., Masdemont, J. J., & García-Gómez, C. 2007, A&A, 472, 63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  103. Rosse, T. E. O. 1850, Philos. Trans. R. Soc. London, Ser. I, 140, 499 [NASA ADS] [CrossRef] [Google Scholar]
  104. Sahu, N., Graham, A. W., & Davis, B. L. 2019, ApJ, 876, 155 [NASA ADS] [CrossRef] [Google Scholar]
  105. Salo, H., & Laurikainen, E. 2000a, MNRAS, 319, 377 [NASA ADS] [CrossRef] [Google Scholar]
  106. Salo, H., & Laurikainen, E. 2000b, MNRAS, 319, 393 [NASA ADS] [Google Scholar]
  107. Salo, H., & Laurikainen, E. 2017, ApJ, 835, 252 [NASA ADS] [CrossRef] [Google Scholar]
  108. Salo, H., Rautiainen, P., Buta, R., et al. 1999, AJ, 117, 792 [NASA ADS] [CrossRef] [Google Scholar]
  109. Salo, H., Laurikainen, E., Buta, R., & Knapen, J. H. 2010, ApJ, 715, L56 [NASA ADS] [CrossRef] [Google Scholar]
  110. Salo, H., Laurikainen, E., Laine, J., et al. 2015, ApJS, 219, 4 [NASA ADS] [CrossRef] [Google Scholar]
  111. Salo, H., Ohtsuki, K., & Lewis, M. C. 2018, in Planetary Ring Systems, eds. M. S. Tiscareno, & C. D. Murray, 434 [CrossRef] [Google Scholar]
  112. Sánchez, S. F., Kennicutt, R. C., & Gil de Paz, A. 2012, A&A, 538, A8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  113. Sánchez-Blázquez, P., Ocvirk, P., Gibson, B. K., Pérez, I., & Peletier, R. F. 2011, MNRAS, 415, 709 [NASA ADS] [CrossRef] [Google Scholar]
  114. Sánchez-Janssen, R., & Gadotti, D. A. 2013, MNRAS, 432, L56 [NASA ADS] [CrossRef] [Google Scholar]
  115. Sánchez-Menguiano, L., Sánchez, S. F., Pérez, I., et al. 2017, A&A, 603, A113 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  116. Sanders, R. H., & Huntley, J. M. 1976, ApJ, 209, 53 [NASA ADS] [CrossRef] [Google Scholar]
  117. Savchenko, S. S., & Reshetnikov, V. P. 2013, MNRAS, 436, 1074 [NASA ADS] [CrossRef] [Google Scholar]
  118. Schödel, R., Ott, T., Genzel, R., et al. 2002, Nature, 419, 694 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  119. Schweizer, F. 1976, ApJS, 31, 313 [NASA ADS] [CrossRef] [Google Scholar]
  120. Seidel, M. K., Falcón-Barroso, J., Martínez-Valpuesta, I., et al. 2015, MNRAS, 451, 936 [NASA ADS] [CrossRef] [Google Scholar]
  121. Seiden, P. E., & Gerola, H. 1979, ApJ, 233, 56 [NASA ADS] [CrossRef] [Google Scholar]
  122. Seigar, M. S., & James, P. A. 1998a, MNRAS, 299, 685 [NASA ADS] [CrossRef] [Google Scholar]
  123. Seigar, M. S., & James, P. A. 1998b, MNRAS, 299, 672 [NASA ADS] [CrossRef] [Google Scholar]
  124. Seigar, M. S., Chorney, N. E., & James, P. A. 2003, MNRAS, 342, 1 [NASA ADS] [CrossRef] [Google Scholar]
  125. Seigar, M. S., Block, D. L., Puerari, I., Chorney, N. E., & James, P. A. 2005, MNRAS, 359, 1065 [NASA ADS] [CrossRef] [Google Scholar]
  126. Seigar, M. S., Bullock, J. S., Barth, A. J., & Ho, L. C. 2006, ApJ, 645, 1012 [NASA ADS] [CrossRef] [Google Scholar]
  127. Seigar, M. S., Kennefick, D., Kennefick, J., & Lacy, C. H. S. 2008, ApJ, 678, L93 [NASA ADS] [CrossRef] [Google Scholar]
  128. Seigar, M. S., Davis, B. L., Berrier, J., & Kennefick, D. 2014, ApJ, 795, 90 [NASA ADS] [CrossRef] [Google Scholar]
  129. Sellwood, J. A. 2011, MNRAS, 410, 1637 [NASA ADS] [Google Scholar]
  130. Sellwood, J. A., & Kahn, F. D. 1991, MNRAS, 250, 278 [NASA ADS] [Google Scholar]
  131. Sellwood, J. A., & Sparke, L. S. 1988, MNRAS, 231, 25P [NASA ADS] [CrossRef] [Google Scholar]
  132. Shen, J., & Sellwood, J. A. 2004, ApJ, 604, 614 [NASA ADS] [CrossRef] [Google Scholar]
  133. Sheth, K., Regan, M., Hinz, J. L., et al. 2010, PASP, 122, 1397 [NASA ADS] [CrossRef] [Google Scholar]
  134. Shields, D. W., Boe, B., & Pfountz, C. 2015a, ArXiv e-prints [arXiv:1511.06365] [Google Scholar]
  135. Shields, D. W., Boe, B., & Pfountz, C. 2015b, Astrophysics Source Code Library [record ascl:1512.015] [Google Scholar]
  136. Toomre, A. 1981, in Structure and Evolution of Normal Galaxies, eds. S. M. Fall, & D. Lynden-Bell, 111 [Google Scholar]
  137. Villa-Vargas, J., Shlosman, I., & Heller, C. 2010, ApJ, 719, 1470 [NASA ADS] [CrossRef] [Google Scholar]
  138. Voglis, N., Tsoutsis, P., & Efthymiopoulos, C. 2006, MNRAS, 373, 280 [NASA ADS] [CrossRef] [Google Scholar]
  139. Watkins, A. E., Mihos, J. C., & Harding, P. 2015, ApJ, 800, L3 [NASA ADS] [CrossRef] [Google Scholar]
  140. Werner, M. W., Roellig, T. L., Low, F. J., et al. 2004, ApJS, 154, 1 [NASA ADS] [CrossRef] [Google Scholar]
  141. Willett, K. W., Lintott, C. J., Bamford, S. P., et al. 2013, MNRAS, 435, 2835 [NASA ADS] [CrossRef] [Google Scholar]
  142. Yoshizawa, M., & Wakamatsu, K. 1975, A&A, 44, 363 [NASA ADS] [Google Scholar]
  143. Yu, S.-Y., & Ho, L. C. 2018, ApJ, 869, 29 [NASA ADS] [CrossRef] [Google Scholar]
  144. Yu, S.-Y., & Ho, L. C. 2019, ApJ, 871, 194 [NASA ADS] [CrossRef] [Google Scholar]
  145. Yu, S.-Y., Ho, L. C., Barth, A. J., & Li, Z.-Y. 2018, ApJ, 862, 13 [NASA ADS] [CrossRef] [Google Scholar]
  146. Zwikker, C. 1963, The Advanced Geometry of Plane Curves and Their Applications (Mineola, NY: Dover Phoenix Editions) [Google Scholar]

Appendix A: Tabulated pitch angles and spiral amplitudes

The following data are listed in Table A.1 (the complete version for the 391 spiral galaxies in our sample can be found at the CDS).

Table A.1.

Average pitch angles and spiral amplitudes.

Appendix B: Bar strength versus spiral strength: extended analysis

In Sect. 7.2 we showed a correlation between the bar strength and the mean strength of the spirals. Here we test this dependence when only the mean of the amplitude the two innermost and the two outermost spiral segments are considered (see Fig. B.1). We only use galaxies with more than four measured spiral segments. The strength of the spiral segments is measured from the m = 2 Fourier density amplitudes (left), tangential-to-radial forces (centre), and also torques corrected for the halo-dilution (right).

thumbnail Fig. B.1.

As in Fig. 14, but taking the mean of the maximum amplitude of the two innermost (upper panels) and outermost (lower panels) spiral logarithmic segments.

In Fig. B.2 we confirm the same trends when using bar-only forces (see Sect. 4.1). Interestingly, the correlation between the amplitude of the bar and that of the spiral modes holds in these extreme radial ranges.

thumbnail Fig. B.2.

As in Fig. B.1, but eliminating the spiral contribution to the bar gravitational torque.

Appendix C: Pitch angle as a function of the relative mass of cold gas

Here we assess a possible dependence between |ϕ| and the relative content of gas (e.g. Davis et al. 2015). Atomic gas masses are estimated as (e.g. Giovanelli & Haynes 1988; Erwin 2018):

M HI = 2.356 × 10 5 · D 2 · 10 0.4 · ( 17.4 m 21 c ) , $$ \begin{aligned} M_{\rm HI}=2.356 \times 10^5 \cdot D^2 \cdot 10^{0.4 \cdot (17.4-m21c)}, \end{aligned} $$(C.1)

where m21c is the corrected 21-cm line flux in magnitude from HyperLEDA, D is the distance to the galaxy (in megaparsecs) adopted by Muñoz-Mateos et al. (2015). H I surface gas densities are obtained as follows:

Σ HI = M HI / ( π · R 25.5 2 ) , $$ \begin{aligned} \Sigma _{\rm HI}=M_{\rm HI}/(\pi \cdot R_{25.5}^2), \end{aligned} $$(C.2)

where R25.5 indicates the isophotal radii at the surface brightness 25.5 mag arcsec−2 in the 3.6 μm images. The pitch angles of the spirals in the S4G are weakly correlated with gas fraction or density (Fig. C.1). A slightly tighter trend is seen when we take only multi-armed and grand-design spiral galaxies.

thumbnail Fig. C.1.

Pitch angle versus H I gas fraction (upper panels) and H I surface density (lower panels) for S4G galaxies with inclination < 65°, including barred and non-barred and all types of spirals (multi-armed, flocculent, and grand-design). The Spearman’s correlation coefficient and significance are indicated on top of the plots. Different spiral types are shown with different colours and symbols, as indicated in the legend. The red line shows a linear fit to the cloud of points.

Appendix D: Pitch angle as a function of galaxy mass

We test a possible dependence of the spiral pitch angle on the total stellar mass of the galaxies, estimated by Muñoz-Mateos et al. (2015), and on dynamical mass (Mdyn) and (dark) halo-to-stellar mass ratio (Mhalo/M*). Díaz-García et al. (2016a) derived a first-order estimate of Mhalo/M* within the optical radius (Ropt ∼ 3.2 ⋅ hR), which we use here11:

M halo / M ( < R opt ) 1.34 · ( ( V HI max ) 2 V 3.6 μ m 2 ( R opt ) 1 ) , $$ \begin{aligned} M_{\rm halo}/M_{*}\,( < R_{\rm opt})\approx 1.34 \cdot \bigg (\frac{(V_{\rm HI}^\mathrm{max})^{2}}{V_{3.6 \mu \mathrm m}^{2}\,(R_{\rm opt})}-1\bigg ), \end{aligned} $$(D.1)

where V HI max $ V_{\mathrm{HI}}^{\mathrm{max}} $ refers to the inclination-corrected H I velocity amplitude from HyperLEDA and V3.6 μm corresponds to the stellar component of the circular velocity (inferred from 3.6 μm imaging). We also estimate the dynamical mass within Ropt (assuming a spherical mass distribution):

M dyn ( < R opt ) ( V HI max ) 2 · R opt G , $$ \begin{aligned} M_{\rm dyn}\,( < R_{\rm opt})\approx \frac{(V_{\rm HI}^\mathrm{max})^{2} \cdot R_{\rm opt}}{G}, \end{aligned} $$(D.2)

where G is the Newtonian constant of gravitation.

No dependence of the global spiral pitch angle on total stellar mass, dynamical mass, or halo-to-stellar mass ratio is found in the S4G (Fig. D.1). This questions the effect of the global mass distribution in determining the winding angle of the arms (e.g. Kennicutt 1981; Kennicutt & Hodge 1982; Seigar & James 1998a; Seigar et al. 2006; Grand et al. 2013).

thumbnail Fig. D.1.

Galaxy mass versus pitch angle. We plot the total stellar mass (upper row), dynamical mass (central row), and halo-to-stellar mass ratio (lower row) versus the mean pitch angle (left column), mean weighted by the arc length of the arms (central column), and mean of the innermost logarithmic segments (right column). Different colours and symbols represent different types of spirals, as indicated in the legend of the upper left panel (as in Fig. 14).

Appendix E: Strength versus pitch angle of spiral segments

In Sect. 7.2 we showed a correlation between the mean pitch and the mean strength of the spirals (see Fig. 17). Here we confirm this dependence when all the individual logarithmic segments (often several per galaxy) are used (Fig. E.1).

thumbnail Fig. E.1.

As in Fig. 17 but using the strength of all the spiral segments, instead of the mean, and showing also the m = 2 Fourier density amplitude of the spirals (right panel).

Appendix F: Pitch angle versus shear

We calculate the shear (Γ) from H I and Hα rotation curves – V(r) – from the SPARC compilation (Lelli et al. 2016) as follows (e.g. Fujii et al. 2018):

Γ = d ln Ω / d ln r , $$ \begin{aligned} \Gamma =-\mathrm{d} \mathrm{ln}\,\Omega /\mathrm{d} \mathrm{ln}\,r, \end{aligned} $$(F.1)

where Ω = V/r is the angular velocity. The shear is obtained for the 17 SPARC galaxies with reliable measurements of the pitch angle in Herrera-Endoqui et al. (2015). Specifically, Γ was calculated in the same radial ranges where |ϕ| was fitted, and thus several measurements per galaxy are provided. We do not find any correlation between Γ and |ϕ|.

thumbnail Fig. F.1.

Pitch angle versus shear for a sample of 17 galaxies. Every point (72) corresponds to a logarithmic segment. In green we show the flocculent spirals, while black squares and filled circles correspond to grand-design and multi-armed galaxies, respectively (see legend). The dashed (solid) lines show the expected value of Γ for a flat (Keplerian) rotation curve.

All Tables

Table 1.

Pitch angle of a subsample of S4G galaxies determined based on 2D Fourier analysis.

Table 2.

Average pitch angles (in degrees) as a function of Hubble stage.

Table 3.

Average pitch angles as a function of radius.

Table 4.

As in Table 3, but defining the radial ranges with respect to the disc scale-length (hR).

Table A.1.

Average pitch angles and spiral amplitudes.

All Figures

thumbnail Fig. 1.

3.6 μm image of NGC 5194 (upper panels) and NGC 3992 (lower panels) in sky plane (left), with fitted spiral arm segments, and logarithmic polar plot of the same galaxies de-projected to the disc plane (right). Images and measurements are taken from Herrera-Endoqui et al. (2015).

In the text
thumbnail Fig. 2.

Comparison between mean and (a) median pitch angle, (b) mean of the innermost logarithmic segments, and (c) mean weighted by the arc length of the arms. The y = x straight line is shown in red, and the dotted green line shows the linear fit to the cloud of points. Error bars indicate the standard deviation of the mean calculated from the segments fitted in each galaxy.

In the text
thumbnail Fig. 3.

Upper panels: logarithmic fit to the spiral arms in the de-projected 3.6 μm image of NGC 5194 (left) (only central parts considered) (18–23 μ3.6μm (AB) magnitude scale) and NGC 3992 (right) (21–25 μ3.6μm (AB) magnitude scale), using 2D Fourier transform spectral analysis (see text). Lower panels: logarithm of the m = 2 Fourier amplitude as a function of the pitch angle for the two galaxies shown above. The vertical line denotes the maximum value, that determines the global pitch angle assigned to the galaxy.

In the text
thumbnail Fig. 4.

For a subsample of 11 S4G galaxies, comparison between the pitch angles obtained in this work and in the literature (see legend) applying Fourier analysis (see text). The y = x straight line is shown in red.

In the text
thumbnail Fig. 5.

Comparison between mean and weighted mean pitch angles calculated in this work, and the pitch angle obtained by Yu et al. (2018) (black) and Davis et al. (2017) (blue diamonds) via Fourier analysis (see the text). The dashed red line corresponds to the straight line y = x, while the dotted green line shows the linear fit to the data cloud.

In the text
thumbnail Fig. 6.

For a subsample of S4G galaxies, comparison between the pitch angles determined with Fourier methods (from this work, see text) and (a) mean pitch angle, (b) inner mean pitch angles, and (c) mean pitch angle weighted by the arc length of the segments. We show measurements with very reliable (=1, filled circle) and acceptable (=2, empty diamonds) flags in Herrera-Endoqui et al. (2015). As in Fig. 5, the red dashed line corresponds to the straight line y = x, while the green dotted line shows the linear fit to the cloud of points.

In the text
thumbnail Fig. 7.

Comparison between the proxies of the amplitude of the spirals used in this work (Eqs. (9)–(11)). We show the mean value per galaxy, and measurements over individual segments. The Spearman’s rank correlation coefficients (significances) are 0.62 (4.14 × 10−43) and 0.82 (0.0) for the upper and lower panels, respectively. With a red dashed line we show the linear fit to the red data points.

In the text
thumbnail Fig. 8.

Upper panel: histograms of the distributions of the galaxies in our sample (showing the barred galaxies separately) and of the different arm classes (grand-design, multi-armed, and flocculent) as a function of revised Hubble stage. Small offsets (≤0.1) have been added in the x-axis for the sake of avoiding line overlapping. Lower panel: bar fraction for the different spiral types. With vertical lines we represent binomial counting errors.

In the text
thumbnail Fig. 9.

Mean pitch angle (left panel), mean of the two innermost spiral segments (middle panel), and weighted mean pitch angle (right panel) as a function of the integer value of the revised numerical Hubble stage, for all the grand-design and multi-armed spirals in our sample. The running mean and standard deviation of the mean are shown for the barred (red) and non-barred (blue) galaxies separately, whenever more than one data point lies within the bin (i.e. when T ≤ 8 for barred galaxies and T ≤ 5 for non-barred ones). We have added small random offsets (≲0.3) to the T values in the x-axis (integers) to avoid point overlapping.

In the text
thumbnail Fig. 10.

As in Fig. 9, but including flocculent spirals (green). The running mean and standard deviation of the mean are shown for the galaxies that are either barred (red) or non-barred (blue, covering only the T-type range 1 ≤ T ≤ 6 where they appear).

In the text
thumbnail Fig. 11.

As in Fig. 10, but for the standard deviation of |ϕ| in each of the galaxies. The lines indicate the running median per T-type for barred and non-barred galaxies, and the vertical lines correspond to the standard deviation of the mean. The horizontal orange dotted line traces the mean σ|ϕ| (≈9.5°). Mean residual differences with respect to this value for each T bin (typically ≲2°) are shown in the inner panel on the upper left corner.

In the text
thumbnail Fig. 12.

Upper panel: radial variation of the pitch angle based on the individual measurements of logarithmic segments by Herrera-Endoqui et al. (2015). In the x-axis we plot the boundaries of the radial interval where the pitch was measured, normalised by the 25.5 mag arcsec−2 isophotal radius. The lines correspond to the running mean (bin size of 0.1), and the error bars are the standards deviation of the mean (the values appear listed in Table 3). Barred and non-barred galaxies are plotted separately. Lower panel: same as above but using the disc scale-length for the normalization, and showing only the running mean (bin size of 1, see Table 4). The colours indicate flocculent, multi-armed, and grand-design galaxies, as indicated in the legend.

In the text
thumbnail Fig. 13.

Weighted mean pitch angle (upper panel) and mean of the innermost spiral segments (middle panel) versus tangential-to-radial forces evaluated at the bar radius. The blue diamonds and line outline the trend in the simulations by Athanassoula et al. (2009b) for QT (L1). With red triangles we show the few galaxies in our sample hosting outer rings of type R1 or R1R2. On the lower panel we show the same plot but using pitch angles obtained via 2D Fourier analysis from this work and from the literature, using the measurements from Yu et al. (2018) (black points), and from Davis et al. (2017, 2019) (green diamonds).

In the text
thumbnail Fig. 14.

Bar strength (from Díaz-García et al. 2016a) as a function of spiral strength, measured from m = 2 Fourier amplitudes (Eq. (11)) (left), from tangential-to-radial forces (Eq. (9)) (middle), also including the correction for the halo dilution (Eq. (10)) (right). Different types of spirals are plotted with different symbols, as indicated in the legend. The red dashed line shows the linear fit to the cloud of points (in green and black we show the fit for flocculent and grand-design+multi-armed, respectively).

In the text
thumbnail Fig. 15.

Bar-only force versus spiral strength (Eq. (9)). The colour palette and symbols are the same as in Fig. 14.

In the text
thumbnail Fig. 16.

Bar size, in physical units (left) and normalised to the disc size (central and right panels), as a function of the spiral strength, measured from the m = 2 Fourier amplitude and from the tangential-to-radial forces. The colour palette and symbols are the same as in Fig. 14.

In the text
thumbnail Fig. 17.

Mean pitch angle as a function of spiral strength, measured from tangential-to-radial forces (Eq. (6)). The colour palette and symbols are the same as in Fig. 14. The red dashed line represents the linear fit to the data cloud, while the green and black lines show the fit for flocculent and grand-design+multi-armed galaxies, respectively.

In the text
thumbnail Fig. 18.

For a sample of S4G disc galaxies with inclinations lower than 65°, we show a comparison of direct measurements of black hole masses from the compilations by Graham & Scott (2013), Sahu et al. (2019), Davis et al. (2019) – and also from Cisternas et al. (2013) (obtained via MBH − σ* relation and using direct measurements) – with the ones estimated using LEDA’s central velocity dispersions. In red, we show the straight line y = x.

In the text
thumbnail Fig. 19.

Mass of supermassive black holes (left) (estimated from central stellar velocity dispersions) and inner slope of the stellar component of the rotation curve (right) (tracer of central stellar mass concentration) as a function of mean pitch angle. The red dashed straight line (left panel) corresponds to the fit from the tight correlation presented by Davis et al. (2017): their values are consistent with ours, but we do not find such a strong dependence of the pitch angle on MBH. Different colours correspond to different spiral types, as indicated in the legend.

In the text
thumbnail Fig. B.1.

As in Fig. 14, but taking the mean of the maximum amplitude of the two innermost (upper panels) and outermost (lower panels) spiral logarithmic segments.

In the text
thumbnail Fig. B.2.

As in Fig. B.1, but eliminating the spiral contribution to the bar gravitational torque.

In the text
thumbnail Fig. C.1.

Pitch angle versus H I gas fraction (upper panels) and H I surface density (lower panels) for S4G galaxies with inclination < 65°, including barred and non-barred and all types of spirals (multi-armed, flocculent, and grand-design). The Spearman’s correlation coefficient and significance are indicated on top of the plots. Different spiral types are shown with different colours and symbols, as indicated in the legend. The red line shows a linear fit to the cloud of points.

In the text
thumbnail Fig. D.1.

Galaxy mass versus pitch angle. We plot the total stellar mass (upper row), dynamical mass (central row), and halo-to-stellar mass ratio (lower row) versus the mean pitch angle (left column), mean weighted by the arc length of the arms (central column), and mean of the innermost logarithmic segments (right column). Different colours and symbols represent different types of spirals, as indicated in the legend of the upper left panel (as in Fig. 14).

In the text
thumbnail Fig. E.1.

As in Fig. 17 but using the strength of all the spiral segments, instead of the mean, and showing also the m = 2 Fourier density amplitude of the spirals (right panel).

In the text
thumbnail Fig. F.1.

Pitch angle versus shear for a sample of 17 galaxies. Every point (72) corresponds to a logarithmic segment. In green we show the flocculent spirals, while black squares and filled circles correspond to grand-design and multi-armed galaxies, respectively (see legend). The dashed (solid) lines show the expected value of Γ for a flat (Keplerian) rotation curve.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.