Open Access
Volume 649, May 2021
Article Number A62
Number of page(s) 7
Section Stellar structure and evolution
Published online 11 May 2021

© V. Prat and S. Mathis 2021

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1. Introduction

The transport of angular momentum and chemical species plays a key role in stellar evolution. In particular, rotation induces transport both through large-scale motions such as meridional circulation and through turbulent, small-scale motions generated by hydrodynamical instabilities (Zahn et al. 1983; Zahn 1992; Maeder & Zahn 1998; Mathis & Zahn 2004; Maeder 2009). These multi-dimensional small-scale transport processes cannot be resolved in one-dimensional (1D) stellar evolution codes, so their secular effects have to be modelled in those codes.

In the last decades, helio- and asteroseismology have provided us with many constraints on the internal rotation of the Sun and distant stars. In the case of the Sun, helioseismic data bring to light an almost flat rotation profile in the solar radiative zone down to 0.2 R, and a possibly faster core (Brown et al. 1989; Thompson et al. 2003; García et al. 2007; Fossat et al. 2017). Benomar et al. (2015) detected only weak differential rotation on a few solar-type stars. For a large number of subgiant and red giant stars, asteroseismic data show a weak core-surface rotation contrast (Beck et al. 2012; Mosser et al. 2012; Deheuvels et al. 2012, 2014, 2015; Triana et al. 2017; Gehan et al. 2018). Weak differential rotation has also been detected in intermediate-mass and massive stars (Kurtz et al. 2014; Saio et al. 2015; Triana et al. 2015; Murphy et al. 2016; Van Reeth et al. 2016; Aerts et al. 2017; Van Reeth et al. 2018; Ouazzani et al. 2019). Finally, a strong transport of angular momentum is needed to explain data obtained for white dwarves (Suijs et al. 2008; Hermes et al. 2017) and neutron stars (Heger et al. 2005; Hirschi & Maeder 2010).

Zahn (1992) introduced a formalism for turbulent transport accounting for both advection by the meridional circulation and the anisotropic turbulent transport induced by differential rotation, which is modelled thanks to vertical and horizontal turbulent diffusion coefficients. The anisotropy of the transport is due to the combined effect of rotation and stable stratification, which limit horizontal and vertical motions through the buoyancy force and the Coriolis acceleration, respectively. Those coefficients are based on phenomenological arguments and thus introduce large uncertainties in stellar evolution models. A lot of work has been done to improve the physics of the vertical diffusion coefficient by adding physical ingredients such as gradients of mean molecular weight (Maeder & Meynet 1996; Talon & Zahn 1997), and local numerical simulations have been performed to test existing models (Prat & Lignières 2013, 2014; Garaud & Kulenthirarajah 2016) and investigate the effect of viscosity (Prat et al. 2016; Garaud et al. 2017; Gagnier & Garaud 2018). Current models predict less transport than what is needed to explain observations of internal rotation (Eggenberger et al. 2012; Ceillier et al. 2013; Marques et al. 2013; Cantiello et al. 2014; Ouazzani et al. 2019).

Concerning the horizontal diffusion coefficient, since the inital model of Zahn (1992) that was based on phenomenological arguments, other models based on energetic arguments (Maeder 2003) or results of an unstratified Taylor-Couette experiment (Mathis et al. 2004) have been proposed, but this coefficient is still poorly constrained, and its impact on stellar evolution is clearly overlooked. The first numerical simulations of the turbulent transport triggered by the instabilities of a horizontal shear in stellar conditions have also been performed (Cope et al. 2020; Garaud 2020). Finally, Park et al. (2020, 2021) studied the linear instabilities of the horizontal differential rotation in stellar radiation zones as a function of stratification, rotation, and thermal diffusivity.

However, the widely used formalism of Zahn (1992) assumes that the turbulent transport of angular momentum is viscous, which is not necessarily the case. Indeed, anisotropic turbulent motions can generate non-zero turbulent fluxes of angular momentum even in the absence of differential rotation, which is a characteristic of non-diffusive transport. The existence of such non-diffusive turbulent transport associated with the Coriolis force, known as the Λ effect (Rüdiger 1989), has been confirmed by direct numerical simulations (see e.g. Käpylä 2019). Originally, this effect was introduced in the context of solar and stellar convective zones to explain the persistance of differential rotation despite the presence of a magnetic field, which was supposed to erase it. Kitchatinov & Brandenburg (2012, hereafter referred to as KB12) later investigated the effect of the stable stratification of stellar radiative zones on the Reynolds stress.

