Free Access
Issue
A&A
Volume 643, November 2020
Article Number A14
Number of page(s) 24
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/202038448
Published online 27 October 2020

© ESO 2020

1. Introduction

A large number of observational and theoretical studies have been providing mounting evidence of the important physical processes driven by bars in massive disc galaxies. For example, in the area of the disc within the bar radius, stars are continuously trapped by the bar and gas is funnelled to the central region, where new stellar structures such as nuclear discs, inner bars, nuclear rings and nuclear spiral arms are thus built (see Sanders & Tubbs 1980; Simkin et al. 1980; Kormendy 1982; Prendergast 1983; Louis & Gerhard 1988; Shlosman et al. 1989; Athanassoula 1992; Knapen et al. 1995, 2002; Piner et al. 1995; Sakamoto et al. 1999; Gadotti & dos Anjos 2001; Sheth et al. 2005; Allard et al. 2006; Wozniak 2007; Coelho & Gadotti 2011; Ellison et al. 2011; de Lorenzo-Cáceres et al. 2012, 2013; Athanassoula et al. 2013; Cole et al. 2014; Fragkoudi et al. 2016; Kim et al. 2016; Kruk et al. 2018; Seo et al. 2019; Donohoe-Keyes et al. 2019). There is increasing evidence that these processes start to play a major role at redshifts z ∼ 1−2 in the most massive disc galaxies, and at z ∼ 0 these processes take part in the evolution of about two-thirds of disc galaxies (see e.g. Eskridge et al. 2000; Menéndez-Delmestre et al. 2007; Sheth et al. 2008, 2012; Kraljic et al. 2012; Simmons et al. 2014; Melvin et al. 2014; Gadotti et al. 2015; Pérez et al. 2017).

A major process in the evolution of barred galaxies is the buckling (or bulging) of the inner part of the bar, which grows vertically from the disc plane, creating the so-called box/peanut/X-shaped bulges, such as the one hosted by our own Milky Way (e.g. Combes & Sanders 1981; de Souza & Dos Anjos 1987; Shaw 1987; Combes et al. 1990; Binney et al. 1991; Kuijken & Merrifield 1995; Merrifield & Kuijken 1999; Bureau & Freeman 1999; Lütticke et al. 2000; Chung & Bureau 2004; Bureau & Athanassoula 2005; Bureau et al. 2006; Martinez-Valpuesta et al. 2006; Erwin & Debattista 2017; Fragkoudi et al. 2017, 2018; Kruk et al. 2019). Naturally, these box/peanuts are easier to identify in galaxies close to or at an edge-on projection, but a number of studies made an effort to produce diagnostics able to uncover box/peanuts in inclined galaxies (e.g. Erwin & Debattista 2013, from photometry), or even in face-on galaxies (e.g. Debattista et al. 2005; Méndez-Abreu et al. 2008; Łokas 2019; Fragkoudi et al. 2020, from stellar kinematics).

More recently, it has been shown that the inner parts of bars not only expand vertically from the plane of the disc, but also radially, parallel to the disc plane, away from the bar major axis (see Laurikainen et al. 2005, 2007, 2011; Athanassoula et al. 2015; Laurikainen & Salo 2016). When bars are seen face-on, this additional structure appears as a less eccentric central component, wider and shorter than the main body of the bar (see Fig. 8 in Gonzalez & Gadotti 2016). This component was named “barlens”, but it is important to note that the terms barlens and “box/peanut” refer to different projections of the same stellar structure, namely the inner part of the bar, which expands due to dynamical processes as the bar evolves. It is also very important to highlight that although box/peanuts have a substantial vertical component (and despite their bulge-like morphology), the formation of this bar structure is a process internal to the bar, unlike and unrelated to the violent formation of kinematically hot spheroids: Box/peanuts and barlenses are simply the same inner part of the bar seen at different projections.

Photometrically, all these bar-built structures produce an excess of light in the central region of the galaxy on top of the inward extrapolation of the exponential profile of the main disc. They would thus satisfy one of the criteria to identify bulges, even though they are unrelated to the classical picture of a bulge as a merger-built, kinematically hot spheroid with stars in radial orbits. To distinguish bar-built structures from classical bulges the term “pseudo-bulge” is commonly used (e.g. Kormendy & Kennicutt 2004; Gadotti 2009; Méndez-Abreu et al. 2014; Fisher & Drory 2016; Neumann et al. 2017, amongst many others). However, as described above, bar-built structures also come in two different flavours according to their physical nature. Nuclear discs, inner bars, nuclear rings and nuclear spiral arms are all thought to be built mostly from gas brought to the central region where star formation takes place, whereas box/peanuts (and barlenses) are composed by stars that gradually move from one bar orbital family to another. Therefore, using the term pseudo-bulge collectively to describe all these structures, these two different flavours of bar-built structures, can be misleading. To minimise confusion, Athanassoula (2005) introduced the term “disc-like bulge” (which later produced variations such as “discy pseudo-bulge” and “discy bulge”) to distinguish nuclear discs and related structures from box/peanuts.

Photometric disc-like bulges typically have exponential surface density profiles, and therefore seem to be simply nuclear discs that in contrast to the main galaxy disc are built by bar-driven processes that develop in the main galaxy disc (e.g. Falcón-Barroso et al. 2004, 2006). These nuclear discs often host inner bars, nuclear spiral arms and nuclear rings. The formation of inner bars and nuclear spiral arms appear to be simply scaled-down versions of the formation of a bar and spiral arms in the main galaxy disc (Méndez-Abreu et al. 2019; de Lorenzo-Cáceres et al. 2019), confirming the theoretical work by Wozniak (2015) and Du et al. (2015), who also found that nuclear discs and inner bars can be long-lived. Crucially, however, excluding the nuclear disc, the central regions of disc galaxies are dynamically hotter than the regions where the main disc dominates, and thus these processes still need to be better understood.

On the other hand, the formation of nuclear rings is more strongly connected to the properties of the main bar. In fact, a number of studies using hydrodynamical simulations suggest that nuclear rings form close to the inner Lindblad resonance (ILR) of the main bar, or in the region where the x2 orbital family of the main bar dominates (e.g. Kim et al. 2012; Li et al. 2015; Sormani et al. 2018, and references therein). These theoretical studies show that the radius of the nuclear ring depends on the size and other properties of the main bar. Observationally, Knapen (2005) and Comerón et al. (2010) find indeed corroborating evidence. The interplay and evolutionary connection between nuclear discs and nuclear rings are still unclear, but Cole et al. (2014) argue that the nuclear ring is part of the nuclear disc, namely, its outer rim.

It is important to point out that the build-up of nuclear discs in processes unrelated to bars, as via the accretion of external gas onto unbarred galaxies and mergers, has also been explored in numerical simulations (see Mayer et al. 2008; Chapon et al. 2013). This may explain the presence of nuclear discs in unbarred galaxies, but we note that the nuclear discs produced in these unbarred merger simulations are an order of magnitude less extended than those built by bars. In a different study, Eliche-Moral et al. (2011) argued that their idealised, collisionless simulations of minor mergers do create nuclear discs, sometimes without forming a noticeable bar, provided that the satellite galaxy is not too dense and that the primary galaxy has a massive classical bulge before the merger. We will discuss these simulations further below and conclude that they do not reproduce the observed properties of nuclear discs. Another possibility is that nuclear discs would form via gas inflow due to oval distortions in the main disc. This is akin to the bar-driven formation process, only with a weaker non-axisymmetric component. To date, it is still unclear how rare are nuclear discs in unbarred galaxies (but see Comerón et al. 2010).

In this paper, we take advantage of the integral-field spectroscopy data from the Time Inference with MUSE in Extragalactic Rings (TIMER) project to study the kinematic properties of barred galaxies and nuclear stellar structures with unprecedented spatial resolution. We also take advantage of the vast ancillary data for the TIMER sample to (i), provide evidence that bar-driven processes appear to be the main mechanism responsible for the formation of nuclear discs and related structures; (ii), demonstrate more rigorously the connection between nuclear discs, detected via their kinematic properties, and central exponential components, found via photometric decompositions, and (iii), provide further evidence that nuclear rings are the outer rims of nuclear discs (see also our accompanying paper, Bittner et al. 2020).

This paper is organised as follows. In the next section we introduce the TIMER project and summarise the main aspects concerning sample selection, observations and data reduction. In Sect. 3, we discuss the derivation of the parameters characterising the stellar kinematics, and describe the detection of kinematic signatures of nuclear discs, bars and box/peanuts, as well as the observed kinematic properties of galaxies showing barlenses. We connect the kinematic and photometric properties of nuclear discs in Sect. 4, and discuss the origin of nuclear discs in Sect. 5. In Sect. 6 we present a more general discussion and summarise our main conclusions.

2. The TIMER project

The TIMER project is a survey with the Very Large Telescope (VLT) MUSE integral-field spectrograph (Bacon et al. 2010) of 24 nearby barred galaxies1 with prominent central structures, such as nuclear rings, nuclear spiral arms, inner bars and nuclear discs (see Gadotti et al. 2019, hereafter Paper I). One of the projects’ main goals is to study the star formation histories of such structures to infer the cosmic epoch of the formation of the bar and the dynamical settling of the main disc of the host galaxy. The methodology was demonstrated with a pilot study of NGC 4371 (Gadotti et al. 2015).

The TIMER sample was drawn from the Spitzer Survey of Stellar Structure in Galaxies (S4G, Sheth et al. 2010), which includes only galaxies at distances below 40 Mpc, brighter than 15.5 B-mag, and larger than 1′. The TIMER galaxies are all barred, with stellar masses above 1010M, inclinations below ≈60°, and nuclear stellar structures. The presence of the bar and nuclear structures was assessed from the morphological classifications of Buta et al. (2015), who used the S4G images for their work. The TIMER sample is thus biased towards conspicuous bars and nuclear structures, and it is important to keep this in mind when considering the results discussed in this paper.

Most of the observations were performed during ESO Period 97 (April to September 2016) with a typical seeing of 0.8−0.9″, mean spectral resolution of 2.65 Å (full width half maximum; FWHM), and spectral coverage from 4750 Å to 9350 Å. MUSE covers an almost square 1′×1′ field of view with a contiguous sampling of 0.2″ × 0.2″, which corresponds to a massive dataset of about 90 000 spectra per pointing. The spectral sampling is 1.25 Å per pixel. The total integration time on source for each galaxy was typically 3840 s.

The MUSE pipeline (version 1.6) was used to reduce the dataset (Weilbacher et al. 2012), correcting for bias and applying flat-fielding and illumination corrections, as well as wavelength calibration. The exposures were flux-calibrated through the observation of a spectrophotometric standard star, which was also used to remove telluric features. Dedicated empty-sky exposures and a PCA methodology were employed to remove signatures from the sky background. Finally, the exposures were also finely registered astrometrically, so that the point spread function (PSF) of the combined cube is similar to that in individual exposures. Typically, the averaged signal-to-noise ratio (S/N) per spectral and spatial pixel at the central spaxels of our fully reduced data cubes is approximately 100. We refer the reader to Paper I for further details on the sample selection, observations and data reduction.

Within the TIMER collaboration, Méndez-Abreu et al. (2019) presented the discovery of the first box/peanut found in an inner bar, and de Lorenzo-Cáceres et al. (2019) found evidence indicating that inner bars are long-lived. Neumann et al. (2020) discovered variations in the stellar population properties across galaxy bars that were predicted in idealised simulations and can be reproduced by state-of-the-art cosmological simulations. Furthermore, Leaman et al. (2019), adding data from ALMA, showed a spectacular example of the effects of bars on the interstellar medium, central star formation and stellar feedback.

3. Stellar kinematics

3.1. Methodology

The technical details behind our derivation of the stellar kinematics in TIMER were extensively presented in Paper I and Gadotti et al. (2015). Here we simply summarise the essential aspects. The stellar line of sight velocity distributions (LOSVDs) were parameterised as Gauss-Hermite functions with four parameters (following van der Marel & Franx 1993): velocity (v), velocity dispersion (σ), and the h3 and h4 higher-order moments. The h3 and h4 parameters can be used to examine in further detail the orbital structure of the stellar systems in question. For instance, near-circular orbits, with a distribution of v/σ, produce h3 values that are anti-correlated with v, whereas elongated orbits result in a correlation between v and h3. High values of h4 suggest the superposition of structures with different LOSVDs (see e.g. Bender et al. 1994, and references therein). For the derivation of 2D maps of v, σ, h3 and h4, the spectra from the data cube of each galaxy were spatially binned to ensure a minimum S/N of approximately 40 per spectral pixel. This was done using the Voronoi binning technique as described in Cappellari & Copin (2003). To derive systemic velocities and central velocity dispersions, we combined all spectra within a circular aperture. The spectra from each Voronoi bin or aperture were fitted with the pPXF code of Cappellari (2017, see also Cappellari & Emsellem 2004; Cappellari et al. 2011). The entire ensemble of procedures was performed using the GIST2 pipeline (which employs pPXF; Bittner et al. 2019).

