On the shape of rapidly rotating stars
J.-P. Zahn1 - C. Ranc1 - P. Morel2
1 - LUTH, Observatoire de Paris, CNRS UMR 8102, Université Paris Diderot, 5 place Jules Janssen, 92195 Meudon, France
2 - Université de Nice-Sophia Antipolis, CNRS UMR 6202, Laboratoire Cassiopée, Observatoire de la Côte d'Azur, BP 4229, 06304 Nice Cedex 04, France
Received 7 December 2009 / Accepted 17 April 2010
Aims. The critical surface of a rapidly rotating star is determined, assuming that the rotation is either uniform or shellular (angular velocity constant on level surfaces, but increasing with depth).
Methods. A step beyond the classical Roche model, where the entire mass is assumed to be gathered at the center of the star, here the quadrupolar moment of the mass distribution is taken into account through a linear perturbation method.
Results. The flattening (defined here as the ratio between the equatorial and the polar radius) can somewhat exceed the 3/2 value of the Roche model, depending on the strength of the interior rotation. The result is applied to a star of 7 solar masses, which is the mass of Achernar, the star with the largest flattening detected so far through optical interferometry.
Key words: stars: imaging - stars: rotation - stars: individual: Achernar ( Eri)2001; Domiciano de Souza et al. 2003, 2004). Because that flattening is due to the centrifugal force, one may then ask whether its measure can constrain the interior rotation, and specifically the type and amount of differential rotation.
The question was raised in particular when Domiciano de Souza et al. (2003) announced that the Be star Achernar ( Eri) had an oblate shape: projected to the sky, the ratio between the major and the minor axes of the stellar disk amounted to 1.56. This meant that the ratio between the equatorial and the polar radius was even higher, thus substantially exceeding the value of 3/2 predicted for an ideal star whose entire mass is gathered at the center and which rotates uniformly at break-up speed. But these high flattening ratios are easily achieved when the star is rotating differentially, with the angular velocity decreasing outward, as it was demonstrated by Jackson et al. (2004).
Afterwards Achernar was found to emit a jet, and the ratio between the major and minor axes of the projected star was reduced to 1.41 (Kervella & Domiciano de Souza 2006). This value - the highest measured so far in any star - is compatible with a uniformly rotating star, even considering that the rotation axis is inclined with respect to the line of sight. One may thus be tempted to conclude that the question is settled. However other stars will be analyzed in the future, and we should be prepared to confront their flattening values with those predicted by models.
For this reason we wish to determine here which is the highest flattening allowed for a realistic star that is in uniform rotation or, more generally, in shellular rotation (where the angular velocity is constant on level surfaces, but varies with depth). By ``realistic'' we mean that we take into account the mass distribution distorted by the centrifugal force, which contributes a quadrupolar moment to the gravitational field exerted by the star.
Here designates the density, the gravitational potential, the angular velocity, and the (vector) distance from the rotation axis; the effective gravity is the sum of the gravity and of the centrifugal force. The isobars coincide with the level surfaces, which by definition are orthogonal to the effective gravity.
If the angular velocity depends only on s (i.e. if
is constant on cylinders), the centrifugal force derives from a
potential, and the total potential (gravitational + centrifugal) is
Then the isobars coincide with the equipotentials of the total field , as do the surfaces of constant density; therefore the pressure is a function of density only: - the star has a barotropic structure. In the particular case of uniform rotation, we retrieve the familiar expression
In the general case, is a function also of the coordinate z along the rotation axis; then the density is no longer constant on isobars, but varies with latitude, as can be seen by taking the curl of Eq. (1):
the star is said to be in a baroclinic state.
Here we are mainly interested in the ``surface'' of the rotating star. In all rigor, it should be defined as the surface on which the optical thickness (for a given wavelength or some average of it) is unity (or 2/3, or , depending on the treatment of radiative transfer). For most purposes, however, one may simply assume that this ``surface'' coincides with a level surface of suitably chosen (low) pressure.
The outermost level surface a rotating star can fill is that on which the effective gravity vanishes in the equatorial plane; there the centrifugal force is thus in exact balance with the gravity.1992), but also on more physical grounds, because differential rotation tends to be smoothed out in latitude through shear turbulence.
One can then still define a generating function inspired by Eq. (3)
here is no longer a potential, but because
the surfaces of constant coincide with the isobars, as we can check by comparing with Eq. (1):
The level surfaces of shellular rotation thus have the same shape as for uniform rotation, and this is true in particular for the ``surface''. That property was first pointed out by Meynet & Maeder (1997).
The expression in factor of
is constant on an isobar, and we can thus write
where is the coordinate of that isobar on the polar axis; note that this expression is just another form of the baroclinic relation (4).
here r2 = s2 + z2, M is the mass of the star, and G is the gravitational constant.
The generating function (5) is then
and it allows us to calculate the flattening of the critical surface. At the equator of that surface, where and , it takes the value
and at its poles
hence, because is constant on an isobar,
a value that applies both to uniform and to shellular rotation.
From the generating function (10) we can deduce the equation of the critical surface: scaling the coordinates by the equatorial radius, i.e.
Finally, we may define the mean critical radius R0 as the average of r over the solid angle seen from the center:
Replacing the integrant in terms of and , we obtain
and thus .
We now take into account that the mass distribution is actually oblate,
due to the centrifugal force; this generates a gravitational potential
outside the star which may be expanded in a succession of multipoles as
Here is the colatitude, the Legendre polynomial of order (an even number for obvious symmetry reason), and is a non-dimensional constant which measures the strength of that multipole. R0 designates the radius of the spherically symmetric reference model, in which the horizontal average of the vertical component of the centrifugal force has been subtracted from the radial gravity force (cf. Kippenhahn et al. 1970); as is well known, this has the effect of increasing the radius.
To lowest order, we keep only the quadrupolar term (), and
the generating function is then
As before, we evaluate it at the equator
and at the poles
Next we replace the last term above by its expression drawn from the hydrostatic balance at the equator
to obtain the flattening ratio to lowest order in J2:
where we have approximated and by their values given in Eq. (16) for the critical Roche surface. Again, this result applies to stars that are either in uniform or in shellular rotation. As expected, the value of the critical flattening is increased beyond 3/2 when the quadrupolar moment is taken into account. 1950). We expand the centrifugal acceleration in spherical functions
and proceed likewise for the perturbations of the pressure p', of the density and of the gravitational potential , so that
Next we introduce these expressions in the hydrostatic equation Eq. (1) and expand it in :
where and g0(r) are the local density and gravity of the unperturbed star. Through the elimination of we obtain an expression for the density perturbation , which is introduced in the Poisson equation to yield
Here we are interested in the special case of shellular rotation, where depends only on r, to first approximation. Then, for ,
Rescaling the radial coordinate by the radius, x= r/R0, the angular velocity profile by its surface value, , and the potential perturbation as , we obtain the following Poisson equation (cf. Mathis & Zahn 2004):
To ensure regularity at x=0 and , the solution of this second order o.d.e. must satisfy the boundary values
The quadrupolar moment is given by the surface value of :
thus when the star rotates at critical speed.
Note that this quadrupolar moment is closely related to the apsidal motion constant k2, when the star is in uniform rotation. The reason is that the centrifugal potential then has the same functional behavior as the tidal potential in a binary star component: both scale as (cf. Zahn 1966). One easily finds that in this case J2 = 2 k2/3.
To illustrate the impact of including the quadrupolar moment, we numerically solved this o.d.e. (29) for a 7
star, which according to its spectral type should be about the mass of
Achernar, the star with the highest flattening detected so far. We
adopted the following interior rotation profile:
x=r/R0 being the radial coordinate scaled by the mean radius. For the present purpose, this bell-shaped profile sufficiently resembles those predicted by models including the transport of angular momentum through rotational mixing (cf. Talon et al. 1997; Meynet & Maeder 2000); it imposes a contrast of between the central and surface values of .
The spherical reference models were built with the stellar evolution code CESAM (Morel 1997) for a standard initial composition (X=0.70, Z=0.02). These models define the mean radius R0 and provide the density and mass distribution that enter in the perturbed Poisson Eq. (29). The hydrostatic equation was modified to include the horizontal average of the vertical component of the centrifugal force, as was done by Kippenhahn et al. (1970). We used two rotation profiles (32): one uniform (a=0), and the other with a center-to-surface ratio of 4 in angular velocity (a=15). Three stages of evolution were considered: one at the ZAMS (defined as the location in the HR diagram where half of the luminosity originates from nuclear reactions and half from the star's contraction), the other at mid-MS (where the central hydrogen concentration has been reduced to ), and the third at the end of the MS ( ).
The results are given in Table 1, which shows how the quadrupolar moment J2 increases with the interior rotation at given age, when the surface rotates at critical speed. For uniform rotation, our values agree with the apsidal motion constants calculated by Claret (1995). From Eq. (23) we then deduce the maximum flattening that such a star can achieve and, combining this with the value of given in Eq. (16), we derive the equatorial and polar radii.
For simplicity, one could be tempted to take as reference model that of the non-rotating star, but at critical speed the swelling and mass redistribution due to the centrifugal force cannot be treated as a mere perturbation. This can be seen in Table 1 by comparing the radius of the non-rotating models with that of the models including the effect of rotation. The discrepancy is even more pronounced for the quadrupolar moments J2, which would be underestimated by a factor 2.5 when using the non-rotating model.
Table 1: Characteristics of a 7 star rotating at critical speed.
The main goal of this paper was to investigate how much the flattening of a star rotating at critical speed can exceed the classical value of 3/2 that characterizes the Roche model. Our resultswere obtained by a linear perturbation method treating the quadrupolar moment as a small quantity, which is justified for massive MS stars, because these are sufficiently centrally condensed. Nevertheless, it would be worthwhile to refine our results through genuine two-dimensional calculations, like those described by Clement (1978) or by Rieutord (2006).
Although their analysis applies strictly only to conservative laws of rotation, Jackson et al. (2004) have convincingly shown that the flattening may assume rather high values when the star is rotating differentially, with the angular velocity decreasing outward. Here we instead considered a star rotating uniformly at critical speed, but we improved upon the Roche model by taking into account the quadrupolar moment of the mass distribution, due the centrifugal force. We found that the flattening of the critical surface is then slightly increased from 1.50 to about 1.52, for a 7 star which was chosen to represent Achernar. This increase would be more pronounced for lesser central condensation, i.e. for MS stars of lower mass.
Then we examined the so-called shellular rotation regime, where the angular velocity is constant on level surfaces, but increases with depth. Due to this ``hidden'' rotation, the flattening of the star can increase to substantially higher values: for example to 1.54 or 1.56, depending on the state of evolution, for a contrast of 4 between its angular velocity at center and surface. This differential shellular rotation was first invoked by Zorec et al. (2005) in their analysis of Achernar; they concluded that the observations were compatible with a center-to-surface ratio for the angular velocity of 2.7, and an inclination of of the rotation axis with respect to the line of sight. Thus, provided the interoferometric determinations become sufficiently precise, the shape of stars may provide a valuable constraint on their internal rotation, complementary to the powerful asteroseismic diagnostic.Acknowledgements
We thank P. Kervella for kindly updating us on the current state of optical interferometry and its application to rotating stars. And we express our gratitude to the referee Y. Eriguchi for his pertinent remarks, which incited us to substantially improve our paper.
- Claret, A. 1995, A&AS, 109, 441 [NASA ADS] [Google Scholar]
- Clement, M. J. 1978, ApJ, 222, 967 [NASA ADS] [CrossRef] [Google Scholar]
- Domiciano de Souza, A., Kervella, P., Jankov, S., et al. 2003, A&A, 407, 47 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kervella, P., & Domiciano de Souza, A. 2006, A&A, 453, 1059 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kippenhahn, R., Meyer-Hofmeister, E., & Thomas, H. C. 1970, A&A, 5, 155 [NASA ADS] [Google Scholar]
- Jackson, S., MacGregor, K. B., & Skumanich, A. 2004, ApJ, 606, 1196 [NASA ADS] [CrossRef] [Google Scholar]
- Maeder, A. 1999, A&A, 347, 185 [NASA ADS] [Google Scholar]
- Mathis, S., & Zahn, J.-P. 2004, A&A, 425, 242 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Meynet, G., & Maeder, A. 1997, A&A, 321, 46 [NASA ADS] [Google Scholar]
- Meynet, G., & Maeder, A. 2000, A&A, 361, 101 [NASA ADS] [Google Scholar]
- Morel, P. 1997, A&ASS, 124, 597 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [PubMed] [Google Scholar]
- Pijpers, F. P. 1998, MNRAS, 297, L76 [NASA ADS] [CrossRef] [Google Scholar]
- Rieutord, M. 2006, in SF2A 2006: Proc. Ann. Meet. French Soc. Astron. Astrophys., 501 [Google Scholar]
- Roxburgh, I. W. 2004, A&A, 428, 171 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Sweet, P. A. 1950, MNRAS, 110, 548 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
- Talon, S., Zahn, J.-P., Maeder, A., & Meynet, G. 1992, A&A, 322, 209 [NASA ADS] [Google Scholar]
- van Belle, G. T., Ciardi, D. R., Thompson, R. R., Akeson, R. L., & Lada, E. A. 2001, ApJ, 559, 1155 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
- Zahn, J.-P. 1966, Ann. Astrophys, 29, 313 [NASA ADS] [Google Scholar]
- Zahn, J.-P. 1992, A&A, 265, 115 [NASA ADS] [Google Scholar]
- Zorec, J., Domiciano de Souza, A., Frémat, Y., & Vakili, F. 2005, in SF2A 2005, ed. F. Casoli, T. Contini, J.-M. Hameury, & L. Pagani (EDP Sciences), 363 [Google Scholar]
- In a subsequent paper, Maeder (1999) took a different definition for the shellular rotation, namely that the angular velocity is constant on spheres: , and then this property no longer holds.
Table 1: Characteristics of a 7 star rotating at critical speed.
Copyright ESO 2010