Based on this work, Mathis et al. (2018, hereafter referred to as M+18) proposed a new model of the horizontal transport induced by the vertical shear instability. In particular, M+18 added a vertical shear (due to radial differential rotation) to the formalism of KB12 and show that it has no significant effect on the anisotropy of turbulence. Besides, they used estimates of mean velocity correlations to propose new prescriptions for the horizontal turbulent diffusion coefficient of chemical elements. The implementation of these prescriptions in stellar evolution computations shows a slightly enhanced transport of angular momentum throughout the main sequence, but not enough to fit helio- and asteroseismic observations. Furthermore, a major caveat of this implementation is that it assumed that the transport of angular momentum had the same (diffusive) nature as the transport of chemical elements, and thus completely ignored the Λ effect.

In the present work, we perform a new generalisation of the formalism of KB12 in the presence of a general (both vertical and horizontal) shear. It allows us to propose new prescriptions for the transport of both angular momentum and chemical elements that include the Λ effect.

In Sect. 2, we derive spectral properties of the flow anisotropised by stratification, rotation, and shear. Then, in Sect. 3, we deduce scalings for velocity correlations as a function of stratification, rotation, and shear. We propose possible interpretations of the new scalings in Sect. 4. We explain how non-viscous turbulent transport of angular momentum should be implemented in stellar evolution codes to account for the Λ effect in Sect. 5. Finally, we conclude on the impact of these results on the modelling of turbulence in stellar evolution codes in Sect. 6.

2. Turbulent spectrum

In the Boussinesq approximation, which is justified here because turbulent motions are at length scales much smaller than the pressure scale height and velocity scales much smaller than the sound speed, the equations governing the flow are




where V is the velocity, P′ and ρ′ are pressure and density fluctuations, respectively, ρ is the background density, g is the gravity vector, ν is the kinematic viscosity, f is a forcing term, s is the specific entropy, and κ is the thermal diffusivity. If one neglects the large-scale meridional circulation, velocity can be split into a mean rotation part and a fluctuating part:


where r, θ, and φ are the standard spherical coordinates, (er, eθ, eφ) is the associated orthonormal basis, Ω is the rotation rate, and u is the velocity fluctuation. The continuity Eq. (1) implies that


The left-hand side of the momentum Eq. (2) can be developed into


where s is the distance to the rotation axis and es the associated unit vector. The second term, which is the centrifugal acceleration, will be neglected in the following. Then, in the frame that locally rotates at the same rate as the fluid, Eq. (2) becomes


where S = r sin θΩ is the shear rate.

Similarly, the entropy can be split into a mean part ⟨s⟩ and a fluctuating one s′:


Equation (3) yields


Non-linear and dissipative terms are usually difficult to model analytically in turbulent flows. To simplify the problem, we use here a relaxation approximation, also known as τ approximation, which assumes that in a stationary steady state, the main effect of these terms is to tend the flow to relax with a given time constant τ. Equations (7) and (9) can then be approximated by



For a complete treatment of thermal diffusion, we refer the reader to Park et al. (2020, 2021).

These two equations can be combined into a single one by expressing entropy fluctuations as a function of density fluctuations:


where cP is the specific heat capacity at constant pressure. This yields


where the Brunt-Väisälä frequency N is defined by


Using the Fourier decomposition


the continuity Eq. (5) becomes


and pressure fluctuations can be eliminated from Eq. (13):


where , , , , , and . This way of writing the forcing term means that we assume that it comes from a pre-existing background turbulent flow. The previous equation can be written as a matrix relation


where is the matrix


and .

When is inversible, Eq. (18) implies


Using the fact that also verifies the continuity equation, the previous equation reduces to




is the determinant of , , 𝒟 is the matrix given by


and δij and εijk are the Kronecker and Levi-Civita tensors, respectively.

The spectrum of the background turbulence is assumed to be characterised by the correlation tensor , where ⟨⟩ denotes a statistical average. In the continuity of Mathis et al. (2018), we assume a purely horizontal background turbulence, for which the Fourier transform of is


where E(k) is the kinetic energy spectrum (see e.g. Davidson 2004) such that


The turbulent spectrum of the sheared turbulence is deduced from that of the background turbulence using the relation


After some algebra, this leads to






