Angular momentum redistribution by mixed modes in evolved low-mass stars. II. Spin-down of the core of red giants induced by mixed modes

The detection of mixed modes in subgiants and red giants by the CoRoT and emph{Kepler} space-borne missions allows us to investigate the internal structure of evolved low-mass stars. In particular, the measurement of the mean core rotation rate as a function of the evolution places stringent constraints on the physical mechanisms responsible for the angular momentum redistribution in stars. It showed that the current stellar evolution codes including the modelling of rotation fail to reproduce the observations. An additional physical process that efficiently extracts angular momentum from the core is thus necessary. Our aim is to assess the ability of mixed modes to do this. To this end, we developed a formalism that provides a modelling of the wave fluxes in both the mean angular momentum and the mean energy equations in a companion paper. In this article, mode amplitudes are modelled based on recent asteroseismic observations, and a quantitative estimate of the angular momentum transfer is obtained. This is performed for a benchmark model of 1.3 $M_{odot}$ at three evolutionary stages, representative of the evolved pulsating stars observed by CoRoT and Kepler. We show that mixed modes extract angular momentum from the innermost regions of subgiants and red giants. However, this transport of angular momentum from the core is unlikely to counterbalance the effect of the core contraction in subgiants and early red giants. In contrast, for more evolved red giants, mixed modes are found efficient enough to balance and exceed the effect of the core contraction, in particular in the hydrogen-burning shell. Our results thus indicate that mixed modes are a promising candidate to explain the observed spin-down of the core of evolved red giants, but that an other mechanism is to be invoked for subgiants and early red giants.

Angular momentum redistribution by mixed modes in evolved low-mass stars. I. Theoretical formalism

Seismic observations by the space-borne mission emph{Kepler} have shown that the core of red giant stars slows down while evolving, requiring an efficient physical mechanism to extract angular momentum from the inner layers. Current stellar evolution codes fail to reproduce the observed rotation rates by several orders of magnitude, and predict a drastic spin-up of red giant cores instead. New efficient mechanisms of angular momentum transport are thus required. In this framework, our aim is to investigate the possibility that mixed modes extract angular momentum from the inner radiative regions of evolved low-mass stars. To this end, we consider the Transformed Eulerian Mean (TEM) formalism, introduced by Andrews & McIntyre (1978), that allows us to consider the combined effect of both the wave momentum flux in the mean angular momentum equation and the wave heat flux in the mean entropy equation as well as their interplay with the meridional circulation. In radiative layers of evolved low-mass stars, the quasi-adiabatic approximation, the limit of slow rotation, and the asymptotic regime can be applied for mixed modes and enable us to establish a prescription for the wave fluxes in the mean equations. The formalism is finally applied to a $1.3 M_odot$ benchmark model, representative of observed CoRoT and emph{Kepler} oscillating evolved stars. We show that the influence of the wave heat flux on the mean angular momentum is not negligible and that the overall effect of mixed modes is to extract angular momentum from the innermost region of the star. A quantitative and accurate estimate requires realistic values of mode amplitudes. This is provided in a companion paper.

Non-perturbative effect of rotation on dipolar mixed modes in red giant stars

The space missions CoRoT and Kepler provide high quality data that allow us to test the transport of angular momentum in stars by the seismic determination of the internal rotation profile. Our aim is to test the validity of the seismic diagnostics for red giants rotation that are based on a perturbative method and to investigate the oscillation spectra when the validity does not hold. We use a non-perturbative approach implemented in the ACOR code (Ouazzani et al. 2012) that accounts for the effect of rotation on pulsations, and solves the pulsation eigenproblem directly for dipolar oscillation modes. We find that the limit of the perturbation to first order can be expressed in terms of the rotational splitting compared to the frequency separation between consecutive dipolar modes. Above this limit, non-perturbative computations are necessary but only one term in the spectral expansion of modes is sufficient as long as the core rotation rate remains significantly smaller than the pulsation frequencies. Each family of modes with different azimuthal symmetry, m, has to be considered separately. In particular, in case of rapid core rotation, the density of the spectrum differs significantly from one m-family of modes to another, so that the differences between the period spacings associated with each m-family can constitute a promising guideline toward a proper seismic diagnostic for rotation.

