Issue 
A&A
Volume 531, July 2011



Article Number  A143  
Number of page(s)  11  
Section  Astronomical instrumentation  
DOI  https://doi.org/10.1051/00046361/201016386  
Published online  01 July 2011 
Accurate stellar rotational velocities using the Fourier transform of the cross correlation maximum
^{1}
Universidad Nacional de San Juan, Av. J. I. de la Roza 590 oeste, 5400 Rivadavia, San Juan, Argentina
email: cgonzadiaz@gmail.com
^{2}
ICATE, CONICET, Av. España 1512 sur, J5402DSP San Juan, Argentina
email: fgonzalez@icateconicet.gob.ar; hlevato@icateconicet.gob.ar
Received: 22 December 2010
Accepted: 29 March 2011
Aims. We propose a method for measuring the projected rotational velocity vsini with high precision even in spectra with blended lines. Though not automatic, our method is designed to be applied systematically to large numbers of objects without excessive computational requirement.
Methods. We calculated the cross correlation function (CCF) of the object spectrum against a zerorotation template and used the Fourier transform (FT) of the CCF central maximum to measure the parameter vsini taking the limb darkening effect and its wavelength dependence into account. The procedure also improves the definition of the CCF base line, resulting in errors related to the continuum position under 1% even for vsini = 280 km s^{1}. Tests with highresolution spectra of Ftype stars indicate that an accuracy well below 1% can be attained even for spectra where most lines are blended.
Results. We have applied the method to measuring vsini in 251 Atype stars. For stars with vsini over 30 km s^{1} (2–3 times our spectra resolution), our measurement errors are below 2.5% with a typical value of 1%. We compare our results with Royer et al. (2002a) using 155 stars in common, finding systematic differences of about 5% for rapidly rotating stars.
Key words: stars: fundamental parameters / stars: rotation / methods: data analysis
© ESO, 2011
1. Introduction
The projected axial rotational velocity of single stars can be measured directly from the broadening of their spectral lines. That is why the rotational velocity is an important observable in statistical studies of stellar astrophysics. Several methods have been developed to measure vsini, but the problem of systematic differences between different authors and methods has always been present. The first workable model for the stellar rotation was established by Carroll (1928, 1933) and Carroll & Ingram (1933), but Shajn & Struve (1929) presented a very simple graphical model called the classical model of a rotating star (CMRS) by Collins & Truax (1995), which was used as a standard for most of the work done in the 20th century (for more details see the excellent review paper by Collins 2004). Shajn & Struve (1929) did not consider limb darkening, which was introduced in the CMRS by Carroll (1933). In a series of papers Slettebak measured vsini for stars in the spectral range OG (Slettebak 1954, 1955, 1956; Slettebak & Howard 1955) by considering a linear limb darkening law. These stars were used in many subsequent papers by different authors to calibrate the full width at half depth of the line used to measure the vsini parameter. Therefore, determination of vsini was conditioned by the calibration made by the author.
In the last quarter of the 20th century Slettebak et al. (1975) established a new standard system of stars distributed over both hemispheres and covering a range in spectral type from O9 to F9. The measure of vsini was again based on a calibration of vsini with the full width at half depth of stellar absorption lines in the star spectrum. Slettebak et al. (1975) established this system by comparing their digital data with numerical models constructed by adopting Roche geometry, uniform angular velocity at the surface, von Zeipel gravity darkening, and numerical integration of angledependent model atmosphere intensities. The standards for the highvelocity specimens were established by Slettebak (1982) by visual comparison of line widths on the spectrograms that were recorded on photographic plates. Systematic effects between the old and new Slettebak systems were studied by García & Levato (1984). The calibrations of the systems described were made for only three lines, He i 4471 Å, Mg ii 4481 Å, and Fe i 4476 Å, and the theoretical profiles were calculated with main sequence models. Besides this, determination of the full width at half depth has shown to be very sensitive to the continuum position, which represents an important error source.
A significant improvement was possible thanks to methods based on Fourier transform (FT) of line profiles providing a deeper analysis. In the first place, this tool is used to avoid an external calibration by suppressing the error related to this stage of the process. In second place, with high signaltonoise ratio (S/N) it is possible to identify another broadening agents present in the line profile, at least in a qualitative way. Finally, it is possible to analyze secondorder effects like differential rotation (Gray 1977, 1982; Bruning 1981; Reiners & Schmitt 2002) and to put a limit on the inclination of the rotational axes (Reiners & Royer 2004). Regardless of the method adopted to determine vsini, it is not an easy task to find the spectral lines that have the minimum conditions required to be used in the measurement of rotational velocity, namely: 1) lines with no blends, 2) lines intense enough to be identified in rapid rotators, and 3) lines mainly broadened by rotation. For instance, in Otype and Btype stars, only the Balmer lines and some helium lines are intense enough, but they all show a significant Stark effect. In stars of spectral type later than A4, almost all lines are blended if vsini is grater than 100 km s^{1}, making it impossible to find an isolated line. This problem is evident in the work of Royer et al. (2002a) where fewer than three lines were measured in Atype stars with vsini > 60 km s^{1}. Because this problem increases with the spectral type, other methods have been developed.
One solution to the blending limitation is the use of a least square deconvolution procedure to derive the broadening function in a selected wavelength region instead of a single line profile. This methodology implies the application of an iterative process to fit the equivalent width of the template’s spectral lines, the broadening function, and the continuum position (for details see Reiners & Schmitt 2003). This is a very powerful technique for a detailed study of the rotational profile. However, for extensive applications like the development of a catalog of rotational velocities, a more direct method that does not involve fitting the intrinsic spectrum or including any atmospheric parameter other than limb darkening, might be more suitable.
Even though CCF had been originally proposed for determining radial velocities, the projected axial rotational velocity can be inferred from the width of the CCF maximum. This tool has often been applied to determine vsini in cool stars (ML spectral type). The standard procedure is to calculate the CCF between an observed spectrum and a template spectrum of the same temperature with vsini = 0 km s^{1}. Then, a fitting profile for the maximum of this function is calculated and the width of the fit can be used to measure vsini through an empirical calibration vsiniwidth. The absence of single lines is somehow solved by means of the CCF. Nevertheless, as in any other method that depends on empirical calibrations, various strategies, not always equivalent, have been adopted in the literature. In some works the central maximum is fitted with a Gaussian (BailerJones 2004), while in others a parabola (Tinney & Reid 1998), or even a Gaussian plus a quadratic function are used (White & Basri 2003). The adoption of different functions to evaluate the rotational broadening could lead to systematic differences in results from different authors.
Templates selection is also very heterogeneous. Even though some authors calculate synthetic template spectra, real stellar spectra have been used in most works based on the CCF (see White & Basri 2003; Mohanty & Basri 2003; BailerJones 2004). Since the CCF contains information from the template and the object spectrum, using an observed template could introduce an external error source from the unknown broadening factors present in the template spectral lines. Nowadays, the best methods dealing with stars that show intense line blending in their spectra require significant computational resources.
Motivated by the need for a precise and expeditious technique to be applied extensively for the construction of a catalog of rotational velocities of bright Atype stars (Levato et al., in prep.), we develop here an alternative method based on the CCF and, at the same time, independent of any external calibration. Our procedure uses the FT to measure the parameter vsini, taking the dependence of the transform with limb darkening into account. In Sect. 2 a full description of the methodology is presented. Limb darkening consideration and other practical details of the procedure are explained in Sect. 3. Section 4 describes specific application to A and latertype stars. The precision of the obtained rotational velocities is discussed in Sect. 5, and our main conclusions are summarized in Sect. 6.
2. Method
As a first approximation, according to the classical model of a rotating star revised by Collins & Truax (1995; see also Brown & Verschueren 1997), we can consider an observed stellar line profile as the convolution of an intrinsic line profile and a rotational broadening function. Then, the observed spectrum D(λ) can be approximated by a convolution between a template spectrum T(λ) and the rotational broadening function G(λ),
where T(λ) includes any other broadening effect different from rotation: natural line broadening, thermal broadening, microturbulence, Stark effect, etc.
Under the previous assumption the CCF between the object spectrum and a template of the same spectral type, but without rotation, results in (1)where ∥ T ∥ and ∥ D ∥ are the norm of each spectrum. Since the rotational velocity is the main broadening agent of metallic lines in stars with vsini > 10 km s^{1}, i.e. in most normal mainsequence stars (Abt & Morrell 1995; Abt et al. 2002), the intrinsic width of a line profile σ is usually σ ≪ vsini.
The autocorrelation function CCF_{TT} presents a narrow peak centered at zero, which for the current application can be considered as a Gaussian whose width is times the width of the template lines. Side lobes are present in this function, but they are usually very small in comparison with the central peak since their intensity is inversely proportional to the number of spectral lines involved in the crosscorrelation. Therefore, regardless the shape of the line profile of the template, the function CCF_{DT} is very similar to the rotation function G and can be used to derive the parameter vsini.
The idea of the proposed method is to calculate the FT of the central maximum of the CCF_{DT} to derive vsini from the position of its first zero. In fact, CCF_{DT} is essentially the function G convolved with a narrow profile that has no significant impact on the position of the first zero of the FT. The central maximum of the CCF has the same information about the rotational velocity as a single line profile. The main benefit of utilizing the CCF is that having to select lines without blending, which represents a problem in stars later than lateA earlyF, is avoided.
The CCF peak, however, might be blended with side lobes, but the more lines used to calculate the CCF, the bigger the intensity difference between the central maximum and the secondary maxima. In the proposed method, we include a procedure to remove the side lobes contribution as described in Sect. 2.1.
Finally, the S/N of the CCF peak is much higher than that of a single line. Therefore, using the CCF calculated from a large spectral region instead of individual line profiles presents advantages over both line blending and S/N. Nevertheless, there are some restrictions on the size of the spectral region to be correlated owing to the variation in the limb darkening coefficient with wavelength, which will be discussed in the next section.
2.1. Extraction of the rotational profile
The CCF_{DT} would be equal to the rotational profile G only in the ideal case in which CCF_{TT} has a single peak of negligible width (a Dirac δ function). In practice, however, besides the finite width of the central peak, the CCF_{DT} contains several small subsidiary peaks coming from the fortuitous coincidence of each spectral line of the object spectrum with different spectral lines in the template as the former is shifted with respect to the latter. These small peaks can be blended with the central peak due to the rotational broadening of the object, affecting the determination of the base line of the CCF central maximum and hence the measurement of the rotational velocity.
To overcome this problem, a key stage in the process is the remotion of secondary maxima. To this aim we consider the function CCF_{TT} as the sum of two functions: CCF1_{TT}, representing only the central peak, and CCF2_{TT}, which includes all the subsidiary peaks. Then, from Eq. (1)
To calculate the contribution of the subsidiary peaks, we remove the central maximum from the CCF_{TT} to obtain CCF2_{TT}, which is convolved with a first approximation of the rotational broadening profile G_{1} extracted from the central peak of the CCF_{DT}. The result is then subtracted from the CCF_{DT}, thereby improving the determination of the base of the rotational profile. Considering the central peak of the CCF1_{TT} to be much narrower than G, we find from the previous equation (2)Before applying the subtraction, two corrections are required: a radial velocity correction for both CCFs to be centered, and a scale correction to account for an eventual global line intensity difference between the object and the template spectra. As illustration, Fig. 1 shows the application of this process to a spectrum of the A2type star HR 892 with a synthetic template for T = 9000 K.
Fig. 1 Calculation of the rotational profile using crosscorrelations. Upper panel: object spectrum and template spectrum. Lower panel, from top to bottom: a) crosscorrelation function templatetemplate CCF_{TT} (thin line) and the same function after removing the main peak and convolving with the provisional rotational profile, i.e. CCF2_{TT} ∗ G (thicker line), b) the same function CCF2_{TT} ∗ G scaled to the object intensity, c) crosscorrelation function objecttemplate CCF_{DT}, and d) rotational profile calculated subtracting CCF2_{TT} ∗ G from CCF_{DT}. In both panels arbitrary vertical shifts have been applied for clarity. 
Finally, once the base line of the central maximum has been improved, we fit the background in the surrounding region of this peak and subtract the fit to have the base of the rotational profile at zero. Through this process the rotational broadening profile G is obtained.
In principle, reliable and precise rotational profile can be reconstructed from the CCF, as long as the rotational broadening is not comparable to other photospheric or instrumental broadening effects. In Fig. 2 we compare the theoretical rotational profile of Eq. (3) with the profile recovered from the CCF of a 100 Å long spectrum (5350–5450 Å). As shown in the lower panel, the only noticeable difference is the smoothing of the sharp cut at the edges, owing to the intrinsic spectral line width (0.08–0.16 Å). The example corresponds to an atmosphere model with T_{eff} = 8000 K and log g = 4.0, convolved with a rotational profile of vsini = 60 km s^{1} and ε = 0.6.
Fig. 2 Rotational profile recovered from the CCF. Upper panel: the functions involved in Eq. (2), from top to bottom: CCF_{DF}, CCF2_{TT}*G_{1}, and the recovered rotational profile G. Lower panel: comparison of the retrieved rotational profile (filled circles) with the input theoretical profile (solid line). 
2.2. Calculation of v sin i
The rotational profile of a spectral line centered on a wavelength λ_{o} for a spherical star rotating as a rigid body and whose limb darkening law is linear with a coefficient ε, is (3)for x < 1, and G(x) = 0 for x > 1, where x = ln(λ/λ_{o})·c/(vsini). Expressing the rotational profile in terms of Δlnλ instead of Δλ has the advantage of being independent of the line wavelength (cf. with Eq. (17.12) in Gray 1992; Eq. (4) in Reiners & Schmitt 2002). In addition, the Doppler formula in logarithmic scale, Δlnλ ≈ v/c, is a better approximation than the classical Δλ/λ ≈ v/c.
The FT of G(x) presents zeros at different positions σ_{n}, which are related to vsini by (4)where k_{n}(ε) are functions of the limb darkening coefficient (Reiners & Schmitt 2002). For the first zero, k_{1} is given implicitly by the following expression:
where J_{1} is the firstorder Bessel function and β = 2πk_{1}. In our procedure we used the following approximate formula
which is accurate enough (better than 0.01%) for the whole range of limb darkening coefficient of normal mainsequence stars: ε = 0.00–1.10. We note that this cubic is not a series expansion around ε = 0 but a fit of the function in this interval. In fact, it is more accurate than the fourthdegree polynomial calculated by Dravins et al. (1990,Eq. (6); Eq. (8) in Reiners & Schmitt 2002 in the range of stellar limb darkening coefficients.
Once k_{1}(ε) has been calculated, the determination of vsini results from Eq. (4). We calculate the power spectrum of the extracted rotational broadening function to evaluate the position of the first zero σ_{1} of the FT and calculate vsini from Eq. (4). As mentioned before, this procedure is valid for stars rotating rigidly. The eventual presence of differential rotation might be detected by studying the shape of the derived rotational profile or using the first two zeros of the FT (Reiners & Schmitt 2002).
Fig. 3 Fourier transform of the CCF of 5 synthetic spectra (Δλ = 120 Å) with T_{eff} = 9000 K, vsini = 40 km s^{1} and convolved with the following Voigt profiles: thick solid line: original rotational profile. Thin solid line: gfwhm = 0.2 Å. Dotted line: gfwhm = 0.1 Å and lfwhm = 0.1 Å. Dashed line: gfwhm = 0.4 Å. Dotdashed line: gfwhm = 0.2 Å and lfwhm = 0.2 Å. 
Finally, it is important to mention that a linear limb darkening law was used in this work, but it is possible to consider other types of limb darkening law like the ones compared by Brown & Verschueren (1997). We estimate that a linear law could produce errors of ~1%; therefore, a more adequate limb darkening law is a potential improvement for the future.
3. Procedure and assumptions
The procedure was programmed as an IRAF task divided in five stages:

1.
calculation of the CCF between an object spectrum and a template spectrum;

2.
extraction and cleaning of the CCF central maximum;

3.
determination of the first zero of the FT of the CCF maximum;

4.
calculation of vsini;

5.
measurement error estimation.
3.1. Influence of nonrotational broadening effects
We evaluated the impact of nonrotational broadening effects on the position of the first zero and on the shape of the main lobe of the FT using a synthetic spectrum for a T_{eff} = 9000 K atmosphere model in the wavelength range 4480–4600 Å. This spectrum was convolved with a rotational profile of vsini = 40 km s^{1} and different Voigt profiles to simulate additional nonrotational broadening. The FT of the various broadening profiles, derived using the CCF, are compared in Fig. 3. We used Voigt profiles corresponding to different combinations of Gaussian (gfwhm) and Lorentz (lfwhm) components’ fullwidthhalfmaximum (FWHM). The change in the general shape of the first lobe with the nonrotational broadening effects is more significant than the shift of the first root, which in the four examples lies within 0.4 km s^{1} from the original rotational profile. This has been already noted by Reiners & Schmitt (2002).
We estimate that acceptable vsini measurements are obtained as long as the rotational profile is at least twice wider than the nonrotational profile. In mainsequence stars, typically the intrinsic FWHM of metallic lines is 0.07–0.12 Å. Therefore, for highresolution spectra, the lower limit for vsini measurements through the zero of the FT is about vsini = 5–8 km s^{1}. Btype stars are more difficult to measure since their most conspicuous spectral features are He i lines, whose intrinsic widths are usually on the order of 0.8 Å. Such lines would not be suitable for rotation measurement with this technique, except for very fast rotators (vsini ≳ 100 km s^{1}). In fact, these lower limits depend on the S/N, since the effect of the intrinsic line broadening is to reduce the intensity of the first subsidiary lobe, making it more difficult to determine the first zero position.
All in all, for the usual values of S/N and instrumental broadening, the variation in the first zero position caused by additional broadening and noise is below 1%. A typical case is shown in Fig. 4, where the FT of 20 spectra with different random noise corresponding to S/N = 100 are plotted. The atmospheric parameters, vsini, and spectral range are the same as for Fig. 3, and the assumed instrumental profile is 6.6 km s^{1}. We measured ⟨ vsini ⟩ = 39.88 km s^{1} with σ = 0.26 km s^{1} (0.7%) for S/N = 100. Using S/N = 50 the typical error is σ = 1.3%. Considering that we used just a 120 Å wide region, we conclude that, with medium quality spectra (e.g. Δλ ≈ 1000 Å and S/N = 50 or Δλ ≈ 100–200 Å and S/N = 100), errors are well below 1%.
Fig. 4 Fourier transform of the CCF of 20 synthetic spectra (Δλ = 120 Å) with T_{eff} = 9000 K, vsini = 40 km s^{1} and random noise simulating S/N = 100. 
In light of these results and under the assumption of rigid rotation, we consider that a global fitting of the FT of the rotational profile does not present much of an advantage for measuring rotational velocity, since the general shape of the function depends on the nonrotational broadening effects much more strongly than the position of the first root does, and therefore, additional free parameters should be included in the fitting in order to model the nonrotational contributions simultaneously. Moreover, even though the S/N is evidently lower around the first root, it is usually high enough for this purpose, since the CCF peak has a much higher S/N than individual spectral lines.
We note that the FT zero position is even less sensitive to nonrotational broadening than the BesselFourier transform proposed by Piters et al. (1996). In fact, according to Fig. 2 of Piters et al. (1996), the maximum of the BesselFourier transform is shifted by about 1% (2%) when an additional Gaussian broadening of gfwhm = 0.14 × vsini (0.28 × vsini) is present. These values approximately correspond to 0.2 and 0.4 Å in the calculations of our Fig. 3, for which we found a shift in the FT zero of only 0.2% (0.4%).
3.2. Considerations on limb darkening
The influence of the limb darkening coefficient on the rotational broadening function G(λ) is well known (see Reiners & Schmitt 2002; Collins & Truax 1995). Although the errors can rise to nearly 16%, most authors assume a fixed value of ε = 0.6 (SimónDiaz & Herrero 2007; SimónDiaz et al. 2006; Reiners & Royer 2004; Royer et al. 2002a,b; Collins & Truax 1995; Ramella et al. 1989), even for lowmass stars (spectral types M or L) and brown dwarfs (White & Basri 2003; Mohanty & Basri 2003; BailerJones 2004; Tinney & Reid 1998).
It is also worth noting that the limb darkening coefficient in the core of the lines may differ substantially from that of the nearby continuum. Collins & Truax (1995) treated this problem and indicate that the difference may be important. For a nonrotating B9V model they calculated ε = 0.57 for the continuum at λ = 4475 Å and ε = 0.29 for the core of the Mg ii line at λ4481.13. Even though the limb darkening coefficient varies within the intrinsic line profile, it is valid to define the rotational profile G(x) as in Eq. (3) using an effective limb darkening coefficient (e.g. weighted average over the line profile), provided the rotational broadening is much larger than the intrinsic line profile. Our method does not use single lines but spectral regions that include absorption lines with ε values that differ from the values of ε for the continuum in different amounts. This can still be considered by defining an effective limb darkening coefficient for each spectral range. However, we do not consider this issue in the present paper. Instead, we split the spectrum to perform the measurements in regions of about 200–400 Å and adopt the limb darkening coefficient for the continuum at the central wavelength of each region as representative of the spectral region.
Moreover, owing to the dependence of ε on temperature and wavelength, the parameter k_{n}(ε) in Eq. (4) varies with spectral type and along the spectrum. To account for this effect, we include a calibration k_{1}(λ,T_{eff}) (see next section) specifically obtained and tested for A type stars. The purpose of having an empirical calibration for a particular range of temperatures or spectral type is to maximize the precision and to keep the procedure as simple as possible. Another effect that has not been considered is the gravity darkening, which is important for vsini larger than 100 km s^{1}. Collins & Truax (1995) also discusses this effect and warns that the error may be as large as 10% for vsini larger than 200 km s^{1}. Consideration of this effect may also be a second step towards improving the method.
3.3. Template spectrum and spectral region selection
The template spectrum must be morphologically similar to the object spectrum but with zero rotational velocity. As the maximum of the CCF is the result of the contribution of all the coincident spectral lines between both spectra, slight differences resulting from small discrepancies in spectral type, metallicity, or the presence of spectral peculiarities have little influence on the CCF maximum. We tested the incidence of using templates that differ morphologically from the target spectrum, measuring a synthetic spectrum for an A7V star (T_{eff} = 8000 K and log g = 4.0) and vsini = 60 km s^{1} with various templates. We measured six spectral regions ≥ 240 Å in the range 3985–5870 Å with 19 main sequence templates ranging from T_{eff} = 5000 to 14 000 K, and we found that the error produced by spectral type mismatch is below 1% for templates in the range A0–G0 (T_{eff} = 6000–10 000 K) as shown in Fig. 5. Although the present work is focused on Atype stars, from Fig. 5 it is also evident that earlier spectral types are more sensitive to this effect, so this analysis should be carried out for the spectral range of interest in future works.
We also used an A5V template to measure A0 and A9 stars in our program and a dwarf template to measure giants and supergiants of the same type. In all cases, the differences in the vsini are always below 1%. Such a small error shows the strength of the CCF, and that spectral type mismatch is not the main source of error in Atype stars. For this error to be comparable to the dispersion of values obtained from different regions in the same spectrum, the error in the template temperature has to be ≳ 2000 K.
Fig. 5 Error in ⟨ vsini ⟩ from the mismatch of spectral type between the object and the template spectrum. The xaxis has the temperature of the template. The temperature of the object is indicated with the dashed line (T_{eff} = 8000 K). 
The selection of the spectral regions to be correlated depends on the quality and the wavelength range of the observed spectra. Since the central wavelength of each region is used to calculate k_{1}, we consider 500 Å as the maximum region width to minimize the influence of the limb darkening coefficient on k_{1}. Typically, in an interval of 500 Å, the limb darkening coefficient varies within ± 0.015, which affects the measured rotational velocity in about ± 0.25%. In practice, the resulting error would be significantly lower since it depends on the difference between coefficient for central wavelength (used for the calculations) and the effective coefficient, which depends on the spectral line distribution in the region under consideration.
In addition, since the method is based on the rotational velocity as the main broadening factor in all spectral lines, any line strongly deviated from a rotational profile must be avoided. Otherwise, it will introduce an unwanted distortion into the CCF. Even though some of them might be known a priori, e.g., He i in early B stars or Ca ii in late A to F stars, it is convenient to evaluate whether other regions introduce a significant distortion into the CCF. The crosscorrelation between two template spectra in different spectral regions can be used to detect regions that introduce significant distortions in the CCF.
3.4. Single measurement error calculation
The error assigned to a single measurement is calculated from the intensity and the FWHM of the CCF central maximum, and the noise present in the CCF. To calibrate the error formulae we used synthetic spectra artificially broadened with six different rotational profiles and random noise simulating S/N in the range S/N = 70–200. Then, ten spectra were generated for each pair of parameters and four regions were measured in each spectrum. For each region the standard deviation of the ten measures was adopted as the true error.
Fig. 6 Percentage error on ⟨ vsini ⟩ from single regions of artificially broadened spectra. For each velocity two values of S/N are plotted. S/N = 100 and 200 for ⟨ vsini ⟩ = 80, 180, and 280 km s^{1}. S/N = 100 and 70 for ⟨ vsini ⟩ = 20, 120, and 220 km s^{1}. Triangles: standard deviation of 10 consecutive measures (true error). Circles: designated error (calibration). 
Finally, we calibrated the error as a function of the intensity I (height of the CCF peak), the FWHM, and the noise of the CCF (rms). As a result, we obtained the following expression to calculate the error of a single region measurement: Error(km s^{1}) = 4.42 FWHM^{0.520}·rms·I^{1.08}, valid for spectra with S/N = 70–200, the range used for the calibration. Figure 6 shows no systematic difference among real error and its calibration.
4. Application
4.1. Application to Atype stars
As part of our current research on rotational velocity, we have been engaged for several years in a program to measure vsini for all southern Atype stars of the Bright Star Catalogue (BSC). This sample includes more than 800 stars.
To test our method, we measured ⟨ vsini ⟩ for 251 stars in the spectral range A1–A5^{1}, the results of which are presented in Table 1. Among them 155 stars were also measured by Royer et al. (2002a). In a forthcoming paper we will publish the complete results of vsini measurements for the whole sample of southern Atype stars of the BSC.
The southern Atype stars of the BSC were observed spectroscopically using the bench echelle spectrograph (EBASIM) fed with fibers coming from the Cassegrain focus of the 2.1 m telescope at CASLEO. The spectrograph has been described by Pintado & Adelman (2003). The resolving power of the spectra varies from ~30 000 in the blue end to about 25 000 in the red end.
The observing material was reduced using IRAF packages. Templates of spectral type A1–A5 with resolving power 500,000 were selected from the Kurucz online database (Kurucz 1991) to measure vsini in objects with luminosity class VIV, independently of the presence of spectroscopic peculiarities. For Am stars we used the template corresponding to the spectral type of the metallic spectrum.
We divided our spectral range into five regions of different sizes according to the density of spectral lines. Near the blue end of the spectra, regions of 250 Å were sufficient, while in the red end 400 Å were necessary to obtain a welldefined central maximum in the CCF. We avoided a rather small region between 4690 and 4730 Å in all spectral types, because it introduces a triangular base in the CCF maximum, which makes it impossible to unambiguously identify the first zero in Fourier space. For the same reason all lines from the Balmer series were also excluded.
To simplify the calculation of the appropriate limb darkening coefficient for each star and wavelength region, we implemented a calibration k_{1}(λ,T_{eff}), where λ is the central wavelength of the spectral region. We used the limb darkening coefficients tabulated by Claret (2000) and DíazCordovés et al. (1995) for various photometric bands. Then, we fitted ε(T_{eff}) for our temperature range of interest, i.e. T_{eff} = 7500 K–11 500 K (spectral types B9–A9), for log g = 4.0 and for the filters B,V, v, b, and y, whose central wavelength are within the range of our spectra (4000 Å–6000 Å). The residuals of the fit for each filter lie within ± 0.012 (rms = 0.006), which represents an error in vsini ~ 0.15%. The final wavelength calibration was made as a function of the central wavelength of each filter, giving ε(T_{eff},λ) as a result. By means of this calibration it is possible to calculate the k_{1} value appropriate for the star’s T_{eff} and the central wavelength of spectral region.
Results for 251 A1–A5 stars from the Bright Star Catalogue in the southern hemisphere.
The number of regions used to determine vsini depends on the quality of the spectra and the wavelength range, and is totally independent of vsini. The observational material analyzed here consist of 32 spectra in the range 4000–6000 Å in which we measured five regions, 116 spectra in the range 4000–5500 Å in which four regions were measured, 99 spectra in the range 3850–5000 Å or in the range 4000–5250 Å in which three regions were measured, and finally four objects in which only two regions were measured owing to excessive noise in part of the spectra.
Once vsini and its error were obtained in all the regions of the same spectrum, the weighted mean value of vsini was calculated, along with two different estimates of its uncertainty, ξ_{1} and ξ_{2}, following the formulae proposed by González & Lapasset (2000, Eq. (1)) for radial velocities.
The error ξ_{1} is the standard deviation of the weighted mean of a sample of n uncorrelated observations with standard deviations e_{i}: . On the other hand, ξ_{2} is a generalization of the standard error of the mean σ/ for a set of measurements weighted according to their individual errors. Therefore, ξ_{2} is calculated from the dispersion of vsini from different spectral regions, while ξ_{1} is computed from the single measurement errors estimated through the calibration described in Sect. 3.4. In general, ξ_{1} would be more appropriate when the number of measurements is small (2–3 spectral regions), and consequently the dispersion of measurements is a less precise estimate of the true error. In Fig. 7 we have plotted the values of ξ_{1} and ξ_{2} computed for each object. Excluding slow rotators (vsini below 30 km s^{1}) whose line widths are close to the limit imposed by the spectral resolution (c/R = 10–12 km s^{1}), the error ξ_{1} has an average value of 1.1%, always smaller than 2.5%. Regarding the dispersion of the values obtained from different regions of the same spectrum, only three objects with ⟨ vsini ⟩ > 30 km s^{1} were measured with ξ_{2} over 5%, and the average error ξ_{2} of these objects is 1.5%.
Fig. 7 Left: percentage error of the 251 measured stars. Asterisks: ξ_{1} (%), individual region error indicator. Triangles: ξ_{2} (%), measurements dispersion indicator. Lines represent constant values of ξ in km s^{1}. Solid: 1 km s^{1}. Dashed: 2 km s^{1}. Dotdashed: 5 km s^{1}. Dotted: 10 km s^{1}. Right: histogram of the percentage error. Solid line: ξ_{1} (%). Dotted line: ξ_{2} (%). 
Fig. 8 Position of the first zero of the FT of the rotational profiles for 8 spectral regions of 8 spectra of the star HD 77370. Dotted lines are theoretical lines corresponding to different values of vsini and have been calculated from the limb darkening coefficients in DíazCordovés et al. (1995). 
4.2. Application to late type stars
One field of application for the proposed method is the measurement of rotation in stars of spectral type F or later, where it is not possible to measure individual lines even for low rotators. As a test, we measured the rotational velocity of HD 77370, an F3 V star with projected rotational velocity of about 60 km s^{1}, high enough for all spectral lines to appear blended. Eight highresolution HARPS spectra of this star were downloaded from the ESO archive database. These spectra cover a spectral range 3800–6900 Å with a S/N ≈ 200. Rotational velocities were measured for eight spectral regions of about 200–400 Å, using as template a synthetic spectrum with T_{eff} = 6750 K and solar chemical abundances. Figure 8 shows the measured vsini values for the eight spectra.
The dependence of the position of the first zero of the FT of the rotational profile with wavelength – due to the variation of the limb darkening coefficient with wavelength – is noticeable. The rotational velocity of each spectrum was calculated as the average of the eight spectral regions. The mean value of the 64 measurements is 59.16 km s^{1}. The standard deviation of the measurements of the 8 spectra in one single spectral region is on average 0.41 km s^{1}. If we take this value as the typical uncertainty of one of the 64 individual measurement, the uncertainty of the mean vsini for each region should be times lower (0.14 km s^{1}). However, if we calculate the average for each of the eight spectral regions, then the dispersion of these mean rotational velocities is 0.37 km s^{1}, indicating that the main source of uncertainty in these measurements is not random, but small systematic differences between the spectral regions. All in all, the uncertainty of the vsini measured from a highquality spectrum of an Ftype star would be about 0.6%, demonstrating that highprecision rotational velocities can be measured even when lines are blended. The key point in this respect is the use of the templatetemplate correlation function to model the subsidiary lobes of the objecttemplate CCF during the calculation of the rotational profile. As an illustration, Fig. 9 shows the reconstruction of the rotational profile for a small spectral region (70 Å) of the spectrum of HD 77370, where all spectral lines are blended.
Fig. 9 Example of the rotational profile calculated from a highS/N spectrum of only 70 Å for an Ftype star. Upper panel: object spectrum and template spectrum. Lower panel: the same as Fig. 1. 
5. Discussion
To evaluate the performance of our method in determining vsini, we compared our results with those from Royer et al. (2002a). They used FT of line profiles to provide accurate vsini of a large sample of Atype stars in the southern hemisphere observed with similar resolving power (~28 000). Figure 10 shows Royer et al. values of vsini and the ones obtained in this work for 155 objects in common with their sample. Even though a good agreement is found in stars of moderate, projected rotational velocities, a significant deviation is noticed for vsini > 150 km s^{1} with an average difference of 5 ± 1%, our values larger than theirs.
Fig. 10 Comparison between our results and Royer et al. (2002a) for objects in common. Squares: normal stars. Circles: peculiar stars. Solid line: Normal stars fit f(⟨ vsini ⟩ ) = 1.43 ⟨ vsini ⟩ ^{0.92}. 
There are several differences between the methodology used by Royer et al. (2002a) and the present work that can account for the effect present in Fig. 10. First, they measured single lines in the range from 4200 Å to 4500 Å, i.e. a region of 300 Å wide which is the mean width of each region used in this work. Second, the blending effect becomes evident through the dependence of the number of measured lines with vsini. They used in average less than 3 lines on objects with vsini > 60 km s^{1}, and over 90 km s^{1} the average number of measured lines is ≤ 2. Even though the dependence on the position of the first zero of the Fourier transform with the limb darkening coefficient was not taken into account by Royer et al., the adopted value ε = 0.6 is typical of the limb darkening coefficient in the range of temperatures of the objects in the sample. However, an error of ± 0.01 in k_{1} represents an error of 1.5% in vsini. Thus, using a wrong k value could account for 1.5% of the deviation detected in the comparison.
Finally, the main source of error is the uncertainty on the continuum position. Royer et al. (2002a) analyzed the magnitude of this effects by broadening synthetic spectra of T_{eff} from 7500 to 10 000 K with vsini = 10, 50, and 100 km s^{1}. They found that in the coldest models with the largest vsini the error can reach 3%. This uncertainty affects spectral lines giving as a result a lower value of vsini than the real value. To evaluate the intensity of this type of systematic error with our method, we measured vsini in a synthetic A5V spectrum broadened by vsini = 20, 80, 120, 180, 220, and 280 km s^{1}. Random noise was added to simulate S/N = 100, generating ten spectra with different noise pattern for each value of vsini. Then, four regions were used to measure vsini, and the average over the ten measures for each region was calculated. As expected, we found a tendency to measure a lower value of vsini than the actual value. This effect is shown in Fig. 11. Nevertheless, with the method proposed in this work the difference between measures and real values stays under 1% () even for vsini = 280 km s^{1}. In conclusion, a considerable fraction of the difference between our results and those from Royer et al. (2002a) can be attributed to the effect of line blending on the continuum position. The influence of the line blending in high rotational velocity stars is significantly less with the proposed method than with Fourier transform of singleline profiles, by virtue of subtracting the secondary maximums from the CCF before calculating the Fourier transform. Similar results would be obtained by removing the blended lines with a similar procedure in the spectrum: modeling the smaller lines by convolving a rotational profile with a synthetic spectrum in which the line of interest has been subtracted.
Fig. 11 Difference between real vsini and average over ten consecutive measures with distinct random noise. Dotted line: f(⟨ vsini ⟩ ) = 0.0032 ⟨ vsini ⟩ + 0.0456. 
The main originality of the proposed method is the way in which the broadening function is built from the CCF, which, on the one hand, allows measurement of stars with blended spectral lines and, on the other, improves the S/N. Once the broadening function has been obtained, we determine the stellar rotation through the first zero of the FT. However, other alternative analyses are possible, such as direct fitting of the broadening function or FourierBessel transformation (Piters et al. 1996).
As an illustration, we used our method in combination with the FourierBessel transformation to measure one of the objects in Piters et al. (1996) for which a HARPS spectrum was publicly available. These authors applied the FourierBessel transformation method to individual spectral lines for measuring rotation in Ftype stars. Owing to line blending, stars with vsini > 100 km s^{1} could not be measured or had errors of about 10–20%. We selected one of their fast rotating objects: HD 12311, for which they measured 110 ± 22 km s^{1}. We defined 5 spectral region between 3980 Å and 6530 Å obtaining an average of 157.4 km s^{1} with σ = 4.6 km s^{1} using our method with the FT first zero, and 155.9 km s^{1} with σ = 1.8 km s^{1} from our broadening function but using the FourierBessel transform. Therefore, the application of FourierBessel transform method to the broadening function derived from the CCF would be an excellent alternative for latetype fastrotating stars. In the case of slowly rotating stars, however, the FourierBessel transform is more sensitive to broadening effects that are different from rotation, so the zero of the FT would give better results (see Sect. 3.1).
6. Summary
We have developed a method for measuring vsini based on the FT of the central maximum of the CCF. This combination provides a simple and precise solution to the line blending problem with mediumresolution spectra and/or spectral types later than midA. As a result, the number of useful spectral regions per object with our method is independent of spectral type and vsini.
The high precision of the proposed method is supported by two key features of the procedure. The first one is to use an empirical calibration to take limb darkening effects into account in the position of the first zero of the FT, which is directly related to vsini. The second key feature is the subtraction of the subsidiary lobes of the CCF during the calculation of the rotational broadening function. This assures a good definition of the rotational profile even in presence of line blending.
The vsini value is derived from the first zero of the FT of broadening function. However, the S/N of the resulting rotational profile is high enough to also do other detailed studies of the shape of the reconstructed rotational profile.
For slowly rotating stars, a lower limit for measurable vsini is imposed by the instrumental profile and the intrinsic line width, which is included twice in our broadening function owing to the use of a template with lines of finite width. For metallic lines this limit is 5–8 km s^{1}. In earlytype stars the measurement is more difficult since metallic lines become weaker and He lines are usually not useful. Morphological differences between template and object are also less significant in latetype stars. For lateA type stars or later, differences in several subtypes have no impact on the rotation measurements.
We consider that, in accurate vsini determinations, one significant contribution to the errors is the adopted value for the limb darkening. Even though in this work we use a wavelengthdependent coefficient, the continuum limb darkening might differ significantly from the line limb darkening. Using many spectral lines in our calculations, this difference might be averaged out, but systematic differences between lines and continuum might result in small but systematic vsini errors. This is an issue that deserves further research.
We applied the proposed method to a sample of 251 Atype stars. Measurement errors are always under 2.5%, being 1.1% the average error for stars with ⟨ vsini ⟩ over 30 km s^{1}. As regards the dispersion of values measured on different regions across a spectral range of ~1150 Å to 2000 Å wide, we found standard errors of the mean below 5% with an average of 1.5%.
The spectral types were taken from the SIMBAD astronomical database: http://simbad.ustrasbg.fr/simbad
Acknowledgments
We thank the night assistants of CASLEO who helped during the observing procedures. Part of this research was supported by a grant from CONICET PIP 1113. We gratefully acknowledge the use of ESO archival data.
References
 Abt, H., & Morrell, N. 1995, ApJS, 99, 135 [Google Scholar]
 Abt, H., Levato, H., & Grosso, M. 2002, ApJ, 573, 359 [NASA ADS] [CrossRef] [Google Scholar]
 Bruning, D. H. 1981, ApJ, 248, 274 [NASA ADS] [CrossRef] [Google Scholar]
 BailerJones, C. A. L. 2004, A&A, 419, 703 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Brown, A. G. A., & Verschueren, W. 1997, A&A, 319, 811 [NASA ADS] [Google Scholar]
 Carroll, J. A. 1928, MNRAS, 88, 548 [NASA ADS] [CrossRef] [Google Scholar]
 Carroll, J. A. 1933, MNRAS, 93, 478 [NASA ADS] [CrossRef] [Google Scholar]
 Carroll, J. A., & Ingram, L. J. 1933, MNRAS, 93, 508 [NASA ADS] [Google Scholar]
 Claret, A. 2000, A&A, 363, 1081 [NASA ADS] [Google Scholar]
 Collins, G. W. 2004, Proc. IAU Symp., 215, on Stellar Rotation, 3 [Google Scholar]
 Collins, G. W., & Truax, R. J. 1995, ApJ, 439, 860 [NASA ADS] [CrossRef] [Google Scholar]
 DíazCordovés, J., Claret, A., & Giménez, A. 1995, A&AS, 110, 329 [NASA ADS] [Google Scholar]
 Dravins, D., Lindergren, L., & Torkelsson, U. 1990, A&A, 237, 137 [NASA ADS] [Google Scholar]
 García, B., & Levato, H. 1984, RMAA, 9, 9 [Google Scholar]
 González, J. F., & Lapasset, E. 2000, AJ, 119, 2296 [NASA ADS] [CrossRef] [Google Scholar]
 Gray, D. F. 1977, ApJ, 211, 198 [NASA ADS] [CrossRef] [Google Scholar]
 Gray, D. F. 1982, ApJ, 258, 201 [NASA ADS] [CrossRef] [Google Scholar]
 Gray, D. F. 1992, The observation and analysis of stellar photospheres, second edition (Cambridge University Press) [Google Scholar]
 Kurucz, R. L. 1991, Harvard Preprint, 3348 [Google Scholar]
 Mohanty, S., & Basri, G. 2003, ApJ, 583, 451 [NASA ADS] [CrossRef] [Google Scholar]
 Pintado, O. I., & Adelman, S. J. 2003, A&A, 406, 987 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Piters, A. J. M., Groot, P. J., & van Paradijs, J. 1996, A&AS, 118, 529 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ramella, M., Gerbaldi, M., Faraggiana, R., & Böhm, C. 1989, A&AS, 209, 233 [NASA ADS] [Google Scholar]
 Reiners, A., & Royer, F. 2004, A&A, 428, 199 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Reiners, A., & Schmitt, J. H. M. M. 2002, A&A, 384, 155 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Reiners, A., & Schmitt, J. H. M. M. 2003, A&A, 412, 813 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Royer, F., Gerbaldi, M., Faraggiana, R., & Gómez, A. E. 2002a, A&A, 381, 105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Royer, F., Grenier, S., Baylac, M. O., Gómez, A. E., & Zorec, J. 2002b, A&A, 393, 897 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 SimónDíaz, S., & Herrero, A. 2007, A&A, 468, 1063 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 SimónDíaz, S., Herrero, A., Esteban, C., & Najarro, F. 2006, A&A, 448, 351 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Shajn, G., & Struve, O. 1929, MNRAS, 89, 222 [NASA ADS] [Google Scholar]
 Slettebak, A. 1949, ApJ, 110, 498 [NASA ADS] [CrossRef] [Google Scholar]
 Slettebak, A. 1954, ApJ, 119, 146 [NASA ADS] [CrossRef] [Google Scholar]
 Slettebak, A. 1955, ApJ, 121, 653 [NASA ADS] [CrossRef] [Google Scholar]
 Slettebak, A. 1956, ApJ, 124, 173 [NASA ADS] [CrossRef] [Google Scholar]
 Slettebak, A. 1982, ApJS, 50, 55 [NASA ADS] [CrossRef] [Google Scholar]
 Slettebak, A., & Horward, R. F. 1955, ApJ, 121, 102 [NASA ADS] [CrossRef] [Google Scholar]
 Slettebak, A., Collins, G. W., Boyce, P., White, M., & Parkinson, T. 1975, ApJS, 29, 137 [NASA ADS] [CrossRef] [Google Scholar]
 Tinney, C. G., & Reid, I. N. 1998, MNRAS, 301, 1031 [NASA ADS] [CrossRef] [Google Scholar]
 White, R., & Basri, G. 2003, ApJ, 582, 1109 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
Results for 251 A1–A5 stars from the Bright Star Catalogue in the southern hemisphere.
All Figures
Fig. 1 Calculation of the rotational profile using crosscorrelations. Upper panel: object spectrum and template spectrum. Lower panel, from top to bottom: a) crosscorrelation function templatetemplate CCF_{TT} (thin line) and the same function after removing the main peak and convolving with the provisional rotational profile, i.e. CCF2_{TT} ∗ G (thicker line), b) the same function CCF2_{TT} ∗ G scaled to the object intensity, c) crosscorrelation function objecttemplate CCF_{DT}, and d) rotational profile calculated subtracting CCF2_{TT} ∗ G from CCF_{DT}. In both panels arbitrary vertical shifts have been applied for clarity. 