3. Mean correlations

Mean velocity correlations can be computed using the relation


Stellar radiative zones are usually strongly stratified. It implies that (for example, the typical value of the ratio N/(2Ω) is around 500 for the Sun; see also André et al. 2018, 2019, for the evolution of this ratio along stellar evolution for low- and intermediate-mass stars), , and . As a first approximation,


where μ = cos α, and (with these notations, σ = cos θ cos α − sin θ sin α cos β). For the same reason,


Using parity arguments and integrating over k (see Appendix A), this leads, when , to


This term corresponds to the Reynolds stress in the meridional plane. Similarly,




The first two terms correspond to the vertical and horizontal transport of angular momentum, respectively. The other two describe the horizontal turbulent velocity scales.

We note that when tends towards zero, the scalings for ⟨uruθ⟩ and ⟨uθuφ⟩ (in Eqs. (34) and (36), respectively) vanish, while the others tend towards finite values. These values, as well as the zero limit of ⟨uθuφ⟩, are consistent with the results of KB12. In contrast, KB12 found a non-zero value for ⟨uruθ⟩:


which is of lower order in than the new scaling of Eq. (34).

The case of is more complicated, since it contains three terms that come from


and a priori have the same order of magnitude. The first term leads to


which is consistent with the scaling derived in KB12 and M+18. The second term vanishes for parity reasons. The third term leads to


When , this reduces to


which is small compared to . In total,


4. Interpretation

Interestingly, none of the terms computed in Sect. 3 explicitely depend on the radial shear , although it was included in the governing equations. In contrast, all share a common feature: they have a denominator in . This implies that when gets close to 2, the turbulence is enhanced. Physically, this might be related to the limit of the Rayleigh-Taylor instability. The Rayleigh-Taylor stability criterion is


Assuming that Ω > 0, this criterion can be rewritten


When 2Ω + sin θSr > 0, the last inequality is verified for all positive values of cos θSθ, but only for negative values verifying


By comparison with , this suggests that


This expression is similar to τ = (2Ω + Sr)−1, proposed by M+18. However, the presence of cos θ at the numerator means that τ would vanish at the equator, in which case the scalings derived in Sect. 3 are not valid there. Besides, the comparison with the Rayleigh-Taylor instability stands only in the special case where 2Ω + sin θSr > 0 and cos Sθ < 0. Thus, the criterion could be linked to some other hydrodynamical instability instead, such as the inertial or the inflectional instability studied by Park et al. (2020, 2021).

The form of the expression of ⟨uθuφ⟩ in Eq. (36), which describes the horizontal transport of angular momentum, suggests that this flux term is viscous, with a turbulent viscosity coefficient equal to


We note that this coefficient has no explicit dependence on the Brunt-Väisälä frequency or on the vertical shear. However, τ most likely depends on Sr, and not necessarily on N. Thus, the viscosity coefficient probably depends on both Sr and Sθ. The dependence on ⟨u2(0) is problematic for the use in stellar evolution codes, since there is no easy way to reliably prescribe this quantity. One way to roughly estimate it is to assume that horizontal motions generated by the background turbulence are global. This can be written


where R is the radial extent of the radiative zone. Introducing this expression in Eq. (48) yields


In the case where the dynamics is dominated by the horizontal shear, assuming τ ∼ 1/Sθ further leads to a horizontal viscosity coefficient that is proportional to the horizontal shear rate, which is consistent with the viscosity coefficient derived by Park et al. (2021). We note that, in contrast, Garaud (2001) found that the horizontal shear instability can lead to an anti-diffusive horizontal transport of angular momentum.

Although the expression in Eq. (48) is not directly usable in stellar evolution codes, it can be used to express other flux terms as a function of νh. Thus,



In the present work, we consider the case of a horizontally isotropic turbulence (i.e. with identical statistical properties in the latitudinal and azimuthal directions). This is rigourously justified for the radiative zone of slowly rotating stars, in which turbulent structures would have the shape of a horizontal pancake (see e.g. Fig. 1 of M+18). Because of that, there is no horizontal Λ effect. In rapidly rotating stars, turbulent structures are elongated along the rotation axis, so the horizontal isotropy is broken (see e.g. Rüdiger 1980) and a non-zero horizontal flux of angular momentum due to the Λ effect is present. However, in stellar radiative zones, provided that the Brunt-Väisälä frequency is still much larger than the rotation rate, the horizontal flux due to the Λ effect remains negligle compared to the leading diffusive term, which reduces to Eq. (48) in the horizontally symmetric case.