The only fundamental difference between the kinematic maps presented in Paper I (where we presented the maps corresponding to NGC 1097 and NGC 4643 only) and the ones derived here is the rest-frame wavelength range employed. In Paper I we opted to use the range between 4750 and 5500 Å after checking that similar results are obtained when the full MUSE wavelength range is employed. That approach has the advantage of necessitating less computation time. In addition, it avoids complications arising from bright emission lines, particularly if they are not masked, but modelled. Some emission lines are nevertheless present in the restricted wavelength range, and those were masked before the fitting procedure. However, in this paper, we used an extended wavelength range (4800−8950 Å) that is close to the full wavelength range provided by MUSE. A comparison between the maps derived employing these two different approaches shows that the results are qualitatively identical and quantitatively very similar. However, the maps derived using the extended wavelength range reveal some results in a more enhanced fashion. For example, regions with elevated absolute values of h3 and h4 appear more conspicuously and sharply defined. This seems to be related to the fact that the instrumental spectral resolution increases with wavelength, and thus lower values of velocity dispersion, as well as deviations from a pure Gaussian function for the LOSVDs (i.e. values of h3 and h4 different from zero), can be more robustly measured. Further, the signal is boosted by the addition of several absorption features, including the calcium triplet. When using the extended wavelength range, for two galaxies (NGC 7140 and NGC 7755) small differences in the continuum shape between the galaxy spectra and the input library of synthetic models required more than a multiplicative low-order Legendre polynomial included in the fit (as done in Paper I). We thus included in the fit an 8th-order multiplicative polynomial plus a 4th-order additive Legendre polynomial to account for such small differences. In addition, in this paper, we also employed a non-constant line spread function (LSF) to account for the wavelength dependence of the instrumental spectral resolution. We adopted the LSF derived in Bacon et al. (2017, their Eq. (8)).

3.2. Systemic velocities and central velocity dispersions

For each galaxy, the heliocentric systemic radial velocity and central velocity dispersion were measured by combining all central spectra within a circular aperture corresponding to one eighth of the galaxy effective radius, re/8. For re we use the measurements derived by Muñoz-Mateos et al. (2015) using S4G 3.6 μm growth curves. These measurements are presented in Table 1, along with error estimates and the corresponding values available at the Lyon Extragalactic Data Archive (LEDA). Our error estimates were derived from 30 Monte Carlo realisations. We note that to properly compare the measurements of central velocity dispersion we adapted our measurements to follow the same aperture corrections of LEDA, which correspond to an aperture of 0.595 h−1 kpc.

Table 1.

Systemic velocities and central velocity dispersions.

In Fig. 1 we compare our measurements with the LEDA values. The agreement in the measurements of systemic velocity is excellent. Concerning the central velocity dispersion, although there are no available measurements in LEDA for nine galaxies, the agreement is very good, with only one significant outlier, NGC 1300. The LEDA values for the central velocity dispersion in NGC 1300 are a factor of two larger than ours. Curiously, most of the LEDA measurements are quoted as from Davies et al. (1987) but an inspection of the article shows that NGC 1300 is not part of the sample of galaxies in that study. In fact, the measurements by Davies et al. (1987) correspond to NGC 2300. Another measurement quoted in LEDA is from an unpublished study that reports 145 ± 22 km s−1, a value closer to our measurement (102 ± 8 km s−1). Part of the differences between our measurements and the values quoted in LEDA can be due to differences in seeing and pixel size amongst the different observations.

thumbnail Fig. 1.

Left: TIMER measurements of systemic velocity against the values available at LEDA. Right: TIMER measurements of central velocity dispersion against the LEDA values. The TIMER velocity dispersions were adapted to follow the same aperture corrections of LEDA, and the error bars in both panels correspond to 1σ errors. NGC 1300 is the only significant outlier but this is due to a mistake in LEDA (see text for further details; the LEDA team has been informed).

3.3. Kinematic maps

In this section, we present and discuss signatures found in the kinematic maps derived with GIST that reveal the presence of different structural components. Here we focus on the general trends and leave to Appendix A a discussion on the more complex cases. The high-level data products derived to produce all kinematic maps discussed in this study are publicly available3.

3.3.1. Kinematic signatures of nuclear discs

The 2D kinematic maps provide a test to understand if the observed nuclear stellar structures indeed arise from bar-driven secular evolution processes. If they were built via such dissipative processes within the main disc, they are foreseen to have dynamical properties similar to main discs (which are also formed through dissipative processes, albeit not induced by bars), namely, elevated rotational support, low velocity dispersion and near-circular orbits4 with an angular momentum vector aligned to that of the main disc. Further, since they should be dynamically distinct from the underlying main galaxy disc, the LOSVDs measured from spectra from such structures are likely to show elevated values of h4 (and elevated absolute values of h3), a signature of the overlapping of stellar structures with different velocity distributions: the nuclear structure that dominates the light in such regions and the underlying main galaxy disc (see Cole et al. 2014). In principle, external gas accretion unrelated to the bar could also build stellar structures with these properties, but such structures are likely to show at some stage an angular momentum vector that is not aligned to that of the main disc. Moreover, in simulations, bars are the most common mechanism to remove angular momentum from the gas in a uniform fashion down to the central few hundred parsecs (see e.g. Athanassoula 1992). On the other hand, if violent processes such as dry mergers build such nuclear structures, then they are expected in most merger configurations to have elongated orbits (showing no v − h3 anti-correlation) with relatively low rotational support and high velocity dispersion. We discuss such scenarios further in Sect. 5.

A careful analysis of Figs. 2 and A.1 reveals that indeed the TIMER galaxies fit remarkably well the picture in which nuclear structures originate from bar-driven processes. For example, IC 1438 clearly hosts a rapidly rotating nuclear component within a radius of about 5″, with a kinematic axis well aligned with that of the main disc, as seen in the corresponding velocity map. The velocity dispersion map shows that the rapidly rotating component is characterised by low velocity dispersion, producing a well defined and well localised dip in the central region of the velocity dispersion map, which would, in the absence of this component, simply show an increasing trend in σ. The h3 map shows also very clearly that this moment of the LOSVD is anti-correlated with v, a robust signature of near-circular orbits with a range of values for v/σ. This is also the case for the outer, main disc of the galaxy, for example at radii around 15″. Finally, the h4 map shows a sharp increase at the region dominated by the rapidly rotating nuclear component, due to the superposition of the low-σ nuclear component on top of the relatively high-σ underlying component.

thumbnail Fig. 2.

Radial velocity, velocity dispersion, h3 and h4 maps for the stellar components in the TIMER galaxies, as indicated. The plotted range of the parameter measured is indicated at the bottom right corner of each panel. For radial velocity and velocity dispersion these are given in km s−1. The isophotes shown are derived from the MUSE data cube reconstructed intensities and are equally spaced in steps of about 0.5 mag. On average, 1″ corresponds to ≈100 pc. North is up, east to the left. See Appendix A for the maps corresponding to the rest of the galaxies in TIMER.

IC 1438 thus satisfies all the kinematic criteria discussed above that need to be satisfied by nuclear stellar discs built from dissipative processes on a pre-existent underlying disc. In addition, it is important to stress that the kinematic maps show that the properties of the nuclear stellar disc are not part of a continuous distribution covering the whole MUSE field. For example, one sees clearly in IC 1438 that v does not vary monotonically from the outer radii inwards. The stellar velocity peaks at the outer parts of the field and steadily decreases until it increases sharply again at the region of the nuclear disc. This is a clear indication that the nuclear disc is a separate component from the main galaxy disc. In fact, at the same galactocentric radius, stars in the nuclear disc are rotating around the galaxy centre faster than the stars in the main disc (see Fig. 2). This implies that – in the region where the nuclear disc resides – the nuclear disc is dynamically colder than the main disc and its stars have orbits closer to near-circular orbits5. As discussed above, this is also the reason behind the high values of h4 in the nuclear disc.

To corroborate these findings we also produced maps of local v/σ (see Fig. 3). The figure shows again and even more clearly that the nuclear structures morphologically identified by means of visual inspection by Buta et al. (2015) are stellar structures with very elevated dynamical rotational support and separate from the main galaxy disc. We note that (Buta et al. 2015) have not necessarily identified these nuclear structures as nuclear discs but rather as nuclear disc features, such as rings, lenses, spiral arms and bars. However, Figs. 2 and 3 show that these nuclear structures are extended, and in our accompanying paper (Bittner et al. 2020) we present evidence from the stellar population properties that these nuclear structures appear to extend all the way to the galaxy centre. We thus have very strong evidence that the TIMER galaxies host nuclear discs with kinematic properties consistent with a bar-driven origin. In Sect. 5 we will elaborate on alternative scenarios for the building of these nuclear discs, and also discuss further the connection between nuclear discs and nuclear rings.

thumbnail Fig. 3.

Maps of v/σ for the stellar components in the TIMER galaxies, as indicated. The plotted range is indicated on the bottom right corner. The isophotes shown are derived from the MUSE data cube reconstructed intensities and are equally spaced in steps of about 0.5 mag. On average, 1″ corresponds to ≈100 pc. North is up, east to the left.

thumbnail Fig. 3.

continued.

Figure 4 shows deprojected radial profiles of the v/σ measurements at each Voronoi bin, after correcting v for inclination. We use the values determined with S4G data by Muñoz-Mateos et al. (2015) for the inclination and position angle of each galaxy. NGC 1365 and NGC 6902 are not included in this figure, since the signatures of a nuclear disc in these galaxies are not as clear as in the remaining of the sample (these cases are discussed in Appendix A). NGC 1291 is also excluded, as determining v/σ is difficult due to the low inclination of the galaxy (11°). These profiles allow us to obtain the maximum value of v/σ in the nuclear discs, assuming that the nuclear disc is in the same plane as the main disc. The region within which the maximum is searched for is delimited by a circumference centred at the galaxy centre and crossing the position of the first minimum in v/σ along the disc major axis, beyond the centre. The measurements of v/σ considered in this region are those for each Voronoi bin. We also define the kinematic radius (rk) of the nuclear disc as the deprojected distance from the galaxy centre of the Voronoi bin showing the maximum value of v/σ in the region dominated by the rapidly rotating nuclear disc. The values of maximum v/σ and rk are shown in Fig. 4 with the horizontal and vertical dashed lines, respectively.

thumbnail Fig. 4.

Deprojected radial profiles of v/σ at each Voronoi bin, with v corrected for inclination. NGC 1291, 1365 and NGC 6902 are not included. The values of maximum v/σ and rk are shown with the horizontal and vertical dashed lines, respectively.

In Fig. 5 we present the distributions of the values of maximum v/σ and rk. Interestingly, although the v/σ distribution is skewed towards higher values, it also shows a broad range of v/σ values, as low as unity. We note that the bar-built nuclear disc in the hydrodynamical simulations of Cole et al. (2014) has a peak v/σ of ≈1.8 (the vertical dotted line in Fig. 5), fitting well within the distribution of observed values. The same is true for the size of the simulated nuclear disc, although it lies on one of the extremes of the observed rk distribution. The simulated nuclear disc is as large as the largest observed nuclear discs. We note that the values of v/σ and rk taken from Cole et al. (2014) correspond to the end of their simulation, at 10 Gyr, when the nuclear disc has reached its maximum size.

thumbnail Fig. 5.

Distribution of the maximum value of v/σ (top) and of the kinematic radius (bottom) of the observed nuclear discs. NGC 1291, 1365 and NGC 6902 are not included. The vertical dotted lines mark the values corresponding to the bar-built nuclear disc in the model of Cole et al. (2014).

Finally, further insights can be gained by employing the λR parameter introduced by Emsellem et al. (2007), which quantifies the projected stellar angular momentum per unit mass. We thus calculated λrk, which is λR integrated within the radius rk, as in Emsellem et al. (2007), again using the position angle and ellipticity derived for the main disc by Muñoz-Mateos et al. (2015) with S4G data. In Fig. 6 we plot λrk against the observed (projected) ellipticity at the radius of the nuclear disc, rk, derived with the S4G radial profiles of ellipticity. The positions of our measurements in this diagram show clearly that the region within rk in the TIMER galaxies has the same angular momentum support as the fast rotators of Emsellem et al. (2011). In addition, most systems are consistent with the elevated edge-on ellipticities (or intrinsic flatness) of disc systems. However, as discussed in Cappellari et al. (2007), this diagram is only rigorously valid for stellar systems with a density stratified on homologous oblate ellipsoids (as in the model of Binney 2005), which is evidently not the case of our sample galaxies. The central region of the TIMER galaxies hosts not only the nuclear disc but at least also the main disc and bar, which means that the intrinsic flatness of the TIMER nuclear discs as seen in this diagram is at best a lower limit.

