Search

Context. The interactions between magnetic fields and differential rotation in stellar radiative interiors could play a major role in achieving an understanding of the magnetism of intermediate-mass and massive stars and of the differential rotation profile observed in red-giant stars. Aims. The present study is aimed at studying the flow and field produced by a stellar radiative zone which is initially made to rotate differentially in the presence of a large-scale poloidal magnetic field threading the whole domain. We focus both on the axisymmetric configurations produced by the initial winding-up of the magnetic field lines and on the possible instabilities of those configurations. We investigate in detail the effects of the stable stratification and thermal diffusion and we aim, in particular, to assess the role of the stratification at stabilising the system. Methods. We performed 2D and 3D global Boussinesq numerical simulations started from an initial radial or cylindrical differential rotation and a large-scale poloidal magnetic field. Under the conditions of a large rotation frequency compared to the Alfvén frequency, we built a magnetic configuration strongly dominated by its toroidal component. We then perturbed this configuration to observe the development of non-axisymmetric instabilities. Results. The parameters of the simulations were chosen to respect the ordering of time scales of a typical stellar radiative zone. In this framework, the axisymmetric evolution is studied by varying the relative effects of the thermal diffusion, the Brunt-Väisälä frequency, the rotation, and the initial poloidal field strength. After a transient time and using a suitable adimensionalisation, we find that the axisymmetric state only depends on tes/tAp the ratio between the Eddington–Sweet circulation time scale and the Alfvén time scale. A scale analysis of the Boussinesq magnetohydrodynamical equations allows us to recover this result. In the cylindrical case, a magneto-rotational instability develops when the thermal diffusivity is sufficiently high to enable the favored wavenumbers to be insensitive to the effects of the stable stratification. In the radial case, the magneto-rotational instability is driven by the latitudinal shear created by the back-reaction of the Lorentz force on the flow. Increasing the level of stratification then leaves the growth rate of the instability mainly unaffected while its horizontal length scale grows. Conclusions. Non-axisymmetric instabilities are likely to exist in stellar radiative zones despite the stable stratification. They could be at the origin of the magnetic dichotomy observed in intermediate-mass and massive stars. They are also unavoidable candidates for the transport of angular momentum in red giant stars. [ABSTRACT FROM AUTHOR]

The advent of space photometry with CoRoT and Kepler has allowed for the gathering of exquisite and extensive time series for a wealth of main-sequence stars, including γ Doradus stars, whose detailed seismology was not achievable from the ground. γ Doradus stars present an incredibly rich pulsation spectra, with gravito-inertial modes, in some cases supplemented with δ Scuti-like pressure modes – for the hybrid stars – and, in many cases, with Rossby modes. The present paper aims to show that in addition to these modes which have been established in the radiative envelope, pure inertial modes that are trapped in the convective core can be detected in Kepler observations of γ Doradus stars thanks to their resonance with the gravito-inertial modes. We started by using a simplified model of perturbations in a full sphere of uniform density. Under these conditions, the spectrum of pure inertial modes is known from analytical solutions of the so-called Poincaré equation. We then computed coupling factors, which helped select the pure inertial modes which interact best with the surrounding dipolar gravito-inertial modes. Using complete calculations of gravito-inertial modes in realistic models of γ Doradus stars, we are able to show that the pure inertial and gravito-inertial resonances appear as "dips" in the gravito-inertial mode period spacing series at spin parameters that are close to those predicted by the simple model. We find the first evidence of such dips in the Keplerγ Doradus star KIC 5608334. Finally, using complete calculations in isolated convective cores, we find that the spin parameters of the pure inertial and gravito-inertial resonances are also sensitive to the density stratification of the convective core. In conclusion, we have discovered that certain dips in gravito-inertial mode period spacings that have been observed in some Kepler stars are, in fact, signatures of resonances with pure-inertial modes that are trapped in the convective core. This holds the promise that it would be possible to finally access the central conditions, namely, the rotation and density stratification, of intermediate-mass stars in the main sequence. [ABSTRACT FROM AUTHOR]

Context: Mode identification has remained a major obstacle in the interpretation of pulsation spectra in rapidly rotating stars. This has motivated recent work on calculating realistic multi-colour mode visibilities in this type of star. Aims: We would like to test mode identification methods and seismic diagnostics in rapidly rotating stars, using oscillation spectra that are based on these new theoretical predictions. Methods: We investigate the auto-correlation function and Fourier transform of theoretically calculated frequency spectra, in which modes are selected according to their visibilities. Given that intrinsic mode amplitudes are determined by non-linear saturation and cannot currently be theoretically predicted, we experimented with various ad-hoc prescriptions for setting the mode amplitudes, including using random values. Furthermore, we analyse the ratios between mode amplitudes observed in different photometric bands to see up to what extent they can identify modes. Results: When non-random intrinsic mode amplitudes are used, our results show that it is possible to extract a mean value for the large frequency separation or half its value and, sometimes, twice the rotation rate, from the auto-correlation of the frequency spectra. Furthermore, the Fourier transforms are mostly sensitive to the large frequency separation or half its value. The combination of the two methods may therefore measure and distinguish the two types of separations. When the intrinsic mode amplitudes include random factors, which seems more representative of real stars, the results are far less favourable. It is only when the large separation or half its value coincides with twice the rotation rate, that it might be possible to detect the signature of a frequency regularity. We also find that amplitude ratios are a good way of grouping together modes with similar characteristics. By analysing the frequencies of these groups, it is possible to constrain mode identification, as well as determine the large frequency separation and the rotation rate. [ABSTRACT FROM AUTHOR]