In contrast with the horizontal transport, the vertical transport may not be viscous, since the term ⟨uruφ⟩ given in Eq. (35) is not explicitly proportional to Sr. Instead, it is proportional to the rotation rate Ω, which is a characteristic of the Λ effect. However, when the dynamics is dominated by the vertical shear, one can assume that τ ∼ 1/Sr, in which case Eq. (35) becomes


This then leads to a vertical turbulent viscosity coefficient


which is negative. This corresponds to an anti-diffusive transport of angular momentum. Cope et al. (2020) and Garaud (2020) estimated the vertical turbulent diffusion coefficient associated with the horizontal shear instability using scaling laws but they did not provide any proof that the transport of chemical elements is indeed diffusive. In addition, they did not discuss the transport of angular momentum. Using the approximation of Eq. (49) leads to


The same approximation can be used in the non-viscous case to express the flux from Eq. (35):


Following KB12 and M+18, we estimate the turbulent diffusion coefficients using a mixing-length approximation and , where . This allows us to estimate the ratio between vertical and horizontal diffusion coefficients:


which has the same scaling with N as found by M+18. Besides, in the limit where , Eq. (57) reduces to


which also has the same scaling with τ as found by M+18.

Moreover, in contrast with M+18, where ⟨uθuφ⟩ was zero and no horizontal turbulent viscosity coefficient could be directly derived, we can now compare horizontal viscosity and diffusion coefficients. We thus find that Dh ≃ 4νh. This suggests that the common assumption that they are equal may be qualitatively valid, but given the uncertainties introduced by the approximation we used, we cannot really conclude on the quantitative validity. In the vertical direction, we cannot directly compute the turbulent viscosity coefficient, so there is no evidence of such an equality.

5. Non-viscous transport in stellar evolution codes

If a prescription of ⟨u2(0) is found, our expressions for the flux terms in Eqs. (35) and (36) can be used in stellar evolution calculations. For non-viscous transport, the evolution of angular momentum is governed by


where u(r, θ) is the large-scale velocity field of the meridional circulation in stellar radiative zones. 1D stellar evolution codes that take rotation into account usually are based on the shellular approximation, which assumes that , where


and (see e.g. Zahn 1992). In this approximation, the meridional circulation reads


where P2 is the Legendre polynomial of degree 2. Averaging Eq. (59) over isobars yields




We note that the term of horizontal transport is no longer present since it has a zero average. We showed in Sect. 4 that the horizontal transport can be seen as viscous. Therefore, the equation for the fluctuations is the same as in Mathis et al. (2004), for example:


where and


6. Conclusion

In the current paper, we investigated the impact of horizontal shear on the anisotropy of the transport in stably stratified, rotating stellar radiative zones. We were able to derive new scalings for mean velocity correlations that confirmed that the main effect is to significantly enhance turbulent motions. Thus, horizontal shear can be considered as an additional source of transport of angular momentum and chemical elements both in the horizontal and vertical directions.

In the regime where the stratification has a much stronger effect than rotation, we found that the horizontal transport of angular momentum is diffusive, while the vertical transport is dominated by the Λ effect. Additionally, the horizontal turbulent viscosity and diffusion coefficients are proportional to each other, which partly justifies the approximation that they are equal, well-spread in stellar evolution codes. The expressions of all the mean Reynolds stresses computed in the current paper depend on a property of the background turbulent flow, which is unknown a priori. Nevertheless, if one of the components is known by another way (for example using a phenomenological model for the horizontal turbulent viscosity coefficient), this allows one to compute fluxes of angular momentum. We finally show how the non-viscous transport of angular momentum predicted by our model should be implemented in stellar evolution codes to go beyond the current diffusive and viscous formalism.

Although the relaxation approximation has been validated in special cases by numerical simulations (see e.g. KB12), its validity in the configuration considered in the present work is uncertain. In particular, the dependence of the non-linear time τ on the physical parameters is undetermined. In M+18, three different physically motivated but phenomenological expressions were proposed. Here, we derived from more rigourous arguments a new expression for the non-linear turbulent time scale which is close to one of those proposed in M+18, τ = (2Ω + Sr)−1, which corresponds to a turbulence dominated by the effects of rotation and vertical shear. Ultimately, direct numerical simulations should be performed to test the formalism and the new prescriptions in the presence of a general shear.

