Issue 
A&A
Volume 602, June 2017



Article Number  A83  
Number of page(s)  21  
Section  Stellar atmospheres  
DOI  https://doi.org/10.1051/00046361/201628761  
Published online  19 June 2017 
Critical study of the distribution of rotational velocities of Be stars
II: Differential rotation and some hidden effects interfering with the interpretation of the V sin i parameter
^{1} Sorbonne Universités, UPMC Université Paris 06 et CNRS UMR 7095, Institut d’Astrophysique de Paris, 75014 Paris, France
^{2} CNRS UMR 7095, Institut d’Astrophysique de Paris, 98bis Bd Arago, 75014 Paris, France
email: zorec@iap.fr
^{3} Royal Observatory of Belgium, 3 Av. Circulaire, 1180 Bruxelles, Belgium
^{4} Université Côte d’Azur, Observatoire de la Côte d’Azur, CNRS UMR 7293, Lagrange, 28 Avenue Valrose, 06108 Nice Cedex 2, France
^{5} GEPI, Observatoire de Paris, PSL Research University, CNRS UMR 8111, Université Paris Diderot, Sorbonne Paris Cité, 5 place Jules Janssen, 92190 Meudon, France
^{6} Facultad de Ciencias Astronómicas y Geofísicas, Universidad Nacional de La Plata, Paseo del Bosque S/N, 1900 La Plata, Argentina
^{7} Instituto de Astrofísica La Plata, CONICET, 1900 La Plata, Argentina
^{8} Geneva Observatory, University of Geneva, Maillettes 51, 1290 Sauverny, Switzerland
^{9} European Organization for Astronomical Research in the Southern Hemisphere, 3107 Alonso de Cordova, Vitacura, Santiago de Chile, Chile
Received: 21 April 2016
Accepted: 9 February 2017
Aims. We assume that stars may undergo surface differential rotation to study its impact on the interpretation of Vsini and on the observed distribution Φ(u) of ratios of true rotational velocities u = V/V_{c} (V_{c} is the equatorial critical velocity). We discuss some phenomena affecting the formation of spectral lines and their broadening, which can obliterate the information carried by Vsini concerning the actual stellar rotation.
Methods. We studied the line broadening produced by several differential rotational laws, but adopted Maunder’s expression Ω(θ) = Ω_{0}(1 + αcos^{2}θ) as an attempt to account for all of these laws with the lowest possible number of free parameters. We studied the effect of the differential rotation parameter α on the measured Vsini parameter and on the distribution Φ(u) of ratios u = V/V_{c}.
Results. We conclude that the inferred Vsini is smaller than implied by the actual equatorial linear rotation velocity V_{eq} if the stars rotate with α < 0, but is larger if the stars have α > 0. For a given  α  the deviations of Vsini are larger when α < 0. If the studied Be stars have on average α < 0, the number of rotators with V_{eq} ≃ 0.9V_{c} is larger than expected from the observed distribution Φ(u); if these stars have on average α > 0, this number is lower than expected. We discuss seven phenomena that contribute either to narrow or broaden spectral lines, which blur the information on the rotation carried by Vsini and, in particular, to decide whether the Be phenomenon mostly rely on the critical rotation. We show that twodimensional radiation transfer calculations are needed in rapid rotators to diagnose the stellar rotation more reliably.
Key words: stars: emissionline, Be / stars: rotation
© ESO, 2017
1. Introduction
In Paper I of this series (Zorec et al. 2016) we have obtained the distribution Φ(V/V_{c}) of true velocity ratios V/V_{c} (V_{c} is the critical linear velocity at the equator) corresponding to a sample of 233 galactic classical Be stars. This distribution was derived from that of apparent parameters Vsini, which were corrected for measurement uncertainties, and assuming that the inclination angles i were distributed at random. Since Be stars are rapid rotators, the parameters Vsini were corrected for the Stoeckley underestimation (see definition in Paper I) in the frame of rigid rotation using the original von Zeipel (1924) theorem. In this way we ensure that the GD (GD) effect is not underestimated. Our aim was to determine an upper limit of rotational velocities of Be stars by taking some effects that can disturb the determination of the true velocity V into account. By overestimating the Stoeckley corrections on purpose, the inferred Φ(V/V_{c}) distribution has the largest possible mode, i.e., the distribution reveals the greatest possible velocity ratios V/V_{c} that might characterize the Be phenomenon. Then, we corrected the distribution of Vsini/V_{c} ratios for overestimations of the Vsini values induced by macroturbulent motions in the stellar atmospheres and distribution of ratios V/V_{c}^{1}.
Fig. 1 Curves of internal constant angular velocity Ω(r,θ) = const. in the stellar envelope (blue lines). a) Curves corresponding to the condition S = S(j^{2}); b) curves obtained with S = S(Ω^{2}); and c) curves calculated for S = S(ϖ^{2}Ω^{2}). The solutions are valid in the convective envelope, i.e., above the shaded central region (S = specific entropy; j = specific angular momentum; and ϖ = distance to the rotation axis). Ordinates coincide with the rotation axis; abscissas are in the equatorial plane. 
In Paper I we have assumed implicitly that the angular velocity Ω on the surface of stars is uniform. In the present approach, we take the surface angular velocity dependent on the colatitude θ, so that the measured Vsini depends on the characteristics of the surface rotation law Ω = Ω(θ).
The first part of the present contribution (Sect. 2) aims at exploring two questions: what is the effect of the surface differential rotation on the measured Vsini parameters? and what could be the deviation of the observed distribution of the true velocities ratios V/V_{c} from the distribution of actual equatorial rotational velocities V_{eq}/V_{c}? In the second part of this work (Sect. 3) we discuss some measuring and conceptual uncertainties that affect the determination of Vsini. These uncertainties can cast doubts on the quality of the information that this parameter may carry on the actual properties of the stellar surface rotation.
2. Effect carried by the surface differential rotation
2.1. Is differential rotation a timely topic for Be stars?
Independent of the internal rotation law of massive stars, either conservative or nonconservative “shellular”, Clement (1979) and Maeder et al. (2008) have shown that the two convection zones in the envelope associated with increased opacity due to He and Fe ionization can be considerably enlarged in depth from rapid rotation. These regions together establish an entire convective zone beneath the surface that spread out over a nonnegligible region: from 1/8 of the stellar radius at the pole to nearly 1/4 at the equator (Maeder et al. 2008). Since the rotation law in these layers can have a direct incidence on the surface rotation, we are interested in its nature. In general, two extreme approximations have currently been used to account for the angular momentum distribution in convective regions: first, redistribution promoted by the turbulent viscosity (Maeder et al. 2008), which ends up establishing rigid rotation; and, second, redistribution of specific angular momentum j = Ωϖ^{2} carried by the convective plumes (Tayler 1973), which favors the existence of rotation laws characterized by j ≃ const. (Deupree 1998, 2001). The solar convective regions, however, which are not only characterized by significant turbulence but also rotate differentially (Schou et al. 1998), do not conform with any intermediate frame between these two extreme possibilities.
Fig. 2 Curves of constant radial velocity C contributing to the rotational Doppler broadening of spectral lines. a) Ccurves for a Maunderlike surface angular velocity law with ∂Ω(θ) /∂θ > 0; b) Ccurves for a law given by Eq. (3) with γ > 0, so that ∂Ω(θ) /∂θ > 0; c) Ccurves for a law given by Eq. (3) with γ < 0, so that ∂Ω(θ) /∂θ < 0. In cases a) and b) the largest Doppler displacement is produced in the limb at the equator, while in case c) is produced in the center of the upper closed curves. Case b) identifies the lines producing the same Doppler displacement, C_{r} = C_{d}, where C_{r} (red dotted line) indicates rigid surface rotation, and C_{d} (red dashed line) a differential surface rotation. In all cases the parent Doppler displacement from rigid rotation is determined by the straight line touching the corresponding C_{d} curve at the equator. 
According to the solar rotational picture, where the layers are unstable to convection and rotate differentially with a nonconservative nonshellular pattern, we may wonder whether some coupling can also exist between convection and rotation beneath the surface in rapidly rotating massive and intermediate mass stars. Solving the baroclinic balance relation obtained with the curl of the timeindependent momentum equation of an inviscid, axisymmetric rotating star without magnetic fields, Zorec et al. (2011) obtained solutions for the angular velocity distribution in the envelope under several conditions: 1) the surfaces of specific entropy S are parallel to the surfaces of specific angular momentum j, S = S(j^{2}); 2) the surfaces of specific entropy and angular velocity coincide, i.e., S = S(Ω^{2}); and 3) the surfaces of specific entropy S are parallel with the surfaces of constant specific rotational kinetic energy, S = S(ϖ^{2}Ω^{2}). In Fig. 1 we show some examples of curves Ω(r,θ) =const. that obey these assumptions. Should this phenomenon take place, its imprint on the outermost stellar layers would certainly translate into an angular velocity dependent on the latitude. The existence of a possible external differential rotation in Be and other massive and intermediate mass stars was speculated so far by a number of authors (cf. Stoeckley 1968; Zorec 1986; Zorec et al. 1987; Stoeckley & Buscombe 1987; Zorec et al. 1990; Cranmer & Collins 1993; Espinosa Lara & Rieutord 2007; Zorec & Royer 2012; Zorec et al., in prep.).
From all of the above, the interest of studying the external angular velocity law resides in the fact that it may carry information on the subphotospheric stellar structure. However, as the characterization of a star also needs to account for its Vsini parameter, differential rotation affects the Doppler broadening of spectral lines in a particular way, and this has to be taken into account. In Sect. 2.2 we show that depending on the characteristics of the external rotation law, the larger contribution to the Doppler broadening can be produced anywhere in the observed stellar hemisphere, and not necessarily in the equator. Finally, rotation laws similar to those depicted in Fig. 1 not only determine the stellar external geometry, but can also induce GD laws that strongly deviate from the classical von Zeipel description (see Sect. 3.3).
2.2. Rotational line broadening
We adopt a Cartesian reference system (x,y,z) centered on the star, where the zaxis is directed positively toward the observer and the x and yaxis are in the plane of the sky. Assuming that the rotation axis is in the yzplane, the loci of points contributing to a same Doppler displacement Δλ is described by the z component of the vector Ω ∧ R,(1)where Ω(θ) is the angular velocity on the stellar surface, θ is the colatitude and R is the position vector, and i is the inclination angle measured from the stellar polar axis and the direction toward the observer. To characterize the Doppler displacements described by Eq. (1) more clearly, we normalize this relation by the displacement produced in the limb at the equator, (2)where Ω_{e} is the angular velocity at the equator and C is a constant independent of sini that represents an isoradial velocity curve. Curves C=±C are symmetric with respect to the yzplane. Considering only the x ≥ 0 side, we have 0 ≤ C ≤ C_{max}. The C= 0 curve is on the plane containing the zdirection and the rotation axis whatever the rotation law. The value C= 1 corresponds to the stellar limb in the equator (i.e., x_{max} = R_{e}) either for a rigid rotation or for Ω(θ) accelerated from the pole toward the equator. Finally, when Ω(θ) is accelerated from the equator toward the pole, it happens that C_{max} = C(x_{max}) ≥ 1 with x_{max} ≤ 1.
Fig. 3 Rotationally broadened line profiles of an isolated Gaussian line, whose intensity and equivalent width depends on the effective temperature and gravity such as a He i 4471 transition. Calculations include gravitational darkening. The parent nonrotating object has T_{eff} = 22 000 K and log g = 4.0. The angular velocity law is given by Eq. (3) with the same parameters (κ,γ,η,i) as in the corresponding block in Fig. 2, except that ξ has a wider range of values. The case ξ = 0 represents the rigid rotation and the corresponding line profile is highlighted with shading. In all blocks the ordinates are the same. 
To discuss more general rotation laws, let us assume the following expression (3)where ϖ = r(θ)sinθ is the distance to the rotation axis; r(θ) = R(θ) /R_{e} is the normalized radiusvector that describes the stellar surface, Ω_{p} is the polar angular velocity, and Ω_{e} = Ω_{p}(1 + ξ)^{γ} is the equatorial angular velocity; ξ, κ and γ are constants. We can readily see that ∂Ω /∂θ ≷ 0 when ξ × γ × κ ≷ 0, respectively. When ∂Ω /∂θ< 0, it is x_{max} = [−1 /ξ(1 + κγ)] ^{1 /κ}, which leads to C_{max} = C(x_{max}) ≥ 1 on the meridian that contains the xaxis. We also have Ω_{max} = Ω(x_{max}) = Ω_{p} [γκ/ (1 + γκ)] ^{γ} ≥ Ω_{e}, which for appropriate values of the constants (ξ,γ,κ) can imply Ω_{max}x_{max} ≥ Ω_{e}R_{e} in the first meridian somewhere above the equator. The condition Ω_{max}x_{max} ≤ Ω_{c}(x_{max})x_{max}, where Ω_{c}(x_{max}) is the local critical angular velocity, imposes constraints on the values of the constants (α,γ,κ).
Adopting κ = 2 and γ = 1, Eq. (3) becomes
(4)where Ω_{1} = Ω_{p} [1 + ξr^{2}(θ)] and ϵ = ξr^{2}(θ)/ [1 + ξr^{2}(θ)]. For slowly rotating stars r(θ) ≃ 1, ϵ and Ω_{1} can be considered constants, so that the relation in Eq. (4) becomes the known Maunder law.
When the angular velocity in the surface is uniform [Ω(θ) = constant ∀ θ], from Eq. (2) we find that the isoradial velocity curves, i.e., C= const., are straight lines parallel to the yaxis. Otherwise, for Ω = Ω(θ) the isoradial velocity curves C resemble those shown in Fig. 2. Figures 2a and b show C= constant for angular velocities accelerated from the pole toward the equator, while Fig. 2c shows a case where Ω(θ) is accelerated from the equator toward the pole. In cases “a” and “b”, the maximum Doppler displacement is produced at the stellar limb in the equator, while in case “c” this maximum is produced in the middle of the closed circles (“owl” eyes) on the visible part of the meridian contained in the plane. In this case, the curve C= 1 is also closed and contains the equatorial limb point, while all the remaining curves for C> 1 are closed.
Figure 2b shows two curves (red dotted and dashed curves) that produce the same Doppler displacement, i.e., C_{r} = C_{d}, where C_{r} denotes the radial velocity for rigid rotation (red dotted line) and C_{d} is for differential rotation (red dashed curve). By definition, both isoradial velocity curves coincide at the equator, nevertheless, the C_{d} curve is longer than for C_{r}. This is the main source of difference detected in residual intensities of absorption lines broadened by rigid and differential rotation laws.
Figure 3 shows rotationally broadened line profiles of a hypothetical Gaussian line, whose intensity and equivalent width depend on the effective temperature and gravity as the actual He i 4471 line. In all depicted cases, the GD was taken into account. The parent nonrotating object corresponds to T_{eff} = 22 000 K and log g = 4.0. The angular velocity laws used for these lines are the same as those in Fig. 2, except for the parameter ξ, which ranges from negative to positive values given in the inlaid in the lower right corners of Fig. 3. The cases with ξ = 0 are for rigid rotation and the corresponding line profiles are highlighted with shading to ease the comparison. Lines calculated with ξ ≤ 0 have full width at half maximum (FWHM) wider than FWHM_{rigid}, while for ξ> 0 the FWHM are smaller than FWHM_{rigid}.
In Fig. 3c, the line profiles with ξ< 0 have a central emissionlike component. This phenomenon is produced by the isoradial velocity C_{d} curves close to C_{d} = 0, which are shorter than those approaching C_{d} = 1. The absorption in the central wavelengths is then less pronounced and the line profile acquires an emissionlike aspect. In case “c”, when the gradient ∂Ω /∂θ< 0 is large enough, the polar regions become dimpled. All calculations in Figs. 2 and 3 were performed for the equatorial centrifugal to gravitational force ratio η = 0.9, where η is (5)and R_{e} and R_{c} are the stellar equatorial and critical radii, respectively. Frequently, we use the parameter η instead of Ω_{p} or Ω_{e} to quantify the rotation rate, however, Eq. (3) can also be written as (6)where R_{e}(η) /R_{c} is calculated as explained in Zorec et al. (2011) and Zorec & Royer (2012). We note that R_{e} = R_{e}(η,ξ,k,γ,M,t), where M is the stellar mass and t its age. To calculate the geometrical shape of rotating stars and the line profiles in Figs. 2 and 3, respectively, we used M/M_{⊙} = 8.2 and t/t_{MS} = 0.72 (t_{MS} is the time spent by a star in the mainsequence evolutionary phase), which compared to a parent nonrotating object corresponds to T_{eff} = 22 000 K and log g = 4.0.
In what follows, we illustrate the effect of the differential rotation on the determination of the Vsini parameter and its incidence on the distribution of the true equatorial rotational velocity.
Fig. 4 Rotationally broadened line profiles of an intrinsic Gaussian absorption line, whose intensity and equivalent width depends on the effective temperature and gravity like the He i 4471 transition. Calculations include gravitational darkening in an object rotating at the equator with u = V/V_{c} = 0.9, whose parent nonrotating counterpart has M = 8.6 M_{⊙} and t/t_{MS} = 0.5. The angular velocity law is given by Eq. (7). The colors indicate the inclination angles, but in this figure only the i = 30° (blue) corresponding to different values of α are explicitly identified. 
2.3. Effect on the V sin i parameter
In Sect. 2.1 we pointed out that owing to our ignorance about the physical circumstances that may produce differential rotation in stars, no specific rotational law can be preferred to calculate the spectral line broadening. The first reasonable step we can undertake to uncover the trends produced by possible components of differential rotation is to adopt a Ω(θ) law with the fewest possible number of free parameters. We have therefore adopted the known Maunder relation, (7)which comes from the generic form given by Eq. (3) by putting k = 2 and γ = 1, and considering that in Eq. (4) Ω_{1} and ϵ are constants, so that Ω_{e} = Ω_{1}, α = −ϵ and (8)In Eq. (7), α is the free differential rotation parameter, whose possible values are limited in this work to the interval −1 < α < +∞. In principle, we could also impose values α < −1, but in our initial research on the effects produced on the broadening of lines, we avoid the case where the pole and equator rotate in opposite senses by imposing the condition α > −1. On the other hand, because the distance to the rotation axis diminishes as θ → 0, there is no limitation to the positive values of α. Obviously, rigid rotation is described by α = 0. From Eq. (7) it follows that ∂Ω /∂θ > 0 if α < 0, which produces δ < 0, and it is ∂Ω /∂θ < 0 when α > 0, which implies δ > 0. Since 0 ≤ v ≤ 1 and 0 ≤ u ≤ 1, from Eq. (10) it follows that − 1 ≤ δ ≤ 1.
Fig. 5 a) FWHM in km s^{1} observed at different angles i of the HeI 4471 line broadened by differential rotation characterized by different values of the parameter α in an object of “pnrc” parameters T_{eff} = 23 000 K and log g = 4.1 rotating at different rates η. In abscissas are the FWHM of the same line corresponding to rigid rotators having the same rates η. b) Idem for the MgII 4481 line. 
Almost the entire variety of line profiles calculated with Eq. (3) can reasonably be accounted for with the Maunder relation [Eq. (7)], as shown in Fig. 4. For these profiles we used an intrinsic Gaussian absorption line assimilated to the He i 4471 transition in an object with M = 8.6 M_{⊙} and t/t_{MS} = 0.5. The line broadening was calculated for the equatorial velocity ratio V/V_{c} = 0.9 and for several values of α and inclination angles.
Two methods are currently used to determine the Vsini parameter:

