Issue 
A&A
Volume 539, March 2012



Article Number  A90  
Number of page(s)  12  
Section  Stellar structure and evolution  
DOI  https://doi.org/10.1051/00046361/201118151  
Published online  28 February 2012 
Seismic modelling of the late Be stars HD 181231 and HD 175869 observed with CoRoT: a laboratory for mixing processes^{⋆}
^{1} LESIA, Observatoire de Paris, CNRS UMR 8109, UPMC, Université Paris Diderot, 5 place Jules Janssen, 92190 Meudon, France
email: coralie.neiner@obspm.fr
^{2} Laboratoire AIM, CEA/DSM – CNRS – Université Paris Diderot, IRFU/Service d’Astrophysique, CEASaclay, 91191 GifsurYvette Cedex, France
^{3} Astronomical Institute, Graduate School of Science, Tohoku University, 9808578 Sendai, Japan
^{4} Los Alamos National Laboratory, Los Alamos, NM 87545, USA
^{5} Observatoire de Genève, Université de Genève, 51 chemin des Maillettes, 1290 Sauverny, Switzerland
Received: 26 September 2011
Accepted: 7 January 2012
Context. HD 181231 and HD 175869 are two late rapidly rotating Be stars, which have been observed using highprecision photometry with the CoRoT satellite during about five consecutive months and 27 consecutive days, respectively. An analysis of their light curves, by Neiner and collaborators and GutiérrezSoto and collaborators respectively, showed that several independent pulsation gmodes are present in these stars. Fundamental parameters have also been determined by these authors using spectroscopy.
Aims. We aim to model these results to infer seismic properties of HD 181231 and HD 175869, and constrain internal transport processes of rapidly rotating massive stars.
Methods. We used an adiabatic (NRO) and a nonadiabatic (Tohoku) oscillation code that accounts for the combined action of Coriolis and centrifugal accelerations on stellar pulsations as needed for rapid rotator modelling. We coupled these codes with a 2D (ROTORC) stellar structure model to take the rotational deformation of the star into account. The action of transport processes was parametrised with the mixing parameter α_{ov}, which represents the “nonstandard” extension of the convective core, and determined by matching observed pulsation frequencies assuming a single star evolution scenario. In a second step, we used (Geneva) evolution models to evaluate the contribution of the secular rotational transport and mixing processes in the radiative envelope. A Monte Carlo analysis of spectropolarimetric data was also performed to examine the role of a potential fossil magnetic field. Finally, based on stateoftheart modelling of penetrative convection and internal waves, we unravelled their respective contribution to the needed “nonstandard” mixing.
Results. We find that extra mixing of α_{ov} = 0.3−0.35H_{p} is needed in HD 181231 and HD 175869 to match the observed frequencies with those of prograde sectoral gmodes. We also detect the possible presence of rmodes. We investigated the respective contributions of several transport processes to this mixing: the hydrodynamical processes in particular the meridional circulation and shearinduced turbulence caused by the radiative envelope differential rotation, the possible magnetic field, the penetrative convection at the top of the convective core, and the transport by internal waves.
Conclusions. We showed that the extension of the convective core needed to match observations and models may be explained by mixing induced by the penetrative movements at the bottom of the radiative envelope and by the secular hydrodynamical transport processes induced by the rotation in the envelope. We showed how asteroseismology opens a new door to probe transport processes in stellar interiors.
Key words: turbulence / asteroseismology / stars: emissionline, Be / stars: rotation / stars: individual: HD 181231 / stars: individual: HD 175869
© ESO, 2012
1. Introduction
Since the beginning of the theoretical modelling of stellar structure and evolution, many hydro and magnetohydrodynamical processes have been studied. In particular, the impact of the internal rotation and the associated shear has been modelled (e.g. Zahn 1992; Mathis & Zahn 2004) and evaluated (e.g. Meynet & Maeder 2000, for massive stars and Talon & Charbonnel 2005, for lowmass stars). For massive stars, especially for fast rotators close to their breakup velocity such as Be stars, the modification of stellar structure and evolution caused by the rotation has to be taken into account. Indeed, rotation impacts many aspects of stellar evolution such as the stellar lifetime, internal structure (for both the thermic and chemical stratification), chemical yields, etc. However, such an hydrodynamical modelling always implies prescription, in particular for turbulence both in convective and radiative regions, and their limits. Therefore, it is absolutely necessary to obtain constraints from observations of these physical processes and their consequences, i.e. we have to probe stellar interiors to acquire their direct and indirect signatures.
In this context, highprecision asteroseismology from space is the best way to study the internal structure in the HertzsprungRussel (HR) diagram using oscillation mode frequencies, which give the properties of cavities where waves propagate. This is the framework of the present paper, in which we study the two fast rotating late Be stars HD 181231 and HD 175869. They have been observed with the CoRoT photometric satellite (for a description of CoRoT, see Baglin et al. 2006, 2009) during about five consecutive months and 27 consecutive days, respectively. An analysis of their CoRoT light curves showed that several independent periods of variations are present in these stars, which we interpreted in terms of pulsation gmodes (Neiner et al. 2009; GutiérrezSoto et al. 2009). Late Be stars indeed exhibit pulsations of gmodes similar to those of slowly pulsating B (SPB) stars. Therefore they have recently been baptised SPBe stars (Walker et al. 2005). Apart from a slightly different effective temperature, HD 181231 and HD 175869 have similar fundamental parameters. In addition, the amplitude of variations detected in HD 175869 is lower than that in HD 181231, as expected because of its later spectral type. Because the pulsations of HD 181231 and HD 175869 have been characterized with great precision thanks to CoRoT by Neiner et al. (2009) and GutiérrezSoto et al. (2009) and their fundamental parameters were determined from spectroscopy by the same authors, these late Be stars are very good candidates for seismic modelling. Furthermore, because we know that these Be stars rotate with an angular velocity close to their breakup velocity, they are relevant targets for probing and obtaining tight constraints on the physical processes related to their internal rotation.
We therefore studied HD 181231 and HD 175869 using stateoftheart modelling tools to unravel rotation processes using seismology. First, we modelled stellar oscillations using two numerical codes (Lee & Baraffe 1995; Clement 1998) that take into account the Coriolis and the centrifugal accelerations, which play a crucial role in pulsation mode dynamics in fast rotators. This allows us to calculate mode frequencies and amplitudes for each given stellar model. We chose to use classical stellar structure models for the target structures, adding a perturbative Chandrasekhar expansion (Lee & Baraffe 1995) to model stars flattened by the centrifugal force, as well as 2D stellar models (Deupree 1990, 1995) that take into account the flattening in a nonperturbative way.
In both cases of models, the transport mechanisms and the associated “nonstandard” mixing are not treated in a coherent way but are parametrised using the classical overshoot parameter α_{ov} that gives the extension of the induced convective core compared to the “standard” model for a singlestar evolution without rotation. Thanks to the confrontation of mode computations to the CoRoT observational results, we were able to determine the best value of α_{ov} for each target. This motivated us to study transport processes that may modify the observed mixing. To this aim, four mechanisms can be considered: the hydrodynamical processes, in particular the meridional circulation and shearinduced turbulence caused by the radiative envelope differential rotation, the possible magnetic field, the penetrative convection at the top of the convective core, and the transport by internal waves. For a complete review on these processes, we refer the reader to Talon (2008) and Mathis (2010).
In Sect. 2 we summarize the pulsational characteristics and fundamental parameters of the two Be stars. We then present the codes and assumptions we used to model the stars (Sect. 3) and the results we obtained for the two stars (Sect. 4). We interpret (Sect. 5) and discuss (Sect. 6) these results in terms of core overshooting and rotational mixing by comparing them with theoretical predictions and models. Finally, we conclude and discuss perspectives of this study.
2. The stars HD 181231 and HD 175869
Effective temperature (in log), gravity, projected rotational velocity, inclination angle, equatorial radius, mass, luminosity, and rotational frequency of the two late Be stars HD 181231 and HD 175869.
2.1. HD 181231
HD 181231 is a B5IVe star. Its stellar parameters, assuming a ratio between the angular velocity and its critical value , are summarized in Table 1. See also Table 3 of Neiner et al. (2009).
HD 181231 is not known as a binary in the literature. Inspection of its spectra revealed no sign of a companion. In addition, the spectral modelling performed by Neiner et al. (2009) using solar abundances fitted the observed spectra well and provided no sign of anomalous abundances and therefore no sign of former accretion.
Neiner et al. (2009) detected at least ten independent significant frequencies of variations among the 54 frequencies detected in the CoRoT data of HD 181231, which they interpreted in terms of nonradial pulsation gmodes and rotation. Two longerterm variations were also detected: one at about 14 days resulting from a beating effect between the two main frequencies of shortterm variations, the other at about 116 days caused either by a beating of frequencies or by a zonal pulsation mode.
In addition, Neiner et al. (2009) observed HD 181231 with highprecision spectropolarimetry: none of the measurements exibits a magnetic signature, and all longitudinal field measurements are compatible with 0. From the error bars on these measurements, an uppper limit of the possible nondetected longitudinal magnetic field has been set to about 650 G.
2.2. HD 175869
HD 175869 is a B8IIIe star. Its stellar parameters, assuming , are summarized in Table 1 and are very similar to those of HD 181231 except for the cooler temperature. See also Table 3 of GutiérrezSoto et al. (2009).
Like HD 181231, HD 175869 is not known as a binary in the literature, its spectra revealed no sign of a companion, and the spectral modelling performed by GutiérrezSoto et al. (2009) using solar abundances provided no sign of anomalous abundances and therefore no sign of former accretion.
In HD 175869 GutiérrezSoto et al. (2009) detected a frequency of variation compatible with the rotational frequency, as well as six of its harmonics. Lowamplitude periodic variations were also detected. The results were interpreted both in terms of stellar activity just above the photosphere and in terms of pulsation gmodes.
Moreover, GutiérrezSoto et al. (2009) observed HD 175869 using highprecision spectropolarimetry: none of the measurements exibits a magnetic signature, and all longitudinal field measurements are compatible with 0. From the error bars on these measurements, an upper limit of the possible nondetected longitudinal magnetic field has been set to about 400 G.
2.3. Evolutionary scenarios
Two evolution scenarios are proposed to explain the existence of Be stars. First, Be stars can be born as single stars with a high initial rotation rate, which evolves as the star progresses on the main sequence (Martayan et al. 2007; Ekström et al. 2008). Second, Be stars can result from the evolution of a B star in a binary system (Pols et al. 1991; Dewi 2007; de Mink et al. 2011). In this case, the high rotation results from an accretion phase from the companion star during the main sequence. It is not straightforward to distinguish a fast rotator that was born as a rapidly rotating single star from a fast rotator that resulted from binary evolution. Indeed, rapidly rotating stars resulting from these interactions will often appear to be single because the companion tends to be a lowmass, lowluminosity star or compact object in a long period orbit. Nevertheless, the lack of obvious surface abundance anomalies in HD 181231 and HD 175869 favours the singlestar scenario. Most binary evolution scenarios indeed lead to N surface abundance enrichments (see Langer et al. 2008). In the remain of this paper, we therefore assume that the two stars are single stars. A complete study of the binary scenario is beyond the scope of this paper but we provide an outline in Sect. 6.
3. Models of stellar structure, pulsation, and evolution
Modelling the pulsations of a star first requires the calculation of the star’s structure (see Sect. 3.1) and then a pulsation stability analysis (see Sect. 3.2). In addition we used the Geneva stellar evolution code (see Sect. 3.3) to constrain the impact of rotational mixing on stellar structure.
3.1. Stellar structure
The calculation of the stellar structure of Be stars should take into account the deformation of the star through their rapid rotation. For the present seismic study, we modelled the stellar structure of the rapidly rotating late Be stars with two methods: a 2D model and a ChandrasekharMilnelike expansion with a method similar to that of Lee & Baraffe (1995) but slightly modified with the 2D model as input. The nonstandard mixing (and transport) processes, which are not treated, are parametrised by an overshoot region at the top of the convective core, whose length is defined as d_{ov} ≡ α_{ov}H_{p}.
3.1.1. 2D rotational deformation: ROTORC
The 2D structure models were calculated with the rotorc code (Deupree 1990, 1995). This code solves the conservation equations for mass, momentum, energy and composition, as well as Poisson’s equation on a 2D spherical grid using a Henyey method. We made the models rotate uniformly on the zero age main sequence (ZAMS), and the interior angular momentum is conserved locally as the model evolves. The independent variables are the fractional surface equatorial radius and the colatitude, and the surface was assumed to be an equipotential. In this way, we needed no assumptions about the shape of the stellar surface, and we obtained the rotationally flattened structure of the star. In this model we used the OPAL opacities (Iglesias & Rogers 1996).
3.1.2. ChandrasekharMilne expansion: the Tohoku code
In Lee & Baraffe (1995) an isobaric surface is expressed as r = a [1 + ϵ(a,θ)] where ϵ(a,θ) is a function of Ω^{2}. Here solidbody rotation is assumed (Ω is uniform). Lee & Baraffe (1995) then define ϵ(a,θ) so that the equipotential surface and physical quantities in the equilibrium state depend only on the mean radius a. This allows a mapping between the structure of the rotating and nonrotating stars in terms of a, and ϵ(a,θ) is determined by the ChandrasekharMilne expansion (Chandrasekhar 1933a,b). This mapping is only valid when the difference between the rotating and nonrotating structures can be approximated by the Ω^{2} order corrections involving the function ϵ(a,θ), i.e. when the rotation velocity is not too fast. In these conditions, ϵ(a,θ) can be rewritten as ϵ(a,θ) = α(a) + β(a)P_{2}(cosθ) where P_{2} = (3cos^{2}θ − 1)/2. Consequently, r = a [1 + α(a) + β(a)P_{2}(cosθ)] . The functions α and β represent the horizontal averaged effects of the centrifugal force and the deformation of the equilibrium state, respectively. See Lee & Baraffe (1995) for more details.
Since the late Be stars we are studying here rotate close to their critical velocity, the ChandrasekharMilne approximation described above may not be adequate. Therefore we chose to use the 2D structure described in the previous section as an input to our 1D models. This is a first step towards nonadiabatic seismic modelling of rapidly rotating stars (see Sect. 6 for the discussion of the associated drawbacks).
Therefore we assume the isobaric surface to be expressed as , because the unperturbed 1D models already include the effect of the spherical symmetric component of the centrifugal force. The value of for a given rotation frequency is obtained by interpolating the coefficients determined from the 2D ROTORC equilibrium structures. In this model we use OP opacity tables (Badnell et al. 2005). However there is almost no difference in the results when using OPAL opacities.
3.2. Oscillation codes
In the rapidly rotating case, a single l value does no longer represent a pulsation mode for a given m value because of the latitudinal couplings between the different spherical harmonics due to the Coriolis and the centrifugal acceleration. For this reason the eigenfunctions for a fixed m value are expanded into a series of terms proportional to spherical harmonics Y_{l,m}(θ,ϕ) with l ≥ m. Indeed, we express the spatial dependence of the displacement vector ξ and a scalar variable f′ as (1)and (2)where with I = 0 for even modes (in which scalar variables are symmetric with respect to the equator) and I = 1 for odd modes. The terms proportional to T^{j} represent a toroidal displacement needed because of the Coriolis force term in the momentum equation. Obviously a single latitudinal degree l of is not a good parameter anymore for a rotating star. We designate the type of the latitudinal dependence of a mode by using ℓ, which corresponds to the latitudinal degree l of the component with highest amplitude at the stellar surface. This can be best compared to the ℓ values derived from the CoRoT observations.
Furthermore, an excited nonradial mode, which has a frequency ν_{0} in the corotating frame, should be observed at the frequency ν_{0} − mΩ, where we adopt the convention that prograde modes have a negative azimuthal order m.
3.2.1. NRO adiabatic oscillation code
Clement (1998) developed an adiabatic code (Normal modes of Oscillation for Rotating stars, hereafter NRO) focusing on gmodes of SPB stars. This code is appropriate for late Be stars, which are located in the same instability strip as SPB stars.
We used a modified version of this NRO code that includes differential rotation terms. The variables used in the models are the radial and latitudinal components of the Lagrangian displacement, the Eulerian pressure, and the gravitational potential perturbations. These variables are evaluated on level surfaces and normalized to have nonzero values at the boundaries. The linearized adiabatic pulsations are solved on a finite difference grid along several radii, i.e. including several spherical harmonics.
Convergence cannot always be reached with the NRO code when the mode frequency is lower than twice the angular velocity of the star (i.e. the inertial frequency), and the modes calculated below 2Ω should therefore be considered with caution. This makes it impossible to compute highorder gmodes (n > 15) in the presence of fast rotation.
3.2.2. Tohoku nonadiabatic oscillation code
To be able to safely model frequencies below 2Ω, we used a nonadiabatic oscillation code, coupled to the Tohoku structure code (see Sect. 3.1.2), to model the nonradial pulsations of the two rapidly rotating Be stars. This code includes the effects of the centrifugal acceleration (which breaks the spherical symmetry of the star as described above) and of the Coriolis acceleration, which have to be taken into account in our seismic modelling since (see Lee 2008; Reese 2010).
With this nonadiabatic oscillation code, we performed a nonradial pulsation stability analysis for rapidly rotating star models, whose evolutionary track passes close to the position of the Be stars in the HR diagram. The evolutionary models were computed in the same way as in Walker et al. (2005) and the analysis allows us to predict the modes that will be excited by the κmechanism as well as their frequencies and growth rates (η). The effects of the Coriolis and the centrifugal acceleration are included, although the effect of the deformation is low for lowfrequency modes. Moreover, we recall that the rotation is assumed to be uniform.
3.3. Geneva stellar evolution code including transport processes
In the Geneva stellar evolution code, the effects of rotation are included according to the shellular rotation prescription in radiation zones, which assumes a strongly anisotropic turbulence leading to an essentially constant angular velocity on the isobars in these regions (Zahn 1992). This assumes a closure turbulent model for the horizontal turbulence for which we must choose a prescription as for any hydrodynamical model that consistently treats largescale motions. Here we used prescriptions by Zahn (1992). Moreover, the rotation is assumed to be uniform in convection zones.
With this hypothesis of shellular rotation, every quantity depends solely on pressure and can be split into a mean value and its latitudinal perturbation, i.e. , where P_{2}(cosθ) is the second Legendre polynomial. Then, the equations of stellar structure need to be modified to account for the hydrostatic effects of rotation because of the centrifugal acceleration (Kippenhahn & Thomas 1970). The usual spherical coordinates are then replaced by new coordinates characterizing the equipotentials. The Kippenhahn & Thomas (1970) method is only applicable in conservative cases, while the problem is nonconservative when the rotation law is shellular. In this case, however, the equations of stellar structure can still be written consistently in terms of a coordinate referring to the mass inside the isobaric surfaces (see Meynet & Maeder 1997). We used the Roche model in the computations. This approximation uses the gravitational potential emerging from the spherical distribution of mass for computing the gravitational force. Note that models including shellular rotation show an internal rotation profile far below the critical limit when they reach the critical velocity at the surface. Consequently, the conditions required for the Roche approximation to be valid are fulfilled in these models. This is not because the Roche approximation is imposed, but because the various physical processes affect the angular momentum distribution in the stellar interior (Meynet et al. 2010).
The transport and mixing processes related to (differential) rotation were treated in the framework of shellular rotation (Zahn 1992; Maeder & Zahn 1998). Then, a largescale meridional circulation was generated by structural adjustments of the stars and the turbulent viscous transport owing to the secular shear instability. Next, because of the heat advection by this circulation, stellar radiation zones relax onto a new thermal state that generates a new differential rotation, as explained in Rieutord (2006) and Decressin et al. (2009). In this framework, the transport of angular momentum obeys an advectiondiffusion equation, while the vertical transport of chemicals through the combined action of vertical advection and strong horizontal turbulence can be described as a pure diffusive process (Chaboyer & Zahn 1992). Because the description of modelling these rotational mechanisms has already been presented in detail in previous papers (see e.g. Meynet & Maeder 2000), we only recall here that the Geneva code includes rotational transport and mixing processes where the meridional circulation is treated as a truly advective process (see Eggenberger et al. 2008). The meridional circulation as well as thermal and chemical fluctuations induced by the 1D shellular rotation are described in 2D using projections on spherical harmonics (see e.g. Meynet & Maeder 2000; Decressin et al. 2009).
Fig. 1 Estimated positions of HD 181231 (log T_{eff} = 4.155) and HD 175869 (log T_{eff} = 4.079) on the log T_{eff} − log g plane (black dots) with evolutionary tracks of M/M_{⊙} = 3,4,5,6,7 with a standard chemical composition of (X,Z) = (0.07, 0.02). Dotted lines are the evolutionary tracks of M/M_{⊙} = 4,4.5,5,5.5,6 with a convectivecore overshooting of 0.35H_{p}. There is no rotation included in these tracks. Open squares and circles indicate the positions of models for HD 181231 and HD 175869, respectively (see below). 
4. Seismic modelling of the two late Be stars
4.1. Evolutionary stage
Unperturbed 1D Tohoku models for the pulsation analysis were obtained from spherical symmetric evolutionary models including averaged centrifugal force. Figure 1 shows evolutionary tracks in the log T_{eff} − log g plane, where solid lines represent evolutionary tracks without convectivecore overshooting and dotted lines are those including a core overshooting of 0.35H_{p}.
In a postmainsequence star, gmodes are damped by strong radiative damping in the dense radiative core where the spatial wavelengths of gmodes become very short. This means that for gmodes to be excited in a Btype star, the star must be in the mainsequence phase with a convective core. Because HD 181231 and HD 175869 show gmode oscillations, they should be in the mainsequence phase. Their positions on the log T_{eff} − log g diagram in Fig. 1 indicate that considerable mixing should be occurring in the radiative layers around the convectivecore boundary.
4.2. HD 181231
4.2.1. Tohoku models without mixing
We computed models of HD 181231 with the stellar structure and Tohoku oscillation codes described above. We first calculated models without mixing. We adopted a mass of 5 M_{⊙}, which is consistent with the spectroscopically estimated mass range (see Table 1) and marginally consistent with the position on the log T_{eff} − log g diagram shown in Fig. 1. Three values of rotation frequencies (f_{rot} = 0.432, 0.778 and 1.037 c d^{1}) were examined. These values lie within the spectroscopically estimated range of f_{rot} = 0.76 ± 0.33 c d^{1}.
Figure 2 compares frequencies of HD 181231 observed with CoRoT (bottom panel) with frequencies of lowdegree (<2) nonradial modes excited in models with log T_{eff} ~ 4.15 (near the end of the main sequence). The azimuthal order m is also shown for one of the models. It shows that in the most slowly rotating model (f_{rot} = 0.432 c d^{1}), frequencies of excited gmodes cover the observed frequency range but do not show the observed frequency groupings.
Fig. 2 Frequencies of HD 181231 observed with CoRoT (bottom panel) and growth rate of the modelled frequencies of low degree (ℓ ≤ 2) (first, third and fourth panels from the top) excited in a 5 M_{⊙} star without mixing, calculated with the Tohoku code. The second panel from the top shows the azimuthal order m of the pulsation mode corresponding to the modelled frequencies in the third panel. Colours indicate the parity of the mode: even modes in solid red and odd modes in dashed blue lines. Stellar parameters used in the three models are indicated at the top of each panel. 
Frequency groupings become visible when most of the tesseral gmodes are damped, which occurs when , where R is the mean radius. The critical rotation corresponds to . Since f_{rot} = 0.432 c d^{1} corresponds to , no frequency groupings occur.
In more rapidly rotating models (f_{rot} = 0.778 and 1.037 c d^{1}), frequency groupings indeed appear because most of the tesseral g modes (ℓ − m ≠ 0) are damped. The grouped frequencies, which consist of sectoral gmodes of m = −1 and − 2, are too high compared to the observed ones however.
Note that among the lowest frequency around 0.1 c d^{1}, some frequencies are combinations of higher frequencies or multiples of the frequency resolution (see Neiner et al. 2009). The others can be modelled with rmodes with m = 1.
Fig. 3 Frequencies of HD 181231 observed with CoRoT (bottom panel) and growth rate and azimuthal order of the modelled frequencies of low degree (ℓ ≤ 2) (top panels) excited in a 5.5 and 6 M_{⊙} star with a core overshooting of 0.35 H_{p}, calculated with the Tohoku code. Colours indicate the parity of the mode: even modes in solid red and odd modes in dashed blue lines. Stellar parameters used in the two models are indicated at the top of each model. 
To obtain frequency groupings at f_{rot} ~ 0.4–0.5 c d^{1}, we used models with a larger radius (or a lower surface gravity). Since the models adopted in Fig. 2 are close to the termination of the conventional main sequence, we adopted models beyond the mainsequence band by including extra mixing in the radiative layers. This is consistent with the positions of the two stars in the log T_{eff} − log g diagram (Fig. 1).
4.2.2. Tohoku models with mixing
We then introduced mixing in the models. We calculated models in the effective temperature range log T_{eff} = 4.155 ± 0.046. Figure 3 shows two models of 5.5 and 6.0 M_{⊙}, which roughly reproduce the frequency groups around 0.6 and 1.2 c d^{1} with prograde sectoral gmodes having m = −1 and − 2, respectively. Both models are near the end of the mainsequence phase with a core overshooting of 0.35 H_{p}. The positions on the log T_{eff} − log g plane are shown by open squares in Fig. 1.
The masses used in the two models lie in the upper part and slightly above the limit of the spectroscopically determined range M = 4.4 ± 1.3 M_{⊙}. While similar fits are obtained with slightly more massive models, no good fits are obtained with M < 5 M_{⊙}, i.e. our seismic modelling allows us to narrow the possible mass range of HD 181231 down to 5 < M/M_{⊙} < 6.5.
The rotation frequencies are 0.475 c d^{1} in the 5.5 M_{⊙} model and 0.518 c d^{1} in the 6.0 M_{⊙} model, both of which are within the spectroscopically determined range f_{rot} = 0.76 ± 0.33 c d^{1} and were chosen to qualitatively fit the observed frequency groups. The normalized rotation frequency is 0.52 and 0.62 for the 5.5 and 6 M_{⊙} models, respectively. These correspond to ~75 and 90% of the critical rotation.
The models do not reproduce the comparatively high frequencies observed at 1.87 and 2.5 c d^{1}. These observed frequencies are indeed aliases of the main frequency at 0.62 c d^{1}, as noted in Neiner et al. (2009).
4.2.3. NRO models
We calculated three NRO models with a full 2D ROTORC stellar structure, with 4, 4.5 and 5 M_{⊙} uniformly rotating on the ZAMS with velocities of 200, 250 and 300 km s^{1}. The advantage of NRO compared to Tohoku models is that they are fully 2D. The drawback is that the NRO code is adiabatic. It is therefore useful to compare the results obtained with both methods.
Fig. 4 Frequency match to HD 181231 calculated using NRO models. Solid lines indicate observed frequencies, arbitrarily attributed to ℓ = 2. Dashed lines indicate the model frequencies and their value of ℓ. In this case, the frequencies are best matched by modes with ℓ = 1 to 5, and only one mode is axisymmetric. 
Each NRO model was evolved with a convective core overshoot of 0, 0.15 and 0.3 H_{p}. The evolution sequence can be seen in Fig. 1 of Lovekin et al. (2010). For each model, we then computed a χ^{2} value by comparing the modelled frequencies to the observed ones. We first fitted axisymmetric modes only and then included also nonaxisymmetric modes. Details about these models are published in Lovekin et al. (2010).
With axisymmetric modes, the best model has M = 5 M_{⊙}, is rotating at 105 km s^{1}, has an overshoot of 0.3 H_{p}, and was evolved to an age of 91 Myr. This model has log T_{eff} = 4.156,R = 5.33 R_{⊙} and f_{rot} = 0.389 c d^{1}. While T_{eff} and R fit the spectroscopically determined parameters well, and M = 5 M_{⊙} is compatible with the results on mass obtained with the Tohoku models above, f_{rot} is lower than the observed f_{rot} = 0.76 ± 0.33 c d^{1}.
Including nonaxisymmetric modes, the bestfitting model as M = 5 M_{⊙}, is rotating at 110 km s^{1}, has an overshoot of 0.15 H_{p}, and was evolved to an age of 84 Myr. This model has log T_{eff} = 4.154,R = 5.10 R_{⊙} and f_{rot} = 0.426 c d^{1}. In this case, f_{rot} becomes marginally compatible with the spectroscopically determined value.
Figure 4 shows the match between the modelled frequencies and the ten independent frequencies observed in HD 181231 (Neiner et al. 2009) and gives an indication of the ℓ values. From the CoRoT observations, only the ℓ value for the frequency 0.695 c d^{1} could be estimated: ℓ ~ 3 (Neiner et al. 2009). Here we obtain ℓ = 4 from the NRO models for a frequency at 0.708 c d^{1}. Because the NRO code is adiabatic, the results show all possible modelled modes without an indication of their excitation or damping. In reality only some of the modelled frequencies would indeed be visible in the CoRoT observations.
4.3. HD 175869
4.3.1. Tohoku models with mixing
We modelled HD 175869, which is slightly cooler than HD 181231, in the same way with the Tohoku code. Reasonable models were searched for in the effective temperature range log T_{eff} = 4.079 ± 0.022. Figure 5 shows 4.5 and 5.0 M_{⊙} models, in which sectoral prograde gmodes of m = −1 and − 2 are excited at frequencies around ~0.55 and 1.3 c d^{1} with a core overshooting of 0.35 H_{p}. The rotation frequencies of the two models, f_{rot} = 0.475 and 0.518 c d^{1}, correspond to about 70% and 80% of the critical rotation.
The mass range in which similar agreements can be obtained in the above effective temperature domain is 4.5 < M/M_{⊙} < 5.5. This range is consistent with, and considerably narrower than, the spectroscopically determined range (4.4 ± 1.3 M_{⊙}).
In HD 175869, only four of the observed frequencies are independent ones: 0.639, 1.451, 2.0 and 2.475 c d^{1}. Other observed frequencies are harmonics or combination of these four frequencies (see GutiérrezSoto et al. 2009). The models do reproduce the lower frequencies, but the highest frequencies at 2.0 and 2.475 c d^{1} cannot be explained by the models.
4.3.2. NRO models
The NRO models were calculated for HD 175869 in the same way as presented above for HD 181231. However, with only four independent frequencies and considering the high density of gmodes in these adiabatic models, many models with different parameters but with similar and reasonable χ^{2} values are possible fits of the four frequencies. The NRO models are therefore not well constrained in this case. Nevertheless, we found that all models that fitted the observations with a low χ^{2} require a high overshoot of 0.3 H_{p}.
4.4. Pulsation g and rmodes
The theoretical pulsation models of the two late Be stars HD 181231 and HD 175869 that take the stellar deformation into account show that introducing a mixing of α_{ov} = 0.3–0.35 H_{p} is necessary to reproduce the observed pulsation frequencies with sectoral prograde gmodes. The value of overshooting α_{ov} we obtain is a signature of nonstandard mixing processes, which are discussed below.
A few higher frequencies are observed in HD 175869 and are not reproduced by the Tohoku models. These modes could be either gmodes that are poorly reproduced by our models because we did not treat rotation perfectly, or gmodes excited by another mechanism than the κ mechanism, e.g. stochastic gmodes.
In addition, very low frequencies ( ≤ 0.2 c d^{1}) are observed in both stars and especially in HD 181231. These frequencies could be attributed to m = 1 r modes of high radial order. As Townsend (2005) and Lee (2006) have shown, r modes are indeed excited by the κ mechanism at the Fe opacity bump in Btype stars (which also excites the gmodes). Therefore they appear in the models presented in Figs. 3 and 5.
Indeed, models show rmodes near the zero frequency (<0.15 c d^{1}) of m = 1 with ℓ′ = ℓ − 1 = 1 (see Figs. 3 and 5), where ℓ′ is the latitudinal degree for the toroidal component. Frequencies of rmodes in the observers frame can be written as (3)This equation indicates that an rmode of m = 1 with ℓ′ = 1 (i.e. ℓ = 2) has a low frequency of in the observer’s frame. These modes are odd modes (blue dashed lines in Figs. 3 and 5) because ℓ − m = 1 (the latitudinal degree for spheroidal component determines the symmetry of scalar variables such as temperature perturbations).
5. A laboratory for unravelling transport and mixing processes
As we pointed out in the introduction, four nonstandard transport mechanisms may be modifying mixing in the studied stars:

the hydrodynamical processes induced by rotation: i.e. the action of differential rotation and the associated shearinduced turbulence and largescale meridional circulation; this combines with the centrifugal flattening due to the centrifugal acceleration;

the action of a possible fossil magnetic field in the radiative envelope;

the action of convective flows of the core, which penetrate the surrounding radiative envelope;

the action of internal waves generated by this penetrative convection.
5.1. The rotational effects
To investigate the effects of (differential) rotation on the properties of the Be stars HD 181231 and HD 175869, and in particular on the size of their convective core, 5 M_{⊙} models were computed with the Geneva stellar evolution code described in Sect. 3.3. These models were computed with a solar chemical composition as given by Grevesse & Noels (1993) and a solar calibrated value for the mixinglength parameter (α_{⊙} = 1.59 with the input physics used for these computations).
Fig. 6 Evolutionary tracks in the HR diagram for 5 M_{⊙} models, calculated with the Geneva code. The continuous and dashed lines indicate models with an initial velocity of 350 km s^{1} and without rotation, respectively. Both models are computed without overshooting. The dotted line corresponds to a nonrotating models with an overshoot parameter of 0.1. 
Figure 6 shows the evolutionary tracks in the HR diagram for 5 M_{⊙} models with and without rotation. The rotating model had an initial velocity on the ZAMS of 350 km s^{1} and was computed without overshooting from the convective core into the surrounding radiatively stable layers. A corresponding nonrotating model was computed with exactly the same initial parameters except for the inclusion of shellular rotation. The effects of an overshooting from the convective core into the surroundings layers (see Sect. 5.3.1) on a distance d_{ov} ≡ α_{ov}min [H_{p},r_{core}] (Maeder & Meynet 1989) were also compared by computing the evolution of the same nonrotating model but with an overshoot parameter α_{ov} of 0.1.
As discussed in detail in Eggenberger et al. (2010) for 3 M_{⊙} models, including rotation changes the properties of the star in various ways. First, the hydrostatic effects of rotation result in a shift of the track towards lower luminosities and effective temperatures (compare the dashed and continuous lines in Fig. 6). This is because the effective gravity of the star is reduced by the inclusion of the centrifugal acceleration term. As the evolution proceeds, rotational mixing also begins to play an important role through two main effects. First, this brings fresh hydrogen fuel into the convective core, which slows down its decrease in mass during the evolution on the main sequence. Second, rotational mixing transports helium and other Hburning products into the radiative zone. Including shellular rotation consequently increases the size of the convective core for rotating models compared to nonrotating models. This can be clearly seen in Fig. 6 since this increase results in a shift of the turnoff point at the end of the mainsequence to lower effective temperatures.
Moreover, Fig. 6 shows that including overshooting leads to a similar shift of the turnoff point to the red part of the HR diagram, which is directly related to the increase of the mass of the convective core. Interestingly, this effect is very similar for the rotating model computed without overshooting and for the nonrotating model computed with an overshoot parameter of α_{ov} = 0.1 (solid and dotted lines respectively in Fig. 6). This indicates that the increase of the mass of the convective core due to rotational mixing is roughly reproduced by a nonrotating model with an overshoot parameter of α_{ov} = 0.1.
Fig. 7 Evolution of the enrichment in N for a 5 M_{⊙} model calculated with the Geneva code. The dashed blue line shows the evolution for a model without rotation, while the solid black line shows the rotating model with an initial velocity of 350 km s^{1}. 
For a 5 M_{⊙} model computed with an inital velocity of 350 km s^{1} we find indeed that rotational mixing typically leads to an increase of the size of the convective core of about 0.10 to 0.13 H_{p} along the main sequence. The exact value of this increase within this range depends on the evolutionary stage of the model because the effects of rotational mixing become more and more visible as the star is close to the end of the main sequence.
Figure 7 shows the enrichment in nitrogen for this rotating model along the main sequence. The enrichment increases as the star evolves, but the final enrichment remains low and probably undetectable in the data mentioned in Sect. 2.
5.2. The contribution of the magnetic field
In the past decade, studies of the magnetic fields in massive stars (e.g. by the MiMeS collaboration, see Wade et al. 2011) have shown that a magnetic field exists in a fraction of these stars. The fields are probably fossil ones in the external radiative envelope of massive stars. The magnetic field is therefore another major mechanism that should be studied in HD 181231 and HD 175869, because it could strongly modify the transport of angular momentum and mixing in these stars (Maeder & Meynet 2004; Mathis & Zahn 2005).
Neiner et al. (2009) and GutiérrezSoto et al. (2009) showed that HD 181231 and HD 175869 do not display magnetic signatures in their Stokes V spectra. Moreover, each longitudinal field value determined from their many (18 per star) Stokes V measurements is compatible with 0. The authors estimated upper limits of a possible nondetected longitudinal magnetic field to the largest error bars of the longitudinal field measurements: B_{l} ~ 650 and ~400 G for HD 181231 and HD 175869, respectively.
However, the relevant value for the present work is the upper limit of the polar magnetic field at the surface rather than the longitudinal field value, which depends on the magnetic configuration as seen by the observer. In addition, since the rotation period of the stars is not known with a very high precision, we do not know if the 18 Stokes V measurements have been obtained at a magnetic phase where the amplitude of the putative Zeeman signature is strong or flat. We therefore need to perform a statistical analysis.
Taking these facts into account, we performed a Monte Carlo simulation of the spectropolarimetic data to investigate the level of the magnetic field that could be hidden in the data noise.
Fig. 8 Chances that a magnetic field in HD 181231 and HD 175869 would have been detected at 3σ and 5σ detection levels in at least one of the 18 Narval measurements according to the strength of the dipolar magnetic field. Dashed horizontal lines indicate the 50% and 90% detection probability. 
First, we reanalysed the Narval spectropolarimetric data of Neiner et al. (2009) and GutiérrezSoto et al. (2009) using the leastsquares deconvolution (LSD) technique (Donati et al. 1997). We used an LSD line mask including 80 and 227 photospheric lines for HD 181231 and HD 175869, respectively. The LSD method assumes that the intrinsic broadening of each line is similar, therefore hydrogen lines are never used in LSD line masks. For each spectral line, the mask contains the wavelength, depth, and Landé factor, to be used by the LSD programme. Landé factors were extracted from Kurucz models. Line depths were adjusted to the observed spectrum. For each spectrum and each star, LSD Stokes I and V profiles were extracted.
For each star and for various values of the polar magnetic field B_{pol}, we then calculated 1000 oblique dipole models of each of the 18 Stokes V profiles with random obliquity angle, random rotational phase, and noise corresponding to the signaltonoise ratio of the Narval data.
To calculate these oblique dipole models, we used Gaussian local intensity profiles with a width calculated according to the resolving power of Narval and a typical macroturbulence. The depth of the intensity profile was determined by fitting the observed LSD I profiles. We then calculated local Stokes V profiles assuming the weakfield case and integrated over the visible hemisphere of the star. We obtained synthetic LSD Stokes V profiles, which we normalized to the intensity continuum. We used the mean Landé factor and wavelength obtained from the Narval observations. The model includes three parameters: the inclination angle i, the dipole obliquity angle β, and the polar field strength B_{pol}.
We then computed among the 1000 models of each 18 Stokes V profiles the chances of detection of the field in at least one of the 18 measurements. This translates into an upper limit for the possible nondetected polar field B_{pol} = 2300 G for HD 181231 and 700 G for HD 175869, for a 50% chance of detection at 3σ. Pushing the limit to 90% chance of detection at 5σ, a field with any oblique dipolar configuration above 5000 G for HD 181231 and 1500 G for HD 175869 would have been detected (see Fig. 8).
Local stronger magnetic fields could exist if the hypothetical magnetic field were more complex than dipolar. However, for massive (nonBp) stars, detected fields usually are dipolar (see e.g. Neiner et al. 2003), a predicted behaviour for fossil fields (Braithwaite & Spruit 2004; Duez & Mathis 2010). A few nondipolar examples are known, such as τ Sco (Donati et al. 2006), but they remain very rare.
The relative action between the magnetic field and (differential) rotation can be calculated by comparing the mean Alfvén pulsation (4)where is the average of the magnetic field amplitude and the mean density, and the characteristic surface angular velocity (5)where R and V_{surf} are the stellar radius and surface velocity, respectively (Moss 1992). Stellar parameters are similar for HD 181231 and HD 175869. The theory of dipolar fossil fields shows that the maximum internal field amplitude is only about 30 times stronger than the amplitude at the surface (see Fig. 8 in Braithwaite 2008; Braithwaite & Nordlund 2006). Therefore, using the conservative estimate G, where B_{pol} = 5000 G is the highest value determined above for HD 181231 and a 90% chance of detection, we find that , i.e. that Ω^{2} is more than six thousand times higher than . Consequently, rotational mechanisms would largely dominate magnetic effects, even if an undetected field were present in the studied stars. For the magnetic field to have an impact, the stars would need to have an averaged magnetic field amplitude G. Such a mean amplitude is clearly rejected by the Narval observations (even with 100% chance of detection at 5σ) and is unrealistic in a mainsequence Be star. We conclude that in Be stars in general, magnetic field impact on internal mixing will be dominated by rotational processes.
Until now, we assumed that the magnetic field in the radiative envelope is fossil. For completion, below we discuss the possibility for dynamo mechanisms to occur in massive stars. First, Spruit (2002) suggested that a dynamo may develop in radiation zones because of the instability of a predominently axisymmetric toroidal field generated by the shearing of the poloidal field by the differential rotation: this is the socalled TaylerSpruit dynamo. The only way to close the dynamo loop in this configuration is if the electromotive force from the coupling of the instability fluctuating velocity and magnetic field becomes sufficiently efficient to trigger the axisymmetric poloidal field against its Ohmic decay (Zahn et al. 2007). If this mechanism were to become active in radiative layers, it might modify the mixing (Maeder & Meynet 2004). However, Zahn et al. (2007) and Ruediger et al. (2011) showed how difficult it is to excite this dynamo in stellar interiors and we therefore exclude it as a probable scenario. Second, convective cores of intermediatemass and massive stars host a dynamo action that generates a magnetic field because of the turbulence and the differential rotation driven by convection (Brun et al. 2005; Augustson et al. 2011) that can also interact with the fossil field at the border of the convective core (Featherstone et al. 2009). The net effect of these processes on the mixing is to modify the convective flows because of the feedback of the Lorentz force on the dynamics and therefore on the penetrative convection, which is discussed hereafter. Finally, the potential smallscale field generated by a dynamo occuring in the subsurface convective layer (Cantiello & Braithwaite 2011) is not relevant for the largescale mixing studied here.
5.3. Penetrative convection and internal waves
Here we discuss the extramixing rate, which could be induced by penetrative convection at the border of the convective core and the associated related phenomenon such as the generation of internal waves in the surrounding stablystratified radiative envelope (we ignore here the action of a potential subsurface convection zone).
5.3.1. Penetrative convection
Penetrative convection, commonly called overshooting in stellar evolution, occurs in stars when convective elements penetrate into adjacent radiative zones. This corresponds here to the turbulent convective flows of the core that penetrate into the radiative envelope due to their inertia (see the pedagogical review by Dintrans 2009). One therefore has to obtain the value of the associated extension of the convective core, which is added to its standard radius determined using the classical Schwarzschild or Ledoux criterion. This cannot be derived using the usual mixinglength phenomenology but requires a physical prescription using theoretical considerations or numerical simulations.
From the theoretical point of view, Roxburgh (1978, 1989, 1992) provided constraints on the size of the convective core and the extent of overshooting. He showed that an estimate of the total size of the convective core can be deduced by looking at the following criterion (see also Zahn 1991) (6)which is deduced from the entropy equation. L_{nuc} and L_{rad} are the total nuclear and the radiative luminosities, respectively; Φ is the viscous dissipation per unit volume, is the volume, T is the temperature, r is the radius, and r_{c} the radius of the convective core including overshooting.
In the nuclear core some of the nuclear energy must be carried by convection, i.e. L_{nuc} > L_{rad}. Consequently, if the viscous dissipation is neglected (Φ ~ 0), there must be an overshooting region where L_{rad} > L_{nuc} whose extent is estimated from Eq. (6). This criterion has been applied for example by Doom (1985), who claimed that evolutionary models that take into account the Roxburgh criterion represent the observational constraints much better than standard models (see also Meynet et al. 1993, for the impact of overshooting on isochrones). However, convective flows are highly turbulent and the effects of viscous friction are enhanced. Indeed, if the microscopic viscosity of stellar interiors (ν_{m}) is low, it has to be replaced by a turbulent viscosity ν_{t} ~ l_{t}v_{t}, where l_{t} and v_{t} are the characteristic sizes of the largest turbulent eddies. Then ν_{t}/ν_{m} ≫ 1 and the adiabatic version of Eq. (6) (Φ ~ 0) gives only an upper limit for the penetration (Zahn 1993).
However, the necessary characterisation of turbulent flows can only be tested using numerical simulations in the domain of parameters accessible today. Moreover, we have to take into account rapid rotation, which is one of the key parameters for our two rapid rotating Be stars.
This is the reason why it is interesting to discuss the most recent global highly nonlinear numerical simulations of the dynamics of the rotating convective cores in intermediatemass and massive stars (see e.g. Browning et al. 2004; Meakin & Arnett 2007; Augustson et al. 2011). If we focus on the detailed work by Browning et al. (2004), many interesting features can be isolated concerning convective cores, which are illustrated in their work devoted to the case of an Astar. First, for a given value of the entropy gradient at the convection/radiation border, the degree of turbulence (i.e. the Reynolds number) influences the penetration length, which decreases as the flow is more and more turbulent because of a smaller filling factor for the convective upflows (the plumes) (see e.g. the 3D nonlinear simulations of penetrative compressible convection in localised planar domains, Brummell et al. 2002). Second, when the rotation rate is increased (for given stratification and diffusivities), the penetration becomes more intense because of the stabilisation of the flow by the Coriolis acceleration. This behaviour has also been predicted in planetary studies of, e.g., Takehiro & Lister (2001). Finally, Browning et al. (2004) showed that the 0.2 H_{p} extension needed here to complete the rotational transport contribution to match our seismic modelling is attainable.
Because of rapid rotation, penetrative convection therefore constitutes the natural candidate to explain the additional H_{p} extension that cannot be explained by the combined action of centrifugal acceleration and secular rotational transport computed in Sect. 5.1. However, to be able to conclude, one has to examine the last transport mechanism that may operate in the studied stars, namely the action of internal waves.
5.3.2. Internal waves
Internal waves are generated by penetrative convective flows at the top of the convective core as observed in the previously cited simulations (Browning et al. 2004; Meakin & Arnett 2007). These internal waves can induce mixing in two ways:

First, internal waves propagate in the radiative envelope and deposit their associated angular momentum at the location where they are damped by radiative damping. There, the differential rotation is modified because of Reynolds stresses, and mixing is modified in an indirect way because of the shear instability (Zahn et al. 1997; Talon & Charbonnel 2005).

Second, internal waves act directly on the chemical distribution as a diffusion mechanism, as proposed by Press (1981) and Schatzman (1996), or because they propagate in their nonlinear regime close to their excitation region, which makes them break (Mocák et al. 2011).

First, internal waves are trapped in an equatorial belt when ω_{c} < 2Ω because of the Coriolis acceleration action (Lee & Saio 1997; Mathis 2009; Ballot et al. 2010). The convective kinetic energy is then transmitted less efficiently to the internal waves in the radiative envelope. The associated transport is consequently decreased (see Pantillon et al. 2007; Mathis et al. 2008).

Next, because of the Coriolis acceleration, internal waves are damped through thermal diffusion closer to the excitation region than in the slowly rotating case (i.e. close to the convective core) and just below the surface (see again Pantillon et al. 2007; Mathis et al. 2008).
6. Discussion
Our seismic models have shown that mixing of 0.3–0.35 H_{p} is necessary to reproduce the pulsation gmodes observed in the CoRoT data of the two late Be stars HD 181231 and HD 175869. We also reproduced rmodes at low frequencies.
Our calculation of stellar evolution models including rotation showed that 0.1 to 0.13 H_{p} of the nonstandard convective core extension is due to the combined action of the centrifugal flattening and the rotational mixing mechanisms that lead to the increase of its size. Our present knowledge of penetrative convection and transport processes induced by internal waves suggests that the other 0.2 H_{p} nonstandard core extension may be caused by the convective overshoot with a possible weak action of excited internal waves. Linearly combining the rotational mixing and overshoot contributions is allowed because of their different characteristic timescales (secular and dynamical, respectively). Finally, we showed that a potential fossil magnetic field does not significantly contribute to mixing in the studied stars and in Be stars in general.
The drawbacks of our present analysis are the following:

First, the current knowledge on internal waves generation and related transport processes in rapidly rotating stars is quite poor. Indeed, no secular evolution models, such as those applied to lowmass stars by Talon & Charbonnel (2005), are available in the literature for massive stars. Moreover, the study of the impact of the Coriolis and centrifugal accelerations on transport is only in its infancy, as is the prediction of the stochastic excitation of gravitoinertial waves.
Second, we combined results obtained with several tools that have been computed in 1, 2, and 3D. To obtain a coherent picture in the near future, a more rigourous approach would be to combine completely 2D stellar evolution codes that include transport processes (Rieutord et al. 2005; Espinosa Lara & Rieutord 2007) with corresponding highprecision 2D stellar oscillation codes (Reese et al. 2009). However, this constitutes an extremely complex task, and the approach presented here, in which stateoftheart models are used, constitutes a good first exploration of highly coupled dynamical processes induced by rapid rotation. Moreover, 3D simulations of penetrative convection that have been discussed in Sect. 5.3.1 do not take into account the stellar flattening by the centrifugal acceleration. However, this effect becomes important only for external layers and can be neglected for the convective core.
Finally, we assumed a singlestar evolution scenario. If the studied Be stars had followed a binary evolution, their structure and their mixing would have been modified by the accretion and tidal interactions (the dynamical tide, see e.g. Zahn 1975, 1977). These effects should then be coupled to previously studied transport processes, which are also efficient in binaries (Zahn 1994; Talon & Kumar 1998; de Mink et al. 2009). This constitutes one of the main challenge for the modelling of nonstandard mechanisms in stellar evolution.
7. Conclusions
Using various stateoftheart codes to model CoRoT observations of pulsations in two rapidlyrotating late Be stars, we have shown that these stars host prograde sectoral gmodes as well as possible rmodes.
With these models and theoretical calculations, we have also shown that strong mixing is present in these stars. Part of the mixing is due to penetrative convection at the top of the convective core with a possible contribution of excited internal waves. The rest of the mixing is due to rotational processes, in particular the meridional circulation and shearinduced turbulence due to the radiative envelope differential rotation and the centrifugal flattening.
Finally, we have shown that if a magnetic field is present in rapidly rotating stars, such as Be stars, its impact on mixing will be dominated by the impact of rotation.
Acknowledgments
The CoRoT space mission, launched on December 27, 2006, has been developed and is operated by CNES, with the contribution of Austria, Belgium, Brazil, ESA, Germany and Spain. We thank J.P. Zahn for valuable discussions and the referee for his/her suggestions. C.N. thanks the Agence National pour la recherche (ANR) programme SIROCO, the Programme National de Physique Stellaire (PNPS) and CNES for their support. C.N. is grateful to E. Alecian for providing her Monte Carlo code.
References
 Augustson, K. C., Brun, A. S., & Toomre, J. 2011, in IAU Symp., 271, 361 [Google Scholar]
 Badnell, N. R., Bautista, M. A., Butler, K., et al. 2005, MNRAS, 360, 458 [NASA ADS] [CrossRef] [Google Scholar]
 Baglin, A., Auvergne, M., Barge, P., et al. 2006, in ESA SP 1306, ed. M. Fridlund, A. Baglin, J. Lochard, & L. Conroy, 33 [Google Scholar]
 Baglin, A., Auvergne, M., Barge, P., et al. 2009, in IAU Symp., 253, 71 [Google Scholar]
 Ballot, J., Lignières, F., Reese, D. R., & Rieutord, M. 2010, A&A, 518, A30 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Braithwaite, J. 2008, MNRAS, 386, 1947 [NASA ADS] [CrossRef] [Google Scholar]
 Braithwaite, J., & Nordlund, Å. 2006, A&A, 450, 1077 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Braithwaite, J., & Spruit, H. C. 2004, Nature, 431, 819 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Browning, M. K., Brun, A. S., & Toomre, J. 2004, ApJ, 601, 512 [NASA ADS] [CrossRef] [Google Scholar]
 Brummell, N. H., Clune, T. L., & Toomre, J. 2002, ApJ, 570, 825 [Google Scholar]
 Brun, A. S., Browning, M. K., & Toomre, J. 2005, ApJ, 629, 461 [NASA ADS] [CrossRef] [Google Scholar]
 Cantiello, M., & Braithwaite, J. 2011, A&A, 534, A140 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Chaboyer, B., & Zahn, J.P. 1992, A&A, 253, 173 [NASA ADS] [Google Scholar]
 Chandrasekhar, S. 1933a, MNRAS, 93, 390 [NASA ADS] [CrossRef] [Google Scholar]
 Chandrasekhar, S. 1933b, MNRAS, 93, 462 [NASA ADS] [Google Scholar]
 Clement, M. J. 1998, ApJS, 116, 57 [NASA ADS] [CrossRef] [Google Scholar]
 de Mink, S. E., Cantiello, M., Langer, N., et al. 2009, A&A, 497, 243 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 de Mink, S. E., Langer, N., & Izzard, R. G. 2011, in IAU Symp. 272, ed. C. Neiner, G. Wade, G. Meynet, & G. Peters, 531 [Google Scholar]
 Decressin, T., Mathis, S., Palacios, A., et al. 2009, A&A, 495, 271 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Deupree, R. G. 1990, ApJ, 357, 175 [NASA ADS] [CrossRef] [Google Scholar]
 Deupree, R. G. 1995, ApJ, 439, 357 [NASA ADS] [CrossRef] [Google Scholar]
 Dewi, J. D. M. 2007, in Massive Stars in Interactive Binaries, ed. N. St.Louis, & A. F. J. Moffat, ASP Conf. Ser., 367, 315 [Google Scholar]
 Dintrans, B. 2009, Commun. Asteroseis., 158, 45 [NASA ADS] [Google Scholar]
 Donati, J.F., Semel, M., Carter, B. D., Rees, D. E., & Collier Cameron, A. 1997, MNRAS, 291, 658 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Donati, J.F., Howarth, I. D., Jardine, M. M., et al. 2006, MNRAS, 370, 629 [NASA ADS] [CrossRef] [Google Scholar]
 Doom, C. 1985, A&A, 142, 143 [NASA ADS] [Google Scholar]
 Duez, V., & Mathis, S. 2010, A&A, 517, A58 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Eggenberger, P., Meynet, G., Maeder, A., et al. 2008, Ap&SS, 316, 43 [NASA ADS] [CrossRef] [Google Scholar]
 Eggenberger, P., Miglio, A., Montalban, J., et al. 2010, A&A, 509, A72 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ekström, S., Meynet, G., Maeder, A., & Barblan, F. 2008, A&A, 478, 467 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Espinosa Lara, F., & Rieutord, M. 2007, A&A, 470, 1013 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Featherstone, N. A., Browning, M. K., Brun, A. S., & Toomre, J. 2009, ApJ, 705, 1000 [NASA ADS] [CrossRef] [Google Scholar]
 Grevesse, N., & Noels, A. 1993, in Origin and evolution of the elements: proceedings of a symposium in honour of H. Reeves, held in Paris, June 22 − 25, 1992, ed. N. Prantzos, E. VangioniFlam, & M. Casse (Cambridge, England: Cambridge University Press), 14 [Google Scholar]
 GutiérrezSoto, J., Floquet, M., Samadi, R., et al. 2009, A&A, 506, 133 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Iglesias, C. A., & Rogers, F. J. 1996, ApJ, 464, 943 [NASA ADS] [CrossRef] [Google Scholar]
 Kippenhahn, R., & Thomas, H. 1970, in Stellar Rotation, IAU Colloq., 4, 20 [Google Scholar]
 Langer, N., Cantiello, M., Yoon, S.C., et al. 2008, in IAU Symp. 250, ed. F. Bresolin, P. A. Crowther, & J. Puls, 167 [Google Scholar]
 Lee, U. 2006, MNRAS, 365, 677 [NASA ADS] [CrossRef] [Google Scholar]
 Lee, U. 2008, Commun. Asteroseis., 157, 203 [Google Scholar]
 Lee, U., & Baraffe, I. 1995, A&A, 301, 419 [NASA ADS] [Google Scholar]
 Lee, U., & Saio, H. 1997, ApJ, 491, 839 [NASA ADS] [CrossRef] [Google Scholar]
 Lovekin, C., Neiner, C., Saio, H., Mathis, S., & GutiérrezSoto, J. 2010, Astron. Nachr., 331, 1061 [NASA ADS] [CrossRef] [Google Scholar]
 Maeder, A., & Meynet, G. 1989, A&A, 210, 155 [NASA ADS] [Google Scholar]
 Maeder, A., & Meynet, G. 2004, A&A, 422, 225 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Maeder, A., & Zahn, J.P. 1998, A&A, 334, 1000 [NASA ADS] [Google Scholar]
 Martayan, C., Frémat, Y., Hubert, A.M., et al. 2007, A&A, 462, 683 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mathis, S. 2009, A&A, 506, 811 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mathis, S. 2010, Astron. Nachr., 331, 883 [NASA ADS] [CrossRef] [Google Scholar]
 Mathis, S., & Zahn, J.P. 2004, A&A, 425, 229 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mathis, S., & Zahn, J.P. 2005, A&A, 440, 653 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mathis, S., Talon, S., Pantillon, F.P., & Zahn, J.P. 2008, Sol. Phys., 251, 101 [NASA ADS] [CrossRef] [Google Scholar]
 Meakin, C. A., & Arnett, D. 2007, ApJ, 667, 448 [Google Scholar]
 Meynet, G., & Maeder, A. 1997, A&A, 321, 465 [NASA ADS] [Google Scholar]
 Meynet, G., & Maeder, A. 2000, A&A, 361, 101 [NASA ADS] [Google Scholar]
 Meynet, G., Mermilliod, J.C., & Maeder, A. 1993, A&AS, 98, 477 [NASA ADS] [Google Scholar]
 Meynet, G., Georgy, C., Revaz, Y., et al. 2010, in Rev. Mex. Astron. Astrofis. Conf. Ser., 38, 113 [Google Scholar]
 Mocák, M., Meakin, C. A., Müller, E., & Siess, L. 2011, ApJ, 743, 55 [NASA ADS] [CrossRef] [Google Scholar]
 Moss, D. 1992, MNRAS, 257, 593 [NASA ADS] [Google Scholar]
 Neiner, C., Geers, V. C., Henrichs, H. F., et al. 2003, A&A, 406, 1019 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Neiner, C., GutiérrezSoto, J., Baudin, F., et al. 2009, A&A, 506, 143 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pantillon, F. P., Talon, S., & Charbonnel, C. 2007, A&A, 474, 155 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pols, O. R., Cote, J., Waters, L. B. F. M., & Heise, J. 1991, A&A, 241, 419 [NASA ADS] [Google Scholar]
 Press, W. H. 1981, ApJ, 245, 286 [NASA ADS] [CrossRef] [Google Scholar]
 Reese, D. R. 2010, Astron. Nachr., 331, 1038 [NASA ADS] [CrossRef] [Google Scholar]
 Reese, D. R., MacGregor, K. B., Jackson, S., Skumanich, A., & Metcalfe, T. S. 2009, A&A, 506, 189 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Rieutord, M. 2006, A&A, 451, 1025 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Rieutord, M., Corbard, T., Pichon, B., Dintrans, B., & Lignières, F. 2005, in SF2A2005: Semaine de l’Astrophysique Francaise, ed. F. Casoli, T. Contini, J. M. Hameury, & L. Pagani, 759 [Google Scholar]
 Roxburgh, I. W. 1978, A&A, 65, 281 [NASA ADS] [Google Scholar]
 Roxburgh, I. W. 1989, A&A, 211, 361 [NASA ADS] [Google Scholar]
 Roxburgh, I. W. 1992, A&A, 266, 291 [NASA ADS] [Google Scholar]
 Ruediger, G., Kitchatinov, L. L., & Elstner, D. 2011, MNRAS, submitted [arXiv:1107.2548] [Google Scholar]
 Schatzman, E. 1996, J. Fluid Mech., 322, 355 [NASA ADS] [CrossRef] [Google Scholar]
 Spruit, H. C. 2002, A&A, 381, 923 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Takehiro, S.I., & Lister, J. R. 2001, Earth Planet. Sci. Lett., 187, 357 [NASA ADS] [CrossRef] [Google Scholar]
 Talon, S. 2008, in EAS Pub. Ser. 32, ed. C. Charbonnel, & J.P. Zahn, 81 [Google Scholar]
 Talon, S., & Charbonnel, C. 2005, A&A, 440, 981 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Talon, S., & Kumar, P. 1998, ApJ, 503, 387 [NASA ADS] [CrossRef] [Google Scholar]
 Townsend, R. H. D. 2005, MNRAS, 364, 573 [NASA ADS] [CrossRef] [Google Scholar]
 Wade, G. A., Alecian, E., Bohlender, D. A., et al. 2011, in IAU Symp. 272, ed. C. Neiner, G. Wade, G. Meynet, & G. Peters, 118 [Google Scholar]
 Walker, G. A. H., Kuschnig, R., Matthews, J. M., et al. 2005, ApJ, 635, L77 [NASA ADS] [CrossRef] [Google Scholar]
 Zahn, J.P. 1975, A&A, 41, 329 [NASA ADS] [Google Scholar]
 Zahn, J.P. 1977, A&A, 57, 383 [NASA ADS] [Google Scholar]
 Zahn, J.P. 1991, A&A, 252, 179 [NASA ADS] [Google Scholar]
 Zahn, J.P. 1992, A&A, 265, 115 [NASA ADS] [Google Scholar]
 Zahn, J.P. 1993, Space Sci. Rev., 66, 285 [NASA ADS] [CrossRef] [Google Scholar]
 Zahn, J.P. 1994, A&A, 288, 829 [NASA ADS] [Google Scholar]
 Zahn, J.P., Talon, S., & Matias, J. 1997, A&A, 322, 320 [Google Scholar]
 Zahn, J.P., Brun, A. S., & Mathis, S. 2007, A&A, 474, 145 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