thumbnail Fig. 6.

λrk plotted against the observed (projected) ellipticity at the radius of the nuclear disc. The dashed line is the upper envelope of slow rotators in Emsellem et al. (2011). The solid line shows the relation for a model of oblate stellar systems viewed edge-on (see Binney 2005; Cappellari et al. 2007; Emsellem et al. 2007). The dotted lines correspond to the location of the same model for varying edge-on ellipticity (intrinsic flatness), from top (0.85) to bottom (0.35) in steps of 0.1, with edge-on systems on the relation and face-on systems towards the origin, as in Emsellem et al. (2011). NGC 1291, 1365 and NGC 6902 are not included.

3.3.2. Kinematic signatures of bars

The kinematic maps in Figs. 2 and A.1 also show kinematic signatures of the presence of bars. Take IC 1438 again as an example. Between the main disc and the rapidly rotating nuclear disc, both showing an anti-correlation between v and h3, one sees a region along the bar in which v and h3 are actually correlated. From the bottom-left corner of the h3 map, towards the top-right, one first sees blue/green bins, then red/yellow bins, then blue bins again and the inverse patterns after crossing the centre. The correlation between v and h3 was shown by Bureau & Athanassoula (2005) and Iannuzzi & Athanassoula (2015) to result from the superposition of a bar and a disc, and is seen where both components contribute more or less equally (see also Li et al. 2018).

A similar signature is also seen in NGC 1300, 1433, 3351, 4643, 5850, 7140 and NGC 7755. Possibly due to projection effects, dust extinction and/or lower physical spatial resolution, this signature is somewhat less clear in NGC 613, 4303, 4981, 4984, 5248, 5728 and NGC 7552. Since all these galaxies are known to have bars from studies of their morphology, this result is not terribly surprising, but it lends support to our theoretical understanding of the stellar dynamics in barred galaxies.

Interestingly, NGC 1433 shows a correlation between v and h3 also within ∼4 arcsec from the centre, which suggests the presence of an inner bar (see also Bittner et al. 2019, for a zoom-in view of the kinematics in this region). However, archival HST images inspected by de Lorenzo-Cáceres et al. (2019) show no clear morphological signature of an inner bar (but see Erwin 2004; Buta et al. 2015).

We also point out the conspicuous drops in velocity dispersion at the ends of the inner bars in NGC 1291 and NGC 5850. These σ-hollows were found by de Lorenzo-Cáceres et al. (2008) to be a characteristic signature of inner bars, and these two galaxies were studied in detail already in Méndez-Abreu et al. (2019) and de Lorenzo-Cáceres et al. (2019).

3.3.3. Kinematic signatures of box/peanuts

Our maps also reveal kinematic signatures of the presence of box/peanuts. Debattista et al. (2005) employed numerical simulations to show that box/peanuts can be detected in face-on galaxies by examining the spatial distribution of h4: Box/peanuts imprint two significant minima along the bar major axis at the positions of the box/peanut vertices (see also Iannuzzi & Athanassoula 2015; Li et al. 2018). Méndez-Abreu et al. (2008) showed that this diagnostic works observationally, detecting box/peanuts in primary bars, and, more recently, Méndez-Abreu et al. (2019) were able to find for the first time a box/peanut in an inner bar, that of NGC 1291.

IC 1438 provides again a good example. The very central region, that is, within a radius of about 2″, shows values of h4 close to zero. This region is surrounded by elevated values of h4 where the rapidly rotating nuclear disc dominates, as discussed above. Just outside this region – at radii between about 5 and 10 arcsec – h4 takes negative values, becoming close to zero again at larger radii. This is the signature just mentioned of the presence of a box/peanut.

We find the same signature clearly in NGC 613, the inner bar of NGC 1291 (as shown already in Méndez-Abreu et al. 2019), NGC 1300, 4303, 4643, 4981, 4984, 5728, 5850, 7140, 7552, and NGC 7755. Simulations show that this signature is weaker for galaxies with inclination angles larger than 30° (see e.g. Fig. 30 in Iannuzzi & Athanassoula 2015). Most of the aforementioned galaxies are below or close to this threshold, with the exception of NGC 4981, 4984, 7140 and NGC 7755, where this interpretation should be taken carefully.

NGC 5728 is an interesting case. The h4 map shows clear minima along the bar major axis, and the isophotal contours show regions that are slightly offset from the bar major axis on opposite sides of the bar at each side from the centre. This S-shape configuration is produced by what Erwin & Debattista (2013) called spurs. Such spurs are photometric signatures of box/peanuts seen at an angle (Athanassoula & Beaton 2006, see also Beaton et al. 2007), and thus this galaxy is a beautiful illustration where both photometric and kinematic signatures of box/peanuts coincide.

Considering the TIMER sample altogether, we find a clear signature of the presence of a box/peanut in 13 of the 21 galaxies. As discussed in Appendix A, at least in some cases, the signature absence may be due to the restricted fields studied here. Therefore, our results indicate a lower limit for the fraction of box/peanuts in our sample of massive barred galaxies of at least 62%. This is consistent with the estimate presented by Erwin & Debattista (2017) also for massive barred galaxies, namely, 79%, even though our analysis relies entirely on kinematics and theirs on photometry. On the other hand, in a study employing the kinematic diagnostic on a sample of 10 galaxies, Méndez-Abreu et al. (2014) found a box/peanut fraction of 50%.

3.3.4. Kinematic signatures of barlenses

As mentioned in the Introduction, barlenses are thought to be the face-on projection of box/peanuts, and thus also part of the bar, even though they extend further from the bar major axis than the remaining of the bar (Laurikainen et al. 2005, 2007, 2011; Athanassoula et al. 2015; Laurikainen & Salo 2016). If that is the case, one expects to see, in the region dominated by the barlens, both the h4 minima that characterise box/peanuts and the v − h3 correlation that characterise bars (see Iannuzzi & Athanassoula 2015, although one should be careful to disregard regions dominated by a nuclear disc). In the S4G images, we see a clear barlens morphology in five of the galaxies studied here: NGC 1300, 3351, 4643, 4984 and NGC 7755. Except for the latter, all are classified by Buta et al. (2015) as indeed showing a barlens. In these five galaxies, but more strongly in NGC 4984 and NGC 7755, the expected signatures are seen. In NGC 4984, the h4 minima form a thick ring with inner and outer radii of about 10″ and 20″, respectively. In the same region one sees the v − h3 correlation. Likewise, NGC 7755 shows the same signatures in radii between 5″ and 10″. These features are very similar to what is seen in the simulations of Iannuzzi & Athanassoula (2015, their Fig. 29).

While this is consistent with our current understanding of barlenses and box/peanuts, we point out that seven other galaxies in this sample are classified by Buta et al. (2015) as having barlenses, but the expected signatures are not clear. Evidently, part of the reason could be attributed to the relatively small spatial coverage in some of these galaxies. On the other hand, the visual classification of barlenses is prone to ambiguities and should be considered carefully, as nuclear discs, and even classical bulges, can be confused as barlenses. As discussed in Gadotti et al. (2015), the case of NGC 4371 is very illustrative. Morphologically, the signature of a barlens appears very clear. However, once kinematic information is included in the analysis, one concludes that the concerned structure is actually a rapidly rotating – and rather large – nuclear disc. This demonstrates the importance of assessing the nature of morphological components through a detailed analysis of the corresponding stellar kinematics.

Apart from NGC 4371, Buta et al. (2015) also include NGC 613, 1097, 1291, 5728, 5850 and NGC 7552 as having barlenses. In the case of NGC 1097 and NGC 1291 the TIMER fields are too restricted to properly study the corresponding regions. Recently, de Lorenzo-Cáceres et al. (2019) found that the component visually classified by Buta et al. (2015) in NGC 5850 as a barlens is actually the nuclear disc, and that this is possibly the case for NGC 1291 also but it is difficult to ascertain that in the case of this galaxy given its low inclination. For the remaining three galaxies, the signatures are weak: The h4 minima, clearly seen along the bar major axis are not prominent in the perpendicular direction. As discussed in Sect. 3.3.2, the presence of a v − h3 correlation is not very significant in NGC 613, 5728 and NGC 7552, possibly due to projection effects, dust extinction and/or lower physical spatial resolution.

3.4. Previous observations with integral-field spectrographs

Some galaxies in the TIMER sample have been observed before with different integral-field spectrographs. It is particularly instructive to qualitatively compare our observations with those performed with SAURON (Bacon et al. 2001; de Zeeuw et al. 2002) and WiFeS (Dopita et al. 2007, 2010), which as MUSE operate in the optical wavelength range. Four TIMER galaxies were studied previously with SAURON: NGC 4371 (Cappellari et al. 2011), NGC 4643 (Cappellari et al. 2011; Seidel et al. 2015a), NGC 5248 (Dumas et al. 2007) and NGC 5850 (de Lorenzo-Cáceres et al. 2008, 2013). NGC 7552 was studied with WiFeS by Seidel et al. (2015b). The spectral resolution obtained with SAURON is about a factor of two lower than in MUSE observations, which in turn is a factor of ≳2 lower than that provided by WiFeS. However, MUSE excels in sensitivity and spatial sampling. Part of the gains in sensitivity comes from the hosting telescopes. The collecting area of the VLT is roughly about four times larger than that of the William Herschel Telescope, which hosts SAURON, and roughly 16 times larger than the 2.3 m telescope at Siding Spring Observatory, which hosts WiFeS. The spatial sampling of MUSE corresponds to 0.2″ per spatial element, whereas in both SAURON and WiFeS this is about 1″.

The remarkable spatial sampling and sensitivity of MUSE weigh significantly in studies of central structures in disc galaxies, such as this study. This can be appreciated when comparing the TIMER observations with those mentioned above. The MUSE kinematic maps reveal a multitude of details that remain hidden in the SAURON and WiFeS datasets (see e.g. the analysis concerning NGC 5850 in de Lorenzo-Cáceres et al. 2019). Surely, indications of a kinematically colder nuclear component are seen in the SAURON and WiFeS data, but the lower spatial sampling makes it more difficult to ascertain that the component is a nuclear disc, separate from the main galaxy disc. Furthermore, a more detailed analysis of the kinematic maps, for example by searching for the kinematic signatures of bars, box/peanuts and barlenses, as performed above, is certainly precluded by the relatively low spatial resolution of the SAURON and WiFeS data.

We can also compare our measurements of velocity and velocity dispersion with those presented by Venturi et al. (2018, for NGC 1365), Shimizu et al. (2019, for NGC 5728), and Bosma et al. (2010, for NGC 1291). The first two studies employed the exact same datasets as employed here, but performed independent analyses, and, reassuringly, the maps presented agree both qualitatively and quantitatively with those we present in Figs. 2 and A.1. The study presented in Bosma et al. (2010) is based on spectra taken with the deployable integral-field units of the FLAMES/GIRAFFE spectrograph at the VLT. The spectral resolution is about five times better than that of MUSE but the spectra are restricted to the calcium triplet wavelength range. The authors reported a central velocity dispersion of ≈195 km s−1, which is slightly above the value we derive, namely 168 km s−1, but given the typical uncertainty of ∼10 km s−1, the two measurements are statistically equivalent. Bosma et al. (2010) also reported measurements off the main bar in NGC 1291 at about 20″ from the centre: one can see a drop in velocity dispersion to ≈110 km s−1 and a radial velocity of only a few km s−1 (the galaxy is very close to face-on), and both features match our own measurements very well (see Fig. 2).

4. A comparison with structural analysis from photometry

As shown above, detailed maps of v, σ, h3 and h4 are a powerful tool to understand the nature of the different stellar structures in a galaxy. These maps show that in all galaxies studied here (with no more than two possible exceptions) there clearly is a fast rotating nuclear disc, separate from the main galaxy disc.

In photometric decompositions, the central component in disc galaxies is often fitted using the Sérsic (1968) function, where the Sérsic index n regulates how centrally concentrated stars are, with large values of n produced by high concentrations of stars. Typically, one would expect that nuclear discs have an exponential surface brightness radial profile (i.e. with n ≈ 1), just as main discs. In practice, it is common to associate disc-like bulges to values of n ≤ 2, whereas classical bulges are associated to higher values.

However, establishing the physical nature of bulges via photometric structural analysis is not straightforward, since a stellar component with n ≈ 1 is, in principle, not necessarily dynamically supported by rotation. The latter can only be directly probed with measurements of the stellar kinematics. Conversely, a stellar structure with n >  2 may not necessarily be dispersion-dominated. Given the large body of ancillary data for the TIMER galaxies, we are in an excellent position to directly compare our stellar kinematic measurements with results from studies on the photometric properties of these galaxies. This allows us to test whether the nuclear discs we find in the TIMER sample are correctly identified via photometric structural analysis, and this is the main goal of this section.