1.
FWHM against Vsini. This simple method is based on a polynomialrelation between the FWHM of spectral lines that are not stronglyaffected by the Stark broadening, with theVsini determined for wellstudied stars or with synthetic line profiles of rigidly rotating model atmospheres (Struve 1930; Slettebak & Howard 1955; Slettebak et al. 1975; Slettebak 1982). Taking the He i 4471 and Mg ii 4481 lines as witnesses of different sensitivities to physical conditions of line formation in stellar atmospheres (see also discussion in Sect. 3.3.2), we calculated their theoretical FWHM broadened using Ω(θ) given in Eq. (7) for several values of α including α = 0. The results are shown in Fig. 5, where the sensitivity of FWHM to the differential parameter α is apparent. This clearly shows that if differential rotation exists, depending on the sign of α, the Vsini parameters determined using the above FWHMbased method barely reflects the true equatorial rotation of stars. As reference to inquire as to what could be the expected effects in stars, it is worth noting that for the Sun it is α = −0.3.

2.
Fourier transform. The rather widespread method based on the Fourier transform (FT) to determine the Vsini assumes that an observed line profile F(λ) can be represented by the convolution of the flux line profile F(λ) of a nonrotating star with a “rotation broadening function” G_{R}(λ). This method requires that the specific intensity I(λ,μ) contributing to the observed line flux F(λ) has the same shape over the entire stellar disc, which is possible only if the temperature and gravity are uniform over the entire stellar disc. The analytical expression of G_{R}(λ), which is independent of the inclination angle, also requires that the star be spherical and behave as a rigid rotator (Gray 1975, 1992). The Vsini is finally obtained by comparing the zeroes of the FT of flux line profile F(λ) with the corresponding zeroes of the FT of the function G_{R}(λ).
However, the surface of a rotating star is neither spherical nor does it have uniform effective temperature and gravity. From Fig. 2 it is clearly seen that in a differential rotator the lines of constant radial velocity contributing to a given rotational Doppler displacement impose a rather complicated analytical representation, which is dependent on the a priori unknown inclination anglei (Huang 1961). We call G_{D}(λ) a hypothetical “rotation broadening function” that we can derive for a given surface rotation law. Even if we wanted to maintain a specific intensity I(λ,μ) behaving as in a nonrotating star, the zeroes of the FT of the function G_{D}(λ) would not correspond to those derived from the FT of G_{R}(λ) currently used in the literature. So, the zeroes of the FT of F(λ) are not compatible with those of FT[G(λ)]. As a consequence, the Vsini that is derived using the classical FT method cannot be considered a reliable parameter to represent the rotation of a given star.
We calculated He i 4471 and Mg ii 4481 lines broadened by differential rotators and calculated their FT. Depending on the spectral line and the value of α, the FT of these lines may have unusual shapes where the zeroes are difficult to interpret. In Fig. 6 are shown some FT of the He i 4471 and Mg ii 4481 lines broadened by differential rotators characterized by the pnrc parameters T_{eff} = 23 000 K, log g = 4.1, different values of α, where all have the same equatorial velocity rate η = 0.8. The resulting Vsini parameters obtained using the classical FT method are compared in Fig. 7 with the actual V_{eq}sini corresponding to several values of η, parameters α, and different inclination angles. Deviations to the y = x also exist for rigid rotators. These deviations are due to the gravitational darkening effect, which in the classical rigidrotation broadening function G_{R}(λ) is not taken into account. If differential rotation actually existed and if it were at least on the order of that observed in the Sun (α ≃ −0.3), the Vsini obtained with the FT method would largely deviate from those which are currently believed to characterize the stellar equatorial rotation.
Fig. 6 a)−c) Fourier transforms (FT) of the HeI 4471 line broadened by differential rotation characterized by α = −0.5, 0.0 and +0.5, respectively, in an object of “pnrc” parameters T_{eff} = 23 000 K and log g = 4.1 rotating at η = 0.8. d)−f) Idem for the MgII 4481 line. 
Fig. 7 a) Measured Vsini parameters in km s^{1} by means of the FT technique of the HeI 4471 line broadened by differential rotation characterized by different values of the parameter α in an object of “pnrc” parameters T_{eff} = 23 000 K and log g = 4.1 rotating at different rates η. In abscissas are the true equatorial Vsini for the same rates η. b) Idem for the MgII 4481 line. 
Figures 5 and 7 suggest that the Vsini obtained using either the FWHM or the FT method can be overestimated when α> 0, or underestimated for α< 0 as compared with its value for rigid rotators. Also, for a given  α  value the deviation is larger when α< 0. This means that if stars were actually differential rotators, the distributions derived thus far of ratios V/V_{c} would obviously not represent the distributions of actual equatorial rotations. In the following section we explore the effects induced by the differential rotation on the distribution of V/V_{c} ratios.
2.4. Effect on the distribution of ratios V/V_{c}
2.4.1. Some preliminary relations
We assume that a given studied star undergoes differential rotation and that V_{D}sini is its projected rotational parameter derived by one of the methods evoked above. This quantity can be written as (9)where V_{e} is the equatorial velocity of the star assumed to rotate as a rigid rotator, and ΔV(sini) is the contribution to the spectral line broadening by the hemisphereaveraged deviation from rigid rotation. In what follows, we deal with linear velocity components normalized to the equatorial critical velocity V_{c} as (10)where δ< 0 implies that the stellar surface rotation is accelerated from the pole toward the equator, δ> 0 indicates a rotation accelerated from the equator toward the pole, and δ = 0 is for rigid rotation.
To estimate the deviation δ = δ(i) we used intrinsic Gaussian absorption lines, as in Sect. 2.2, and calculated rotationally broadened flux line profiles with the GD relation put forward by Espinosa Lara & Rieutord (2011). This GD relation is not adapted for surfaces undergoing differential rotation. However even a less detailed treatment of these effects is far out of the scope of the present statistical correction of rotational velocity distributions. These effects will be presented in detail in a series of forthcoming papers.
Using the classical method of Vsini determination with the Fourier transform, we calculated V_{D}sini and V_{r}sini for a large space of parameters (u,α,sini).
The Fourier transformation of the gravitydarkened line profiles broadened by a surface differential rotation and rigid rotation was carried out with a rotationalbroadening function produced by a linear limbdarkening law in a gray atmosphere of a rigid rotator, as is still currently used in the literature, where the first FT root σ_{1} is related with the Vsini by the relation λ(Vsini/c)σ_{1} = 0.66. This was carried out to render similar in nature the model Vsini parameters to the stellar ones used in Paper I.
Because the function Φ(u) is meant to represent the distribution of truevelocity ratios u = V/V_{c} obtained from inclinationangle averaged apparent velocity ratios, in what follows we have to use inclinationangle averaged deviations ⟨ δ ⟩ that we define as follows: (11)where P(i) di = sini di. In Fig. 8a are shown the deviations δ(α,u,i) as a function of α for an object with mass M/M_{⊙} = 8.6, fractional mainsequence age t/t_{MS} = 0.5 (t_{MS} is the time spent by a star in the mainsequence evolutionary phase) rotating with u = V_{e}/V_{c} = 1.0, and for several inclination angles 5° < i < 90°. The iaveraged deviation ⟨ δ(α,u) ⟩ against α for masses M/M_{⊙} = 3.1, 8.6, 13.2, which are all at the fractional age t/t_{MS} = 0.5 and rotating at ratios u = V_{e}/V_{c} = 0.4, 0.7, 1.0, is shown in Fig. 8b. Figure 8c shows the iaveraged ⟨ δ(α,u) ⟩ deviations obtained for u = 1 in M/M_{⊙} = 8.6 at three fractional ages t/t_{MS} = 0.5, 0.5, and 0.9.
Owing to the rather significant dependence of ⟨ δ(α,u) ⟩ with the stellar mass, we ought to study Φ(u) distributions obtained for several representative intervals of stellar masses. Unfortunately, the total number of objects in the entire sample is too low to define statistically reliable subgroups. Having then a unique distribution Φ(u) for masses that are mostly in the interval 3.0 ≲ M/M_{⊙} ≲ 30.0, we obtained ⟨ δ(α,u,M) ⟩ _{t} curves, similar to those shown in Fig. 8b, for the 11 massage stellar groups given in Table 1. The adopted deviation ⟨ δ(α,u) ⟩, shown in Fig. 8d and used in the final calculations, is the weighted average by the fraction of stars in each massage group.
Fig. 8 a) Deviation parameter δ(α,u,i) as a function of α for u = V_{e}/V_{c} = 1.0, mass M/M_{⊙} = 8.6, and fractional age t/t_{MS} = 0.5, and for several inclination angles 5° < i < 90°. b) Inclinationaveraged deviation ⟨ α(α,u) ⟩ for velocity ratios u = V_{e}/V_{c} = 0.4, 0.7, 1.0 in objects with masses M/M_{⊙} = 3.1, 8.6, 13.2 at fractional age t/t_{MS} = 0.5. c) Inclinationaveraged deviation ⟨ α(α,t) ⟩ for M/M_{⊙} = 8.6 at fractional agest/t_{MS} = 0.2, 0.5, 0.9 rotating at u = 1.0. d) The adopted inclinationmassageaveraged deviation ⟨ δ ⟩ as a function of α for different values of the ratio u. 
To simplify notations, in what follows we continue to use the notation δ for the aspect angle averaged deviation ⟨ δ(α,u) ⟩ and write (12)Since by definition 0 ≤ v ≤ 1 and 0 ≤ u ≤ 1, from Eq. (12), the theoretical range of possible values of δ is restricted to − 1 ≤ δ ≤ 1. Nevertheless, from the calculated line profiles it follows that lim_{α → + ∞}δ(u,α) =δ_{L}(u) < 1 for all values of u.
2.4.2. Distributions of apparent and actual equatorial velocities ratios V/V_{c}
Following the principles for convolutions established in Appendix D of Paper I, the relation between the observed distribution Ψ(v) and the corrected distribution Φ(u), takes the form (13)The function Ψ(v) is the distribution of ratios v = V/V_{c}, where V does not represent a actual equatorial velocity, but a kind of linear velocity averaged over the apparent stellar hemisphere that accounts for the observed line broadening. Conversely, the variable u represents the ratio of actual true equatorial rotational velocities V_{eq}/V_{c}. In Eq. (14) the function φ(α) is the probability for the occurrence of the differential rotation parameter α. The lower limit of integration over u in Eq. (14) is determined by the condition u_{1}(v) = v−δ_{L}(u_{1}) as noted after the relation in Eq. (12), while u_{2} = 1 is the natural upper limit. Obviously, this convolution can also be understood as an integration over the interval of possible values of δ of the functional product Φ(v−δ)π(δ).
In Eq. (14) only Ψ(v) is available, while neither of the functions Φ(u) and φ(α) is known. Moreover, since δ(u,α) depends on u and α, the functions Φ(u) and φ(α)) are mutually dependent. The problem is then to apply an iterative blind deconvolution method to obtain a solution for Φ(u) (e.g., Ayers & Dainty 1988; Tsumuraya et al. 1994; Li et al. 2016). The application of such a method would require us to derive the probability function φ(α), which obeys the constraining condition given by the relation δ = δ(u,α), and to deal with a distribution function π(v−u) whose mode and dispersion change with v. Instead of discussing in detail the convergence of the method and the uniqueness of the possible inferred solutions for Φ(u) and φ(α), which could lead us beyond the scope of the present section, we preferred to estimate only the order of magnitude of the deviations that might be expected to be induced by the differential rotation on the distributions of rotational velocity ratios V/V_{c}, by using in Eq. (14) analytical expressions for Φ(u) and φ(α) characterized by free parameters.
In Fig. 7 we see that the deviations produced on the Vsini values can probably be unrealistic if we impose differential rotation parameters  α  > 1. Hence, we assume simple functional forms for φ(α) restricted to the interval − 1 ≤ α ≤ + 1 and characterized by a dispersionlike parameter a_{α} () and the modelike parameter α_{0} as follows: (14)where φ_{0} is the normalization constant that depends on a_{α} and α_{0}. We adopted a_{α} = 4.5.0 and 250.0 and several values for α_{0} from − 0.5 to + 0.5.
Massage groups formed by the sample of Be stars entering the “observed” distribution Φ(u).
Fig. 9 a) Assumed occurrence probability functions φ(α) of the differential rotation parameter α with a_{α} = 4.5; α_{0}< 0 (blue curves), α_{0} = 0 (black curve), α_{0}> 0 (red curves). b) Idem as in a), but for a_{α} = 250.0. a’) Occurrence probability functions π(δ) of the differential rotation parameter δ with a_{α} = 4.5 for different values of α_{0} and true velocity ratios u. b’) Idem as in a’), but for a_{α} = 250.0. a”) Distributions Ψ(v) of the equatorial velocity ratios v = V/V_{c}, where v is now affected by the differential rotation according to the distributions φ(α) of the parameter α given in a). b”) Idem as in a”), but for φ(α) shown in b). In a”) and b”) is also shown the assumed distribution Φ(u_{eq}) (magenta crosses). 
The distribution of equatorial velocity ratios Φ(u_{eq}) [u_{eq} = V_{eq}/V_{c}] is unknown by definition. The function Φ(u) obtained in Paper I for Be stars, after corrections for GD effect and microturbulent motions, must be considered to be, in principle, affected by the differential rotation of stars. It is then taken here as representing Ψ(v). To remain, however, not far from the physical circumstances that determine the function Φ(u), for the present exercise we assume Φ(u_{eq}) given by an analytical expression with roughly the same mode and FWHM as Φ(V/V_{c}) of Paper I mentioned above. So, we have (15)with Φ_{0} as the normalization constant, a = 0.649, b = 29.94 and u_{0} = 0.671.
Figures 9a, b show the assumed functions φ(α) used in this exercise according to the combination of parameters a_{α} and α_{0} given in the inserts. Figures. 9a, b’ reproduce some functions π(δ) = π(v−u) [see Eq. (14)] for three values of α_{0} in φ(α) and corresponding to three different values of v. Figures 9a”, b” show the respective distributions Ψ(v) calculated from Eq. (14). In these figures the function Φ(u_{eq}) (red crosses), given by Eq. (15), is also shown. To give a quantitative overview of the effects carried by the differential rotation on the distribution of ratios u = V/V_{c}, Table 2 reproduces the modes of the functions φ(α), the calculated distributions Ψ(v), and the FWHM of the resulting “apparent” functions Ψ(v).
2.5. Comments on the results obtained
From the results presented in Sect. 2.3 we conclude that if the stars had their atmosphere in differential rotation, the parameters Vsini obtained with the currently used measurement methods would depend on the degree of differential rotation and sensitively differ from the values expected if they were rigid rotators. The difference is due to, on one hand, the geometry of the constant radialvelocity curves contributing to the individual monochromatic Doppler displacements that contribute to the broadening of a spectral line, and on the other hand, to the inadequacy of the measurement methods, which are not adapted to the physical circumstances imposed by the differential rotation. In stars rotating with an angular velocity accelerated from the pole toward the equator (α< 0), the measured parameter Vsini implies that V<V_{eq}, where V_{eq} is the actual linear rotational velocity of the equator. Conversely, a differential rotation accelerated from the equator toward pole (α> 0) leads to a Vsini parameter where V>V_{eq}. It also follows from these results that for the rotation law given in Eq. (7) and for a given absolute value  α , the differences  V−V_{eq}  are larger for α< 0 than for α> 0. The precise value of the difference ΔV = V−V_{eq} depends on the specific differential rotation law, on the value of V_{eq}, and on the spectral line used to carry out the measurements due to the different sensitivities to the physical conditions in their formation regions (see Sect. 3.3.2, item 1)).
Modes of functions φ(α) and modes and FWHM of distributions Ψ(v).
Regarding the distribution of velocity ratios v = V/V_{c}, where V is not the stellar equatorial velocity but accounts for the spectral line broadening as discussed in Sect. 2.4.2, Table 2 reveals that if the surface rotation laws are sketched with the relation in Eq. (7), for α< 0 the modes of the distributions Ψ(v) are , and when α > 0. It also follows that for a given  α , the is larger when α> 0. On average, the FWHM of Ψ(v) is an increasing function of α. For a given value of α, the larger the FWHM of Ψ(v) is the smaller a_{α} is [larger “dispersion” of the probability function φ(α)]. When a_{α} ≳ 4.5 the FWHM is roughly uniform over a large interval of α_{0} values.
To avoid confusion due to the notation used, we note Ψ(v)_{obs}, the resulting distribution Φ(u) in Paper I, after corrections for GD and macroturbulent motions. The distributions Ψ(v) issued from Eq. (14) are thus meant to be different cases of Ψ(v)_{obs}. So, because of the displacement of the distribution of apparent velocity ratios caused by the differential rotation with respect to the distribution of the actual equatorial velocity ratios, we can expect that the number of Be stars with equatorial velocities attaining V_{eq} ≳ 0.9V_{c} can be larger than suggested by Ψ(v)_{obs} if on average the stars have rotation laws with α < 0. On the contrary, the number of rotators with V_{eq} ≳ 0.9V_{c} could be lower if they have rotation laws with α > 0. Finally, if our program Be stars rotate with α ≃ −0.3 as the Sun, or α ≃ −0.2 as predicted for stars with masses 2 ≲ M/M_{⊙} ≲ 4 by Espinosa Lara & Rieutord (2013), the Ψ(v)_{obs} will have its mode displaced toward lower values by Δv ≃ −0.07 ± 0.01 or Δv ≃ −0.11 ± 0.01, respectively, as compared with the mode of the distribution of their actual equatorial ratios u_{eq} = V_{eq}/V_{c}; the noted uncertainty accounts for the two values of a_{α} used here.
3. Some hidden effects on the V sin i parameter
A critical study of the distribution of the Vsini parameter and the conclusions that could be drawn from it cannot be completed if, nevertheless, we do not mention a number of conceptual and measurement uncertainties related to the determination of this quantity, other than differential rotation, which can also render the study of the true rotational properties of stars difficult.
From the existing compilations of Bestar Vsini parameters, we realize that the typical dispersions of independent estimations for the same objects carried out by different authors are typically 20 ≲ 1σ_{Vsini} ≲ 50 km s^{1} (Slettebak 1982; Chauville et al. 2001; Yudin 2001; Głȩbocki & Gnaciński 2005). On the other hand, more or less regular and homogeneous followup observations of Be stars reveal variations of Vsini on the order of 30 km s^{1} (Chauville et al. 2001; Rivinius et al. 2013a). Deviations can sometimes be as high as 100 to 150 km s^{1} as for HD 45910, HD 52721, and HD 135734 (Yudin 2001), or they can even attain 200 km s^{1} as in γ Cas (Slettebak 1982; Chauville et al. 2001; Harmanec 2002). Such discrepancies are possibly not only related to the method used to determine the Vsini, but probably also with physical phenomena that affect the spectral line formation.
Apart from macroturbulent motions evoked in Sect. 4 of Paper I, which pertain to low Vsini (Vsini ≲ 150 km s^{1}), the Vsini of Be stars and active stars, in general, are sensitive to spectral line variation and broadening owing to a sum of physical phenomena that can contribute individually with systematic and/or irregular components of a few tenths of km s^{1}. In some cases these stars can significantly blur or obliterate the information on the stellar rotation carried by the Vsini. Among such phenomena and/or uncertainties, the following deserve some attention:

i)
the angular momentum content in stars;

ii)
doublevalued Vsini;

iii)
GD effect on the radiation source function;

iv)
effects carried by intrinsically asymmetric rotational broadening functions;

v)
deviations produced by expanding layers;

vi)
deformations of line profiles due to tidal interactions in binary systems;

vii)
presence of circumstellar envelopes or discs (CD).
To this list we should add:

viii)
wavelengthdependent limbdarkening within a spectralline (Underhill 1968; Stoeckley& Mihalas 1973; Collins &Truax 1995; Levenhagen 2014);

ix)
uncertainties due to nonradial pulsations (Aerts et al. 2009, 2014).
These uncertainties were discussed by the cited authors and the specific publications provide details.
3.1. Angular momentum content in stars
The most frequent descriptions of rotating objects are based on stellar models that start evolving in the ZAMS as rigid rotators. Evolution then rapidly brings about the angular momentum redistribution so that stars become differential rotators (Maeder 2009). Resuming the angular momentum content in a star with the ratio τ_{Ω} = K/  W , where K is its total rotational kinetic energy and W represents its gravitational potential energy, these models are characterized by ratios (16)where τ_{Ω} ≈ 0.01 identifies the rigid critical rotators. At this limit, owing to the masscompensation effect, the bolometric luminosity produced in the stellar core can decrease by no more than some 5% with respect to its value in stars at rest having masses M/M_{⊙} ≃ 60, and to about 10% for M/M_{⊙} ≃ 1.5 (Bodenheimer 1971; Maeder & Meynet 2000; Ekström et al. 2008). However, secularly stable axisymmetric models can be constructed with energy ratios as large as (17)which necessarily correspond to differential rotators in depth and where the reduction of the core bolometric luminosity at τ_{Ω} ~ 0.10 attains from 40% in stars with M/M_{⊙} ≃ 60 to 80% in those with M/M_{⊙} ≃ 1.5 (Bodenheimer 1971; Clement 1979; Eriguchi & Mueller 1985; Deupree 2001; Jackson et al. 2005).
Except for the Sun and perhaps KIC 10526294 (Triana et al. 2015), where an upper limit could be attempted for τ_{Ω}, at the moment it is impossible to produce systematic estimates of this ratio in other stars. However, if actual Be stars and other massive rapid rotators had ratios τ_{Ω} as high as in Eq. (17), significant difficulties could appear to match the right distributions of surface effective temperatures and gravities to predict correctly profiles of spectral lines that enable us to obtain consistent values of Vsini. Needless to say, the diagnostic of stellar masses and ages would also be uncertain if models of rotating stars did not correspond to the right rotationenergy content in the analyzed stars (Zorec 1986; Collins & Smith 1985; Lovekin et al. 2006; Gillich et al. 2008; Castañeda & Deupree 2014; Castañeda et al. 2016).
Fig. 10 Angular velocity ratio Ω/Ω_{c} as a function of η, where the uncertainty bars indicate differences for masses and ages in the mainsequence phase. It is shown that Ω/Ω_{c} ≃ 0.95 (red line) identifies the force ratio η ≃ 0.6 at the equator, which implies a 40% undercritical rotation. The ratio V/V_{c} ≃ 0.95 (red line) corresponds to η ≃ 0.86. 
3.2. Doublevalued V sin i
The geometry of the stellar surface depends strongly on its surface rotation law. Objects with laws accelerated toward the pole may present polar dimples and, thus, the largest contribution to the Doppler broadening of spectral lines comes from a point of intermediate latitude situated on the stellar meridian facing the observer. For large inclination angles, the contribution to the line broadening produced by the closed curves corresponding to C> 1 (see Fig. 2c) can be partially hidden in the troughs, which in addition can be more or less dimmed by the limb darkening. The relation between the FWHM of spectral lines and the Vsini then becomes double valued, as shown in Zorec (1986, 1994), Zorec et al. (2004).
3.3. Comments on the gravity darkening
3.3.1. Uncertainties related to the gravity darkening exponent
The upper layers of the envelope of massive and intermediatemass stars are globally in hydrostatic and radiative equilibrium. Accordingly, a relation between the emerging bolometric radiation flux F_{bol} at a given colatitude θ and the corresponding effective gravity g_{eff}(θ) takes the form (Espinosa Lara & Rieutord 2011; Rieutord 2016; Zorec et al., in prep.). (18)where f(η,θ) is a function of the rotation law in the stellar atmosphere and of the equatorial rotational rate η defined in Eq. (5). However, the GD relation is still frequently written as (19)where the coefficient κ(η) and the exponent β_{1} are considered independent of the stellar colatitude θ. Actually, if radiative equilibrium is imposed in the atmosphere and conservative rotation laws, in particular the rigid rotation, it can be shown that β_{1} = 1 (von Zeipel 1924; Tassoul 1978). Otherwise, conserving the expression given in Eq. (19) and the above assumption of constancy for κ(η) and β_{1}, it follows that for nonconservative laws the GD exponent becomes β_{1} ≤ 1 (Kippenhahn 1977; Maeder 1999; Claret 2012). In general, the twodimensional nature of the radiative transfer in geometrically deformed stars (or threedimensional as for rapidly rotating components in binary systems) induces a horizontal diffusion of light (Osaki 1966). Accordingly, the GD effect decreases the higher the stellar oblateness, so that β_{1} ≤ 1 or even β_{1} → 0 (Pustyl’Nik 1970; Smith & Worley 1974; Hadrava 1992).
Fig. 11 a) Effective temperature ratio T_{eff}/T_{0} against the colatitude θ derived by Zorec et al. (in prep.) using Eq. (19) for η = 0.8 and several values of the differential rotation parameter α. The paramater T_{0} is the effective temperature of the stellar counterpart without rotation. b) Ratio against η and for the same values of α as in a). 
Fig. 12 Rotationally broadened Gaussian line (assimilated to He i 4471) by a model star in rigid rotation with η = 0.9 and pnrc parameters T_{eff} = 16 000 K and log g = 4.0. The dotted lines indicate no GD effect; the dashed lines indicate a GD effect using the von Zeipel formulation with β_{1} = 1.0; and the solid line indicates the GD effect using EspinosaLara & Rieutord’s (2011) formulation. The colors indicate the inclination angle i. 
The interferometric imaging of stellar surfaces with models of rapid rotators based on Eq. (19), where κ(η) and β_{1} are kept constant, produces β_{1} ≲ 1.0 with β_{1} decreasing as the flattening increases (Monnier et al. 2007; Zhao et al. 2009; Che et al. 2011; Monnier et al. 2012; Delaa et al. 2013; Domiciano de Souza et al. 2014). However, Espinosa Lara & Rieutord (2011) showed that even in rigid rotators β_{1} decreases as the η increases and it is also a function of the colatitude θ. The dependence with θ implies that the GD exponent determined empirically is necessarily inclination angle idependent (Rieutord 2016; Zorec et al., in prep.). The information that β_{1} might carry on the rotation of stars is then somewhat obliterated, which is similar to what happens with the Vsini parameter due to the unknown sini.
According to relation in Eq. (18) the variation of the effective temperature with θ is somewhat different as obtained from Eq. (19). This difference is even stronger if the stellar surfaces undergo differential rotation as shown in Zorec et al. (in prep.). Figure 11a shows the variation of the ratio of effective temperatures T_{eff}/T_{0} obtained in Zorec et al. (in prep.) (T_{0} is the effective temperature of the nonrotating stellar counterpart) against the colatitude θ in differentially rotating models with equatorial force ratio η = 0.8 and several values of the parameter α in Ω(θ,α) given by Eq. (7). Figure 11b shows the ratio against η and for the same values of α as in Fig. 11a. In these figures, the closest variations of the effective temperature with θ and η to the case described by Eq. (19) with β_{1} = 1.0 are those depicted by the green curves, which correspond to α = 0.
Differences in the predicted line profiles, in particular their equivalent widths, can appear according to the theory used to describe the GD effect. Figure 12 shows rotationally broadened flux profiles of a Gaussian line whose central intensity and equivalent width respond to the effective temperature and gravity as does the He i 4471 line. These profiles were calculated using in turn Eq. (19) with κ = const. and β_{1} = 1 of the classical formulation of the GD effect, and the expression in EspinosaLara & Rieutord (2011) given by Eq. (18) for α = 0. As noted above, when α = 0, the θdependent effective temperature does not differ greatly on which of these two GD theories is used. According to Fig. 11 much larger differences are expected when α ≠ 0.
3.3.2. Radiation source function in a gravitydarkened atmosphere
To present the effects we want to stress in this section clearly, two types of phenomena deserve a comment: 1) the sensitivity of line source functions to local physical conditions; and 2) the thermal structure on depth of stellar atmospheres as a function of the colatitude θ and its incidence on the line source function.
1) Sensitivity of line source functions.From the classical theory of stellar atmospheres, it is well known that due to nonLTE effects the source function of spectral lines have selective sensitivities to collisional and radiative processes, which dominate the population of atomic levels; these depend on the particular structure of ions (Thomas 1965, 1983; Mihalas 1978). Since a gravitydarkened atmosphere in a rapidly rotating star displays a wide range of electron temperatures and densities over the hemisphere projected toward the observer, spectral lines with different “reactivities” to local formation conditions do not reflect in the same way, or with the same strength, the physical properties in a given region of this hemisphere. To illustrate this effect, we calculated the T_{eff},log g_{eff}diagrams show in Fig. 13, where the color scaling corresponds to the relative equivalent width of the He i 4471 and Mg ii 4481 lines in a rigid critical rotator. We assumed 24 000 K and the color scale ranges from 0, for no contribution, to 1, or maximum contribution to the central flux. The polar regions are situated in the upper right corner, while the equatorial regions are in the lower left corner. From these diagrams we can conclude that a gravitydarkened atmosphere contributes to the rotational broadening of the He i 4471 line likely in the upper latitudes, while the rotational broadening of the Mg ii 4481 line is produced by the entire hemisphere. Owing to these differences, spectra of rotating stars interpreted with models for nonrotating atmospheres can produce not only different values of apparent (T_{eff},log g_{eff})parameters, but also they lead to different estimates of the Vsini, as shown in the next subsection. The contribution of local specific intensities to the observed flux in a spectral line is still modulated by the local continuum flux, which in a gravity darkened atmosphere can be lower in the equatorial region than in other latitudes.
2) The source function at τ_{λ} = 2/3.It was recognized very early that rotation and simultaneous hydrostatic and radiative equilibrium contradict each other, producing nonvanishing θdependent divergence of radiation flux (Osaki 1966). Moreover, because stars are probably baroclinic (nonconservative rotation laws) their isobars and isopycnic surfaces are neither parallel nor isothermal, which makes even more difficult to know the actual behavior of S_{λ} with θ at τ_{λ} = 2/3 , which is considered as the formation region of the radiation observed at λ. Unfortunately, up to now there are no detailed predictions made of atmospheric thermal structures as a function of θ in rapidly rotating stars. The current approximation used to model the spectra emitted by rotating objects is to assume that at each colatitude θ the stars have an internal atmospheric structure similar to that of a classical planparallel modelatmosphere characterized by the local parameters [T_{eff}(θ),log g_{eff}(θ)], where T_{eff}(θ) is inferred with some GD law among those cited in Sect. 3.3.1. Such a thermal structure has never been demonstrated to be consistent. So, at the moment, we can only simulate the effects produced by several extreme physical circumstances to see what the expected changes in the spectra are:

(a)
Uniform source function: the source functionS_{λ} in the radiation formation region is assumed to be the same for all θ and determined by the classic atmospheric temperature structure, whose effective temperature T_{eff} is the θaveraged effective temperature ⟨ T_{eff} ⟩, where each local T_{eff}(θ) is determined with Eq. (18). Each elementary atmosphere responds to the θdependent surface effective gravity g =g_{eff}(θ). However, the amount of the emerging λdependent specific intensity at each θ is assumed to be modulated according to the bolometric flux given in Eq. (18).

(b)
No gravitational darkening: following the predictions by Pustyl’Nik (1970) and Hadrava (1992) that in rapid rotators β_{1} → 0, we assume that the studied spectral characteristics are produced by a nonspherical rotating star with no GD, where the internal structure of each local elementary atmosphere conforms with a uniform T_{eff} =, but has a θdependent log g. The main difference with case (a) is that the local emerging specific intensity is not modulated by the bolometric flux given in Eq. (18) and that .

(c)
Classical gravitational darkening: we adopt the commonly used approach for the GD effect, i.e., Eq. (19) with κ = const. and β_{1} = 1, which according to the present line of thinking represents another extreme case. In this approach S_{λ} varies in depth at each θ as predicted by the local [T_{eff}(θ),log g_{eff}(θ)]dependent planparallel models of stellar atmospheres.
To account for the combined effects 1) and 2), we have calculated rotationally broadened He i 4471 and Mg ii 4481 line profiles with FASTROT (Frémat et al. 2005) adapted for the above listed atmospheric structures in model stars with pnrc fundamental parameters T_{eff} = 12 500 K and 22 000 K (pnrc = parentnonrotatingcounterpart), both having log g = 4.0, rotating at η = 0.99, and seen at i = 30°, 60°and 90°. In FASTROT the sensitivity to the local formation conditions of spectral lines evoked in 1) are taken into account automatically. Although the He i 4471 and Mg ii 4481 lines always appear blended in stellar spectra, here we purposefully study these lines in isolation. In this way we can view these lines as representative of deviations that other isolated spectral lines with differentiated sensitivities to local formation conditions may carry on the fundamental parameter determination. Figure 14 shows the line profiles of isolated He i 4471 and Mg ii 4481 lines according to the above assumptions a)−c).
Fig. 13 Upper panel: diagram showing the degree (in terms of relative equivalent width) of local contribution to the formation of the He i 4471 line in a rigid critical rotator. Lower panel: degree of local contribution to the formation of the Mg i 4481 line. The white stars identify the pole and the equator. 
Fig. 14 Rotationally broadened He i 4471 and Mg ii 44781 line profiles produced by geometrically deformed model stars characterized by and and rotating at η = 0.9. Red lines indicate hybrid GD with local emerging bolometric fluxes modulated according to Eq. (18), θdependent effective gravity and local source functions determined by a uniform electron temperature T_{e} = ⟨ T_{e} ⟩ (assumption “a”). Black lines indicate no GD and θdependent effective gravity (assumption “b”). Blue lines indicate classical GD with β = 1 and the local source functions are determined by T_{e} [T_{eff}(θ),log g(θ)] (assumption “c”). 
Apparent fundamental parameters deduced according to several assumptions about the structure of gravity darkened atmospheres.
Apart from the obvious differences that can be seen in Fig. 14 on the line profile predictions according to the type of GD adopted and the behavior of the line source function, for the intermediate effective temperatures there are significant differences in the predictions based on assumptions a) and c) when inclination angles are i ≳ 60°, under which most stars are seen. Accordingly, conflicting determinations of apparent fundamental parameters can result.
To sketch these conflicting determinations, we used the equivalent widths of the He i 4471 and Mg ii 4481 calculated in Fig. 14 and interpolated the apparent and in the model relations between these lines and the effective temperature and gravity in nonrotating stars shown in Frémat et al. (2005, see their Fig. 4). The obtained results are shown in Table 3, which includes the corresponding Vsini determined using the classical Fourier transform method.
Several conclusions can be drawn from these results:

1.
The apparent and rather strongly depend on the assumption made on the behavior of the radiation source function in GD atmospheres. As expected, the line profiles issued from model c) lead to the lower estimates of .

2.
The Vsini derived with the Mg ii 4481 line are systematically larger than those obtained with the He i 4471 line, which is a consequence of the differentiated sensitivity of these lines to the local formation conditions.

3.
As expected, and already well known, the gravitydarkened models produce lower values of Vsini than models without GD.

4.
Models with GD where the line source functions vary with θ produce the lowest Vsini parameters.

5.
The differences ΔVsini = Vsini(Mg)−Vsini(He)  are larger in those gravitydarkened models where the line source functions vary with θ.
Some of the assumptions made above on the physical structure of stellar atmosphere in rapidly rotating stars can unleash horizontal diffusion of light that has to be treated properly to determine the stable atmospheric thermal structure. Detailed twodimensional radiation transfer calculations in rotationally deformed stars, similar to those started by Pustyl’Nik (1970) and Hadrava (1992), might be able to tackle these questions and predict reliable dependencies of S_{λ} with θ as well as GD relations consistent with the right thermal structures of atmospheres. From these we may then expect to have reliable stellar fundamental parameters and line profiles required to make some progress in determining the properties of the stellar surface rotation.
3.4. Asymmetric rotational broadening function
The lines due to a variety of elements and their different ionization states in a given analyzed wavelength range are formed in atmospheric layers characterized by different physical conditions and, in particular by velocity fields entertained by nonradial pulsations, convection movements, more or less expanding layers of nascent/starting winds, etc. When the rotational broadening function (RBF) is determined by deconvolution of all spectral lines present in that spectral range, there is a small chance that its shape is symmetric (Reiners & Schmitt 2003b,a). It can be shown that, independent of the sign of the skewness of RBF, we always obtain underestimated Vsini parameters. To quantify this effect, we distorted the classic symmetric RBF of rigid rotators G_{S}(Δλ) using a linear function of the normalized wavelength displacement Δλ/λ_{M} as follows: (20)where − 1 ≤ a ≤ 1 is the free distortionparameter and Δλ_{M} is the wavelength displacement induced by the equatorial rotation velocity. Figure 15a shows a set of such asymmetric RBF. Figure 15b shows the He i 4471 line (T_{eff} = 20 000 K, log g = 4.0) rotationally broadened with Eq. (20) with a = −1 for several apparent equatorial rotational velocities Vsini. The lines are not only broadened, but they become increasingly shifted the higher the rotational velocity. Figure 15c shows the FT of the asymmetric RBF G_{A}(Δλ), where the maxima of lobes are the same for all Vsini, but the zeros of the symmetric RBF G_{S}(Δλ) are transformed into functional minima that are slightly displaced as soon as a ≠ 0. Figure 15d shows the relation between the (Vsini)_{A} implied by the distorted spectral lines and the true (Vsini)_{0} value, where the underestimation of the projected rotational velocity appears as an increasing function of the distortion parameter a and (Vsini)_{0}.
Fig. 15 a) Tested asymmetric rotation broadening functions (ARBF). b) Broadened He i 4471 line profiles with ARBF having a = −1 and for several input Vsini parameters. c) FT of the ARBF. d) Resulting (Vsini)_{A} parameters produced by ARBF against the model input (Vsini)_{0} (the colors correspond to the distortion parameters a in a)). 
Fig. 16 a) Tested truncated rotationexpansion broadening functions (REBF) for several velocity ratios . b) Loci of points with the same projected velocity in the interval of Doppler displacements that contribute to the broadening of line profiles when v_{exp} = 0.24. c) He i 4471 line profiles broadened for several input rotation velocities V_{rot}sini and velocity ratio V_{exp}/V_{rot} = 1.0. d) FT of the REBF. e) Resulting Vsini_{exp + rot} parameters produced in the frame of an expanding and rotating atmosphere against the model input V_{rot}sini for several ratios v_{exp}. 
We can then conclude that when a RBF is obtained by deconvolution of a rather wide spectral range, the supposedly more precise RBF introduces deviations that lead to underestimated parameters Vsini. Generally, the common practice is to symmetrize the RBF to derive Vsini from its TF. However, the symmetrization of a RBF changes its nature; this does not necessarily lead to the expected right Vsini value, as follows from the comparisons of Vsini values made in Collins & Truax (1995), where the authors adopt as examples the line broadening produced by elliptical and parabolic RBF functions.
Fig. 17 a) Model line profiles of the primary star broadened by rotation and perturbed by tidal interactions in a binary system with subsynchronous rotation (Ω_{∗} stellar angular velocity; Ω_{b} circular orbital angular velocity). b) FT of spectral lines shown in a). c) Model line profiles of the primary star broadened by rotation and perturbed by tidal interactions in a binary system with supersynchronous rotation. d) FT of spectral lines shown in c). The red lines are either for the averaged line profiles, which represent the tidally unperturbed rotationally broadened line profile, or for the FT of the averaged line profile. 
3.5. Effects carried by expanding layers
During the seventies and eighties a large amount of literature dealt with the spectral line formation problem in moving stellar atmospheres (see references in Sobolev 1960; Mihalas 1978; Underhill & Doazan 1982; Thomas 1983; Kalkofen 1984; Mihalas & Weibel Mihalas 1984; Kunasz 1985; Kalkofen 1987; Conti et al. 1988; Sen & Wilson 1998; Stee 2002; Hubeny et al. 2003; Cannon 2012; Hubeny & Mihalas 2014). These velocity fields were put into evidence mainly through spectroscopic observations carried out with satellites in the farUV spectral range. The farUV line asymmetries and the multicomponents of certain transitions were associated with mass ejections through winds, whose velocities are accelerated by the radiation pressure only once some layers start expanding. The presence of stellar winds thus implies the existence of expanding layers somewhere deeper in the atmosphere. In the literature, there are also claims of evidence of velocity fields detected in the visual spectral range of single Btype stars (Furenlid & Young 1980). The typical outbursts and fadings in Be stars, which are thought to be the consequence of huge discrete ejections of mass (Cook et al. 1995; Hubert & Floquet 1998; Keller et al. 2002; Mennickent et al. 2002), also imply that there are periods with outward accelerated layers in the stellar atmospheres (Hubert et al. 2000; de Wit et al. 2006).
The combination of expanding velocities with rotation produces broadened and blueshifted spectral lines that can be schematically described with rotationexpansion broadening functions (REBF). These functions are asymmetric, and depending on the ratio of characteristic velocities V_{exp}/V_{rot}, they can become truncated. This phenomenon was studied analytically in the frame of uniform expansions by Duval & Karp (1978). Detailed models of stellar atmospheres in rotation and expansion fields with velocity gradients can be found in Mihalas (1979).
We used the relationships obtained by Duval & Karp (1978) to infer the order of magnitude of effects produced on the Vsini determination by expanding velocity fields. Figure 16a shows a series of REBF obtained for different velocity ratios (V_{rot} = Vsini where i = π/ 2). Thus, v_{exp} = 1 corresponds to pure expansion, v_{exp} = 0 is for pure rotation, and identifies V_{exp}/V_{rot} = 1. Figure 16b depicts the loci of points of isoradial velocity contours for v_{exp} = 1/4, where the associated Doppler displacements belong to the interval . Figure 16c shows He i 4471 model line profiles (T_{eff} = 20 000 K, log g = 4.0) broadened with several values of V_{rot} and V_{exp}/V_{rot} = 1, where the blueshift produced by the expansion velocity component is clearly apparent. In Fig. 16d are shown the particular shapes of the FTs of He i 4471 lines broadened and deformed according to the indicated velocity ratios. The relations between the measured Vsini parameters, using the blueshifted and rationally broadened He i 4471 lines against the true V_{rot}sini for several ratios v_{exp}, are shown in Fig. 16e.
The results shown in Fig. 16 suggest at least two conclusions: 1) the underestimation of Vsini can be quite significant if the studied Be stars have been observed at phases when their atmospheres were driven by expanding velocity fields; 2) the blueshifts of lines can be more or less chaotic; they could partially explain irregular radial velocity drifts, which sometimes are suspected due to undisclosed binaries. Unfortunately, there is not enough systematic information to establish a distribution function of expanding velocities in atmospheres of Be stars to estimate their effect on the distribution Φ(u) and put to test the function φ(α) [or φ(δ)] studied in Sect. 2.4.2.
In the present discussion we have assumed that the entire atmosphere of stars undergoes uniform expansion. Nonetheless, this movement may possibly concern only an equatorial strip, as suggested by the correlation against the Vsini of light outbursts and fadings in Be stars obtained by Hubert & Floquet (1998). In such a case the underestimation of the Vsini parameter can be smaller than predicted here.
3.6. Tidal interactions in binary systems
A large percentage of massive and intermediate mass stars are in binary systems. Two interactions in binary systems can then perturb the determination of the Vsini of individual components: mutual irradiation and tidal interaction.
Owing to the mutual irradiation and radiation pressure that become significant in close binary systems (Cranmer 1993; Drechsel et al. 1995; Phillips & Podsiadlowski 2002; Palate et al. 2013), line spectra can change somewhat as a function of the orbital phase and inclination angle of the system. In this discussion we omit the mutual irradiation and radiation pressure.
Moreno & Koenigsberger (1999), Moreno et al. (2005) have studied the variation of photospheric spectral line profiles induced by the tidal deformations in binary systems, which can also induce differential rotation. Their results show that there are traveling bumps produced mainly by the azimuthal components of the velocity perturbations, which are superimposed on the rotationally broadened profiles. The number and strength of these bumps depend on the orbital phase and they are particularly strong in nonsynchronous systems. Dr. G. Koenisberger has graciously given us the entire series of calculations of line profiles published in Moreno et al. (2005), which we have shifted to a common central rest wavelength and applied the FT method to determine the apparent rotational velocity Vsini for each of them. Figure 17a shows the centered rotationally broadened and tidally perturbed line profiles in the case of a subsynchronous rotation with Ω_{∗}/ Ω_{b} = 0.6 (Ω_{∗} is the stellar angular velocity; Ω_{b} is the circular orbital angular velocity) and where the primary star has Vsini = 130 km s^{1}. The respective FTs of these lines are shown in Fig. 17b. Figures 17c, d are equivalent to Figs. 17a, b, but depict the case of a supersynchronous rotation with Ω_{∗}/ Ω_{b} = 1.5. In these figures, the red lines represent the averaged line profiles, which very closely correspond to the unperturbed rotationally broadened line, and their respective FT. In Table 4 are reproduced the Vsini parameters obtained from each individual line profile, which enables us to appreciate the apparent phasedependent change of Vsini.
Vsini parameters in km s^{1} determined using line profiles perturbed by tidal interactions in binary systems.
3.7. Presence of circumstellar envelopes or discs
One of the most outstanding hallmarks of Be stars is the presence of circumstellar envelopes or discs (CD), where their characteristic emissions and shell absorptions are raised. A large amount of literature describes the formation of these spectral features (cf. Underhill & Doazan 1982; Rivinius et al. 2013b). The Vsini determination can then be affected by variable emissions and/or absorptions superimposed on the lines used to determine this parameter. These effects on the He i 4471 line were discussed in Ballereau et al. (1995), where the detected emission can carry overestimations on the order of 30 ± 20 km s^{1} in Be stars with effective temperatures 15 000 ≲ T_{eff} ≲ 28 000 K. The Mg ii 4481 line, which is also of frequent use to measure the stellar rotation, can be affected either by circumstellar emission or shell absorption. The latter deepens the line absorption leading to underestimated values of Vsini. In the visible spectral region of Be stars, no line can be considered entirely free from circumstellar effects, but depending on the subspectral type some lines, such as Si i and Si ii lines, are apparently less damaged. However, when some lines seem to be good photospheric witnesses because of the GD effect their formation can be favored either on polar or equatorial regions as shown in Sect. 3.3.2. Only a careful modeling of such transitions can then lead to more reliable estimates of the Vsini parameters in Be stars.
Be stars are the epitome of stellar rapid rotators. However, their spectra can be marred by numerous disturbances produced by circumstellar envelopes or discs, which make their interpretation difficult and uncertain. These interpretations concern all insights we can draw from spectra on the surface and/or internal distribution of the angular velocity. As much it may concern the study of the angular momentum distribution in rapidly rotating stars, we should not neglect Bn stars because they are very rapid rotators and do not display line emission components and/or shell absorptions in their spectra. Moreover, there are statistical insights suggesting that these stars may represent a precursor stage to the Be phase of late Btype stars (Zorec 2000; Zorec et al. 2005). Finally, most Bn stars are of late B spectral type, which are more numerous and closer to us than earlytype Be stars.
4. Comments and conclusions
Many effects that obliterate the actual information carried by the Vsini and/or V parameters on the rotation cannot be studied individually for each Be star. We have then considered them statistically by studying the distribution of ratios Vsini/V_{c} and V/V_{c}.
In Paper I (Zorec et al. 2016) we studied a sample of 233 Galactic classical Be stars and obtained the distribution of their V/V_{c} by considering that the inclination angles i are distributed at random. We redetermined the fundamental parameters of the studied objects to correct the Vsini from its underestimation due to the GD effect. We also considered the overestimation of Vsini due to macroturbulent velocities.
In the first part of the present contribution (Sect. 2), we assumed that the atmospheres of Be stars undergo differential rotation. Sections 2.1 and 2.2 describe some properties of lines broadened by surface rotations laws characterized by angular velocities accelerated either toward the pole or the equator. We have shown that the main characteristics of the line profiles broadened by these laws can be reasonably accounted for with the Maunder relation, which depends on a single free parameter α.
In Sect. 2.3 we studied the effect produced by the differential rotation on the value of the Vsini parameter and concluded that rotation laws in which the angular velocity is accelerated from the pole toward the equator (α< 0) tend to produce Vsini<V_{eq}sini, where V_{eq} is the linear rotational velocity of the equator, while laws with accelerated from the equator toward the pole (α> 0) lead to Vsini>V_{eq}sini. Moreover, we have shown that methods like the Fourier transform produce unreliable values of Vsini because the zeroes of a putative rotationally broadening function cannot be associated with those inherent to a broadening function derived for rigid rotators. We noted that for stars with differential rotation the line broadening function cannot be defined, since neither the rotation law nor the inclination of the star are known.
To complete the discussion on the effects carried by the differential rotation, in Sect. 2.4 we estimated the effect it can produce on the distribution of ratios V/V_{c} of true rotational velocities. Owing to the displacement produced by the differential rotation of the distributions of apparent velocity ratios with respect to the distribution of the actual equatorial velocity ratios, we concluded that if a dominant number of Be stars rotate with laws characterized by α< 0, the number of rotators with V_{eq} ≳ 0.9V_{c} will be larger than expected from the observed distribution Φ(V/V_{c}). However, they could be lower if the Be stars had on average α> 0.
In the second part of the present contribution (Sect. 3) we examined a number of conceptual and measurement uncertainties other than differential rotation that affect the determination of the Vsini parameter, which are not currently evoked in the studies of stellar rotation. These include