In the text 
Fig. 2 Rotational profile recovered from the CCF. Upper panel: the functions involved in Eq. (2), from top to bottom: CCF_{DF}, CCF2_{TT}*G_{1}, and the recovered rotational profile G. Lower panel: comparison of the retrieved rotational profile (filled circles) with the input theoretical profile (solid line). 

In the text 
Fig. 3 Fourier transform of the CCF of 5 synthetic spectra (Δλ = 120 Å) with T_{eff} = 9000 K, vsini = 40 km s^{1} and convolved with the following Voigt profiles: thick solid line: original rotational profile. Thin solid line: gfwhm = 0.2 Å. Dotted line: gfwhm = 0.1 Å and lfwhm = 0.1 Å. Dashed line: gfwhm = 0.4 Å. Dotdashed line: gfwhm = 0.2 Å and lfwhm = 0.2 Å. 

In the text 
Fig. 4 Fourier transform of the CCF of 20 synthetic spectra (Δλ = 120 Å) with T_{eff} = 9000 K, vsini = 40 km s^{1} and random noise simulating S/N = 100. 

In the text 
Fig. 5 Error in ⟨ vsini ⟩ from the mismatch of spectral type between the object and the template spectrum. The xaxis has the temperature of the template. The temperature of the object is indicated with the dashed line (T_{eff} = 8000 K). 