It is important to stress that an accurate image decomposition depends on the physical spatial resolution, depth and nature of the data (dust effects at short wavelengths are known to bias the results, see e.g. de Jong 1996; Gadotti et al. 2010; Pastrav et al. 2013). It also depends on the procedures employed to account for the exceedingly complex stellar structure in disc galaxies. In fact, many authors have argued that the distinction between disc-like and classical bulges using the Sérsic index alone is prone to uncertainties, and that this analysis is only more robust when several criteria are used together, including criteria based on stellar kinematics and intrinsic shape (see e.g. Fisher & Drory 2016; Neumann et al. 2017; Costantin et al. 2018a,b).

Despite these caveats, to keep the comparison straightforward, we will simply verify what values of the Sérsic index n are obtained for the photometric bulges in decompositions of the TIMER galaxies, and use the criterion n ≤ 2 for a successful identification of a nuclear disc, as opposed to a classical bulge. We will base this comparison on two recent studies separately in the following subsections. These studies were chosen here given the large overlap with the TIMER sample and the effort made to produce accurate results, with sophisticated models and careful procedures (which typically include individual inspection of the fits and numerous fits per galaxy to understand the uncertainties and avoid local χ2 minima).

4.1. Kim et al. (2014)

This work made use of S4G 3.6 μm images of 144 barred galaxies to derive structural parameters of bulges, bars and discs using BUDDA (de Souza et al. 2004; Gadotti 2008), taking advantage of the exquisite depth of the S4G data and the minimised impact of dust at these wavelengths. Nuclear point sources and outer disc breaks were accounted for to avoid biasing the bulge parameters. The PSF was determined for each image separately and modelled as a circular Moffat (1969) function. The surface brightness radial profiles of bars were modelled using a Sérsic function, and a bar was included in the models of all 13 TIMER galaxies that also belong to their sample.

For eight of the 13 galaxies, the photometric bulge Sérsic index is less than two, so in the majority of the cases the photometric bulge is identified as a disc-like bulge, in agreement with our kinematic analysis. For the remaining five galaxies, the bulges are found to have n >  2, which would in principle indicate that in these cases the Sérsic index fails to correctly identify the nature of the central component. However, one important caveat is that some galaxies may host composite bulges, that is, a small classical bulge at the very centre surrounded by a nuclear disc (see e.g. Erwin et al. 2015, who reported effective radii for small classical bulges ranging from 25 to 430 pc). For example, NGC 1291 clearly has a nuclear disc and an inner bar, and other studies have shown that it has also a small classical bulge with an effective radius of 416 pc within the nuclear disc (see Méndez-Abreu et al. 2019; de Lorenzo-Cáceres et al. 2019). Thus, the bulge component in the fits of Kim et al. (2014) accounts for both the nuclear disc and the small classical bulge, which pushes the Sérsic index of the central component to n = 2.7. The other bulges with n >  2 are those in IC 1438, NGC 1433, NGC 4303 and NGC 7140. In IC 1438 there is a clear and substantial peak in the stellar velocity dispersion in the inner 1−2″, combined with a drop in h4, that suggests that this galaxy too has a small dispersion-dominated component within the nuclear disc. Nevertheless, for the remaining three galaxies we see no evidence of a small classical bulge in the kinematic maps.

In this context, it is important to note that the S4G images, with a typical PSF FWHM of about 2″, are inadequate to separate central components on those spatial scales. Therefore, the issue of composite bulges highlights the importance of high physical spatial resolution in imaging data, and the accounting of the different structural components in the models used to fit galaxy images.

We also point out that no bulge was found to have n ≈ 4, which is expected for central regions dominated by classical bulges and for massive elliptical galaxies (e.g. Gadotti 2009, and references therein). In fact, the largest value of n found is below 3, which indicates that there is no galaxy in the TIMER sample with a dominant classical bulge, a conclusion that is corroborated in our accompanying paper (Bittner et al. 2020) through an analysis of stellar populations properties. This may result from the selection of the TIMER sample, which favours galaxies hosting nuclear components with a disc-like morphology. However, it is not clear if this criterion rejects galaxies hosting large classical bulges.

4.2. Salo et al. (2015)

All galaxies in the TIMER sample were studied in Salo et al. (2015), which also employed S4G 3.6 μm images but did not account for disc breaks. A single oversampled PSF image that accurately reproduces the particular PSF 2D shape in the S4G images was used in all decompositions, and bars were modelled using a modified Ferrers profile. The decompositions were performed using GALFIT (Peng et al. 2002, 2010). Interestingly, bars were not included in the models for NGC 4303, 4981, 5248 and NGC 6902.

Only two of the photometric bulges were modelled with n >  2. One is again in NGC 1291, and the other is in NGC 7552. As discussed above, NGC 1291 hosts a small classical bulge, but we find no evidence for a similar component in NGC 7552. As in the study by Kim et al. no bulge was found with n >  3.

A comparison between the values obtained by the two studies for the Sérsic index of the photometric bulges shows some noticeable discrepancies, which again highlight how large is the uncertainty in the measurement of the Sérsic index. While for many galaxies the agreement between the two studies is remarkable, with an absolute difference in n of only 0.1−0.2, in some cases the difference reaches values above unity, and the measurements of Kim et al. are systematically above those of Salo et al. (2015). The median and mean values of the absolute difference in the Sérsic index measurements are 0.9 and 0.7, respectively. These discrepancies also show that even using the same dataset, a different outcome may result if models and/or techniques employed are different.

4.3. Exponential photometric bulges are nuclear discs

In the previous two subsections we have seen that photometric decompositions can retrieve reasonably well the nature of photometric bulges, and do so in the majority of the cases studied here. It is important to stress, however (and in addition to the caveats on photometric decompositions mentioned above), that the TIMER sample consists of nearby galaxies, and thus the physical spatial resolution of the images employed in these studies is relatively high. Decompositions employing images of more distant galaxies, with lower physical spatial resolution, will presumably not perform as well, but to quantify this is beyond the scope of this paper.

A powerful comparison to further test whether the photometric bulges modelled by Kim et al. (2014) and Salo et al. (2015) are indeed the rapidly rotating nuclear discs we identify with the MUSE TIMER data cubes – or are at least dominated by them, in the case of composite bulges – concerns comparing their sizes, as measured by the different techniques. The photometric decompositions readily provide the effective radii (re) of the fitted component, but a different approach is needed to derive the dimensions of the nuclear discs from the kinematic measurements. To this end, we use the kinematic radius (rk) defined above as the distance from the galaxy centre of the spatial bin showing the maximum value of v/σ in the region dominated by the rapidly rotating nuclear disc. If the photometric and kinematic components are in fact the same physical entity then the ratio rk/re must be clustered around a single value, following a normal distribution. Figure 7 shows that this is indeed the case: rk/re is clustered around 0.86. Fitting a normal distribution to our measured values of rk/re yields a mean value of 0.86 and σ = 0.19. In addition, we applied to our data the statistical tests presented in D’Agostino et al. (1990) and Anscombe & Glynn (1983) to verify that the distribution in Fig. 7 has indeed both skewness and kurtosis that are statistically compatible with a normal distribution. Furthermore, the value of σ derived above is compatible with the typical uncertainty in the measurement of re, which is ≈20%.

thumbnail Fig. 7.

Distribution of the ratio between nuclear disc size, defined kinematically as the radius with maximum v/σ, and the effective radius of the photometric bulge component fitted through careful 2D multicomponent photometric decompositions. The values are clustered around a single value, following a normal distribution (rk/re = 0.86 ± 0.19), which shows that the photometric component defined as the bulge is (or is dominated by) the rapidly rotating nuclear disc.

One important aspect to consider is that correlations between rk and re with other parameters could produce the results presented in Fig. 7, even if there is no physical connection between the nuclear discs and photometric bulges in our sample. For example, one could argue that a galaxy with a large disc would naturally have a large photometric bulge and a large nuclear disc, and thus the correlation between rk and re would be trivial. We test this possibility using as a proxy for disc size (R25.5) the 25.5 mag arcsec−2 isophotal radius at 3.6 μm, as derived by Muñoz-Mateos et al. (2015) with S4G images (see Table 2). We find that the Pearson correlation coefficient R for the three relations, rk versus re, rk versus R25.5, and re versus R25.5, is, respectively: 0.87, 0.63 and 0.75. Thus, the correlation between rk and re is stronger, and therefore not trivial. Considering these results altogether, we can safely state that exponential “bulges” in photometric decompositions are nuclear discs. It is also important to bear in mind that rk and re are measured through completely different techniques, which lends support to this conclusion.

Table 2.

Properties of nuclear discs, bars and galaxies in the TIMER sample.

5. The origin of nuclear discs

While the results above are consistent with the scenario in which bars drive gas to the central region and ultimately build the observed nuclear discs, they do not rule out that, at least in some cases, accretion of gas could have been promoted by an interaction for example, before the formation of the bar, and therefore be unrelated to it. Thus, to shed more light in this discussion (and given that bars are expected to grow longer and stronger with time; e.g. Athanassoula 2003), in this section, we explore the properties of nuclear discs as derived from the TIMER kinematic maps in conjunction with other relevant bar properties. Table 2 presents the physical parameters we explore. In particular, we use the values determined by Kim et al. (2014) from their photometric decompositions for the bar radius (Rbar), ellipticity (ϵ) and bar-to-total luminosity ratio (Bar/T). Further, Díaz-García et al. (2016), presented a number of measurements of bar strength, from which we use QB and A2. QB indicates how strong the non-axisymmetric potential of the bar is as compared to the axisymmetric component of the total stellar mass distribution. It thus provides an indication of the impact of the bar on the dynamics of gas and stars in the host galaxy, by taking into consideration the effects of any central spheroid and the main disc itself. On the other hand, A2 is a measure of the m  =  2 Fourier component that is directly connected to the bar non-axisymmetry. More specifically, A2 is calculated radially, normalised by the value of A0 (the axisymmetric Fourier component) at each radius, and the value we employ here is the peak value of A2/A0 in the bar. Therefore, QB and A2 provide complementary information on the bar and associated secular evolution processes.

Figure 8 shows that the radii of nuclear discs correlate significantly with bar radii, which is expected in the theoretical framework of the orbital structure in bars and how it evolves. In fact, the gas brought to the central region by bars is expected to accumulate in the region where the x2 orbits of the bar dominate over the bar x1 orbits6, where a nuclear stellar structure is thus formed (Athanassoula 1992; Kim et al. 2012; Li et al. 2015; Sormani et al. 2018; Seo et al. 2019, see also Romeo & Fathi 2015, 2016 who suggested that the disc instability radial profile can be used to predict the sizes of nuclear structures). As the bar evolves and grows, the region where the x2 orbits dominate grows as well, hence the expected correlation7. This correlation between rk and Rbar is therefore consistent with the picture in which the nuclear stellar discs studied here were formed from gas brought to the central region by the bar, and not built before the formation of the bar. In addition, in this picture, gas falling onto the outer boundary of the region dominated by x2 orbits is prone to form stars and produce star-forming nuclear rings. Indeed, Seo et al. (2019) found in their simulations that nuclear rings are larger when bars evolve for a longer period of time. As expected, Comerón et al. (2010) found that the relation between the sizes of bars and nuclear rings is such that the upper limit of the distribution of nuclear ring sizes correlates with bar size. It is unclear if the larger scatter seen in their work is due to uncertainties in the measurements, or whether it is real and brought out by the larger sample. Nevertheless, the connection between nuclear ring size and bar radius in that work is also clear.

thumbnail Fig. 8.

Relations between the kinematic sizes of nuclear discs and selected bar properties. Top left panel: rk is correlated with bar size (semi-major axis). The solid line is a simple linear regression fit to the data, whereas the dotted line corresponds to rk = 0.1 Rbar, the relation suggested in Shlosman et al. (1989). Top middle and right panels: we present how the kinematic size of nuclear discs relates to bar ellipticity and bar-to-total luminosity ratio, respectively. Bottom panels: relations between rk and different measures of bar strength (from Díaz-García et al. 2016, see text for details). Bar radii and ellipticities are from Kim et al. (2014) and are deprojected. The Pearson correlation coefficient R is indicated in each panel. In Appendix B we present an analysis of the statistical properties of these relations.

Interestingly, in the pioneering work of Shlosman et al. (1989), a model is put forward in which the bar sweeps gas from the region of the main disc within Rbar, building a nuclear disc that is limited by the ILR radius. In their model, this radius is of the order of 0.1 Rbar, and the relation between rk and Rbar we show in Fig. 8 is close to that. This is an indication of the correctness of their model.