i)
The angular momentum content, which can introduce increasedgeometrical deformations and subsequent GD contrasts.

ii)
Bivalued relations between the line broadening and Vsini, which can be due to surface rotation laws accelerated from the equator toward the pole.

iii)
The GD effect, which cannot be reliably characterized with a constant GD exponent as it is the current practice, and we question whether the radiation emitted by a rotationally deformed star can be reliably accounted for with classical models of stellar atmospheres responding to the local effective temperatures and gravities.

iv)
Asymmetric rotational broadening functions, which are determined by deconvolution of wide spectral ranges;

Effects caused on spectral lines by expansion velocities in stellar atmospheres, which can lead to underestimated Vsini.

vi)
Deformations of line profiles by tidal interactions in binary systems, which introduce changes of Vsini correlated with the orbital phase.

vii)
Overestimations or underestimations of the Vsini due to the presence of circumstellar envelopes or discs.
We note that from the numerical values of errors affecting the estimates of the Vsini of Be stars due to phenomena v) and vii), we could expect they can partially compensate each other as they act in opposite senses.
Although there is not enough data to estimate the effects enumerated in Sect. 3 to put to test the probability distribution φ(α), progress regarding the Vsini determination can be made as follows: a) by studying highly resolved spectra over large wavelength intervals with spectral lines that are good tracers of formation conditions (temperatures, densities, velocity fields); b) by obtaining models of the observed spectral region with twodimensional radiation transfer calculations in rotationally deformed atmospheres, including differential rotation, which solely can predict consistent GD laws; and c) by observing Bn stars to avoid spectral perturbations due to circumstellar envelopes and possibly huge expansion velocity fields in the atmospheres.
Thanks to the superresolution capabilities of modern spectrointerferometry, combined with high resolution spectroscopy, it could be possible to obtain direct information on the inclination angle i and the differential parameter α (Domiciano de Souza et al. 2004; Vinicius et al. 2006; Zorec et al. 2011; Delaa et al. 2013; Domiciano de Souza et al. 2014).
The knowledge of the degree of surface differential rotation and taking the effects mentioned in the list into account above will help, in particular, to decide whether or not the Be phenomenon mostly rely on the critical rotation.
Since differential rotation can be present in rapidly rotating objects other than Be stars, increased attention should be put to the observational aspects of this phenomenon. In fact, every observational indication on this phenomenon can be critical to deepen our understanding of the mechanisms involved in the angular momentum redistribution in stars and of the concomitant mixing phenomena of chemical elements, which help to predict the frequency of several stellar populations currently used to scrutinize the evolution of galaxies and/or the massive star formation rates (Maeder & Meynet 1995, 2004).
The correction for a presumable effect associated with the orbital motion of Be stars in binary systems presented in Paper I must not be taken into account (see Appendix A).
Acknowledgments
We thank Dr. Gloria Koenigsberger for having provided us with line profiles perturbed by tidal interaction in binary systems. We are thankful to the referee for his(her) useful and constructive criticisms and suggestions that significantly helped to correct and put in better focus some subjects discussed in this work.
References
 Aerts, C., Puls, J., Godart, M., & Dupret, M.A. 2009, A&A, 508, 409 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Aerts, C., SimónDíaz, S., Groot, P. J., & Degroote, P. 2014, A&A, 569, A118 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ayers, G. F., & Dainty, J. C. 1988, Opt. Lett., 13, 547 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Ballereau, D., Chauville, J., & Zorec, J. 1995, A&AS, 111, 423 [NASA ADS] [Google Scholar]
 Bodenheimer, P. 1971, ApJ, 167, 153 [NASA ADS] [CrossRef] [Google Scholar]
 Cannon, C. J. 2012, The Transfer of Spectral Line Radiation (Cambridge University Press) [Google Scholar]
 Castañeda, D., & Deupree, R. G. 2014, ApJ, 794, 13 [NASA ADS] [CrossRef] [Google Scholar]
 Castañeda, D., Deupree, R. G., & Aufdenberg, J. 2016, in Bright Emissaries: Be Stars as Messengers of StarDisk Physics, eds. T. A. A. Sigut, & C. E. Jones, ASP Conf. Ser., 506, 25 [Google Scholar]
 Chauville, J., Zorec, J., Ballereau, D., et al. 2001, A&A, 378, 861 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Che, X., Monnier, J. D., Zhao, M., et al. 2011, ApJ, 732, 68 [NASA ADS] [CrossRef] [Google Scholar]
 Claret, A. 2012, A&A, 538, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Clement, M. J. 1979, ApJ, 230, 230 [NASA ADS] [CrossRef] [Google Scholar]
 Collins, II, G. W., & Smith, R. C. 1985, MNRAS, 213, 519 [NASA ADS] [CrossRef] [Google Scholar]
 Collins, II, G. W., & Truax, R. J. 1995, ApJ, 439, 860 [NASA ADS] [CrossRef] [Google Scholar]
 Conti, P. S., Underhill, A. B., Jordan, S., & Thomas, R. 1988, NASA Spec. Publ., 497 [Google Scholar]
 Cook, K. H., Alcock, C., Allsman, H. A., et al. 1995, in Astrophysical Applications of Stellar Pulsation, eds. R. S. Stobie, & P. A. Whitelock, IAU Colloq., 155, ASP Conf. Ser., 83, 221 [Google Scholar]
 Cranmer, S. R. 1993, MNRAS, 263, 989 [NASA ADS] [Google Scholar]
 Cranmer, S. R., & Collins, II, G. W. 1993, ApJ, 412, 720 [NASA ADS] [CrossRef] [Google Scholar]
 Delaa, O., Zorec, J., Domiciano de Souza, A., et al. 2013, A&A, 555, A100 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Deupree, R. G. 1998, ApJ, 499, 340 [NASA ADS] [CrossRef] [Google Scholar]
 Deupree, R. G. 2001, ApJ, 552, 268 [NASA ADS] [CrossRef] [Google Scholar]
 de Wit, W. J., Lamers, H. J. G. L. M., Marquette, J. B., & Beaulieu, J. P. 2006, A&A, 456, 1027 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Domiciano de Souza, A., Zorec, J., Jankov, S., et al. 2004, A&A, 418, 781 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Domiciano de Souza, A., Kervella, P., Moser Faes, D., et al. 2014, A&A, 569, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Drechsel, H., Haas, S., Lorenz, R., & Gayler, S. 1995, A&A, 294, 723 [NASA ADS] [Google Scholar]
 Duval, P., & Karp, A. H. 1978, ApJ, 222, 220 [NASA ADS] [CrossRef] [Google Scholar]
 Ekström, S., Meynet, G., Maeder, A., & Barblan, F. 2008, A&A, 478, 467 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Eriguchi, Y., & Mueller, E. 1985, A&A, 146, 260 [NASA ADS] [Google Scholar]
 Espinosa Lara, F., & Rieutord, M. 2007, A&A, 470, 1013 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Espinosa Lara, F., & Rieutord, M. 2011, A&A, 533, A43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Espinosa Lara, F., & Rieutord, M. 2013, A&A, 552, A35 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Frémat, Y., Zorec, J., Hubert, A.M., & Floquet, M. 2005, A&A, 440, 305 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Furenlid, I., & Young, A. 1980, ApJ, 240, L59 [NASA ADS] [CrossRef] [Google Scholar]
 Gillich, A., Deupree, R. G., Lovekin, C. C., Short, C. I., & Toqué, N. 2008, ApJ, 683, 441 [NASA ADS] [CrossRef] [Google Scholar]
 Głȩbocki, R., & Gnaciński, P. 2005, in 13th Cambridge Workshop on Cool Stars, Stellar Systems and the Sun, eds. F. Favata, G. A. J. Hussain, & B. Battrick, ESA Spec. Publ., 560, 571 [Google Scholar]
 Gray, D. F. 1975, ApJ, 202, 148 [NASA ADS] [CrossRef] [Google Scholar]
 Gray, D. F. 1992, The observation and analysis of stellar photospheres, Camb. Astrophys. Ser., Vol. 20 [Google Scholar]
 Hadrava, P. 1992, A&A, 256, 519 [NASA ADS] [Google Scholar]
 Harmanec, P. 2002, in Exotic Stars as Challenges to Evolution, eds. C. A. Tout, & W. van Hamme, ASP Conf. Ser., 279, 221 [Google Scholar]
 Huang, S.S. 1961, ApJ, 133, 130 [NASA ADS] [CrossRef] [Google Scholar]
 Hubeny, I., & Mihalas, D. 2014, Theory of Stellar Atmospheres (Princeton University Press) [Google Scholar]
 Hubeny, I., Mihalas, D., & Werner, K. 2003, Stellar Atmosphere Modeling, ASP Conf. Ser., 288 [Google Scholar]
 Hubert, A. M., & Floquet, M. 1998, A&A, 335, 565 [NASA ADS] [Google Scholar]
 Hubert, A. M., Floquet, M., & Zorec, J. 2000, in The Be Phenomenon in EarlyType Stars, eds. M. A. Smith, H. F. Henrichs, & J. Fabregat, IAU Colloq., 175, ASP Conf. Ser., 214, 348 [Google Scholar]
 Jackson, S., MacGregor, K. B., & Skumanich, A. 2005, ApJS, 156, 245 [NASA ADS] [CrossRef] [Google Scholar]
 Kalkofen, W. 1984, Methods in radiative transfer (Cambridge University Press) [Google Scholar]
 Kalkofen, W. 1987, Numerical radiative transfer (Cambridge University Press) [Google Scholar]
 Keller, S. C., Bessell, M. S., Cook, K. H., Geha, M., & Syphers, D. 2002, AJ, 124, 2039 [NASA ADS] [CrossRef] [Google Scholar]
 Kippenhahn, R. 1977, A&A, 58, 267 [NASA ADS] [Google Scholar]
 Kunasz, P. B. 1985, in The theory of line transfer in expanding atmospheres, eds. J. E. Beckman, & L. Crivellari, NATO ASI, Ser. C, 152, 319 [Google Scholar]
 Levenhagen, R. S. 2014, ApJ, 797, 29 [NASA ADS] [CrossRef] [Google Scholar]
 Li, X., Ling, S., Strohmer, T., & Wei, K. 2016, ArXiv eprints [arXiv:1606.04933] [Google Scholar]
 Lovekin, C. C., Deupree, R. G., & Short, C. I. 2006, ApJ, 643, 460 [NASA ADS] [CrossRef] [Google Scholar]
 Maeder, A. 1999, A&A, 347, 185 [NASA ADS] [Google Scholar]
 Maeder, A. 2009, Physics, Formation and Evolution of Rotating Stars (SpringerVerlag) [Google Scholar]
 Maeder, A., & Meynet, G. 1995, The Messenger, 80, 19 [Google Scholar]
 Maeder, A., & Meynet, G. 2000, ARA&A, 38, 143 [Google Scholar]
 Maeder, A., & Meynet, G. 2004, in Stellar Rotation, eds. A. Maeder, & P. Eenens, IAU Symp., 215, 500 [Google Scholar]
 Maeder, A., Georgy, C., & Meynet, G. 2008, A&A, 479, L37 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mennickent, R. E., Pietrzyński, G., Gieren, W., & Szewczyk, O. 2002, A&A, 393, 887 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mihalas, D. 1978, Stellar atmospheres, 2 edn. (W. H. Freeman and Co.) [Google Scholar]
 Mihalas, D. 1979, MNRAS, 189, 671 [NASA ADS] [CrossRef] [Google Scholar]
 Mihalas, D., & Weibel Mihalas, B. 1984, Foundations of radiation hydrodynamics (Oxford University Press) [Google Scholar]
 Monnier, J. D., Zhao, M., Pedretti, E., et al. 2007, Science, 317, 342 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Monnier, J. D., Che, X., Zhao, M., et al. 2012, ApJ, 761, L3 [NASA ADS] [CrossRef] [Google Scholar]
 Moreno, E., & Koenigsberger, G. 1999, Rev. Mexicana Astron. Astrofis., 35, 157 [Google Scholar]
 Moreno, E., Koenigsberger, G., & Toledano, O. 2005, A&A, 437, 641 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Osaki, Y. 1966, PASJ, 18, 7 [NASA ADS] [Google Scholar]
 Palate, M., Rauw, G., Koenigsberger, G., & Moreno, E. 2013, A&A, 552, A39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Phillips, S. N., & Podsiadlowski, P. 2002, MNRAS, 337, 431 [NASA ADS] [CrossRef] [Google Scholar]
 Pustyl’Nik, I. 1970, Izv. Akad. Nauk Ehstonskoj SSR, Fiz., Mat., Tartu Astron. Obs. Teated, No. 29, 19, 428 [Google Scholar]
 Reiners, A., & Schmitt, J. H. M. M. 2003a, A&A, 412, 813 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Reiners, A., & Schmitt, J. H. M. M. 2003b, A&A, 398, 647 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Rieutord, M. 2016, in Lecture Notes in Physics, eds. J.P. Rozelot, & C. Neiner (Berlin Springer Verlag), 914, 101 [Google Scholar]
 Rivinius, T., Baade, D., Townsend, R. H. D., Carciofi, A. C., & Štefl, S. 2013a, A&A, 559, L4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Rivinius, T., Carciofi, A. C., & Martayan, C. 2013b, A&ARv, 21, 69 [Google Scholar]
 Schou, J., Antia, H. M., Basu, S., et al. 1998, ApJ, 505, 390 [Google Scholar]
 Sen, K. K., & Wilson, S. J. 1998, Radiative transfer in moving media: basic mathematical methods for radiative transfer in spherically symmetrical moving media (SpringerVerlag) [Google Scholar]
 Slettebak, A. 1982, ApJS, 50, 55 [NASA ADS] [CrossRef] [Google Scholar]
 Slettebak, A., Collins, II, G. W., Parkinson, T. D., Boyce, P. B., & White, N. M. 1975, ApJS, 29, 137 [NASA ADS] [CrossRef] [Google Scholar]
 Slettebak, A., & Howard, R. F. 1955, ApJ, 121, 102 [NASA ADS] [CrossRef] [Google Scholar]
 Smith, R. C., & Worley, R. 1974, MNRAS, 167, 199 [NASA ADS] [CrossRef] [Google Scholar]
 Sobolev, V. V. 1960, Moving envelopes of stars (Harvard University Press) [Google Scholar]
 Stee, P. 2002, Radiative transfer and hydrodynamics in astrophysics, EAS Publ. Ser., 5 [Google Scholar]
 Stoeckley, T. R. 1968, MNRAS, 140, 121 [NASA ADS] [Google Scholar]
 Stoeckley, T. R., & Buscombe, W. 1987, MNRAS, 227, 801 [NASA ADS] [Google Scholar]
 Stoeckley, T. R., & Mihalas, D. 1973, Limb Darkening and Rotation Broadening of Neutral Helium and Ionized Magnesium Line Profiles in EarlyType Stars, NCARTN/STR84 [Google Scholar]
 Struve, O. 1930, ApJ, 72, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Tassoul, J.L. 1978, Theory of rotating stars (Princeton University Press) [Google Scholar]
 Tayler, R. J. 1973, MNRAS, 165, 39 [NASA ADS] [Google Scholar]
 Thomas, R. N. 1965, Some aspects of nonequilibrium thermodynamics in the presence of a radiation field (University of Colorado Press) [Google Scholar]
 Thomas, R. N. 1983, NASACNRS Monograph Series on Nonthermal Phenomena in Stellar Atmospheres, NASA SP471 [Google Scholar]
 Triana, S. A., Moravveji, E., Pápics, P. I., et al. 2015, ApJ, 810, 16 [NASA ADS] [CrossRef] [Google Scholar]
 Tsumuraya, F., Miura, N., & Baba, N. 1994, A&A, 282, 699 [NASA ADS] [Google Scholar]
 Underhill, A., & Doazan, V. 1982, B Stars with and without emission lines, parts 1 and 2, NASA SP456 [Google Scholar]
 Underhill, A. B. 1968, Bull. Astron. Inst. Netherlands, 19, 526 [NASA ADS] [Google Scholar]
 Vinicius, M. M. F., Zorec, J., Leister, N. V., & Levenhagen, R. S. 2006, A&A, 446, 643 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 von Zeipel, H. 1924, MNRAS, 84, 665 [NASA ADS] [CrossRef] [Google Scholar]
 Yudin, R. V. 2001, A&A, 368, 912 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Zhao, M., Monnier, J. D., Pedretti, E., et al. 2009, ApJ, 701, 209 [NASA ADS] [CrossRef] [Google Scholar]
 Zorec, J. 1986, Thèse d’État, Université Paris VII, France [Google Scholar]
 Zorec, J. 1994, in Pulsation, Rotation, and Mass Loss in EarlyType Stars, eds. L. A. Balona, H. F. Henrichs, & J. M. Le Contel, IAU Symp., 162, 257 [Google Scholar]
 Zorec, J. 2000, in The Be Phenomenon in EarlyType Stars, eds. M. A. Smith, H. F. Henrichs, & J. Fabregat, IAU Colloq., 175, ASP Conf. Ser., 214, 51 [Google Scholar]
 Zorec, J., Divan, L., Mochkovitch, R., & Garcia, A. 1987, in Physics of Be Stars, eds. A. Slettebak, & T. P. Snow, IAU Colloq., 92, 68 [Google Scholar]
 Zorec, J., & Royer, F. 2012, A&A, 537, A120 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Zorec, J., Mochkovitch, R. A., & Garcia, A. 1990, in Angular Momentum and Mass Loss for Hot Stars, eds. L. A. Willson, & R. Stalio, NATO ASIC Proc., 316, 239 [Google Scholar]
 Zorec, J., Domiciano de Souza, A., & Frémat, Y. 2004, in Stellar Rotation, eds. A. Maeder, & P. Eenens, IAU Symp., 215, 21 [Google Scholar]
 Zorec, J., Frémat, Y., & Cidale, L. 2005, A&A, 441, 235 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Zorec, J., Frémat, Y., Domiciano de Souza, A., et al. 2011, A&A, 526, A87 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Zorec, J., Frémat, Y., Domiciano de Souza, A., et al. 2016, A&A, 595, A132 (Paper I) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