In the text 
Fig. 6 Percentage error on ⟨ vsini ⟩ from single regions of artificially broadened spectra. For each velocity two values of S/N are plotted. S/N = 100 and 200 for ⟨ vsini ⟩ = 80, 180, and 280 km s^{1}. S/N = 100 and 70 for ⟨ vsini ⟩ = 20, 120, and 220 km s^{1}. Triangles: standard deviation of 10 consecutive measures (true error). Circles: designated error (calibration). 

In the text 
Fig. 7 Left: percentage error of the 251 measured stars. Asterisks: ξ_{1} (%), individual region error indicator. Triangles: ξ_{2} (%), measurements dispersion indicator. Lines represent constant values of ξ in km s^{1}. Solid: 1 km s^{1}. Dashed: 2 km s^{1}. Dotdashed: 5 km s^{1}. Dotted: 10 km s^{1}. Right: histogram of the percentage error. Solid line: ξ_{1} (%). Dotted line: ξ_{2} (%). 

In the text 
Fig. 8 Position of the first zero of the FT of the rotational profiles for 8 spectral regions of 8 spectra of the star HD 77370. Dotted lines are theoretical lines corresponding to different values of vsini and have been calculated from the limb darkening coefficients in DíazCordovés et al. (1995). 

In the text 
Fig. 9 Example of the rotational profile calculated from a highS/N spectrum of only 70 Å for an Ftype star. Upper panel: object spectrum and template spectrum. Lower panel: the same as Fig. 1. 

In the text 
Fig. 10 Comparison between our results and Royer et al. (2002a) for objects in common. Squares: normal stars. Circles: peculiar stars. Solid line: Normal stars fit f(⟨ vsini ⟩ ) = 1.43 ⟨ vsini ⟩ ^{0.92}. 

In the text 
Fig. 11 Difference between real vsini and average over ten consecutive measures with distinct random noise. Dotted line: f(⟨ vsini ⟩ ) = 0.0032 ⟨ vsini ⟩ + 0.0456. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext 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 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.