Figure 8 also shows that rk tends to be larger for lower values of bar ellipticity and higher values of Bar/T. Again this can be naturally understood considering the connection between nuclear discs and the extent of the bar x2 orbits. As x2 orbits cannot extend past the bar edges, more elongated bars will tend to have less extended x2 orbits, and therefore smaller nuclear discs, as we observe. As bars evolve, they capture stars from the disc, grow longer and become more massive, which increases Bar/T. In addition, with the buckling of the inner parts of the bar, and the formation of the box/peanut and barlens structure, the inner region of the bar becomes less elongated, decreasing the overall bar ellipticity. The less elongated inner region promotes the expansion of the x2 orbits, creating the trends between rk with bar ellipticity and Bar/T. Nevertheless, it is important to point out that in this study some of the observed trends are not strong correlations, and therefore more work is needed to confirm these conjectures (see Appendix B for an analysis of the statistical properties of these relations).

It is interesting to note that the range of bar radii in Fig. 8 spans a factor of ≈4, although in simulations bars grow by factors no larger than ≈2 (see e.g. Martinez-Valpuesta et al. 2006). This implies that the trends discussed here are not only the result of the way bar properties evolve, but also of the initial bar properties. On the other hand, the bar radius has been shown to be an increasing function of the bulge-to-total luminosity ratio (B/T) when normalised by R25.5 (Kim et al. 2014). Moreover, the function appears to be the same whether the photometric bulge component is a classical bulge or a nuclear disc. Kim et al. (2014) suggested that this relation is expected since more prominent classical bulges can absorb substantial angular momentum from the bar, which leads to more significant bar growth (see Athanassoula & Misiriotis 2002; Athanassoula 2003). On the other hand, Kim et al. argue that a similar relation is expected also in the case of nuclear discs, as longer bars will push more gas inwards to the central region, leading to more prominent nuclear discs (Athanassoula 1992). However, understanding why the relation is the same for classical bulges and nuclear discs remains a puzzle.

One also sees in Fig. 8 that the relation between rk and QB is such that for low values of QB nuclear discs show a wide range of sizes, but only small nuclear discs are seen when QB is large. This can in part be explained by the reduced extent of x2 orbits in more eccentric bars. Again, a similar result was found by Comerón et al. (2010) for the sizes of nuclear rings. Interestingly, we find a trend (albeit with some scatter) between rk and A2 that to the best of our knowledge has not been reported elsewhere. While A2 is a purely photometric parameter, rk is a purely kinematic one, and therefore this trend too is consistent with a picture in which the formation of the nuclear disc is connected to the bar.

Interestingly, the relation between rk and A2 does not show the same properties as the relation between rk and QB. While the connection between rk and QB is similar to the connection between rk and bar ellipticity, the correlation between rk and A2 resembles that between rk and Bar/T. This can be understood in the way A2 and QB are defined. Although both A2 and QB are measurements of bar strength, QB takes into account the overall axisymmetric galactic potential, whereas A2 accounts for it only at the radius where A2 is measured, typically close to the end of the bar. As described above, this is done with the normalisation of A2 by the value of A0 at that radius. Higher values of A2 are produced if the density of stars in the bar increases while it decreases in the region of the disc outside the bar but within the bar radius. Kim et al. (2016) showed that, as bars evolve, they capture stars in the disc region within the bar radius, reducing the density of stars in the region. Further, they showed that evolved bars thus tend to have elevated values of A2, a result that was confirmed in Buta (2017). This is corroborated by the finding from Díaz-García et al. (2016) that longer bars are stronger and that this correlation is particularly tight if bar strength is measured as A2, since bars are expected to grow longer as they evolve. Therefore, the trend between rk and A2 suggests that galaxies with more evolved bars tend to have larger nuclear discs. We already discussed above that more evolved bars should naturally host larger nuclear discs due to the increase in the extension of the bar x2 orbits as bars evolve. Furthermore, we speculate that this is also a result of the bar being able to induce gas inflows for longer, thus producing more massive nuclear discs. In fact, Cole et al. (2014) found that their simulations produce smaller nuclear discs if the inflow of gas is halted. Nevertheless, we stress that there is substantial scatter in the relation between rk and A2 seen in Fig. 8, and thus further studies are needed to confirm this relation (see Appendix B). In addition, more work is necessary to understand if the scatter is caused not only by difficulties in the measurements, but if other physical properties or processes play a role in the evolution of rk and A2.

Further work with larger samples would also be helpful to rule out that some of these correlations could result from underlying correlations, such as between bar size and disc size (see Erwin 2005; Gadotti 2011). However, we highlight that the anti-correlation we observe between rk and bar ellipticity cannot be explained by the underlying correlation between bar size and ellipticity, simply because these correlations have opposite signs. Since the anti-correlation between rk and bar ellipticity is predicted by the bar-driven model for the formation of nuclear discs, it suggests that the nuclear discs in the TIMER sample are indeed built by the bar.

We point out that while all galaxies plotted in Fig. 8 have nuclear discs, not all of them have nuclear rings, according to the evaluation by Comerón et al. (2010). These galaxies are NGC 4643, 4981, 4984, 5850, 7140 and NGC 7755. Nuclear rings are found at the edge of nuclear discs when both structures are present, close to or at the radius where v/σ peaks, that is, at rk. We verified that these galaxies are not systematically off in any of the relations shown in this section. Therefore, both nuclear rings and nuclear discs follow the same relations, suggesting that the formation of both structures stems from the same fundamental process. The results discussed above suggest that the main processes providing the gaseous content for the formation of both nuclear rings and nuclear discs are bar-driven, and the connection between nuclear rings and nuclear discs suggests that a single set of joint mechanisms is in charge of transforming the gas component into these stellar structures. In our accompanying paper (Bittner et al. 2020) we study this connection in more detail through the analysis of the stellar population properties of nuclear discs and nuclear rings.

Processes unrelated to bars have also been studied as possible drivers of the building of nuclear discs. Recent simulations have shown that nuclear discs can be produced at the late stage of gas-rich mergers (see Querejeta et al. 2015; Wang et al. 2015; Sauvaget et al. 2018). However, we note that most of these nuclear discs are substantially larger than the largest nuclear disc we find. Also, it is important to point out that even in this case bar-driven secular evolution processes can still play a major role in the building of the nuclear disc. The merger may simply serve as a mechanism to bring gas to the disc region within the bar radius, and from this point the bar is responsible to stream the gas further inwards. The merger can also trigger the bar instability if it was not spontaneously triggered before the interaction. Bars can be short-lived in major mergers, so it is unclear whether bars have played no role in the aforementioned simulations. While we present above suggestive evidence that in the TIMER sample of barred galaxies the nuclear discs are built by bars, such mergers can result in unbarred galaxies hosting nuclear discs, even if a bar played a major role but is later dissolved due to the merger.

On the other hand, gas-poor mergers are not expected to produce such regular and rapidly rotating structures. However, as we described in the Introduction, Eliche-Moral et al. (2011) argue that nuclear discs are built in their collisionless simulations when specific conditions are met. Nevertheless, an inspection of the corresponding kinematical maps (see their Figs. 12–15) reveals that these simulations do not reproduce the signatures of nuclear discs we discussed in Sect. 3.3.1. Specifically, our nuclear discs are characterised by conspicuous drops in σ and an anti-correlation between v and h3 (indicating near-circular orbits). In contrast, in the aforementioned simulations, the central region shows an inward increase in σ, and, in most cases, a correlation between v and h3. Therefore, we conclude that, to date, there is no evidence that the nuclear discs in TIMER could have been built in gas-poor mergers.

This is not to say that the building of nuclear discs through bar-driven processes is perfectly well understood in simulations. In fact, the bar-built nuclear disc in the model of Cole et al. (2014) is too large compared to its bar and does not fit in the observed relation we present in Fig. 8. While the observed ratio between nuclear disc size and bar size is of the order of 10%, the nuclear disc size in the model of Cole et al. (2014) is about 30% of the bar size. More theory work is necessary, and it would be particularly helpful to understand how simulations can reproduce the observed trends in Fig. 8.

6. General discussion and conclusions

The results presented in Sects. 3.3.1 and 5 show a good agreement with the bar-driven secular evolution scenario for the building of nuclear stellar components in disc galaxies, at least for massive galaxies with conspicuous bars and nuclear structures. The nuclear discs found with the TIMER data have all the kinematic properties expected in a stellar structure built via the collapse of molecular clouds brought to the central region along the leading edges of the bar and put into near-circular orbits in the region where the bar x2 orbits dominate: Large rotational support, low velocity dispersion and an anti-correlation between v and h3 consistent with near-circular orbits (see Figs. 2 and A.1). Further, these nuclear discs are clearly an additional component on top of the original main galaxy disc. This can be seen from the elevated values of h4, resulting from the fact that the stars in the nuclear discs have orbits closer to circular than those in the main discs at the radii where the nuclear disc dominates (Figs. 2, 3 and A.1). These properties are also consistent with a formation scenario with gas accretion that is unrelated to bars. However, Fig. 8 shows how the size of the nuclear disc depends on the bar radius, ellipticity, bar-to-total ratio and bar strength, which is qualitatively understood in the theoretical framework of bar evolution and the impact of bars on the gas component. These connections between the nuclear disc size and bar properties are only expected in a scenario where the accretion of gas that builds the nuclear disc is due to the bar. Figure 8 also suggests that more evolved bars tend to build larger nuclear discs, which may help putting constraints on the ages of bars, although more theoretical work on the formation and evolution of nuclear discs is necessary. Finally, in our accompanying paper (Bittner et al. 2020), we present stellar population properties of the TIMER nuclear discs, and discuss how these properties too are consistent with the picture in which such nuclear discs are built via bar-driven processes.

As mentioned above, it is still unclear how rare are nuclear discs in unbarred galaxies (but see Comerón et al. 2010). Weak bars and oval distortions in the disc can be difficult to identify morphologically but may as well produce nuclear discs in the same way as prominent bars. Nevertheless, a similar study as we present here but with a sample of unbarred galaxies showing nuclear discs would be very beneficial. A comparison between the physical properties of nuclear discs in barred and unbarred galaxies (such as size and angular momentum) would shed light on their formation process.

As we have seen in Sect. 4, in photometric decompositions of galaxy images, one often finds an exponential or near-exponential central component, that is, a photometric bulge with low Sérsic index (n ≲ 2), which is thought to originate from bar-driven secular processes. In agreement with this expectation, the results presented in Sect. 4 show that the nuclear discs we find spectroscopically in TIMER, through an assessment of the stellar kinematics, are indeed often recognised as (near-)exponential bulges in the photometric studies we considered. This is encouraging, despite the caveats discussed above on using the Sérsic index alone to separate nuclear discs from classical bulges. However, we stress again that the TIMER sample was built to include galaxies for which a visual assessment alone is already capable of identifying nuclear components that appear to have disc-like properties. This is a result of the relatively high physical spatial resolution of the images employed. Arguably, the recovery of nuclear discs via photometry is more prone to errors when the spatial resolution is not suitable, for example for more distant galaxies.

Even in these optimal circumstances, Sect. 4 shows that some difficulties are encountered, for example when performing photometric decompositions of galaxies with composite bulges. A small classical bulge within a nuclear disc can dominate the light emission in such a way that the presence of a bar-built nuclear component goes unnoticed in the photometry. This problem is of course exacerbated if the spatial resolution is not high enough to separate the two components, and/or a single model is used to fit the central region, erroneously encompassing all structures therein. It is as yet not clear how often disc galaxies host small classical bulges embedded in nuclear discs so the severity of this problem is still unknown. However, our results show that composite bulges may produce values for the Sérsic index between 2 and 3, suggesting that pure classical bulges can be identified only by putting a more stringent threshold at n >  3 (especially considering the corresponding uncertainties). Nevertheless, this is to be considered with caution and, particularly, in a statistical sense (see Méndez-Abreu et al. 2018).

The TIMER maps shown in Figs. 2, 3 and A.1 show several kinematic signatures corroborating theoretical work on the dynamical properties of bars and inner bars. This is also the case for box/peanuts (Sect. 3.3.3). This allows us to put a lower limit in the fraction of massive barred galaxies with box/peanuts at 62%, in broad agreement with previous results (with the caveat that our sample selection favours conspicuous bars and nuclear components). The case of NGC 5728 is remarkable, in that it shows an agreement between kinematic and photometric considerations on the properties of box/peanuts in inclined disc galaxies. Concerning barlenses (Sect. 3.3.4), we find evidence corroborating previous studies that find that barlenses are simply the face-on projection of box/peanuts, that is, of the inner parts of bars, since we see in the barlens region kinematic signatures of both a bar and a box/peanut, as expected from numerical simulations. This is the first time kinematic evidence is presented to support this picture. Barlenses are often difficult to identify photometrically and have morphologies similar to classical bulges, and therefore identifying them through their kinematic properties can help putting more accurate constraints on the impact of mergers in the evolution of disc galaxies. Altogether, these results show the power of spatially resolved kinematics in producing a straightforward understanding of the physical properties and nature of stellar structures in the inner regions of galaxies, and suggest that high-quality integral-field spectroscopy data is necessary to accurately decompose these complex inner regions.


1