Appendix A: Binarity
The total radial velocity of a given point in the surface of a stellar component in a binary system has two components: one is due to the stellar rotation proper and the other to the orbital motion. In Paper I the second component was supposed to induce an additional broadening of lines through a differential Doppler shift. Actually, two points in the stellar surface situated at the same latitudefrom the apparent pole in an axially symmetric object and symmetrically, at both sides of the meridian plane that contains the rotation axis and the line of sight, have the same radial velocity in the direction toward the observer. Thus, the suggested differential Doppler shift does not exist and the effect due to binarity discussed in Paper I on Vsini has to be dismissed (J.Z.).
All Tables
Massage groups formed by the sample of Be stars entering the “observed” distribution Φ(u).
Apparent fundamental parameters deduced according to several assumptions about the structure of gravity darkened atmospheres.
Vsini parameters in km s^{1} determined using line profiles perturbed by tidal interactions in binary systems.
All Figures
Fig. 1 Curves of internal constant angular velocity Ω(r,θ) = const. in the stellar envelope (blue lines). a) Curves corresponding to the condition S = S(j^{2}); b) curves obtained with S = S(Ω^{2}); and c) curves calculated for S = S(ϖ^{2}Ω^{2}). The solutions are valid in the convective envelope, i.e., above the shaded central region (S = specific entropy; j = specific angular momentum; and ϖ = distance to the rotation axis). Ordinates coincide with the rotation axis; abscissas are in the equatorial plane. 