All Tables
Effective temperature (in log), gravity, projected rotational velocity, inclination angle, equatorial radius, mass, luminosity, and rotational frequency of the two late Be stars HD 181231 and HD 175869.
All Figures
Fig. 1 Estimated positions of HD 181231 (log T_{eff} = 4.155) and HD 175869 (log T_{eff} = 4.079) on the log T_{eff} − log g plane (black dots) with evolutionary tracks of M/M_{⊙} = 3,4,5,6,7 with a standard chemical composition of (X,Z) = (0.07, 0.02). Dotted lines are the evolutionary tracks of M/M_{⊙} = 4,4.5,5,5.5,6 with a convectivecore overshooting of 0.35H_{p}. There is no rotation included in these tracks. Open squares and circles indicate the positions of models for HD 181231 and HD 175869, respectively (see below). 

In the text 
Fig. 2 Frequencies of HD 181231 observed with CoRoT (bottom panel) and growth rate of the modelled frequencies of low degree (ℓ ≤ 2) (first, third and fourth panels from the top) excited in a 5 M_{⊙} star without mixing, calculated with the Tohoku code. The second panel from the top shows the azimuthal order m of the pulsation mode corresponding to the modelled frequencies in the third panel. Colours indicate the parity of the mode: even modes in solid red and odd modes in dashed blue lines. Stellar parameters used in the three models are indicated at the top of each panel. 