Observations are still lacking for three galaxies.

2

Available at ascl.net/1907.025

4

In the presence of an inner bar one expects to see signatures of elongated orbits where the inner bar dominates.

5

Here we assume that both discs are in the same plane, which is justified by the alignment of their kinematic axes seen in Figs. 2, 3 and A.1.

6

The x1 orbits are eccentric and parallel to the bar major axis, and are present throughout the bar. They are thus considered the “backbone” of bars. The x2 orbits are less eccentric and perpendicular to the bar major axis, and are only found in the bar inner region.

7

We point out that, with all things being equal, more eccentric bars will have smaller nuclear discs, since the extent of the x2 orbits is then shorter (see Athanassoula 1992). A nuclear disc will not grow in a bar that gets longer but keeps its semi-minor axis constant. In fact, Fig. 8 shows that more eccentric bars tend to have smaller nuclear discs.

8

We note, however, that elevated values of h4 can be produced by radial orbits in non-rotating or slowly rotating systems (see Fig. 2 in van der Marel & Franx 1993).

Acknowledgments

We thank the anonymous referee for timely and helpful reports. Based on observations collected at the European Southern Observatory under ESO programmes 097.B-0640(A), 095.B-0532(A), 094.B-0321(A) and 060.A-9313(A). J. F-B, AdLC, and PSB acknowledge support through the RAVET project by the grants AYA2016-77237-C3-1-P, AYA2016-77237-C3-2-P and PID2019-107427GB-C31 from the Spanish Ministry of Science, Innovation and Universities (MCIU). J.F.-B. and AdLC acknowledge support through the IAC project TRACES which is partially supported through the state budget and the regional budget of the Consejería de Economía, Industria, Comercio y Conocimiento of the Canary Islands Autonomous Community. JMA acknowledges support from the Spanish Ministry of Economy and Competitiveness (MINECO) by grant AYA2017-83204-P. TK was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (No. 2019R1A6A3A01092024). The Science, Technology and Facilities Council is acknowledged by JN for support through the Consolidated Grant Cosmology and Astrophysics at Portsmouth, ST/S000550/1. GvdV acknowledges funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme under grant agreement No 724857 (Consolidator Grant ArcheoDyn). This research has made use of NASA’s Astrophysics Data System Bibliographic Services. We acknowledge the usage of the HyperLeda database (http://leda.univ-lyon1.fr). This research has also made use of the NASA/IPAC Extragalactic Database (NED), which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration.

References

  1. Allard, E. L., Knapen, J. H., Peletier, R. F., & Sarzi, M. 2006, MNRAS, 371, 1087 [NASA ADS] [CrossRef] [Google Scholar]
  2. Anscombe, F. J., & Glynn, W. J. 1983, Biometrika, 70, 227 [Google Scholar]
  3. Athanassoula, E. 1992, MNRAS, 259, 345 [NASA ADS] [CrossRef] [Google Scholar]
  4. Athanassoula, E. 2003, MNRAS, 341, 1179 [NASA ADS] [CrossRef] [Google Scholar]
  5. Athanassoula, E. 2005, MNRAS, 358, 1477 [NASA ADS] [CrossRef] [Google Scholar]
  6. Athanassoula, E., & Beaton, R. L. 2006, MNRAS, 370, 1499 [NASA ADS] [CrossRef] [Google Scholar]
  7. Athanassoula, E., & Misiriotis, A. 2002, MNRAS, 330, 35 [NASA ADS] [CrossRef] [Google Scholar]
  8. Athanassoula, E., Machado, R. E. G., & Rodionov, S. A. 2013, MNRAS, 429, 1949 [NASA ADS] [CrossRef] [Google Scholar]
  9. Athanassoula, E., Laurikainen, E., Salo, H., & Bosma, A. 2015, MNRAS, 454, 3843 [NASA ADS] [CrossRef] [Google Scholar]
  10. Bacon, R., Copin, Y., Monnet, G., et al. 2001, MNRAS, 326, 23 [NASA ADS] [CrossRef] [Google Scholar]
  11. Bacon, R., Accardo, M., & Adjali, L. 2010, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, 7735, 8 [Google Scholar]
  12. Bacon, R., Conseil, S., Mary, D., et al. 2017, A&A, 608, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Beaton, R. L., Majewski, S. R., Guhathakurta, P., et al. 2007, ApJ, 658, L91 [NASA ADS] [CrossRef] [Google Scholar]
  14. Bender, R., Saglia, R. P., & Gerhard, O. E. 1994, MNRAS, 269, 785 [NASA ADS] [CrossRef] [Google Scholar]
  15. Binney, J. 2005, MNRAS, 363, 937 [NASA ADS] [CrossRef] [Google Scholar]
  16. Binney, J., Gerhard, O. E., Stark, A. A., Bally, J., & Uchida, K. I. 1991, MNRAS, 252, 210 [NASA ADS] [CrossRef] [Google Scholar]
  17. Bittner, A., Falcón-Barroso, J., Nedelchev, B., et al. 2019, A&A, 628, A117 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Bittner, A., Sánchez-Blázquez, P., Gadotti, D. A., et al. 2020, A&A, in press, https://doi.org/10.1051/0004-6361/202038450 [Google Scholar]
  19. Bosma, A., Gadotti, D. A., de Blok, W. J. G., & Athanassoula, E. 2010, in Galaxies in Isolation: Exploring Nature Versus Nurture, eds. L. Verdes-Montenegro, A. Del Olmo, & J. Sulentic, ASP Conf. Ser., 421, 53 [Google Scholar]
  20. Bureau, M., & Athanassoula, E. 2005, ApJ, 626, 159 [NASA ADS] [CrossRef] [Google Scholar]
  21. Bureau, M., & Freeman, K. C. 1999, AJ, 118, 126 [NASA ADS] [CrossRef] [Google Scholar]
  22. Bureau, M., Aronica, G., Athanassoula, E., et al. 2006, MNRAS, 370, 753 [NASA ADS] [CrossRef] [Google Scholar]
  23. Buta, R. J. 2017, MNRAS, 470, 3819 [NASA ADS] [CrossRef] [Google Scholar]
  24. Buta, R. J., Sheth, K., Athanassoula, E., et al. 2015, ApJS, 217, 32 [NASA ADS] [CrossRef] [Google Scholar]
  25. Cappellari, M. 2017, MNRAS, 466, 798 [NASA ADS] [CrossRef] [Google Scholar]
  26. Cappellari, M., & Copin, Y. 2003, MNRAS, 342, 345 [NASA ADS] [CrossRef] [Google Scholar]
  27. Cappellari, M., & Emsellem, E. 2004, PASP, 116, 138 [NASA ADS] [CrossRef] [Google Scholar]
  28. Cappellari, M., Emsellem, E., Bacon, R., et al. 2007, MNRAS, 379, 418 [NASA ADS] [CrossRef] [Google Scholar]
  29. Cappellari, M., Emsellem, E., Krajnović, D., et al. 2011, MNRAS, 413, 813 [NASA ADS] [CrossRef] [Google Scholar]
  30. Chapon, D., Mayer, L., & Teyssier, R. 2013, MNRAS, 429, 3114 [CrossRef] [Google Scholar]
  31. Chung, A., & Bureau, M. 2004, AJ, 127, 3192 [NASA ADS] [CrossRef] [Google Scholar]
  32. Coelho, P., & Gadotti, D. A. 2011, ApJ, 743, L13 [NASA ADS] [CrossRef] [Google Scholar]
  33. Cole, D. R., Debattista, V. P., Erwin, P., Earp, S. W. F., & Roškar, R. 2014, MNRAS, 445, 3352 [NASA ADS] [CrossRef] [Google Scholar]
  34. Combes, F., & Sanders, R. H. 1981, A&A, 96, 164 [NASA ADS] [Google Scholar]
  35. Combes, F., Debbasch, F., Friedli, D., & Pfenniger, D. 1990, A&A, 233, 82 [NASA ADS] [Google Scholar]
  36. Comerón, S., Knapen, J. H., Beckman, J. E., et al. 2010, MNRAS, 402, 2462 [NASA ADS] [CrossRef] [Google Scholar]
  37. Costantin, L., Corsini, E. M., Méndez-Abreu, J., et al. 2018a, MNRAS, 481, 3623 [NASA ADS] [CrossRef] [Google Scholar]
  38. Costantin, L., Méndez-Abreu, J., Corsini, E. M., et al. 2018b, A&A, 609, A132 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. D’Agostino, R. B., Belanger, A. J., & D’Agostino, R. B., Jr. 1990, Am. Stat., 44, 316 [Google Scholar]
  40. Davies, R. L., Burstein, D., Dressler, A., et al. 1987, ApJS, 64, 581 [NASA ADS] [CrossRef] [Google Scholar]
  41. Debattista, V. P., Carollo, C. M., Mayer, L., & Moore, B. 2005, ApJ, 628, 678 [NASA ADS] [CrossRef] [Google Scholar]
  42. de Jong, R. S. 1996, A&A, 313, 377 [NASA ADS] [Google Scholar]
  43. de Lorenzo-Cáceres, A., Falcón-Barroso, J., Vazdekis, A., & Martínez-Valpuesta, I. 2008, ApJ, 684, L83 [NASA ADS] [CrossRef] [Google Scholar]
  44. de Lorenzo-Cáceres, A., Vazdekis, A., Aguerri, J. A. L., Corsini, E. M., & Debattista, V. P. 2012, MNRAS, 420, 1092 [NASA ADS] [CrossRef] [Google Scholar]
  45. de Lorenzo-Cáceres, A., Falcón-Barroso, J., & Vazdekis, A. 2013, MNRAS, 431, 2397 [NASA ADS] [CrossRef] [Google Scholar]
  46. de Lorenzo-Cáceres, A., Sánchez-Blázquez, P., Méndez-Abreu, J., et al. 2019, MNRAS, 484, 5296 [CrossRef] [Google Scholar]
  47. de Souza, R. E., & Dos Anjos, S. 1987, A&AS, 70, 465 [NASA ADS] [Google Scholar]
  48. de Souza, R. E., Gadotti, D. A., & dos Anjos, S. 2004, ApJS, 153, 411 [NASA ADS] [CrossRef] [Google Scholar]
  49. de Zeeuw, P. T., Bureau, M., Emsellem, E., et al. 2002, MNRAS, 329, 513 [NASA ADS] [CrossRef] [Google Scholar]
  50. Díaz-García, S., Salo, H., Laurikainen, E., & Herrera-Endoqui, M. 2016, A&A, 587, A160 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  51. Donohoe-Keyes, C. E., Martig, M., James, P. A., & Kraljic, K. 2019, MNRAS, 489, 4992 [CrossRef] [Google Scholar]
  52. Dopita, M., Hart, J., McGregor, P., et al. 2007, Ap&SS, 310, 255 [NASA ADS] [CrossRef] [Google Scholar]
  53. Dopita, M., Rhee, J., Farage, C., et al. 2010, Ap&SS, 327, 245 [NASA ADS] [CrossRef] [Google Scholar]
  54. Du, M., Shen, J., & Debattista, V. P. 2015, ApJ, 804, 139 [CrossRef] [Google Scholar]
  55. Dumas, G., Mundell, C. G., Emsellem, E., & Nagar, N. M. 2007, MNRAS, 379, 1249 [NASA ADS] [CrossRef] [Google Scholar]
  56. Eliche-Moral, M. C., González-García, A. C., Balcells, M., et al. 2011, A&A, 533, A104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  57. Ellison, S. L., Nair, P., Patton, D. R., et al. 2011, MNRAS, 416, 2182 [NASA ADS] [CrossRef] [Google Scholar]
  58. Emsellem, E., Cappellari, M., Krajnović, D., et al. 2007, MNRAS, 379, 401 [Google Scholar]
  59. Emsellem, E., Cappellari, M., Krajnović, D., et al. 2011, MNRAS, 414, 888 [Google Scholar]
  60. Erwin, P. 2004, A&A, 415, 941 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Erwin, P. 2005, MNRAS, 364, 283 [NASA ADS] [CrossRef] [Google Scholar]
  62. Erwin, P., & Debattista, V. P. 2013, MNRAS, 431, 3060 [NASA ADS] [CrossRef] [Google Scholar]
  63. Erwin, P., & Debattista, V. P. 2017, MNRAS, 468, 2058 [NASA ADS] [CrossRef] [Google Scholar]
  64. Erwin, P., Saglia, R. P., Fabricius, M., et al. 2015, MNRAS, 446, 4039 [NASA ADS] [CrossRef] [Google Scholar]
  65. Eskridge, P. B., Frogel, J. A., Pogge, R. W., et al. 2000, AJ, 119, 536 [NASA ADS] [CrossRef] [Google Scholar]
  66. Falcón-Barroso, J., Peletier, R. F., Emsellem, E., et al. 2004, MNRAS, 350, 35 [NASA ADS] [CrossRef] [Google Scholar]
  67. Falcón-Barroso, J., Bacon, R., Bureau, M., et al. 2006, MNRAS, 369, 529 [Google Scholar]
  68. Fisher, D. B., & Drory, N. 2016, in Galactic Bulges, eds. E. Laurikainen, R. Peletier, & D. Gadotti, Astrophys. Space Sci. Lib., 418, 41 [NASA ADS] [CrossRef] [Google Scholar]
  69. Fragkoudi, F., Athanassoula, E., & Bosma, A. 2016, MNRAS, 462, L41 [NASA ADS] [CrossRef] [Google Scholar]
  70. Fragkoudi, F., Di Matteo, P., Haywood, M., et al. 2017, A&A, 606, A47 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  71. Fragkoudi, F., Di Matteo, P., Haywood, M., et al. 2018, A&A, 616, A180 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  72. Fragkoudi, F., Grand, R. J. J., Pakmor, R., et al. 2020, MNRAS, 494, 5936 [CrossRef] [Google Scholar]
  73. Gadotti, D. A. 2008, MNRAS, 384, 420 [NASA ADS] [CrossRef] [Google Scholar]
  74. Gadotti, D. A. 2009, MNRAS, 393, 1531 [NASA ADS] [CrossRef] [Google Scholar]
  75. Gadotti, D. A. 2011, MNRAS, 415, 3308 [NASA ADS] [CrossRef] [Google Scholar]
  76. Gadotti, D. A., & dos Anjos, S. 2001, AJ, 122, 1298 [NASA ADS] [CrossRef] [Google Scholar]
  77. Gadotti, D. A., Athanassoula, E., Carrasco, L., et al. 2007, MNRAS, 381, 943 [NASA ADS] [CrossRef] [Google Scholar]
  78. Gadotti, D. A., Baes, M., & Falony, S. 2010, MNRAS, 403, 2053 [NASA ADS] [CrossRef] [Google Scholar]
  79. 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]
  80. Gadotti, D. A., Sánchez-Blázquez, P., Falcón-Barroso, J., et al. 2019, MNRAS, 482, 506 [Google Scholar]
  81. Gonzalez, O. A., & Gadotti, D. 2016, in Galactic Bulges, eds. E. Laurikainen, R. Peletier, & D. Gadotti, Astrophys. Space Sci. Lib., 418, 199 [NASA ADS] [CrossRef] [Google Scholar]
  82. Iannuzzi, F., & Athanassoula, E. 2015, MNRAS, 450, 2514 [NASA ADS] [CrossRef] [Google Scholar]
  83. Isobe, T., Feigelson, E. D., Akritas, M. G., & Babu, G. J. 1990, ApJ, 364, 104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  84. Kim, W.-T., Seo, W.-Y., Stone, J. M., Yoon, D., & Teuben, P. J. 2012, ApJ, 747, 60 [NASA ADS] [CrossRef] [Google Scholar]
  85. Kim, T., Gadotti, D. A., Sheth, K., et al. 2014, ApJ, 782, 64 [NASA ADS] [CrossRef] [Google Scholar]
  86. Kim, T., Gadotti, D. A., Athanassoula, E., et al. 2016, MNRAS, 462, 3430 [NASA ADS] [CrossRef] [Google Scholar]
  87. Knapen, J. H. 2005, A&A, 429, 141 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  88. Knapen, J. H., Beckman, J. E., Heller, C. H., Shlosman, I., & de Jong, R. S. 1995, ApJ, 454, 623 [NASA ADS] [CrossRef] [Google Scholar]
  89. Knapen, J. H., Pérez-Ramírez, D., & Laine, S. 2002, MNRAS, 337, 808 [NASA ADS] [CrossRef] [Google Scholar]
  90. Kormendy, J. 1982, ApJ, 257, 75 [NASA ADS] [CrossRef] [Google Scholar]
  91. Kormendy, J., & Kennicutt, R. C., Jr. 2004, ARA&A, 42, 603 [NASA ADS] [CrossRef] [Google Scholar]
  92. Kraljic, K., Bournaud, F., & Martig, M. 2012, ApJ, 757, 60 [NASA ADS] [CrossRef] [Google Scholar]
  93. Kruk, S. J., Lintott, C. J., Bamford, S. P., et al. 2018, MNRAS, 473, 4731 [NASA ADS] [CrossRef] [Google Scholar]
  94. Kruk, S. J., Erwin, P., Debattista, V. P., & Lintott, C. 2019, MNRAS, 490, 4721 [NASA ADS] [CrossRef] [Google Scholar]
  95. Kuijken, K., & Merrifield, M. R. 1995, ApJ, 443, L13 [NASA ADS] [CrossRef] [Google Scholar]
  96. Laurikainen, E., & Salo, H. 2016, in Galactic Bulges, eds. E. Laurikainen, R. Peletier, & D. Gadotti, Astrophys. Space Sci. Lib., 418, 77 [NASA ADS] [CrossRef] [Google Scholar]
  97. Laurikainen, E., Salo, H., & Buta, R. 2005, MNRAS, 362, 1319 [NASA ADS] [CrossRef] [Google Scholar]
  98. Laurikainen, E., Salo, H., Buta, R., & Knapen, J. H. 2007, MNRAS, 381, 401 [NASA ADS] [CrossRef] [Google Scholar]
  99. Laurikainen, E., Salo, H., Buta, R., & Knapen, J. H. 2011, MNRAS, 418, 1452 [NASA ADS] [CrossRef] [Google Scholar]
  100. Leaman, R., Fragkoudi, F., Querejeta, M., et al. 2019, MNRAS, 488, 3904 [CrossRef] [Google Scholar]
  101. Li, Z., Shen, J., & Kim, W.-T. 2015, ApJ, 806, 150 [NASA ADS] [CrossRef] [Google Scholar]
  102. Li, Z.-Y., Shen, J., Bureau, M., et al. 2018, ApJ, 854, 65 [CrossRef] [Google Scholar]
  103. Łokas, E. L. 2019, A&A, 629, A52 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  104. Louis, P. D., & Gerhard, O. E. 1988, MNRAS, 233, 337 [NASA ADS] [CrossRef] [Google Scholar]
  105. Lütticke, R., Dettmar, R.-J., & Pohlen, M. 2000, A&AS, 145, 405 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  106. Martinez-Valpuesta, I., Shlosman, I., & Heller, C. 2006, ApJ, 637, 214 [NASA ADS] [CrossRef] [Google Scholar]
  107. Mayer, L., Kazantzidis, S., & Escala, A. 2008, Mem. Soc. Astron. It., 79, 1284 [Google Scholar]
  108. Melvin, T., Masters, K., Lintott, C., et al. 2014, MNRAS, 438, 2882 [NASA ADS] [CrossRef] [Google Scholar]
  109. Méndez-Abreu, J., Corsini, E. M., Debattista, V. P., et al. 2008, ApJ, 679, L73 [NASA ADS] [CrossRef] [Google Scholar]
  110. Méndez-Abreu, J., Debattista, V. P., Corsini, E. M., & Aguerri, J. A. L. 2014, A&A, 572, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  111. Méndez-Abreu, J., Aguerri, J. A. L., Falcón-Barroso, J., et al. 2018, MNRAS, 474, 1307 [NASA ADS] [Google Scholar]
  112. Méndez-Abreu, J., de Lorenzo-Cáceres, A., Gadotti, D. A., et al. 2019, MNRAS, 482, L118 [NASA ADS] [CrossRef] [Google Scholar]
  113. Menéndez-Delmestre, K., Sheth, K., Schinnerer, E., Jarrett, T. H., & Scoville, N. Z. 2007, ApJ, 657, 790 [NASA ADS] [CrossRef] [Google Scholar]
  114. Merrifield, M. R., & Kuijken, K. 1999, A&A, 345, L47 [NASA ADS] [Google Scholar]
  115. Moffat, A. F. J. 1969, A&A, 3, 455 [NASA ADS] [Google Scholar]
  116. Muñoz-Mateos, J. C., Sheth, K., Regan, M., et al. 2015, ApJS, 219, 3 [Google Scholar]
  117. Neumann, J., Wisotzki, L., Choudhury, O. S., et al. 2017, A&A, 604, A30 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  118. Neumann, J., Fragkoudi, F., Pérez, I., et al. 2020, A&A, 637, A56 [CrossRef] [EDP Sciences] [Google Scholar]
  119. Pastrav, B. A., Popescu, C. C., Tuffs, R. J., & Sansom, A. E. 2013, A&A, 557, A137 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  120. Peng, C. Y., Ho, L. C., Impey, C. D., & Rix, H.-W. 2002, AJ, 124, 266 [Google Scholar]
  121. Peng, C. Y., Ho, L. C., Impey, C. D., & Rix, H.-W. 2010, AJ, 139, 2097 [NASA ADS] [CrossRef] [Google Scholar]
  122. Pérez, I., Martínez-Valpuesta, I., Ruiz-Lara, T., et al. 2017, MNRAS, 470, L122 [NASA ADS] [CrossRef] [Google Scholar]
  123. Piner, B. G., Stone, J. M., & Teuben, P. J. 1995, ApJ, 449, 508 [NASA ADS] [CrossRef] [Google Scholar]
  124. Prendergast, K. H. 1983, in Internal Kinematics and Dynamics of Galaxies, ed. E. Athanassoula, IAU Symp., 100, 215 [CrossRef] [Google Scholar]
  125. Querejeta, M., Eliche-Moral, M. C., Tapia, T., et al. 2015, A&A, 573, A78 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  126. Romeo, A. B., & Fathi, K. 2015, MNRAS, 451, 3107 [NASA ADS] [CrossRef] [Google Scholar]
  127. Romeo, A. B., & Fathi, K. 2016, MNRAS, 460, 2360 [NASA ADS] [CrossRef] [Google Scholar]
  128. Rousseeuw, P. J. 1984, J. Am. Stat. Assoc., 79, 871 [CrossRef] [Google Scholar]
  129. Rousseeuw, J. P., & Leroy, A. M. 1987, Robust Regression and Outlier Detection (New York: John Wiley) [CrossRef] [Google Scholar]
  130. Sakamoto, K., Okumura, S. K., Ishizuki, S., & Scoville, N. Z. 1999, ApJ, 525, 691 [NASA ADS] [CrossRef] [Google Scholar]
  131. Salo, H., Laurikainen, E., Laine, J., et al. 2015, ApJS, 219, 4 [NASA ADS] [CrossRef] [Google Scholar]
  132. Sanders, R. H., & Tubbs, A. D. 1980, ApJ, 235, 803 [NASA ADS] [CrossRef] [Google Scholar]
  133. Sauvaget, T., Hammer, F., Puech, M., et al. 2018, MNRAS, 473, 2521 [NASA ADS] [CrossRef] [Google Scholar]
  134. Seidel, M. K., Falcón-Barroso, J., Martínez-Valpuesta, I., et al. 2015a, MNRAS, 451, 936 [NASA ADS] [CrossRef] [Google Scholar]
  135. Seidel, M. K., Cacho, R., Ruiz-Lara, T., et al. 2015b, MNRAS, 446, 2837 [NASA ADS] [CrossRef] [Google Scholar]
  136. Seo, W.-Y., Kim, W.-T., Kwak, S., et al. 2019, ApJ, 872, 5 [CrossRef] [Google Scholar]
  137. Sérsic, J. L. 1968, Atlas de Galaxias Australes (Cordoba, Argentina: Observatorio Astronomico) [Google Scholar]
  138. Shaw, M. A. 1987, MNRAS, 229, 691 [NASA ADS] [CrossRef] [Google Scholar]
  139. Sheth, K., Vogel, S. N., Regan, M. W., Thornley, M. D., & Teuben, P. J. 2005, ApJ, 632, 217 [NASA ADS] [CrossRef] [Google Scholar]
  140. Sheth, K., Elmegreen, D. M., Elmegreen, B. G., et al. 2008, ApJ, 675, 1141 [NASA ADS] [CrossRef] [Google Scholar]
  141. Sheth, K., Regan, M., Hinz, J. L., et al. 2010, PASP, 122, 1397 [NASA ADS] [CrossRef] [Google Scholar]
  142. Sheth, K., Melbourne, J., Elmegreen, D. M., et al. 2012, ApJ, 758, 136 [NASA ADS] [CrossRef] [Google Scholar]
  143. Shimizu, T. T., Davies, R. I., Lutz, D., et al. 2019, MNRAS, 490, 5860 [NASA ADS] [CrossRef] [Google Scholar]
  144. Shlosman, I., Frank, J., & Begelman, M. C. 1989, Nature, 338, 45 [NASA ADS] [CrossRef] [Google Scholar]
  145. Simkin, S. M., Su, H. J., & Schwarz, M. P. 1980, ApJ, 237, 404 [NASA ADS] [CrossRef] [Google Scholar]
  146. Simmons, B. D., Melvin, T., Lintott, C., et al. 2014, MNRAS, 445, 3466 [NASA ADS] [CrossRef] [Google Scholar]
  147. Sormani, M. C., Sobacchi, E., Fragkoudi, F., et al. 2018, MNRAS, 481, 2 [NASA ADS] [CrossRef] [Google Scholar]
  148. van der Marel, R. P., & Franx, M. 1993, ApJ, 407, 525 [NASA ADS] [CrossRef] [Google Scholar]
  149. Venturi, G., Nardini, E., Marconi, A., et al. 2018, A&A, 619, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  150. Wang, J., Hammer, F., Puech, M., Yang, Y., & Flores, H. 2015, MNRAS, 452, 3551 [NASA ADS] [CrossRef] [Google Scholar]
  151. Weilbacher, P. M., Streicher, O., Urrutia, T., et al. 2012, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, 8451, 84510B [Google Scholar]
  152. Wozniak, H. 2007, A&A, 465, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  153. Wozniak, H. 2015, A&A, 575, A7 [CrossRef] [EDP Sciences] [Google Scholar]