Seismic diagnostics for transport of angular momentum in stars 1. Rotational splittings from the PMS to the RGB

Rotational splittings are currently measured for several main sequence stars and a large number of red giants with the space mission Kepler. This will provide stringent constraints on rotation profiles. Our aim is to obtain seismic constraints on the internal transport and surface loss of angular momentum of oscillating solar-like stars. To this end, we study the evolution of rotational splittings from the pre-main sequence to the red-giant branch for stochastically excited oscillation modes. We modified the evolutionary code CESAM2K to take rotationally induced transport in radiative zones into account. Linear rotational splittings were computed for a sequence of $1.3 M_{odot}$ models. Rotation profiles were derived from our evolutionary models and eigenfunctions from linear adiabatic oscillation calculations. We find that transport by meridional circulation and shear turbulence yields far too high a core rotation rate for red-giant models compared with recent seismic observations. We discuss several uncertainties in the physical description of stars that could have an impact on the rotation profiles. For instance, we find that the Goldreich-Schubert-Fricke instability does not extract enough angular momentum from the core to account for the discrepancy. In contrast, an increase of the horizontal turbulent viscosity by 2 orders of magnitude is able to significantly decrease the central rotation rate on the red-giant branch. Our results indicate that it is possible that the prescription for the horizontal turbulent viscosity largely underestimates its actual value or else a mechanism not included in current stellar models of low mass stars is needed to slow down the rotation in the radiative core of red-giant stars.

The underlying physical meaning of the $ u_{rm max}- u_{rm c}$ relation

Asteroseismology of stars that exhibit solar-like oscillations are enjoying a growing interest with the wealth of observational results obtained with the CoRoT and Kepler missions. In this framework, scaling laws between asteroseismic quantities and stellar parameters are becoming essential tools to study a rich variety of stars. However, the physical underlying mechanisms of those scaling laws are still poorly known. Our objective is to provide a theoretical basis for the scaling between the frequency of the maximum in the power spectrum ($ u_{rm max}$) of solar-like oscillations and the cut-off frequency ($ u_{rm c}$). Using the SoHO GOLF observations together with theoretical considerations, we first confirm that the maximum of the height in oscillation power spectrum is determined by the so-called emph{plateau} of the damping rates. The physical origin of the plateau can be traced to the destabilizing effect of the Lagrangian perturbation of entropy in the upper-most layers which becomes important when the modal period and the local thermal relaxation time-scale are comparable. Based on this analysis, we then find a linear relation between $ u_{rm max}$ and $ u_{rm c}$, with a coefficient that depends on the ratio of the Mach number of the exciting turbulence to the third power to the mixing-length parameter.

Modelling the excitation of acoustic modes in Alpha Cen A

We infer from different seismic observations the energy supplied per unit of time by turbulent convection to the acoustic modes of Alpha Cen A (HD 128620), a star which is similar but not identical to the Sun. The inferred rates of energy supplied to the modes (i.e. mode excitation rates) are found to be significantly larger than in the Sun. They are compared with those computed with an excitation model that includes two sources of driving, the Reynolds stress contribution and the advection of entropy fluctuations. The model also uses a closure model, the Closure Model with Plumes (CMP hereafter), that takes the asymmetry between the up- and down-flows (i.e. the granules and plumes, respectively) into account. Different prescriptions for the eddy-time correlation function are also confronted to observational data. Calculations based on a Gaussian eddy-time correlation underestimate excitation rates compared with the values derived from observations for Alpha Cen A. On the other hand, calculations based on a Lorentzian eddy-time correlation lie within the observational error bars. This confirms results obtained in the solar case. With respect to the helioseismic data, those obtained for Alpha Cen A constitute an additional support for our model of excitation. We show that mode masses must be computed taking turbulent pressure into account. Finally, we emphasize the need for more accurate seismic measurements in order to discriminate, in the case of Alpha Cen A, between the CMP closure model and the quasi-Normal Approximation as well as to confirm or not the need to include the excitation by the entropy fluctuations.