In the text 
Fig. 3 Frequencies of HD 181231 observed with CoRoT (bottom panel) and growth rate and azimuthal order of the modelled frequencies of low degree (ℓ ≤ 2) (top panels) excited in a 5.5 and 6 M_{⊙} star with a core overshooting of 0.35 H_{p}, calculated with the Tohoku code. Colours indicate the parity of the mode: even modes in solid red and odd modes in dashed blue lines. Stellar parameters used in the two models are indicated at the top of each model. 

In the text 
Fig. 4 Frequency match to HD 181231 calculated using NRO models. Solid lines indicate observed frequencies, arbitrarily attributed to ℓ = 2. Dashed lines indicate the model frequencies and their value of ℓ. In this case, the frequencies are best matched by modes with ℓ = 1 to 5, and only one mode is axisymmetric. 

In the text 
Fig. 5 Same as Fig. 3 but for the star HD 175869. 

In the text 
Fig. 6 Evolutionary tracks in the HR diagram for 5 M_{⊙} models, calculated with the Geneva code. The continuous and dashed lines indicate models with an initial velocity of 350 km s^{1} and without rotation, respectively. Both models are computed without overshooting. The dotted line corresponds to a nonrotating models with an overshoot parameter of 0.1. 

In the text 
Fig. 7 Evolution of the enrichment in N for a 5 M_{⊙} model calculated with the Geneva code. The dashed blue line shows the evolution for a model without rotation, while the solid black line shows the rotating model with an initial velocity of 350 km s^{1}. 

In the text 
Fig. 8 Chances that a magnetic field in HD 181231 and HD 175869 would have been detected at 3σ and 5σ detection levels in at least one of the 18 Narval measurements according to the strength of the dipolar magnetic field. Dashed horizontal lines indicate the 50% and 90% detection probability. 

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.