Appendix A: Further kinematic maps and discussion on exceptional cases

In this appendix we present in Fig. A.1 the kinematic maps for the remainder of the sample as in Fig. 2. We also elaborate on exceptions to the general trends discussed in Sects. 3.3.13.3.3, where we presented kinematic signatures of nuclear discs, bars and box/peanuts, respectively.

thumbnail Fig. A.1.

Same as Fig. 2 but for the remainder of our sample.

thumbnail Fig. A.1.

continued.

thumbnail Fig. A.1.

continued.

A.1. Nuclear discs

As discussed in Sect. 3.3.1, we find clear kinematic signatures of the presence of nuclear discs in most galaxies in this study. Apart from NGC 1291 – where the fact that the galaxy is very close to face-on prevents us from performing this assessment meaningfully – there are only two possible exceptions: NGC 1365 and NGC 6902. In NGC 1365, the maps of v and v/σ do not show a discontinuous behaviour and the clear presence of an additional rapidly rotating nuclear component, even though a nuclear star-forming ring is evident in S4G images, as well as in our MUSE reconstructed images (see Paper I). The galaxy is quite inclined (52°) and our MUSE field is heavily affected by dust, which complicates the interpretation of the kinematic maps. Moreover, it hosts an active galactic nucleus with associated large scale outflow and shock-induced emission lines. Nevertheless, it shows elevated values of h4 in the central region, covering about half of the MUSE field8. In addition, it shows a peculiar spiral-shaped region of elevated absolute values of h3 anti-correlated with v. Furthermore, an extended region of low σ is seen at an angle with the bar, which is consistent with a nuclear, kinematically cold stellar structure seen in projection. It is therefore plausible that the MUSE field is dominated by the nuclear component, which precludes one from seeing the main galaxy disc. The spiral pattern in the h3 map is likely a result from the overly strong dust lanes seen along the leading edges of the bar.