Context. Recent spectropolarimetric surveys of main-sequence intermediate-mass stars have exhibited a dichotomy in the distribution of the observed magnetic field between the kG dipoles of Ap/Bp stars and the sub-Gauss magnetism of Vega and Sirius. Aims. We would like to test whether this dichotomy is linked to the stability versus instability of large-scale magnetic configurations in differentially rotating radiative zones. Methods. We computed the axisymmetric magnetic field obtained from the evolution of a dipolar field threading a differentially rotating shell. A full parameter study including various density profiles and initial and boundary conditions was performed with a 2D numerical code. We then focused on the ratio between the toroidal and poloidal components of the magnetic field and discuss the stability of the configurations dominated by the toroidal component using local stability criteria and insights from recent 3D numerical simulations. Results. The numerical results and a simple model show that the ratio between the toroidal and the poloidal magnetic fields is highest after an Alfvén crossing time of the initial poloidal field. For high density contrasts, this ratio converges towards an asymptotic value that can thus be extrapolated to realistic stellar cases. We then consider the stability of the magnetic configurations to nonaxisymmetric perturbations and find that configurations dominated by the toroidal component are likely to be unstable if the shear strength is significantly higher than the poloidal Alfvén frequency. An expression for the critical poloidal field below which magnetic fields are likely to be unstable is found and is compared to the lower bound of Ap/Bp magnetic fields. [ABSTRACT FROM AUTHOR]

Context. Turbulent transport of chemical elements in radiative zones of stars is considered in current stellar evolution codes thanks to phenomenologically derived diffusion coefficients. Recent local numerical simulations suggest that the coefficient for radial turbulent diffusion due to radial differential rotation satisfies Dt ≃ 0.058κ/Ri, in qualitative agreement with the model of Zahn (1992, A&A, 265, 115). However, this model does not apply (i) when differential rotation is strong with respect to stable thermal stratification or (ii) when chemical stratification has a significant dynamical effect, a situation encountered at the outer boundary of nuclear-burning convective cores. Aims. We extend our numerical study to consider the effects of chemical stratification and of strong shear, and compare the results with prescriptions used in stellar evolution codes. Methods. We performed local, direct numerical simulations of stably stratified, homogeneous, sheared turbulence in the Boussinesq approximation. The regime of high thermal diffusivities, typical of stellar radiative zones, is reached thanks to the so-called small-Péclet-number approximation, which is an asymptotic development of the Boussinesq equations in this regime. The dependence of the diffusion coefficient on chemical stratification was explored in this approximation. Results. Maeder’s extension of Zahn’s model in the strong-shear regime (Maeder 1995, A&A, 299, 84) is not supported by our results, which are better described by a model found in the geophysical literature. As regards the effect of chemical stratification, our quantitative estimate of the diffusion coefficient as a function of the mean gradient of mean molecular weight leads to the formula Dt ≃ 0.45κ(0.12−Riμ) /Ri, which is compatible in the weak-shear regime with the model of Maeder & Meynet (1996, A&A, 313, 140) but not with Maeder’s (1997, A&A, 321, 134). [ABSTRACT FROM AUTHOR]

Context. In stellar interiors, rotation is able to drive turbulent motions, and the related transport processes have a significant influence on the evolution of stars. Turbulent mixing in the radiative zones is currently taken into account in stellar evolution models through a set of diffusion coefficients that are generally poorly constrained. Aims. We want to constrain the form of one of them, the radial diffusion coefficient of chemical elements due to the turbulence driven by radial differential rotation, derived by Zahn (1974, IAU Symp., 59, 185 and 1992, A&A, 265, 115) on phenomenological grounds and largely used since. Methods. We performed local, direct numerical simulations of stably stratified homogeneous sheared turbulence using the Boussinesq approximation. The domain of low Péclet numbers found in stellar interiors is currently inaccessible to numerical simulations without approximation. It is explored here thanks to a suitable asymptotic form of the Boussinesq equations. The turbulent transport of a passive scalar is determined in statistical steady states. Results. We provide a first quantitative determination of the turbulent diffusion coefficient and find that the form proposed by Zahn is in good agreement with the results of the numerical simulations. [ABSTRACT FROM AUTHOR]