FILOU oscillation code

The present paper provides a description of the oscillation code FILOU, its main features, type of applications it can be used for, and some representative solutions. The code is actively involved in CoRoT/ESTA exercises (this volume) for the preparation for the proper interpretation of space data from the CoRoT mission. Although CoRoT/ESTA exercises have been limited to the oscillations computations for non-rotating models, the main characteristic of FILOU is, however, the computation of radial and non-radial oscillation frequencies in presence of rotation. In particular, FILOU calculates (in a perturbative approach) adiabatic oscillation frequencies corrected for the effects of rotation (up to the second order in the rotation rate) including near degeneracy effects. Furthermore, FILOU works with either a uniform rotation or a radial differential rotation profile (shellular rotation), feature which makes the code singular in the field.

Stochastic excitation of non-radial modes I. High-angular-degree p modes

Turbulent motions in stellar convection zones generate acoustic energy, part of which is then supplied to normal modes of the star. Their amplitudes result from a balance between the efficiencies of excitation and damping processes in the convection zones. We develop a formalism that provides the excitation rates of non-radial global modes excited by turbulent convection. As a first application, we estimate the impact of non-radial effects on excitation rates and amplitudes of high-angular-degree modes which are observed on the Sun. A model of stochastic excitation by turbulent convection has been developed to compute the excitation rates, and it has been successfully applied to solar radial modes (Samadi & Goupil 2001, Belkacem et al. 2006b). We generalize this approach to the case of non-radial global modes. This enables us to estimate the energy supplied to high-($ell$) acoustic modes. Qualitative arguments as well as numerical calculations are used to illustrate the results. We find that non-radial effects for $p$ modes are non-negligible: - for high-$n$ modes (i.e. typically $n > 3$) and for high values of $ell$; the power supplied to the oscillations depends on the mode inertia. - for low-$n$ modes, independent of the value of $ell$, the excitation is dominated by the non-diagonal components of the Reynolds stress term. We carried out a numerical investigation of high-$ell$ $p$ modes and we find that the validity of the present formalism is limited to $ell < 500$ due to the spatial separation of scale assumption. Thus, a model for very high-$ell$ $p$-mode excitation rates calls for further theoretical developments, however the formalism is valid for solar $g$ modes, which will be investigated in a paper in preparation.

Inferring information about rotation from stellar oscillations

The first part of this paper aims at illustrating the intense scientific activity in the field of stellar rotation although, for sake of shortness, we cannot be exhaustive nor give any details. The second part is devoted to the rotation as a pertubation effect upon oscillation frequencies. The discussion focuses on one specific example: the p-modes frequency small separation which provides information about properties of the stellar inner layers. It is shown that the small separation can be affected by rotation at the level of 0.1-0.2 microHz for a 1.4 Mo model rotating with an equatorial velocity of 20 km/s at the surface. This is of the same order of magnitude as the expected precision on frequencies with a 3 months observation and must therefore be taken into account. We show however that it is possible to recover the small separation free of these contaminating effects of rotation, provided enough high quality data are available as will be with space seismic missions such as Eddington.

Seismic diagnostics on stellar convection treatment from oscillation amplitudes of p-modes

The excitation rate P of solar p-modes is computed with a model of stochatic excitation which involves constraints on the averaged properties of the solar turbulence. These constraints are obtained from a 3D simulation. Resulting values for P are found ~9 times larger than when the calculation assumes properties of turbulent convection which are derived from an 1D solar model based on Gough (1977)s formulation of the mixing-length theory (GMLT). This difference is mainly due to the assumed values for the mean anisotropy of the velocity field in each case.Calculations based on 3D constraints bring the P maximum closer to the observational one.We also compute P for several models of intermediate mass stars (1Mo < M < 2Mo). Differences in the values of P_max between models computed with the classical mixing-length theory and GMLT models are found large enough for main sequence stars to suggest that measurements of P in this mass range will be able to discriminate between different models of turbulent convection.