Stellar kinematics across the Hubble sequence in the CALIFA survey: General properties and aperture corrections
Key Words.:Galaxies: kinematics and dynamics – Galaxies: elliptical and lenticular, cD – Galaxies: spiral – Galaxies: structure – Galaxies: evolution – Galaxies: formation
We present the stellar kinematic maps of a large sample of galaxies from the integral-field spectroscopic survey CALIFA. The sample comprises 300 galaxies displaying a wide range of morphologies across the Hubble sequence, from ellipticals to late-type spirals. This dataset allows us to homogeneously extract stellar kinematics up to several effective radii. In this paper, we describe the level of completeness of this subset of galaxies with respect to the full CALIFA sample, as well as the virtues and limitations of the kinematic extraction compared to other well-known integral-field surveys. In addition, we provide averaged integrated velocity dispersion radial profiles for different galaxy types, which are particularly useful to apply aperture corrections for single aperture measurements or poorly resolved stellar kinematics of high-redshift sources. The work presented in this paper sets the basis for the study of more general properties of galaxies that will be explored in subsequent papers of the survey.
The motion of stars within galaxies is a fundamental property set very early on in their life. Ever since the detection of rotation of stars in the Milky Way and nearby systems (e.g., Lindblad 1927; Mayall 1951; Münch & Münch 1960), the study of stellar motions has been a fruitful avenue to pose important constraints on our knowledge about galaxy formation and evolution. The analysis of rotational over random motions in early-type galaxies, for instance, has led to the realization that bright early-type galaxies are likely triaxial objects supported by orbital anisotropy (e.g., Bertola & Capaccioli 1975; Illingworth 1977; Binney 1978), rather than rotation.
The coupling of long-slit spectrographs with telescopes 2 to 4 m in size has provided, over the last three decades, a wealth of spatially resolved observations that has greatly improved our understanding of the overall stellar motion and level of kinematic substructure in external galaxies (e.g., Davies et al. 1983; Bertola et al. 1984; Bender et al. 1994; Fisher 1997; Simien & Prugniel 1997; Rubin et al. 1999; Vega Beltrán et al. 2001; Aguerri et al. 2003; Falcón-Barroso et al. 2003; Pizzella et al. 2004; van den Bosch et al. 2015).
While the first integral-field units (IFUs) were already in place in the mid-90’s (e.g., Bacon et al. 1995), the first serious efforts to measure stellar kinematics on large samples of galaxies using these kinds of instruments did not occur until year 2001. One of the pioneer projects in this respect was the SAURON survey (Bacon et al. 2001; de Zeeuw et al. 2002). With a representative sample of 72 galaxies (24 ellipticals, 24 lenticulars, and 24 early-type spirals, later extended with observations of 18 late-type spirals), this survey has set the reference for stellar kinematic IFU studies (e.g., Emsellem et al. 2004; Falcón-Barroso et al. 2006; Ganda et al. 2006). The discovery of the slow and fast rotator families in early-type galaxies (Emsellem et al. 2007) served as the trigger for a larger project: the ATLAS survey (Cappellari et al. 2011), in which a volume complete sample of 260 early-type galaxies revisited many kinematic aspects, from the amount of global angular momentum (Emsellem et al. 2011) to a detailed account of kinemetric features (Krajnović et al. 2006, 2011). In parallel, the DiskMass survey mapped, the stellar kinematic properties of nearby late-type spirals with the aid of the PPak IFU (Roth et al. 2005; Kelz et al. 2006).
The CALIFA survey (Sánchez et al. 2012) was born to fill in existing gaps in other IFU surveys and to provide a morphologically unbiased view of the stellar kinematics in galaxies based on a large ( 600 galaxies) and homogeneous integral-field spectroscopic dataset. The main advantage of CALIFA over existing surveys resides in a sample selection that includes all morphological types, as well as a field-of-view (FoV) that extends up to several effective radii (R). While CALIFA is no longer the IFU survey with the largest number of observed objects in the nearby Universe, it still provides the best compromise between spatial coverage (1.8–3.7 R) and sampling (1 kpc). Currently ongoing IFU surveys are hampered in one way or another by these factors, for example, SAMI covers areas within 1.1–2.9 R with a spatial sampling 1.7 kpc (Croom et al. 2012; Bryant et al. 2015), while MaNGA primary sample targets have a spatial sampling of 3 kpc within 1.5 R (Bundy et al. 2015). The real revolution in this respect will take place when MUSE at the Very Large Telescope (Bacon et al. 2010) is used in survey mode, as anticipated by the very spectacular stellar kinematic cases presented in the first few years of operations (e.g., Emsellem et al. 2014; Krajnović et al. 2015; Gadotti et al. 2015; Iodice et al. 2015).
The goal of this paper is to present the first stellar kinematic maps extracted from the CALIFA survey, describe all the technical details of the extraction, and provide basic stellar velocity dispersion aperture corrections for elliptical and spiral galaxies. The maps presented here have already been used within the survey to establish the effect of galaxy interactions on the stellar kinematics of galaxies (Barrera-Ballesteros et al. 2014, 2015), constrain the pattern speed of barred galaxies across the Hubble sequence (Aguerri et al. 2015), to present a volume-complete Tully-Fisher relation (Bekeraitė et al. 2016a), and the velocity function of galaxies as a benchmark for numerical simulations (Bekeraitė et al. 2016b). Forthcoming papers of the survey will make use of this information, for example, to revisit the distribution of global angular momentum in nearby galaxies and determine their dark matter content. Falcón-Barroso et al. (2015) provides a preview of some highlights. For results on the kinematics of the ionized gas in CALIFA, see García-Lorenzo et al. (2015).
The paper is organized as follows. Section §2 describes the sample of 300 galaxies used in our study and how this sample compares with the full CALIFA sample. Section 3 summarizes the instrumental setup employed during the observations. In section 4 we provide details of our kinematic extraction and comparisons with other major IFU surveys. Section 5 explains the limit set by our instrumental setup in the measurement of stellar velocity dispersions. In section 6 we provide velocity dispersion aperture corrections for elliptical and spiral galaxies. Finally, we summarize our work and conclusions in section 7.
2 The CALIFA sample
This study is based on observations of 300 galaxies drawn from the CALIFA mother and extended samples111Extended sample galaxies, two objects in this study, have CALIFA IDs larger than 1000. See Table 1 and Sánchez et al. (2016) as well., which are part of the photometric catalog of the seventh data release (Abazajian et al. 2009) of the Sloan Digital Sky Survey (SDSS). The main selection criteria in the survey is an angular isophotal diameter (), which is followed by a limited range in redshift, . These constraints ensure an efficient use of the PPak IFU and excludes, together with an apparent magnitude cut at band Petrosian magnitude of 20 mag, the presence of too many dwarf galaxies in our sample. Walcher et al. (2014) provides more details about the sample selection criteria and an in-depth discussion of the effects they introduce in the survey.
The CALIFA sample contains a large number of galaxies with diverse kinematic properties: from slow rotating ellipticals, to disk-dominated fast rotating galaxies, and perturbed interacting systems. This paper is based on the V1200 data (see §3) available until June 2014. We removed from our original sample of 375 galaxies those cases where the quality of the resulting stellar kinematic maps was not sufficient (e.g., poor spatial sampling due to low-quality data) to guarantee a meaningful analysis. We also selected out those cases whose stellar kinematics appeared highly disturbed by the presence of large nearby companions or had clear indications of being in final stages of a merging process. While this criteria excluded cases like ARP 220 (shown in Figure 3), it did not remove cases like the Mice galaxies (see Wild et al. 2014 for a detailed CALIFA study of this system), where the interaction has not drastically affected the observed kinematics. Barrera-Ballesteros et al. (2014, 2015) carefully examine the stellar kinematics of merging systems in the CALIFA survey. Our final sample thus consists of 300 galaxies.
In Fig. 1 we show the distribution of CALIFA galaxies presented in this paper as a function of Hubble type, stellar mass, and total absolute magnitude in the band. Hubble type classification was determined after a careful visual inspection by several members of the team. Stellar masses and total absolute magnitudes were derived following the prescriptions described in Walcher et al. (2014). Stellar masses assume a Chabrier initial-mass function (Chabrier 2003). While the number of galaxies represents a major improvement over other integral-field surveys, the selection criteria adopted in the CALIFA survey introduce an important shortcoming: the lack of low-mass, low-luminosity early-type systems and high-mass, high-luminosity late-type galaxies. Another important aspect is that our selection criteria favors edge-on orientations for the lowest mass and fainter systems (i.e., Sd galaxies). The advantage of this selection, however, is that it allows us to volume-correct averaged quantities and thus provide kinematic results that are representative of the general population of galaxies. Table 1 contains the basic properties of the subset of galaxies of our study.
We illustrate how representative our subsample is with respect to the mother sample in Fig. 2. The top and middle panels of the figure show the distribution of both the mother sample and our subsample in redshift, isophotal diameter (A) and petrosian -band magnitude (M)222Total absolute magnitudes are used throughout this paper, except in Fig. 2 where petrosian magnitudes are employed instead for consistency with Walcher et al. (2014).. The vertical lines indicate the limits in absolute magnitude in which the CALIFA mother sample is representative. In this space of parameters, the distribution of our subsample is consistent with that shown by the mother sample. Furthermore our galaxies cover all areas sampled by the mother sample. The bottom panel compares the luminosity function of SDSS (Blanton et al. 2003), the CALIFA mother sample, and the subset of 300 galaxies of the kinematic sample. We have applied a Kolmogorov-Smirnov test to the different parameters and confirm that the kinematic and mother samples are statistically consistent. Therefore the set of 300 galaxies studied in this paper are a good representation of the overall population of galaxies of all morphological types in the nearby Universe, within the luminosity and size constraints imposed by the CALIFA target selection.
3 Instrumental setup
The data presented in this paper is part of the CALIFA survey and as such were observed with the PMAS instrument (Roth et al. 2005) in the PPak mode (Verheijen et al. 2004; Kelz et al. 2006), mounted at the 3.5 m telescope of the Calar Alto observatory. For each galaxy, our observations cover the central 74″ 64″ using a hexagonal fiber bundle. For a detailed description of the observations and data reduction, see the CALIFA presentation article (Sánchez et al. 2012) and the CALIFA Data Release papers 1 and 2 (Husemann et al. 2013; García-Benito et al. 2015). The stellar kinematics presented in this paper is based on data from the v1.4 data reduction pipeline. Here we give a brief overview of the features of the observational setup that are relevant to our scientific interests.
The CALIFA survey is conducted in two instrumental setups: a low resolution mode (V500) with 850 at 5000 Å and a medium resolution mode (V1200) with 1650 at 4500 Å. The V500 grating covers a broad spectral range (3700–7300 Å) and includes a number of absorption and emission features, from the Ca H+K and [O II]3727 to H and [S II]6731 lines. The V1200 grating covers a smaller spectral window (3400–4750 Å). After careful evaluation of the spectral resolutions of the two gratings we established a value of 6 Å (FWHM 327 km s) for the V500 and 2.3 Å (FWHM 169 km s) for the V1200 gratings, respectively (see Husemann et al. 2013).
4 Stellar kinematics extraction
We extracted the stellar kinematics from every galaxy in a uniform way using both instrumental setups, i.e., V500 and V1200. Before accomplishing this, we applied spatial masks to remove spurious effects such as bad pixels, nearby objects, and/or foreground stars. We then logarithmically rebinned the spectra in each data cube to conserve a linear step in velocity space. We trimmed the data to contain only a useful spectral range: 3800–7000 Å for the V500 and 3750–4550 Å for the V1200 setup. We then selected for future use all spaxels within the isophote level where the average signal-to-noise ratio333We define our S/N as the average within the spectral range used in the fitting process. (S/N) was larger than 3. This cut ensured the removal of low-quality spaxels, which could introduce undesired systematic effects in our data at low surface brightness regimes. The next step was to spatially bin the data cubes to achieve an approximately constant S/N of 20 (per pixel). This value allows us to conserve a good spatial resolution while still being able to reliably estimate the first two moments of the line-of-sight velocity distribution (LOSVD). For this step we used the Voronoi 2D binning method of Cappellari & Copin (2003). Special care was taken in the S/N calculation to account for the correlation in the error spectrum of nearby spaxels (see Husemann et al. 2013 for details).
We measured the stellar kinematics of all galaxies in our sample using the pPXF code of Cappellari & Emsellem (2004). We used as templates the Indo-U.S. spectral library (Valdes et al. 2004) from which we selected 330 stars that uniformly cover the parameter space in gravity, metallicity, and effective temperature. The careful choice of stellar spectra is crucial to minimize template mismatch effects. We confirmed that, using our subset of 300 stars, we could reproduce the same results obtained using the full library. A non-negative linear combination of those stellar templates, convolved with a Gaussian LOSVD, was fitted to the spectrum of each Voronoi bin. The best-fitting parameters were determined by minimization in pixel space. In the wavelength region covered by CALIFA, there are several emission lines that needed to be masked during the fitting process, for example, [O II], [Ne III], H, H, [S II], H, [Fe II], H, [O III], He II, [Ar IV], H, [N I], He I, [O I], [N II], and H. We used a generous band width of 500 km s around those lines during the fitting process. This window was enough to mask the emission in all our galaxies. We also masked the regions affected by sky line residuals and the sodium doublet at 5890 Å. Additionally, a low-order additive Legendre polynomial was included in the fit to account for small differences in the continuum shape between the galaxy spectra and the input library. An order 6 polynomial was the minimum that ensured no large-scale wiggles in the residual spectra. In the end, the best-fitting values ( and ), and their associated uncertainties, were determined as the bi-weight mean and standard deviations of a set of 100 Monte Carlo realizations of the fitting. As expected, the distribution of best-fitting parameters from the Monte Carlo iterations are well-behaved and follow a Gaussian distribution. The bi-weight values measured from those distributions agree very well with those obtained from the direct fitting of the original spectra.
In Fig. 3 we show a few representative stellar velocity and velocity dispersion maps obtained with the V1200 grating. The remaining maps are presented in Appendix A of the Online Material. The four examples shown in the figure illustrate the diversity in the kinematics observed in the survey. NGC 6125, is a slow-rotator (e.g., low velocity amplitude and overall large velocity dispersion). NGC 1167 is an early-type spiral galaxy with large velocity and central velocity dispersion amplitudes. NGC 4210 is a disk-dominated galaxy (e.g., high velocity amplitude and overall small velocity dispersion). ARP 220 is an interacting system (e.g., with complex stellar velocity and velocity dispersions maps).
4.1 Stellar kinematics coverage
The CALIFA data presented in this paper allow us to produce stellar kinematic maps up to a typical surface brightness level of 19 mag arcsec (and as faint as 20 mag arcsec) in band. We also quantified how far, in terms of R, our maps extend. This is shown in Fig. 4, where we plot the maximum radius reached by the measurements in our V1200 maps over R. More than 90% of the sample covers at least up to 1 R and 39% extends beyond 2 R, with 50% of the galaxies reaching at least 1.8 R. This is a significant improvement over previous IFU surveys (e.g., SAURON, ATLAS), which aimed to probe different properties up to 1 R. The strength of those surveys, however, resides in the study of nearby systems at a higher spatial resolution, which allows them to detect small-scale inner kinematic subcomponents (e.g., McDermid et al. 2006; Krajnović et al. 2011).
4.2 Comparison between V500 and V1200 kinematics
The two instrumental setups used for the CALIFA survey give us the interesting opportunity of measuring the stellar kinematics of galaxies from independent datasets. As described in §3, one of these setups (V500) offers a much lower spectral resolution than the other, which turns out to be not enough to measure the lowest velocity dispersions present in our sample. This issue is clearly shown in Fig. 5, which presents the difference in velocity dispersion for each setup, measured within a 3″diameter aperture centered in each galaxy. In this panel systematic differences appear at dispersion values below 100 km s. We also compared the line-of-sight velocities from each setup (not shown here) and, as expected, found that they are well within the uncertainties of our measurements. Given this limitation, from now on we only report about results coming from the V1200 grating.
4.3 Comparison with other surveys
As an additional test to check the reliability and accuracy of our kinematic extraction, we compared our central velocity dispersion values with those provided by the SDSS DR7 survey (Abazajian et al. 2009) for those galaxies in our sample with SDSS spectroscopy available. We mimicked the SDSS aperture and extracted our velocity dispersions within a 3″ diameter aperture centered in each galaxy. The result of this comparison is presented in Fig. 6. The agreement between the two sets of measurements is good in general showing only a small systematic offset of km s, which is likely due to differences in the extraction method, set of templates, point-spread function (PSF)/seeing effects, and inaccuracies in the determination of the spectral resolution of both the data and templates. Similar levels of discrepancy and systematic differences have been identified in the past with SDSS DR7 measurements, even using the same SDSS dataset (see Fig. 6 in Oh et al. 2011).
An even more stringent test is the direct comparison of our stellar velocity dispersion maps to those of other surveys. We found up to six objects in common with one of the reference IFU surveys today, which is ATLAS. We focus our test on the velocity dispersion maps, as the velocity maps (not shown here) are in good agreement. The results of this comparison are presented in Fig. 7. By construction, all the ATLAS galaxies are early-type systems, which are predominantly red objects with fairly high central velocity dispersions. The figure shows a color-composite SDSS image of each galaxy in common, as well as the dispersion maps of both surveys and radial velocity dispersion profiles (extracted in circular annuli). The overall agreement between the two surveys is very good, despite differences in S/N thresholds applied in each survey. While the ATLAS data was Voronoi binned to a S/N of 40, we deemed it necessary to adopt a threshold S/N of 20 to find a good balance between spatial resolution (i.e., Voronoi bin sizes) and spatial coverage. The bigger footprint of the PPak IFU allows us to reach well beyond 1 R for most of the sample (see Fig. 2), which is a significant improvement over ATLAS. The large bins are also responsible for the smoother trends observed in the CALIFA radial profiles.
The only major difference between the two datasets, however, is on the central dispersion values. In general, ATLAS values are larger. There are two main reasons that could explain this behavior. The PPak IFU is a fiber bundle made of 2.7″ wide fibers, as opposed to ATLAS with 1″ lenslets. While our dithering strategy during observations (see Sánchez et al. 2012 for details) allowed us to resample our final data cube to 1″ per spaxel, the original fiber size could result in lower velocity dispersion values due to beam smearing. We tested this scenario by comparing the central ATLAS values with those obtained by collapsing the ATLAS spectra within a 3″aperture (similar to a CALIFA fiber). Our results show that velocity dispersion values can decrease by up to 15%. This effect can therefore explain part of the discrepancy between the two surveys.
In addition, the effective PSF measured for the CALIFA survey (García-Benito et al. 2015) could also affect these values. While reported seeing conditions between the two surveys appear similar, if the CALIFA PSF was worse than the ATLAS PSF, this could also explain part of the decrease in the central velocity dispersion. Based on some simulations carried out in the context of another CALIFA paper (Méndez-Abreu et al. 2016, submitted), we estimated that the PSF can account for up to 5% difference in the observed values. On top of that, the level of Voronoi binning could play a similar role, although this seems unlikely in our case as the central spaxels remain mostly unbinned.
5 Reliability of velocity dispersion measurements below the instrumental resolution
An important aspect to consider when extracting stellar kinematics of galaxies is to understand the limiting velocity dispersions one can reach given the spectral resolution provided by the instrument used. The safest option is to use an instrumental setup where the spectral resolution is better than the expected values. Under certain circumstances, however, it is possible to push this limit and measure velocity dispersions below the nominal threshold imposed by the instrument. As shown in Ryś et al. (2013), but see also González (1993) and Pedraz et al. (2002), a combination of high signal-to-noise and spectral sampling of the line spread function (e.g., Koleva et al. 2009) makes it possible to overcome, to some extent, this limitation.
5.1 Comparison between DiskMass and CALIFA datasets
The nominal spectral dispersion of the CALIFA V1200 data is 72 km s. We determined the velocity dispersion limit of our data using as a reference three galaxies in the DiskMass survey (DMS; Bershady et al. 2010). This dataset was designed to measure velocity dispersions in face-on, disk galaxies and the spectral resolution of the instrument (FWHM=0.69Å) was chosen to safely reach values around 17 km s (Martinsson et al. 2013). The PPak was custom built for the DMS and subsequently employed in the CALIFA survey, which can help suppress systematic effects inherent in the analysis. The DMS team has kindly provided their data for three galaxies. One was already in common with the CALIFA survey (NGC 234). We observed two more, specifically for these tests, in 21–23 February 2014 (UGC 4256, UGC 4458), using the same V1200 instrumental configuration of the main survey.
Before carrying out our tests, and to account for potential systematic effects, we checked that neither the method (cross-correlation technique versus pPXF) used to measure the stellar kinematics had a strong impact on the resulting velocity dispersions. Our own extraction, using pPXF, of velocity dispersions from DMS data provided fully consistent results. The choice of templates, whether a single star (as the DMS team used) or a full stellar library (like in our case), did not cause any systematic difference in this particular exercise. Therefore template mismatch is not an issue in these tests. The successful comparison of the two methods using the same data was also reported by Westfall et al. (2011).
In addition to the difference in spectral resolution and template mismatch, there are some further differences with respect to the DMS team analysis that can cause systematic effects in the velocity dispersion values. Spatial binning is desirable in general to reach a threshold S/N, but it can also have the negative effect of artificially broadening the line-of-sight velocity distribution. This is more acute in the outer regions of galaxies, where the S/N drops quickly and the combination of a larger number of spectra is required. Despite this drawback, we chose to Voronoi bin the data to ensure a minimum quality of the spectra used to derive the velocity dispersion. The DMS team preferred to extract their values on single spaxels of 3″ diameter, and remove values with uncertainties larger than 8 km s(see § 7.3.2 of Martinsson et al. 2013).
Another important issue is the wavelength range used to derive the velocity dispersion. The DMS values rely on measurements in the short spectral range between 4980–5370 Å. Our CALIFA values are based on fits between 3750–4550 Å. While a longer baseline is in principle preferred, different spectral features may have slightly distinct intrinsic broadening (at the spectral resolutions we are considering here). We believe this may be the case in the CALIFA spectral range with the Ca H+K lines. We attempted to derive our stellar kinematics ignoring those lines, but results were noisier and uncertainties larger as the fits rely on a few weak spectral features, for example, Fe ( 4383 Å), H, and H. As shown in Kobulnicky & Gebhardt (2000), the Ca H+K lines are reliable features to obtain stellar kinematics in all kinds of galaxies, although their results appeared to be more uncertain for late-type systems. This may be the culprit of some of the differences we see with the DiskMass survey (see below). The detailed characterization of all these effects is a complex task, and even if we could measure the systematic deviations introduced by each effect, it is not totally obvious they would affect different kinds of galaxies in the same manner (e.g., emission-free early-type galaxies versus highly star-forming spiral disks).
Figure 8 shows the comparison of our own CALIFA data with the stellar velocity dispersions measured by the DMS team on the three galaxies in common. As in Fig. 7, we plot a SDSS color image of each galaxy, velocity dispersion maps of the two surveys, and also radial profiles (extracted in circular annuli) for a more direct comparison. When displaying the maps and radial profiles, we plot the individual spaxel measurements provided by the DMS team and our Voronoi binned values. The agreement between the two datasets is good overall. We do see discrepancies in some measurements (most noticeable in NGC 234 and UGC 4256). It appears that some of our CALIFA measurements are much larger than those reported by the DMS team at a given radius. We explored the reasons for these discrepant values and concluded that they occur in low surface brightness regions ( 22 mag arcsec) that are often affected by dust or unmasked foreground stars. They often have S/N slightly below the nominal S/N=20 threshold, which is permitted by the Voronoi binning routine within some tolerance (see Cappellari & Copin 2003). These values are also naturally associated with large Voronoi bins, which can also artificially help to increase the broadening. However, the pPXF fits in those regions are not particularly worse than in other areas with similar level of binning, S/N, or surface brightness levels. Given that there might be some physical insight as to why those values are high (e.g., dust obscuration, multiple kinematic components, and kinematic flaring in the outer parts of galaxies), we prefer to keep them in our data release and let the user, based on diagnostic parameters we provide, decide whether to include or exclude them depending on their science case. This effect is not evident in the DMS values owing to the partial field-of-view coverage of their data.
5.2 Limiting velocity dispersion and relative uncertainties
In order to establish the lowest reliable velocity dispersion we can measure, we directly compared the DMS and ATLAS3D values to our CALIFA measurements. This is shown in the top panel of Fig. 9, where we present the ratio of the CALIFA over the DMS and ATLAS3D dispersion values as a function of the DMS or ATLAS3D measurements. This exercise determines at which velocity dispersion values our CALIFA measurements depart systematically from the one-to-one relation. For a fair comparison, we used the Voronoi values of our CALIFA maps at the locations of DMS measurements. This is a better approach than interpolating our maps at those locations, which may produce artifacts. The drawback of this approach is that there is some instrinsic scatter produced by the sampling of our points in locations that could be far from the Voronoi centroids in our data. It is also sensitive to the different levels of scatter of the data points in the surveys (e.g., the scatter of ATLAS3D points is larger than CALIFA, see Fig. 7). While the number of points is not too large for the DMS survey ( 360 measurements), it is enough to compute some statistics. Besides the individual datapoints, we indicate the limiting 16% and 84% percentiles of the distribution with dotted lines. The median of the distribution is marked with a solid line. The panel shows that velocity dispersion values are consistent within the uncertainties down to 40 km s. Below that point, CALIFA measurements are systematically larger up to a factor 3 on average for values around 20 km s. On the high velocity dispersion end, values converge asymptotically to unity, as expected, except for the most massive systems where we suffer the PSF and beam smearing effects discussed in §4.3 for ATLAS3D.
In addition to the potential bias in our measurements, it is also interesting to determine the relative uncertainties of our measurements at different velocity dispersion regimes. This is presented in the bottom panel of Fig. 9. We produced this figure using all the individual Voronoi bin measurements and uncertainties for the 300 CALIFA galaxies presented here. The shaded region and lines as in the top panel. The figure shows that uncertainties are rather small around 5% for 150 km s. Below that value, relative uncertainties increase up to 50% for velocity dispersions 20 km s. The median uncertainty at 40 km s, where our measurements start deviating systematically from the DMS values, is 20%.
6 Aperture profiles
The large number of galaxies across the Hubble sequence in our study allows us to estimate velocity dispersion aperture corrections for different groups of galaxies. These corrections are useful to homogenize dispersion values measured with fiber-fed spectrographs (e.g., SDSS) for galaxies at different distances, and they are particularly important for high-redshift studies.
We studied the behavior of the integrated velocity dispersion profiles of our galaxies, extracted in elliptical apertures with a fixed position angle and ellipticity (as indicated in Table 1). We chose elliptical rather than circular apertures to properly account for inclination effects. We used the position angle and ellipticity measured in the outer parts of the galaxy (as listed in Table 1). While this choice ignores potential radial variations in these two parameters (e.g., due to bars), the velocity dispersion maps do not appear to be clearly influenced by those photometric deviations. This is true in particular for barred galaxies, as shown in Seidel et al. (2015) or Gadotti et al. (2015).
We found three types of radial profiles: (Class 1) those that decrease steadily, (Class 2) those decreasing up to a certain radius and then increasing again, and (Class 3) those that increase steadily with radius. We analyzed the type of galaxies conforming each class and found that class 1 objects are typically early types (e.g., ellipticals, lenticulars and Sa galaxies). Class 2 is made of a rather small set of objects (20) that are mostly early-type galaxies too. They show fairly high dispersions in the center but also reasonably high rotation in the outer parts, which drive the increase of integrated velocity dispersion at large radii. This effect is even more pronounced in class 3 objects, that are predominantly late-type systems (e.g., Sb, Sc, and Sd galaxies).
Given the small number of objects in class 2, we only provide aperture corrections for the other two groups (class 1 and 3). We followed previous works in the literature and fit the individual profiles of each galaxy in these two classes using a power-law function in the form
where the effective radius (R) is used as a normalization factor for both the radius and velocity dispersion. An important aspect to consider during the fitting process was the effect of the PSF in our measurements. We account for this effect by convolving our models for each galaxy with the CALIFA PSF during the fitting process. As illustrated in the comparison with the ATLAS survey data (see § 4.3), our velocity dispersions are probably smaller than they should at the very center of galaxies. Ignoring this effect artificially lowers the parameter in the power-law function.
Figure 10 shows the individual profiles for classes 1 and 3. Class 1, in the top panel, is made of predominantly early-type systems with an average stellar mass of M and absolute magnitude M mag. We determined the average profile for the class by weighting with the volume correction factor (V) of each galaxy. That provides a good representation of the average profile for early-type galaxies with those properties. The average fit and uncertainty is indicated with the solid red line and black dashed lines, respectively. The average value of is in good agreement with corrections reported in previous works (e.g., , Jorgensen et al. 1995; , Mehlert et al. 2003; , Cappellari et al. 2006), but see Ziegler & Bender 1997 for a steeper correction using a different prescription.
The family of late-type systems in class 3 is much more heterogeneous. We decided to divide the sample into three intervals of mass and absolute magnitude. As for early-type galaxies, PSF effects and volume corrections were taken into account for the fitting. As illustrated in Fig. 10 (bottom panel) there are significant differences in the slopes as a function of mass and magnitude. Our results indicate that low-mass and/or low-luminosity spiral galaxies display larger values than high-mass and/or bright systems.
7 Summary & conclusions
In this paper we present stellar kinematic maps for a sample of 300 galaxies that are part of the CALIFA survey. The sample covers a wide range of Hubble types, from ellipticals to late-spiral galaxies. This subset is a good representation of the CALIFA mother sample in terms of redshift, isophotal diameter, and absolute magnitude. The large footprint of the PPak IFU, together with the average distance of the survey, allow us to measure stellar kinematics well beyond 1.8 R for 50% of the galaxies, reaching out to 4–5 R in a few exceptional cases. The penalty, caused by the combination of spatial sampling and distance, is the inability to detect kinematically decoupled components at the centers of galaxies. Still our data is well suited for the study of large-scale kinematic twists or long-axis rotation, which occurs in a handful of objects.
The measurements presented in this paper are in good agreement with those of other well-known IFU surveys (e.g., ATLAS and DiskMass). The detailed comparison with the DiskMass survey allowed us to establish that we can measure reliable velocity dispersion values down to 40 km s (i.e., 30 km s below the instrumental resolution). We also characterized the relative uncertainties of our measurements, which are around 5% for 150 km s. Below that value, relative uncertainties increase up to 50% for velocity dispersions all the way down to 20 km s.
We also took advantage of our large sample to compute integrated stellar velocity dispersion aperture corrections for different sets of galaxies across the Hubble sequence. These corrections are particularly useful to homogenize dispersion values of galaxies at different distances. We find two main classes of integrated aperture radial profiles: steadily decreasing profiles representative of early-type galaxies, and a second class of systematically increasing profiles typical of late-type spiral galaxies. We provide aperture corrections for each class for different stellar masses and absolute magnitudes.
The main properties of the sample and the stellar velocity and velocity dispersion maps introduced in this paper are available as part of the Online Material in Table 1 and Appendix A. The values of the maps themselves, together with many diagnostic parameters to assess the quality of the measurements, will be made available to the community at the CALIFA website (http://califa.caha.es).
Acknowledgements.We would like to thank the anonymous referee for constructive comments that helped improve some aspects of the original manuscript. We are also grateful to the DiskMass survey team for sharing their data with us for the spectral resolution tests, and to Marc Verheijen and Kyle Westfall in particular for in-depth discussions on the topic. This study makes use of the data provided by the Calar Alto Legacy Integral Field Area (CALIFA) survey (http://califa.caha.es). Based on observations collected at the Centro Astronómico Hispano Alemán (CAHA) at Calar Alto, operated jointly by the Max-Planck-Institut für Astronomie and the Instituto de Astrofísica de Andalucía (CSIC). CALIFA is the first legacy survey being performed at Calar Alto. The CALIFA collaboration would like to thank the IAA-CSIC and MPIA-MPG as major partners of the observatory, and CAHA itself, for the unique access to telescope time and support in manpower and infrastructures. The CALIFA collaboration thanks also the CAHA staff for the dedication to this project.
Funding and financial support acknowledgements: J. F-B. from grant AYA2013-48226-C3-1-P from the Spanish Ministry of Economy and Competitiveness (MINECO). J. F-B. and GvdV from the FP7 Marie Curie Actions of the European Commission, via the Initial Training Network DAGAL under REA grant agreement number 289313. J. M-A. and V. W. acknowledge support from the European Research Council Starting Grant (SEDMorph P.I. V. Wild). P.S-B acknowledge financial support from the BASAL CATA Center for Astrophysics and Associated Technologies through grant PFB-06. R. M. GD. from grant AYA2014-57490-P. R. G-B, R. M. GD. and EP acknowledge support from the project JA-FQM-2828. C. J. W. acknowledges support through the Marie Curie Career Integration Grant 303912. L. G. from the Ministry of Economy, Development, and Tourism’s Millennium Science Initiative through grant IC120009 awarded to The Millennium Institute of Astrophysics (MAS), and CONICYT through FONDECYT grant 3140566. I. M. from grant AYA2013-42227-P.
- Abazajian et al. (2009) Abazajian, K. N., Adelman-McCarthy, J. K., Agüeros, M. A., et al. 2009, ApJS, 182, 543
- Aguerri et al. (2003) Aguerri, J. A. L., Debattista, V. P., & Corsini, E. M. 2003, MNRAS, 338, 465
- Aguerri et al. (2015) Aguerri, J. A. L., Méndez-Abreu, J., Falcón-Barroso, J., et al. 2015, A&A, 576, A102
- Bacon et al. (2010) Bacon, R., Accardo, M., Adjali, L., et al. 2010, in Proc. SPIE, Vol. 7735, Ground-based and Airborne Instrumentation for Astronomy III, 773508
- Bacon et al. (1995) Bacon, R., Adam, G., Baranne, A., et al. 1995, A&AS, 113, 347
- Bacon et al. (2001) Bacon, R., Copin, Y., Monnet, G., et al. 2001, MNRAS, 326, 23
- Barrera-Ballesteros et al. (2014) Barrera-Ballesteros, J. K., Falcón-Barroso, J., García-Lorenzo, B., et al. 2014, A&A, 568, A70
- Barrera-Ballesteros et al. (2015) Barrera-Ballesteros, J. K., García-Lorenzo, B., Falcón-Barroso, J., et al. 2015, A&A, 582, A21
- Bekeraitė et al. (2016a) Bekeraitė, S., Walcher, C. J., Falcón-Barroso, J., et al. 2016a, ArXiv e-prints [\eprint[arXiv]1605.07632]
- Bekeraitė et al. (2016b) Bekeraitė, S., Walcher, C. J., Wisotzki, L., et al. 2016b, ArXiv e-prints [\eprint[arXiv]1608.01677]
- Bender et al. (1994) Bender, R., Saglia, R. P., & Gerhard, O. E. 1994, MNRAS, 269, 785
- Bershady et al. (2010) Bershady, M. A., Verheijen, M. A. W., Swaters, R. A., et al. 2010, ApJ, 716, 198
- Bertola et al. (1984) Bertola, F., Bettoni, D., Rusconi, L., & Sedmak, G. 1984, AJ, 89, 356
- Bertola & Capaccioli (1975) Bertola, F. & Capaccioli, M. 1975, ApJ, 200, 439
- Binney (1978) Binney, J. 1978, MNRAS, 183, 501
- Blanton et al. (2003) Blanton, M. R., Hogg, D. W., Bahcall, N. A., et al. 2003, ApJ, 592, 819
- Bryant et al. (2015) Bryant, J. J., Owers, M. S., Robotham, A. S. G., et al. 2015, MNRAS, 447, 2857
- Bundy et al. (2015) Bundy, K., Bershady, M. A., Law, D. R., et al. 2015, ApJ, 798, 7
- Cappellari et al. (2006) Cappellari, M., Bacon, R., Bureau, M., et al. 2006, MNRAS, 366, 1126
- Cappellari & Copin (2003) Cappellari, M. & Copin, Y. 2003, MNRAS, 342, 345
- Cappellari & Emsellem (2004) Cappellari, M. & Emsellem, E. 2004, PASP, 116, 138
- Cappellari et al. (2011) Cappellari, M., Emsellem, E., Krajnović, D., et al. 2011, MNRAS, 413, 813
- Chabrier (2003) Chabrier, G. 2003, PASP, 115, 763
- Croom et al. (2012) Croom, S. M., Lawrence, J. S., Bland-Hawthorn, J., et al. 2012, MNRAS, 421, 872
- Davies et al. (1983) Davies, R. L., Efstathiou, G., Fall, S. M., Illingworth, G., & Schechter, P. L. 1983, ApJ, 266, 41
- de Zeeuw et al. (2002) de Zeeuw, P. T., Bureau, M., Emsellem, E., et al. 2002, MNRAS, 329, 513
- Emsellem et al. (2011) Emsellem, E., Cappellari, M., Krajnović, D., et al. 2011, MNRAS, 414, 888
- Emsellem et al. (2007) Emsellem, E., Cappellari, M., Krajnović, D., et al. 2007, MNRAS, 379, 401
- Emsellem et al. (2004) Emsellem, E., Cappellari, M., Peletier, R. F., et al. 2004, MNRAS, 352, 721
- Emsellem et al. (2014) Emsellem, E., Krajnović, D., & Sarzi, M. 2014, MNRAS, 445, L79
- Falcón-Barroso et al. (2006) Falcón-Barroso, J., Bacon, R., Bureau, M., et al. 2006, MNRAS, 369, 529
- Falcón-Barroso et al. (2003) Falcón-Barroso, J., Balcells, M., Peletier, R. F., & Vazdekis, A. 2003, A&A, 405, 455
- Falcón-Barroso et al. (2015) Falcón-Barroso, J., Lyubenova, M., & van de Ven, G. 2015, in IAU Symposium, Vol. 311, Galaxy Masses as Constraints of Formation Models, ed. M. Cappellari & S. Courteau, 78–81
- Fisher (1997) Fisher, D. 1997, AJ, 113, 950
- Gadotti et al. (2015) Gadotti, D. A., Seidel, M. K., Sánchez-Blázquez, P., et al. 2015, A&A, 584, A90
- Ganda et al. (2006) Ganda, K., Falcón-Barroso, J., Peletier, R. F., et al. 2006, MNRAS, 367, 46
- García-Benito et al. (2015) García-Benito, R., Zibetti, S., Sánchez, S. F., et al. 2015, A&A, 576, A135
- García-Lorenzo et al. (2015) García-Lorenzo, B., Márquez, I., Barrera-Ballesteros, J. K., et al. 2015, A&A, 573, A59
- González (1993) González, J. J. 1993, PhD thesis, Thesis (PH.D.)–University of California, Santa Cruz, 1993.Source: Dissertation Abstracts International, Volume: 54-05, Section: B, page: 2551.
- Husemann et al. (2013) Husemann, B., Jahnke, K., Sánchez, S. F., et al. 2013, A&A, 549, A87
- Illingworth (1977) Illingworth, G. 1977, ApJ, 218, L43
- Iodice et al. (2015) Iodice, E., Coccato, L., Combes, F., et al. 2015, A&A, 583, A48
- Jorgensen et al. (1995) Jorgensen, I., Franx, M., & Kjaergaard, P. 1995, MNRAS, 276, 1341
- Kelz et al. (2006) Kelz, A., Verheijen, M. A. W., Roth, M. M., et al. 2006, PASP, 118, 129
- Kobulnicky & Gebhardt (2000) Kobulnicky, H. A. & Gebhardt, K. 2000, AJ, 119, 1608
- Koleva et al. (2009) Koleva, M., Prugniel, P., Bouchard, A., & Wu, Y. 2009, A&A, 501, 1269
- Krajnović et al. (2006) Krajnović, D., Cappellari, M., de Zeeuw, P. T., & Copin, Y. 2006, MNRAS, 366, 787
- Krajnović et al. (2011) Krajnović, D., Emsellem, E., Cappellari, M., et al. 2011, MNRAS, 414, 2923
- Krajnović et al. (2015) Krajnović, D., Weilbacher, P. M., Urrutia, T., et al. 2015, MNRAS, 452, 2
- Lindblad (1927) Lindblad, B. 1927, MNRAS, 87, 553
- Martinsson et al. (2013) Martinsson, T. P. K., Verheijen, M. A. W., Westfall, K. B., et al. 2013, A&A, 557, A130
- Mayall (1951) Mayall, N. U. 1951, Publications of Michigan Observatory, 10, 19
- McDermid et al. (2006) McDermid, R. M., Emsellem, E., Shapiro, K. L., et al. 2006, MNRAS, 373, 906
- Mehlert et al. (2003) Mehlert, D., Thomas, D., Saglia, R. P., Bender, R., & Wegner, G. 2003, A&A, 407, 423
- Méndez-Abreu et al. (2016) Méndez-Abreu, J., Aguerri, J. A. L., Falcón-Barroso, J., et al. 2016, A&A
- Münch & Münch (1960) Münch, G. & Münch, L. 1960, ApJ, 131, 253
- Oh et al. (2011) Oh, K., Sarzi, M., Schawinski, K., & Yi, S. K. 2011, ApJS, 195, 13
- Pedraz et al. (2002) Pedraz, S., Gorgas, J., Cardiel, N., Sánchez-Blázquez, P., & Guzmán, R. 2002, MNRAS, 332, L59
- Pizzella et al. (2004) Pizzella, A., Corsini, E. M., Vega Beltrán, J. C., & Bertola, F. 2004, A&A, 424, 447
- Roth et al. (2005) Roth, M. M., Kelz, A., Fechner, T., et al. 2005, PASP, 117, 620
- Rubin et al. (1999) Rubin, V. C., Waterman, A. H., & Kenney, J. D. P. 1999, AJ, 118, 236
- Ryś et al. (2013) Ryś, A., Falcón-Barroso, J., & van de Ven, G. 2013, MNRAS, 428, 2980
- Sánchez et al. (2016) Sánchez, S. F., García-Benito, R., Zibetti, S., et al. 2016, ArXiv e-prints [\eprint[arXiv]1604.02289]
- Sánchez et al. (2012) Sánchez, S. F., Kennicutt, R. C., Gil de Paz, A., et al. 2012, A&A, 538, A8
- Seidel et al. (2015) Seidel, M. K., Falcón-Barroso, J., Martínez-Valpuesta, I., et al. 2015, MNRAS, 451, 936
- Simien & Prugniel (1997) Simien, F. & Prugniel, P. 1997, A&AS, 126
- Valdes et al. (2004) Valdes, F., Gupta, R., Rose, J. A., Singh, H. P., & Bell, D. J. 2004, ApJS, 152, 251
- van den Bosch et al. (2015) van den Bosch, R. C. E., Gebhardt, K., Gültekin, K., Yıldırım, A., & Walsh, J. L. 2015, ApJS, 218, 10
- Vega Beltrán et al. (2001) Vega Beltrán, J. C., Pizzella, A., Corsini, E. M., et al. 2001, A&A, 374, 394
- Verheijen et al. (2004) Verheijen, M. A. W., Bershady, M. A., Andersen, D. R., et al. 2004, Astronomische Nachrichten, 325, 151
- Walcher et al. (2014) Walcher, C. J., Wisotzki, L., Bekeraité, S., et al. 2014, A&A, 569, A1
- Westfall et al. (2011) Westfall, K. B., Bershady, M. A., & Verheijen, M. A. W. 2011, ApJS, 193, 21
- Wild et al. (2014) Wild, V., Rosales-Ortega, F., Falcón-Barroso, J., et al. 2014, A&A, 567, A132
- Ziegler & Bender (1997) Ziegler, B. L. & Bender, R. 1997, MNRAS, 291, 527
Note. – Col. 1: galaxy name. Col. 2: CALIFA identification number for each galaxy. Col. 3 redshift of the galaxy from SDSS (Abazajian et al. 2009). Col. 4: position angle of the galaxy measured in the outer parts, using SDSS images. Col. 5: average ellipticity measured in the outer parts of the galaxy, using SDSS images. Col. 6: Hubble type of the galaxy from Walcher et al. (2014). Col. 7: total stellar mass of the galaxy, measured as described in Walcher et al. (2014). Col. 8: total absolute magnitude in band from SDSS (Abazajian et al. 2009). Col. 9: effective radii (in arcsec) of the galaxy, measured as described in Walcher et al. (2014). Col. 10: maximum radial extent of our kinematic maps (in arcsec).
Appendix A Stellar kinematic maps
This online material presents all the stellar velocity (Figs. A1–A17) and velocity dispersion (Figs. A18–A34) maps extracted from the V1200 grating used in this paper. The complete sample comprises 300 galaxies of Hubble morphological types ranging from ellipticals to late-type spirals. Velocity maps are in km s and use a fixed range in the interval [,150] km s. Velocity dispersion maps are also expressed in km s and use a fixed range from 20 to 300 km s. Color schemes as in Fig. 3. Overlaid contours come from SDSS band images and have been limited to the isophote reaching 2 R. All panels cover an area of 80″100″.