On the other hand, NGC 6902 is not heavily affected by dust and has a mild inclination angle (37°), but it appears to host a rather modest nuclear stellar structure. The v and h3 maps show no clearly distinguished and rapidly rotating nuclear component, although this is better seen in the v/σ map. Nevertheless, the central region of the MUSE field, where the nuclear component resides, is dominated by high values of σ rather than low values. In addition, the h4 map shows a ring-shaped region of elevated values around the central region, not so much in the region dominated by the nuclear component. Therefore, the interpretation of the kinematic maps for NGC 6902 is not as straightforward as for most of the TIMER sample. We note that NGC 6902 has the weakest bar of the TIMER sample, in the morphological classification of Buta et al. (2015).

A.2. Bars

We do not find the v − h3 correlation foreseen in bars in six barred galaxies: NGC 1097, 1291, 1365, 4371, 5236, and NGC 6902. Apart from the latter, all these galaxies show prominent bars, but in most of these cases it is likely that the TIMER data do not cover enough of the bar to show this signature. In fact, NGC 5236, 1291, 1097 and NGC 1365 are, in this order, the largest projected bars in the TIMER sample, and our MUSE fields cover less than a third of the bar in these cases. The bar of NGC 4371 is well covered by our MUSE pointing but this is the most inclined galaxy in this study and the bar is close to perpendicular to the line of nodes, and thus projection effects may be preventing us from seeing the v − h3 correlation. As mentioned above, NGC 6902 is the weakest bar in the sample and again appears as an exception.

A.3. Box/peanuts

Since our MUSE fields seem to not go much further than any box/peanut vertices, the h4 drops expected for box/peanuts are not clearly present in NGC 1097, 1365, 1433, and NGC 3351. We also do not find this box/peanut kinematic signature in NGC 4371, but in this case this is possibly due to projection effects. The bar is seen inclined around an axis close to its minor axis, and the galaxy inclination angle is relatively large (at 59°, this galaxy is the most inclined galaxy in this study).

In NGC 5236 the signature is not so clear. One sees regions of low h4 on average along the bar major axis but the MUSE field seems to not be large enough to show where h4 would rise again. NGC 5248 shows h4 minima but it is not clear how they are associated to the rather weak bar, which in addition becomes harder to distinguish due to projection effects.

NGC 6902 is again an interesting case. As discussed above, it is the weakest bar in TIMER and there are no clear kinematic signatures of the presence of a nuclear disc. Its weak bar also does not show in the kinematic maps. Tentatively, there is weak evidence for a box/peanut, since inside the ring of elevated values of h4 mentioned above one sees dips to values close to zero along the bar major axis.

Finally, NGC 4643 shows a curious behaviour. Along the bar major axis outwards, h4 drops to a minimum after the central region of elevated values, further out rises again reaching values close to zero, and then drops again to very low values towards the end of the bar. This behaviour would be produced if the galaxy has an inner bar with its own box/peanut (as does NGC 1291) plus the box/peanut of the primary bar. However, this cannot be the case here as the inner bar should be within the nuclear disc radius, but the first pair of h4 minima is outside that radius. In addition, the two pairs of h4 minima are very well aligned, suggesting that the properties they indicate concern only the (main) bar. To understand the presence of two pairs of h4 minima in NGC 4643 is beyond the scope of this paper.

Appendix B: Statistical analysis of the trends between rk and selected bar properties (Fig. 8)

In order to better understand the statistical significance of the relations presented in Fig. 8 and discussed in Sect. 5, we present in Table B.1 the values of the slope and correlation coefficient as determined through different methods for each relation. We employed the SIXLIN IDL implementation of the formulae provided by Isobe et al. (1990) to calculate the slopes through ordinary least squares regression (OLS) and orthogonal regression (OR). OLS can be calculated with the Y variable (in this case, rk) as the dependent variable (Y|X) or as the independent variable (X|Y). The slope of the bisector of the region between the two OLS lines is also shown.

Table B.1.

Slope (with standard deviation) and correlation coefficient corresponding to the relations shown in Fig. 8, as determined through various statistical methods as indicated (see main text for details).

Isobe et al. (1990) discussed how, by definition, OLS(Y|X) and OLS(X|Y) can often lead to different results, why the former is typically preferred over the latter, and why OR should only be used with scale-free variables. The authors then recommend that the OLS bisector estimate is to be preferred, in particular when the goal is to probe the underlying relation between two variables.

Table B.1 also shows the slopes derived using the reweighted least squares regression (RLS) of Rousseeuw (1984), which is particularly robust against contamination from outliers. We used the code PROGRESS (Rousseeuw & Leroy 1987) to calculate the RLS slopes, as well as three correlation coefficients also shown in Table B.1: the Pearson correlation coefficient, the Spearman rank correlation coefficient and the RLS correlation coefficient.

One sees that the correlation between rk and Rbar is strong: The RLS correlation coefficient is particularly strong and the slope derived via the different methods is relatively stable (within the expected differences). The rk × A2 relation also sees a significant increase in the Spearman and RLS correlation coefficients, as compared to the Pearson coefficient. However, the three correlation coefficients are similar for the relations between rk and ϵbar and Bar/T. The relatively large variation between the slopes derived with OLS(Y|X) and the OLS bisector indicates that the correlations between rk and A2, ϵbar and Bar/T are only moderately significant. As already mentioned in Sect. 5, there is no significant correlation between rk and QB, and this is shown again by the results in Table B.1.

We stress again that further work is necessary to fully probe and understand these relations, in particular with larger samples. While the correlation between rk and Rbar is strong, and the trends between rk and ϵbar, Bar/T and A2 do not seem fortuitous, these are not proofs that nuclear discs are built by bars. However, these relations are consistent with that picture, and there is currently no obvious reason to think that a scenario in which bars are irrelevant to the formation of nuclear discs would produce such relations.

All Tables

Table 1.

Systemic velocities and central velocity dispersions.

Table 2.

Properties of nuclear discs, bars and galaxies in the TIMER sample.

Table B.1.

Slope (with standard deviation) and correlation coefficient corresponding to the relations shown in Fig. 8, as determined through various statistical methods as indicated (see main text for details).

All Figures

thumbnail Fig. 1.

Left: TIMER measurements of systemic velocity against the values available at LEDA. Right: TIMER measurements of central velocity dispersion against the LEDA values. The TIMER velocity dispersions were adapted to follow the same aperture corrections of LEDA, and the error bars in both panels correspond to 1σ errors. NGC 1300 is the only significant outlier but this is due to a mistake in LEDA (see text for further details; the LEDA team has been informed).

In the text
thumbnail Fig. 2.

Radial velocity, velocity dispersion, h3 and h4 maps for the stellar components in the TIMER galaxies, as indicated. The plotted range of the parameter measured is indicated at the bottom right corner of each panel. For radial velocity and velocity dispersion these are given in km s−1. The isophotes shown are derived from the MUSE data cube reconstructed intensities and are equally spaced in steps of about 0.5 mag. On average, 1″ corresponds to ≈100 pc. North is up, east to the left. See Appendix A for the maps corresponding to the rest of the galaxies in TIMER.

In the text
thumbnail Fig. 3.

Maps of v/σ for the stellar components in the TIMER galaxies, as indicated. The plotted range is indicated on the bottom right corner. The isophotes shown are derived from the MUSE data cube reconstructed intensities and are equally spaced in steps of about 0.5 mag. On average, 1″ corresponds to ≈100 pc. North is up, east to the left.

In the text
thumbnail Fig. 3.

continued.

In the text
thumbnail Fig. 4.

Deprojected radial profiles of v/σ at each Voronoi bin, with v corrected for inclination. NGC 1291, 1365 and NGC 6902 are not included. The values of maximum v/σ and rk are shown with the horizontal and vertical dashed lines, respectively.

In the text
thumbnail Fig. 5.

Distribution of the maximum value of v/σ (top) and of the kinematic radius (bottom) of the observed nuclear discs. NGC 1291, 1365 and NGC 6902 are not included. The vertical dotted lines mark the values corresponding to the bar-built nuclear disc in the model of Cole et al. (2014).

In the text
thumbnail Fig. 6.

λrk plotted against the observed (projected) ellipticity at the radius of the nuclear disc. The dashed line is the upper envelope of slow rotators in Emsellem et al. (2011). The solid line shows the relation for a model of oblate stellar systems viewed edge-on (see Binney 2005; Cappellari et al. 2007; Emsellem et al. 2007). The dotted lines correspond to the location of the same model for varying edge-on ellipticity (intrinsic flatness), from top (0.85) to bottom (0.35) in steps of 0.1, with edge-on systems on the relation and face-on systems towards the origin, as in Emsellem et al. (2011). NGC 1291, 1365 and NGC 6902 are not included.

In the text
thumbnail Fig. 7.

Distribution of the ratio between nuclear disc size, defined kinematically as the radius with maximum v/σ, and the effective radius of the photometric bulge component fitted through careful 2D multicomponent photometric decompositions. The values are clustered around a single value, following a normal distribution (rk/re = 0.86 ± 0.19), which shows that the photometric component defined as the bulge is (or is dominated by) the rapidly rotating nuclear disc.

In the text
thumbnail Fig. 8.

Relations between the kinematic sizes of nuclear discs and selected bar properties. Top left panel: rk is correlated with bar size (semi-major axis). The solid line is a simple linear regression fit to the data, whereas the dotted line corresponds to rk = 0.1 Rbar, the relation suggested in Shlosman et al. (1989). Top middle and right panels: we present how the kinematic size of nuclear discs relates to bar ellipticity and bar-to-total luminosity ratio, respectively. Bottom panels: relations between rk and different measures of bar strength (from Díaz-García et al. 2016, see text for details). Bar radii and ellipticities are from Kim et al. (2014) and are deprojected. The Pearson correlation coefficient R is indicated in each panel. In Appendix B we present an analysis of the statistical properties of these relations.

In the text
thumbnail Fig. A.1.

Same as Fig. 2 but for the remainder of our sample.

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.