In the text 
Fig. 2 Curves of constant radial velocity C contributing to the rotational Doppler broadening of spectral lines. a) Ccurves for a Maunderlike surface angular velocity law with ∂Ω(θ) /∂θ > 0; b) Ccurves for a law given by Eq. (3) with γ > 0, so that ∂Ω(θ) /∂θ > 0; c) Ccurves for a law given by Eq. (3) with γ < 0, so that ∂Ω(θ) /∂θ < 0. In cases a) and b) the largest Doppler displacement is produced in the limb at the equator, while in case c) is produced in the center of the upper closed curves. Case b) identifies the lines producing the same Doppler displacement, C_{r} = C_{d}, where C_{r} (red dotted line) indicates rigid surface rotation, and C_{d} (red dashed line) a differential surface rotation. In all cases the parent Doppler displacement from rigid rotation is determined by the straight line touching the corresponding C_{d} curve at the equator. 

In the text 
Fig. 3 Rotationally broadened line profiles of an isolated Gaussian line, whose intensity and equivalent width depends on the effective temperature and gravity such as a He i 4471 transition. Calculations include gravitational darkening. The parent nonrotating object has T_{eff} = 22 000 K and log g = 4.0. The angular velocity law is given by Eq. (3) with the same parameters (κ,γ,η,i) as in the corresponding block in Fig. 2, except that ξ has a wider range of values. The case ξ = 0 represents the rigid rotation and the corresponding line profile is highlighted with shading. In all blocks the ordinates are the same. 

In the text 
Fig. 4 Rotationally broadened line profiles of an intrinsic Gaussian absorption line, whose intensity and equivalent width depends on the effective temperature and gravity like the He i 4471 transition. Calculations include gravitational darkening in an object rotating at the equator with u = V/V_{c} = 0.9, whose parent nonrotating counterpart has M = 8.6 M_{⊙} and t/t_{MS} = 0.5. The angular velocity law is given by Eq. (7). The colors indicate the inclination angles, but in this figure only the i = 30° (blue) corresponding to different values of α are explicitly identified. 

In the text 
Fig. 5 a) FWHM in km s^{1} observed at different angles i of the HeI 4471 line broadened by differential rotation characterized by different values of the parameter α in an object of “pnrc” parameters T_{eff} = 23 000 K and log g = 4.1 rotating at different rates η. In abscissas are the FWHM of the same line corresponding to rigid rotators having the same rates η. b) Idem for the MgII 4481 line. 

In the text 
Fig. 6 a)−c) Fourier transforms (FT) of the HeI 4471 line broadened by differential rotation characterized by α = −0.5, 0.0 and +0.5, respectively, in an object of “pnrc” parameters T_{eff} = 23 000 K and log g = 4.1 rotating at η = 0.8. d)−f) Idem for the MgII 4481 line. 

In the text 
Fig. 7 a) Measured Vsini parameters in km s^{1} by means of the FT technique of the HeI 4471 line broadened by differential rotation characterized by different values of the parameter α in an object of “pnrc” parameters T_{eff} = 23 000 K and log g = 4.1 rotating at different rates η. In abscissas are the true equatorial Vsini for the same rates η. b) Idem for the MgII 4481 line. 

In the text 
Fig. 8 a) Deviation parameter δ(α,u,i) as a function of α for u = V_{e}/V_{c} = 1.0, mass M/M_{⊙} = 8.6, and fractional age t/t_{MS} = 0.5, and for several inclination angles 5° < i < 90°. b) Inclinationaveraged deviation ⟨ α(α,u) ⟩ for velocity ratios u = V_{e}/V_{c} = 0.4, 0.7, 1.0 in objects with masses M/M_{⊙} = 3.1, 8.6, 13.2 at fractional age t/t_{MS} = 0.5. c) Inclinationaveraged deviation ⟨ α(α,t) ⟩ for M/M_{⊙} = 8.6 at fractional agest/t_{MS} = 0.2, 0.5, 0.9 rotating at u = 1.0. d) The adopted inclinationmassageaveraged deviation ⟨ δ ⟩ as a function of α for different values of the ratio u. 

In the text 
Fig. 9 a) Assumed occurrence probability functions φ(α) of the differential rotation parameter α with a_{α} = 4.5; α_{0}< 0 (blue curves), α_{0} = 0 (black curve), α_{0}> 0 (red curves). b) Idem as in a), but for a_{α} = 250.0. a’) Occurrence probability functions π(δ) of the differential rotation parameter δ with a_{α} = 4.5 for different values of α_{0} and true velocity ratios u. b’) Idem as in a’), but for a_{α} = 250.0. a”) Distributions Ψ(v) of the equatorial velocity ratios v = V/V_{c}, where v is now affected by the differential rotation according to the distributions φ(α) of the parameter α given in a). b”) Idem as in a”), but for φ(α) shown in b). In a”) and b”) is also shown the assumed distribution Φ(u_{eq}) (magenta crosses). 

In the text 
Fig. 10 Angular velocity ratio Ω/Ω_{c} as a function of η, where the uncertainty bars indicate differences for masses and ages in the mainsequence phase. It is shown that Ω/Ω_{c} ≃ 0.95 (red line) identifies the force ratio η ≃ 0.6 at the equator, which implies a 40% undercritical rotation. The ratio V/V_{c} ≃ 0.95 (red line) corresponds to η ≃ 0.86. 

In the text 
Fig. 11 a) Effective temperature ratio T_{eff}/T_{0} against the colatitude θ derived by Zorec et al. (in prep.) using Eq. (19) for η = 0.8 and several values of the differential rotation parameter α. The paramater T_{0} is the effective temperature of the stellar counterpart without rotation. b) Ratio against η and for the same values of α as in a). 

In the text 
Fig. 12 Rotationally broadened Gaussian line (assimilated to He i 4471) by a model star in rigid rotation with η = 0.9 and pnrc parameters T_{eff} = 16 000 K and log g = 4.0. The dotted lines indicate no GD effect; the dashed lines indicate a GD effect using the von Zeipel formulation with β_{1} = 1.0; and the solid line indicates the GD effect using EspinosaLara & Rieutord’s (2011) formulation. The colors indicate the inclination angle i. 

In the text 
Fig. 13 Upper panel: diagram showing the degree (in terms of relative equivalent width) of local contribution to the formation of the He i 4471 line in a rigid critical rotator. Lower panel: degree of local contribution to the formation of the Mg i 4481 line. The white stars identify the pole and the equator. 

In the text 
Fig. 14 Rotationally broadened He i 4471 and Mg ii 44781 line profiles produced by geometrically deformed model stars characterized by and and rotating at η = 0.9. Red lines indicate hybrid GD with local emerging bolometric fluxes modulated according to Eq. (18), θdependent effective gravity and local source functions determined by a uniform electron temperature T_{e} = ⟨ T_{e} ⟩ (assumption “a”). Black lines indicate no GD and θdependent effective gravity (assumption “b”). Blue lines indicate classical GD with β = 1 and the local source functions are determined by T_{e} [T_{eff}(θ),log g(θ)] (assumption “c”). 

In the text 
Fig. 15 a) Tested asymmetric rotation broadening functions (ARBF). b) Broadened He i 4471 line profiles with ARBF having a = −1 and for several input Vsini parameters. c) FT of the ARBF. d) Resulting (Vsini)_{A} parameters produced by ARBF against the model input (Vsini)_{0} (the colors correspond to the distortion parameters a in a)). 

In the text 
Fig. 16 a) Tested truncated rotationexpansion broadening functions (REBF) for several velocity ratios . b) Loci of points with the same projected velocity in the interval of Doppler displacements that contribute to the broadening of line profiles when v_{exp} = 0.24. c) He i 4471 line profiles broadened for several input rotation velocities V_{rot}sini and velocity ratio V_{exp}/V_{rot} = 1.0. d) FT of the REBF. e) Resulting Vsini_{exp + rot} parameters produced in the frame of an expanding and rotating atmosphere against the model input V_{rot}sini for several ratios v_{exp}. 

In the text 
Fig. 17 a) Model line profiles of the primary star broadened by rotation and perturbed by tidal interactions in a binary system with subsynchronous rotation (Ω_{∗} stellar angular velocity; Ω_{b} circular orbital angular velocity). b) FT of spectral lines shown in a). c) Model line profiles of the primary star broadened by rotation and perturbed by tidal interactions in a binary system with supersynchronous rotation. d) FT of spectral lines shown in c). The red lines are either for the averaged line profiles, which represent the tidally unperturbed rotationally broadened line profile, or for the FT of the averaged line profile. 

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.