The results presented in this work have been obtained in the Boussinesq approximation, which neglects density fluctuations except in the buoyancy term. This treatment is valid as long as turbulent structures are small compared to the density scale height. When this condition is not satisfied, for example in regions with strong density gradients, similar calculations should be performed in the more general framework of the anelastic approximation, which accounts for density gradients, but filters pressure waves out.


V.P. and S.M. acknowledge support from the European Research Council through ERC grant SPIRE 647383 and from the CNES PLATO/GOLF grant at CEA-Saclay. The authors thank the anonymous referee for their constructive comments, which have helped us improve the manuscript.


  1. Aerts, C., Van Reeth, T., & Tkachenko, A. 2017, ApJ, 847, L7 [NASA ADS] [CrossRef] [Google Scholar]
  2. André, Q., Mathis, S., & Amard, L. 2018, SF2A-2018: Proceedings of the Annual meeting of the French Society of Astronomy and Astrophysics, Di [Google Scholar]
  3. André, Q., Mathis, S., & Amard, L. 2019, [Google Scholar]
  4. Beck, P. G., Montalban, J., Kallinger, T., et al. 2012, Nature, 481, 55 [NASA ADS] [CrossRef] [Google Scholar]
  5. Benomar, O., Takata, M., Shibahashi, H., Ceillier, T., & García, R. A. 2015, MNRAS, 452, 2654 [NASA ADS] [CrossRef] [Google Scholar]
  6. Brown, T. M., Christensen-Dalsgaard, J., Dziembowski, W. A., et al. 1989, ApJ, 343, 526 [NASA ADS] [CrossRef] [Google Scholar]
  7. Cantiello, M., Mankovich, C., Bildsten, L., Christensen-Dalsgaard, J., & Paxton, B. 2014, ApJ, 788, 93 [Google Scholar]
  8. Ceillier, T., Eggenberger, P., García, R. A., & Mathis, S. 2013, A&A, 555, A54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Cope, L., Garaud, P., & Caulfield, C. P. 2020, J. Fluid Mech., 903, A1 [Google Scholar]
  10. Davidson, P. A. 2004, Turbulence: an Introduction for Scientists and Engineers (Oxford University Press) [Google Scholar]
  11. Deheuvels, S., García, R. A., Chaplin, W. J., et al. 2012, ApJ, 756, 19 [NASA ADS] [CrossRef] [Google Scholar]
  12. Deheuvels, S., Doğan, G., Goupil, M. J., et al. 2014, A&A, 564, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Deheuvels, S., Ballot, J., Beck, P. G., et al. 2015, A&A, 580, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Eggenberger, P., Montalbán, J., & Miglio, A. 2012, A&A, 544, L4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Fossat, E., Boumier, P., Corbard, T., et al. 2017, A&A, 604, A40 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Gagnier, D., & Garaud, P. 2018, ApJ, 862, 36 [NASA ADS] [CrossRef] [Google Scholar]
  17. Garaud, P. 2001, MNRAS, 324, 68 [Google Scholar]
  18. Garaud, P. 2020, ApJ, 901, 146 [Google Scholar]
  19. Garaud, P., & Kulenthirarajah, L. 2016, ApJ, 821, 49 [Google Scholar]
  20. Garaud, P., Gagnier, D., & Verhoeven, J. 2017, ApJ, 837, 133 [NASA ADS] [CrossRef] [Google Scholar]
  21. García, R. A., Turck-Chièze, S., Jiménez-Reyes, S. J., et al. 2007, Science, 316, 1591 [Google Scholar]
  22. Gehan, C., Mosser, B., Michel, E., Samadi, R., & Kallinger, T. 2018, A&A, 616, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Heger, A., Woosley, S. E., & Spruit, H. C. 2005, ApJ, 626, 350 [NASA ADS] [CrossRef] [Google Scholar]
  24. Hermes, J. J., Gänsicke, B. T., Kawaler, S. D., et al. 2017, ApJS, 232, 23 [Google Scholar]
  25. Hirschi, R., & Maeder, A. 2010, A&A, 519, A16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Käpylä, P. J. 2019, A&A, 622, A195 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Kitchatinov, L. L., & Brandenburg, A. 2012, Astron. Nachr., 333, 230 [NASA ADS] [CrossRef] [Google Scholar]
  28. Kurtz, D. W., Saio, H., Takata, M., et al. 2014, MNRAS, 444, 102 [Google Scholar]
  29. Maeder, A. 2003, A&A, 399, 263 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Maeder, A. 2009, Physics, Formation and Evolution of Rotating Stars (Berlin Heidelberg: Springer) [Google Scholar]
  31. Maeder, A., & Meynet, G. 1996, A&A, 313, 140 [NASA ADS] [Google Scholar]
  32. Maeder, A., & Zahn, J.-P. 1998, A&A, 334, 1000 [NASA ADS] [Google Scholar]
  33. Marques, J. P., Goupil, M. J., Lebreton, Y., et al. 2013, A&A, 549, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Mathis, S., & Zahn, J. P. 2004, A&A, 425, 229 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Mathis, S., Palacios, A., & Zahn, J.-P. 2004, A&A, 425, 243 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Mathis, S., Prat, V., Amard, L., et al. 2018, A&A, 620, A22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Mosser, B., Goupil, M. J., Belkacem, K., et al. 2012, A&A, 548, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Murphy, S. J., Fossati, L., Bedding, T. R., et al. 2016, MNRAS, 459, 1201 [Google Scholar]
  39. Ouazzani, R. M., Marques, J. P., Goupil, M. J., et al. 2019, A&A, 626, A121 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Park, J., Prat, V., & Mathis, S. 2020, A&A, 635, A133 [EDP Sciences] [Google Scholar]
  41. Park, J., Prat, V., Mathis, S., & Bugnet, L. 2021, A&A, 646, A64 [EDP Sciences] [Google Scholar]
  42. Prat, V., & Lignières, F. 2013, A&A, 551, L3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Prat, V., & Lignières, F. 2014, A&A, 566, A110 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. Prat, V., Guilet, J., Viallet, M., & Müller, E. 2016, A&A, 592, A59 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. Rüdiger, G. 1980, Geophys. Astro. Fluid, 16, 239 [Google Scholar]
  46. Rüdiger, G. 1989, Differential Rotation and Stellar Convection (Gordon & Breach Science Publishers) [Google Scholar]
  47. Saio, H., Kurtz, D. W., Takata, M., et al. 2015, MNRAS, 447, 3264 [Google Scholar]
  48. Suijs, M. P. L., Langer, N., Poelarends, A.-J., et al. 2008, A&A, 481, L87 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  49. Talon, S., & Zahn, J.-P. 1997, A&A, 317, 749 [Google Scholar]
  50. Thompson, M. J., Christensen-Dalsgaard, J., Miesch, M. S., & Toomre, J. 2003, ARA&A, 41, 599 [NASA ADS] [CrossRef] [Google Scholar]
  51. Triana, S. A., Moravveji, E., Pápics, P. I., et al. 2015, ApJ, 810, 16 [NASA ADS] [CrossRef] [Google Scholar]
  52. Triana, S. A., Corsaro, E., De Ridder, J., et al. 2017, A&A, 602, A62 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  53. Van Reeth, T., Tkachenko, A., & Aerts, C. 2016, A&A, 593, A120 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  54. Van Reeth, T., Mombarg, J. S. G., Mathis, S., et al. 2018, A&A, 618, A24 [NASA ADS] [EDP Sciences] [Google Scholar]
  55. Zahn, J. P. 1983, in Saas-Fee Advanced Course 13: Astrophysical Processes in Upper Main Sequence Stars, eds. A. N. Cox, S. Vauclair, & J. P. Zahn, 253 [Google Scholar]
  56. Zahn, J.-P. 1992, A&A, 265, 115 [NASA ADS] [Google Scholar]

Appendix A: Integrals

The derivation of Eq. (34) from Eq. (33) requires the computation of the double integral


The first and third terms of the integrand are antisymmetric with respect to α = π/2, so the corresponding contribution to the integral is zero. The remaining term is then proportional to


Using the periodicity of the integrand and the variables t = tanβ and , one obtains


These integrals are convergent only when . Using , the integral now reads


Similarly, the derivation of Eq. (35) requires the computation of


and the derivation of Eq. (37) requires the computation of


It follows that


which is required for the derivation of Eq. (40). Finally, the derivation of Eq. (41) requires the computation of


Using the partial fraction decomposition


one obtains


Current usage metrics show cumulative count of Article Views (full-text 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 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.