Free Access
Issue
A&A
Volume 622, February 2019
Article Number A187
Number of page(s) 10
Section Stellar structure and evolution
DOI https://doi.org/10.1051/0004-6361/201834330
Published online 18 February 2019

© ESO 2019

1. Introduction

The inclusion of rotation in stellar evolution codes is complicated, due to the limited understanding of how it affects the stellar structure (Maeder 2009). Some observational constraints are available from studying surface abundances and surface rotation values. These help with the calibration of the implementation of rotation in stellar evolution codes affecting both the mixing of chemical elements and the transport of angular momentum (Heger et al. 2000; Meynet et al. 2013). In this paper, we use a set of core and envelope rotation rates from core helium burning stars as an extra set of observational constraints to study the transport of angular momentum.

A big step forward in the determination of internal rotation rates has been provided by the Kepler spacecraft (see Borucki et al. 2010), as mixed modes were found in the spectra. Mixed modes are able to propagate through both convective and radiative zones, and hold information about both. The rotational frequency splittings of mixed models were first measured by Beck et al. (2012) for the red giant KIC 8366239. They were then able to rule out solid-body rotation for this star. Another conclusion of this paper was that the core rotates about ten times faster than the envelope. More papers on internal rotation rates followed, also mainly focussing on observations of low-mass, evolved stars, and supported these conclusions: Deheuvels et al. (2012, 2014, 2015, 2017) analyse in great detail small data sets to obtain both core and surface rotation rates, while Mosser et al. (2012) and Gehan et al. (2018) analyse the core rotation rates of a data set of several hundred giants (and core helium burning) stars.

Stellar evolution codes have been unable to match the observed ratio of core and envelope rotation rates. Predicted core rotation rates are about two orders of magnitude higher than observed rates (see Eggenberger et al. 2012; Marques et al. 2013). This is a major issue within the stellar evolution community as it shows that a process for transport of angular momentum is missing from the current implementation of rotation.

Broadly speaking, there are three types of processes which could improve the transport of angular momentum in stars: hydrodynamical, wave-driven, and magnetic. There is no consensus on which process dominates. Eggenberger et al. (2005), Suijs et al. (2008), Cantiello et al. (2014) show that magnetic fields, and the Tayler-Spruit (TS) dynamo in particular (Spruit 1999, 2002) are effective in coupling the core and envelope to increase the transport of angular momentum. While the TS dynamo is able to reproduce the solar rotation profile, it cannot provide the coupling needed to match the asteroseismic observations of red giants. Also, the inclusion of the TS dynamo has not been tested for the full mass range for which we have observational constraints on the internal stellar rotation rates.

The transport of angular momentum by low-frequency internal gravity waves has been studied with multi-D simulations (Fuller et al. 2014; Rogers 2015; Rogers & McElwaine 2017). They conclude that internal gravity waves are a promising method for transporting angular momentum. However its behaviour is complex and difficult to translate to a 1D parametrisation that can be included in 1D stellar evolution codes. Pinçon et al. (2017) showed that, for red giants, the gravity waves alone do not transport enough angular momentum to match observations. However, they could provide the needed transport during the subgiant phase. Belkacem et al. (2015) show that the transport by mixed models could play a role for evolved red giants.

For the magnetic fields and for the wave-driven option, there is still work to be done on the physical process itself. Therefore, it is important to derive the efficiency of the missing process of angular momentum so that its physical character can be revealed. To this end, Eggenberger et al. (2012, 2017) have performed stellar evolution calculations including both hydrodynamical processes like shear and meridional circulation, and a constant additional artificial viscosity (νadd), that only influences the angular momentum transport and not the chemical composition mixing. This constant has no physical meaning; it is added to investigate the level of efficiency the missing process of angular momentum transport should have, and is used to investigate whether we can determine the evolution of the missing process. Eggenberger et al. (2012) focus on the 1.5 M red giant KIC 8366239 observed by Beck et al. (2012) and the efficiency of the missing process needed to match those observations. The authors found a νadd of 3 × 104 cm2 s−1 as a mean value for the efficiency of the transport process, constrained strongly by the asteroseismic observations. Spada et al. (2016) followed a similar approach to constrain the missing process by including a diffusive process to the transport of angular momentum that varies with the ratio of core to envelope rotation rate, inspired by the azimuthal magneto-rotational instability (AMRI, see Rüdiger et al. 2007). They compare their stellar evolution models to observations of core rotation rate from Deheuvels et al. (2014) and Mosser et al. (2012) and conclude that the missing process of angular momentum transport has to change throughout the evolution of a star to be able to match the post-main sequence rotational evolution of low-mass stars. Eggenberger et al. (2017), following the same strategy as Eggenberger et al. (2012), focussed on KIC 7341231, a 0.84 M red giant, for which Deheuvels et al. (2012) deduced a core rotation rate and an upper limit for the surface rotation rate. For this red giant an additional transport process was again needed, and the authors determined the efficiency to be 1 × 103 <  νadd cm2 s−1 <  1.4 × 104. This value is lower than that found for the efficiency of the missing transport process for the more massive KIC 8366239, so Eggenberger et al. (2017) concluded that the missing process of angular momentum transport is sensitive to both evolutionary phase and the initial mass.

This paper focusses on further constraining the missing process of angular momentum transport by adding a νadd to stellar evolution calculations. The asteroseismic observations used are the seven core helium burning stars from Deheuvels et al. (2015). These stars have a mass around 2.5 M, which means they have not experienced helium flashes during their evolution, and therefore the evolution of their rotational properties post-main sequence are different from KIC 7341231 and KIC 8366239. For each of the seven stars, the surface and core rotation rates are both published. This allows us to put strong constraints on the efficiency of the missing process of angular momentum transport. The structure of the paper is as follows. In Sect. 2 we introduce our methodology, and in Sect. 3 we describe the evolution of our 2.5 M models up to the core helium burning phase. In Sect. 3 we investigate whether the TS dynamo is also unable in this case to provide the coupling needed to match observations. In Sect. 4 we determine the efficiency of the missing process of angular momentum transport needed to match the rotation rates obtained by Deheuvels et al. (2015). Section 5 follows with the comparison of our models to the observed white dwarf rotation rates. We end the paper with our conclusions in Sect. 6.

2. Physics of the models

The stellar evolution calculations presented in this paper were performed with the Module for Experiments in Stellar Astrophysics (MESA). The MESA code is described in the code papers (see Paxton et al. 2011, 2013) and will not be repeated here. Most initial parameters match the papers of the Nugrid collaboration (Pignatari et al. 2016; Battino et al. 2016). In summary, the Schwarzschild criterion is used for the convective boundary placement, while exponentially decaying convective boundary mixing as introduced by Herwig et al. (1997) is used for the boundaries. Mass loss on the red-giant branch (RGB) is set according to Reimers (1975), on the AGB by Blöcker (1995). OPAL Type 2 opacities are used (Rogers et al. 1996), and for the lower temperatures Ferguson et al. (2005).

2.1. Rotation

The implementation of rotation in MESA follows Heger et al. (2000) and we use the default settings as defined in that paper. The transport processes included in calculations in this paper are the Eddington-Sweet circulation, dynamical and secular shear instabilities, the Solberg-Høiland criterion and the TS dynamo. We excluded the Goldreich-Schubert-Fricke instability because both Hirschi & Maeder (2010) and Caleo et al. (2016) show that the GSF instability is not likely to contribute to the transport of angular momentum and might not be present at all in stars. Rotation is included at the zero age main sequence (ZAMS) as rigid body rotation. Although the implementation of rotation in MESA is different from that in GENEC, the code used by Eggenberger et al. (2012,2017), we confirm that we find similar νadd to explain the observations as in those papers (see Appendix B for a comparison between models of the two codes).

Several techniques exist in MESA to smooth the diffusion profiles of the instabilities, individually or their sum. We did not use any of these techniques, apart from a technique that smoothes the diffusion profile of the TS dynamo over time, exactly as was included by Cantiello et al. (2014). Without this smoothing technique, the stellar evolution calculation with and without the TS dynamo show rotation rates in the same order of magnitude.

In MESA, to calculate the transport of angular momentum, we use

( Ω t ) m = 1 j ( m ) t [ ( 4 π r 2 ρ ) 2 j D am ( Ω m ) ] 2 Ω r ( r t ) m ( 1 2 d ln j d ln r ) , $$ \begin{aligned} \left( \frac{\partial \Omega }{\partial t}\right)_m=\frac{1}{j}\left( \frac{\partial }{\partial m}\right)_t \left[ (4\pi r^2 \rho )^2jD_{\rm {am}} \left( \frac{\partial \Omega }{\partial m}\right) \right]-\frac{2\Omega }{r}\left(\frac{\partial r}{\partial t}\right)_m\left( \frac{1}{2}\frac{d\mathrm {ln}j}{d\mathrm {ln}r}\right) ,\end{aligned} $$(1)

where Ω is the angular velocity, j the specific angular momentum, and total diffusion coefficient Dam takes into account all processes that transport angular momentum, including the hydrodynamical ones listed previously. The νadd is also added to this variable Dam. This implementation is identical to the implementation of the νadd by Eggenberger et al. (2012,2017). At the inner and outer boundary ∂Ω/∂t is set to 0, following Heger et al. (2000).

The general evolution of the structure and angular momentum profiles of the models are presented in Appendix A. Uncertainties on the models are discussed in Appendix B.

2.2. The seven KIC stars

In Table 1 we summarise the important parameters of the seven KIC stars used as comparison sample, which are taken from Deheuvels et al. (2015). We include the core and surface rotation rates, the ratio between them, and the surface gravity (log g), all with their error margins. The metallicities of the seven stars are around solar, according to the APOGEE Data Release 14 (Abolfathi et al. 2018), which includes all seven stars. We used a metallicity of Z = 0.014 and the metal abundance mixture of Grevesse et al. (1993), and therefore focus on matching the global trends of the seven stars as a group instead of trying to find best-fit models for each star individually. This allows us to constrain the missing process of angular momentum for core helium burning stars.

Table 1.

Properties of the seven KIC stars from Deheuvels et al. (2015).

The initial mass of our models is chosen to be 2.5 M because this is very close to the mean observed mass of the seven KICs. In Appendix B we will see that the influence of the stellar evolution code used on the rotational properties is negligible.

Other observations of rotation rates in evolved stars in the same mass range have been published in Massarotti et al. (2007), Mosser et al. (2012), Tayar et al. (2015), and Ceillier et al. (2017) and analysed in Tayar & Pinsonneault (2018). These data sets, however, only include either the surface or the core rotation rates. To date, the data set of Deheuvels et al. (2015) is the only data set in the 2 to 3 M mass range that provides both rotation rates. This allows us to constrain our models better than when we only have one of the rates, so we only use the data set of Deheuvels et al. (2015) in this study.

It is important to note that the rotation rates labelled as “core” rotation rates are actually “near core” values, as shown in Fig. 5 of Deheuvels et al. (2015). The comparison of the calculations to the region where the observations of the core rotation originate from is explained in Appendix C.

3. Can the TS dynamo provide enough coupling to explain asteroseismic derived rotation properties of core helium burning stars?

The first goal of this paper is to investigate whether the TS dynamo provides enough coupling between core and envelope to match the observations of the core helium burning stars analysed by Deheuvels et al. (2015). Cantiello et al. (2014) show, for stars with an initial mass of 1.5 M, that during the early RGB inclusion of the TS dynamo provides more coupling between core and envelope but not enough to match the RGB rotation rates provided by Mosser et al. (2012). Thus, they concluded that the RGB phase is the evolutionary phase where more coupling is needed. However, the evolution of 1.5 M and 2.5 M are very different, in particular during the RGB phase. Stars with an initial mass below about 2 M undergo helium flashes in the core after it has become degenerate, and cores of stars with a higher initial mass ignite core helium burning before becoming degenerate. As a consequence, the times between the end of core hydrogen and the start of core helium burning are different; our calculations show a difference of one order of magnitude. For this reason, testing the conclusions of Cantiello et al. (2014) for 2.5 M stars is a valuable task, especially when comparing them with observations of stars that are already past the RGB phase.

Figure 1 shows the core (solid line) and envelope (dashed line) rotation rates of our 2.5 M models as a function of the surface gravity with different initial rotational velocities: 25, 50, and 150 km s−1. The start of the main sequence (MS) is where the core and envelope rotation rates are equal (top left) and the end of the core helium burning phase is where core and surface rotation rates are the furthest apart (middle and bottom right). The core H and core He burning phases are both shown in thick line widths, while the RGB phase is shown in thinner line width. Starting with the comparison of the surface rotation rates (dashed lines), we see that the 50 km s−1 models, with the TS dynamo (wTS) and without (nTS), reach five of the seven data points, while the 25 km s−1 model reaches one of the seven and the 150 km s−1 model reaches none. We therefore set the initial rotation rate of all the models to 50 km s−1. The two other data points can be reached by reducing the initial mass of the models, see Appendix B.

thumbnail Fig. 1.

Effect of inclusion TS dynamo. Rotation rates of the core (solid line) and envelope (dashed line) of the models with the TS dynamo (wTS) and without (nTS). The initial rotation rates of the models are included in the legend. Data points are from Deheuvels et al. (2015).

When focussing on core rotation rates during the core helium burning phase, we see that all models including the TS dynamo (Ωc ≃ 104 nHz) are two orders of magnitude away from the data points. Including the TS dynamo improves the match to the observations as the difference between observations and the model without the TS dynamo (Ωc ≃ 106 − 7 nHz) is more than 3 orders of magnitude worse. We thus conclude that also for the 2.5 M stars, the TS dynamo does not provide enough coupling between core and envelope to reduce the core rotation rates enough to match asteroseismic observations of the core helium burning stars.

4. Additional viscosity needed to reproduce observations of helium burning stars

Now that we have shown that the models with and without the TS dynamo cannot reproduce the asteroseismic observations of the seven secondary clump stars from Deheuvels et al. (2015), we continue by determining the strength of the missing process of angular momentum transport as a first step to revealing its physical nature. To do so, a constant νadd is added to the transport of angular momentum. We stress, however, that we do not believe the missing process of angular momentum transport is constant.

4.1. Determination of the additional viscosity needed to reproduce the Deheuvels et al. (2015) data

From Eggenberger et al. (2017) we know that the efficiency of the unknown transport process for angular momentum increases with stellar mass. Therefore, in this study a stronger process is expected than employed by Eggenberger et al. (2012), studying a 1.5 M star, and Eggenberger et al. (2017), studying a 0.84 M star.

As mentioned before, we did not attempt to fit all stars separately, but we look for global trends instead. Using Fig. 2, we determined the global efficiency of the missing process of angular momentum in the seven KIC stars. Figure 2 shows the ratio of core to envelope rotation rate, which, as mentioned by Eggenberger et al. (2017), allows us to determine νadd independently of the initial rotation rate. The best match in Fig. 2 is νadd = 107 cm2 s−1, which matches five of the seven data points. The other two models included reach none (νadd = 106 cm2 s−1) or two (νadd = 108 cm2 s−1) of the data points. More importantly, the general trend shown by the data points is best matched by the model that includes a νadd of 107 cm2 s−1. Again, the two data points with the highest surface gravities cannot be reached (see Appendix B for how to reach these points).

thumbnail Fig. 2.

Determining the νadd. The ratio of core to surface rotation rate as a function of surface gravity for three models calculated with an initial rotational velocity of 50 km s−1, while the νadd is varied. The data points are from Deheuvels et al. (2015).

When comparing the lines in Fig. 3 to the lines in Figs. A.2 and A.3, we can determine the start of the core He burning phase in Fig. 3. This is at the lowest surface gravity, in the bottom right corner of the figure. Then, both surface gravity g and the core rotation rate Ωc increase in a short amount of time until steady core He burning sets in and a slow decrease in both surface gravity and the core rotation rate characterises the rest of this phase. All data points are positioned around the turning point of the trend in surface gravity. From Fig. A.3 it follows that these seven stars are thus in the early phases of core He burning.

thumbnail Fig. 3.

Determining the νadd. Core and surface (solid and dashed line, respectively) rotation rates as a function of surface gravity for three models in Fig. 2. The data points are from Deheuvels et al. (2015).

In Fig. 3 the core and surface rotation rates are shown for the same three models as in Fig. 2. This figure confirms the choice for the initial rotation rate because the data points for surface rotation are matched. Also in this comparison, the general trend shown by the data points is best matched by the model with a νadd of 107 cm2 s−1.

4.2. Time dependence of the additional viscosity

In Sect. 4.1 we showed that the mean efficiency of the missing transport mechanism in the seven stars of Deheuvels et al. (2015) is around 107 cm2 s−1 when adding the νadd at the start of the main sequence. In this section we investigate whether this is dependent on the evolutionary phase during which νadd is added to the calculation. By doing this we are able to determine whether there is a phase in which the transport of angular momentum dominates the rest of the evolution. In this section we focus on the evolution up to the core helium burning phase and in Sect. 5.2 we focus on the later phases to investigate the influence of the inclusion of νadd on the final white dwarf spin.

We calculated models that include the νadd only from the end of the main sequence and from the start of the core helium burning phase. For the first, we find that adding the same νadd is sufficient to reach the data points, see Fig. 4, and that this model is comparable to the model in which we included νadd from the start of the main sequence. Therefore, we conclude that the main sequence is not a dominant phase for angular momentum transport in our models and we have no arguments to exclude νadd during the main sequence either.

thumbnail Fig. 4.

Dependence on inclusion time. The models presented here have been calculated with the best fit parameters (50 km s−1, 107 cm2 s−1), apart from the model labelled “2.5 × 106 cm2 s−1”, while varying the moment of including the νadd. The labels reflect the phase when the νadd is included. For the 2.5 × 106 cm2 s−1 model, the moment of inclusion is at the start of the core helium burning phase.

The inclusion of νadd only at the start of the core helium burning phase changes the evolution of the rotation rates, see again Fig. 4. Without the νadd earlier in the calculation, the core rotation rate is higher at the start of the core helium burning phase in this model than in the models that do include νadd earlier in the evolution. This is why the line of this model starts at a different point (top-left corner) in the figure. However, this difference has disappeared around log10(g cm−1 s−2) ≃ 2.8. The location of the curve in this model is dependent on νadd, shown by the model labelled “2.5 × 106 cm2 s−1”, this number being the νadd added at the start of the core helium burning phase. Thus, when we add the νadd at the start of the core helium burning phase, we are still able to reach all data points. However, we then have to use a value in the range of 2.5 × 106 <  νadd cm2 s−1 <  107. While the data cannot rule out the models that include νadd at the start of the core helium burning phase, the data does favour earlier inclusion of νadd because no data points are found with an angular velocity of the core above 200 nHz.

When we suppress the νadd from the start of the core helium burning phase onwards, we are unable to reach any data points. The reason for this is that the molecular weight gradient is too strong and without any νadd there is no transport of angular moment over this gradient. Therefore, we conclude that the crucial phase for the transport of angular momentum is the start of the core helium burning phase.

5. White dwarf rotation rates

After the core helium burning phase, we continued the models until they reached the white dwarf phase. In between these two phases, the stars pass through the asymptotic giant branch phase (AGB). During this phase, the energy production comes from the hydrogen and helium burning shell, located between the core and envelope. The helium shell becomes unstable, resulting in thermal pulses (TP-AGB phase). Around 25 to 30 thermal pulses take place in this phase in our models, and between each TP a third dredge-up (TDU) can occur. During the TP-AGB phase mass loss is enhanced, leading to removal of the envelope. Via the planetary nebulae phase, the star moves to the white dwarf track.

5.1. Calculation of the AGB phase

We calculated the full AGB phase as we would have done when studying the s-process nucleosynthesis (see Pignatari et al. 2016; Battino et al. 2016, for details). For instance, for the mass-loss treatment during the AGB phase we used Blöcker (1995) with an efficiency of 0.01 at the start of the AGB phase, 0.04 from when the envelope is carbon rich, and to 0.5 when the convergence issues appear (see below). We also used calibrated parameters for convective boundary mixing specifically for the AGB phase.

This is an improvement compared to the works of Suijs et al. (2008, no AGB specific mass loss, manually stopped models somewhere in AGB phase), Tayar & Pinsonneault (2013, no details given apart from initial mass and rotational velocity), and Cantiello et al. (2014, unphysical large mass loss efficiencies in the AGB phase which shorten this phase). By calculating the whole AGB phase, we can investigate the effects of the νadd on the thermal pulse cycle by investigating both the transport of angular momentum and the s-process nucleosynthesis, and compare them to the standards models without νadd. Details will be presented in a forthcoming paper. We report that the models with νadd included during the TP-AGB phase are able to transport angular momentum during the TDUs. This is due to the TDU reducing the molecular weight gradient and therefore the (local) barrier that has to be overcome to transport angular momentum.

It is common for convergence issues to arise during the final TPs in calculations like these and we report that these issues also occur in all models presented in this paper. There are two options for how to proceed, the first being the continuation of the AGB phase with a higher mass loss rate and the second the ejection of the whole remaining envelope (see Wood & Faulkner 1986; Herwig 2001; Sweigart 1999; Lau et al. 2012). We proceed with the models by increasing the mass loss parameter from 0.04 to 0.5, which allows for a smooth continuation of the models into the white dwarf phase.

5.2. Final spins of best fit models

In this section we show the comparison between the calculated rates and the observed white dwarf rotation rates by Hermes et al. (2017) and the compilation by Kawaler (2015). Most pulsating white dwarfs (WDs) in these two papers are DAVs, variable WDs with spectral type DA having only hydrogen absorption lines in their spectra. These pulsating WDs can be found in a specific temperature regime where their surface hydrogen has to become partially ionised. This regime for about 0.6 M white dwarfs is between 12 600 and 10 600 K, and we show the rotational periods of our models when passing through that same temperature regime in Fig. 5. The observational points from other pulsating white dwarfs are depicted as black crosses, while the DAVs are shown as black diamonds. The number of observed white dwarf periods is still low (36, we removed EPIC 201730811 because it is in a post-common envelope close binary according to Hermes et al. 2015), no statistical comparison is provided. White dwarf spins are also available for magnetic white dwarfs (see Kawaler 2015, for a summary). All of our models are non-magnetic, with only one exception, so we do not include these data points in our comparison.

thumbnail Fig. 5.

WD periods as a function of mass. The data points are from Kawaler (2015) and Hermes et al. (2017). The black diamonds are the DAVs, the black crosses are other pulsating white dwarfs. All coloured symbols are our predicted WD periods: the sphere and hexagon are the models without νadd; the triangles are the models with different values of νadd; the square is the model that excludes νadd at the end of core helium burning; and the star, cross, and plus signs are the models that exclude νadd at different times during the core helium burning phase.

All coloured symbols in Fig. 5 are WDs from our models. The two blue symbols correspond to the models introduced in Sect. 3, where we tested the impact of the TS dynamo. These models are the only ones without νadd in Fig. 5. As already shown by Suijs et al. (2008) and Cantiello et al. (2014), the model without the TS dynamo (nTS, dark blue circle) is orders of magnitude lower than the observed white dwarf periods. The model that does include the TS dynamo (wTS, light blue hexagon) reaches the lower limit of observed white dwarf periods, but as we saw before this model does not reach the observed periods of core He burning stars.

All models that include νadd in Fig. 5, have a spin period that is larger than all observed white dwarf rotation rates. There are three models with νadd of 106, 107, and 108 cm2 s−1 included during the whole calculation (three triangles), and one model where we excluded the νadd of 107 cm2 s−1 from the end of the core He burning phase (square). All these models are introduced in Sect. 4.1, except for the last one. From the previous section, we know that only the models labelled “107 cm2 s−1” and “end core He b” match the core He burning observations. However, they all transport too much angular momentum in the later phases of the evolution to match the white dwarf observations. Even the model that does not include νadd after the core He burning phase is finished does not reach the observed white dwarfs periods. Therefore, the efficiency of the missing process of angular momentum is negligible after the end of the core He burning phase according to our models, and the efficiency of the missing process also has to change during the core helium burning phase itself.

To investigate this last conclusion in more detail, we calculated models where we include νadd at the ZAMS and exclude it at different moments during the core helium burning phase. The whole core helium burning phase lasts for 183 Myr in these models and νadd has been excluded from times that correspond to 1/4, 2/4, and 3/4 of that time span. After excluding νadd we continue the calculation into the white dwarf phase. These three new models have also been included in Fig. 5. Again the rotational period within the DAV temperature range is used1. All three models are located within the range of observed white dwarf periods, and all three therefore match both the core helium burning and white dwarf observed rotation rates.

6. Conclusions

In this paper, we investigated the efficiency of the missing process of angular momentum by calculating 1D stellar evolution models with an initial mass of 2.5 M. As observational tests, we used the observed core and surface rotation rates of core helium burning stars as published by Deheuvels et al. (2015) and white dwarf rotational periods published by Kawaler (2015) and Hermes et al. (2017). The main conclusions of this paper are the following:

  • As for the 1.5 M of Cantiello et al. (2014), the 2.5 M models including the TS dynamo do not provide enough coupling between core and envelope to match asteroseismic observations of core rotation rates.

  • We have added a constant additional viscosity to our model as a first step towards revealing the physical nature of the missing process of angular momentum transport.

  • We are able to match the core rotation rates published by Deheuvels et al. (2015) by adding νadd = 107 cm2 s−1 and using an initial rotational velocity of 50 km s−1. This order of magnitude for νadd is independent of stellar evolution code, and initial mass (see Appendix B).

  • The trends identified by Eggenberger et al. (2017) concerning the increase in νadd with both initial mass and evolutionary phase are confirmed here. See Table 2 for an overview of all published studies on νadd. The strong increase in νadd from the two lower mass studies to this 2.5 M study suggests that when increasing the initial mass of the star, the change from radiative to convective core has less effect on the efficiency of the missing process of angular momentum than the absence of helium flashes in the more massive stars.

    Table 2.

    Summary of all published values for νadd to date.

  • We show that the dynamical instabilities (DSI and SH) are not attributed to the transport of angular momentum from ZAMS to the end of core helium burning in our models (see Appendix B).

  • We show that the extra transport of angular momentum that fits the observations of the core helium burning phase leads to rotation periods in the WD phase that are too high. Our results show that the efficiency of the missing process needs to change during the core helium burning phase, and must be strongly decreased before the end of the core helium burning phase.

  • When excluding νadd at 1/4, 2/4, or 3/4 of the whole duration of the core helium burning phase, our models match the observed rotation rates of both the set of core helium burning stars and the set of white dwarfs.

  • This implies that transport processes for which the efficiency only depends on the amount of differential rotation (such as the diffusive mixing introduced in Spada et al. 2016, based on the AMRI by Rüdiger et al. 2007) are incompatible with the result that the missing process has to be strongly decreased by the end of the core helium burning phase, unless an inhibiting effect is included to facilitate the decrease. A consequence of this work is that we have all initial parameters for the follow-up study, which will focus on the s-process production in rotating AGB stars. For this study, having a core rotation rate in the AGB phase that is consistent with asteroseismic observations of earlier and later evolutionary phases is crucial.


1

Apart from model “1/4” because this model undergoes a very late thermal pulse (VLTP) during the WD phase and is rebrightened before the DAV temperature range is reached. Convergence issues prevent the model from returning to the WD phase. We therefore calculated the rotational period of this WD just before the VLTP.

2

This principle is not a physical law, but an empirical observation confirmed by numerical simulations. It states that when a region within a burning shell contracts the region outside the shell will expand, and vice versa.

Acknowledgments

We thank the referee, Prof. C. Tout, for his careful reading of the manuscript and comments that greatly improved the clarity and flow of the paper. This work has been supported by the European Research Council (ERC-2012-St Grant 306901, ERC-2016-CO Grant no. 724560), and the EU COST Action CA16117 (ChETEC). RH acknowledges support from the World Premier International Research Centre Initiative (WPI Initiative), MEXT, Japan. This work is part of the BRIDGCE UK Network (www.bridgce.ac.uk) a UK-wide network established to Bridge Research in the different disciplines related to the galactic chemical evolution and nuclear astrophysics. This paper includes data collected by the Kepler mission. Funding for the Kepler mission is provided by the NASA Science Mission directorate.

References

  1. Abolfathi, B., Aguado, D. S., Aguilar, G., et al. 2018, ApJS, 235, 42 [NASA ADS] [CrossRef] [Google Scholar]
  2. Aerts, C., Reeth, T. V., & Tkachenko, A. 2017, ApJ, 847, L7 [NASA ADS] [CrossRef] [Google Scholar]
  3. Battino, U., Pignatari, M., Ritter, C., et al. 2016, ApJ, 827, 30 [NASA ADS] [CrossRef] [Google Scholar]
  4. Beck, P. G., Montalban, J., Kallinger, T., et al. 2012, Nature, 481, 55 [NASA ADS] [CrossRef] [Google Scholar]
  5. Belkacem, K., Marques, J. P., Goupil, M. J., et al. 2015, A&A, 579, A31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Blöcker, T. 1995, A&A, 297, 727 [NASA ADS] [Google Scholar]
  7. Borucki, W. J., Koch, D., Basri, G., et al. 2010, Science, 327, 977 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  8. Caleo, A., Balbus, S. A., & Tognelli, E. 2016, MNRAS, 460, 338 [NASA ADS] [CrossRef] [Google Scholar]
  9. Cantiello, M., Mankovich, C., Bildsten, L., Christensen-Dalsgaard, J., & Paxton, B. 2014, ApJ, 788, 93 [NASA ADS] [CrossRef] [Google Scholar]
  10. Ceillier, T., Tayar, J., Mathur, S., et al. 2017, A&A, 605, A111 [NASA ADS] [CrossRef] [EDP Sciences] [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., & Mosser, B. 2015, A&A, 580 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Deheuvels, S., Ouazzani, R. M., & Basu, S. 2017, A&A, 605, A75 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Edelmann, P. V. F., Röpke, F. K., Hirschi, R., Georgy, C., & Jones, S. 2017, A&A, 604, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Eggenberger, P., Maeder, A., & Meynet, G. 2005, A&A, 440, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Eggenberger, P., Meynet, G., Maeder, A., et al. 2008, Ap&SS, 316, 43 [NASA ADS] [CrossRef] [Google Scholar]
  18. Eggenberger, P., Montalbán, J., & Miglio, A. 2012, A&A, 544, L4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Eggenberger, P., Lagarde, N., Miglio, A., et al. 2017, A&A, 599, A18 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Ferguson, J. W., Alexander, D. R., Allard, F., et al. 2005, ApJ, 623, 585 [NASA ADS] [CrossRef] [Google Scholar]
  21. Fuller, J., Lecoanet, D., Cantiello, M., & Brown, B. 2014, ApJ, 796, 17 [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. Grevesse, N., Noels, A., Hauck, B., Paltani, S., & Raboud, D. 1993, Perfectionnement de l’Association Vaudoise des Chercheurs en Physique, 205 [Google Scholar]
  24. Heger, A., Langer, N., & Woosley, S. E. 2000, ApJ, 528, 368 [NASA ADS] [CrossRef] [Google Scholar]
  25. Hermes, J. J., Gänsicke, B. T., Bischoff-Kim, A., et al. 2015, MNRAS, 451, 1701 [NASA ADS] [CrossRef] [Google Scholar]
  26. Hermes, J. J., Gänsicke, B. T., Kawaler, S. D., et al. 2017, ApJS, 232, 23 [NASA ADS] [CrossRef] [Google Scholar]
  27. Herwig, F. 2001, Astrophys. Space Sci., 275, 15 [NASA ADS] [CrossRef] [Google Scholar]
  28. Herwig, F., Blöcker, T., Schoenberner, D., & El Eid, M. 1997, A&A, 324, L81 [NASA ADS] [Google Scholar]
  29. Hirschi, R., & Maeder, A. 2010, A&A, 519, A16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Kawaler, S. D. 2015, in 19th European Workshop on White Dwarfs, eds. P. Dufour, P. Bergeron, & G. Fontaine, ASP Conf. Ser., 493, 65 [NASA ADS] [Google Scholar]
  31. Kippenhahn, R., Weigert, A., & Weiss, A. 2013, Stellar Structure and Evolution (Berlin Heidelberg: Springer-Verlag) [Google Scholar]
  32. Lau, H. H. B., Gil-Pons, P., Doherty, C., & Lattanzio, J. 2012, A&A, 542, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Maeder, A. 2009, Physics, Formation and Evolution of Rotating Stars (Berlin Heidelberg: Springer) [Google Scholar]
  34. Marques, J. P., Goupil, M. J., Lebreton, Y., et al. 2013, A&A, 549, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Massarotti, A., Latham, D. W., Stefanik, R. P., & Fogel, J. 2007, ApJ, 135, 209 [Google Scholar]
  36. Meynet, G., Ekstrom, S., Maeder, A., et al. 2013, in Lect. Notes Phys., eds. M. Goupil, K. Belkacem, C. Neiner, F. Lignières, & J. J. Green (Berlin: Springer Verlag), 865, 3 [NASA ADS] [CrossRef] [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. Paxton, B., Bildsten, L., Dotter, A., et al. 2011, ApJS, 192, 3 [Google Scholar]
  39. Paxton, B., Cantiello, M., Arras, P., et al. 2013, ApJS, 208, 4 [NASA ADS] [CrossRef] [Google Scholar]
  40. Pignatari, M., Herwig, F., Hirschi, R., et al. 2016, ApJS, 225, 24 [Google Scholar]
  41. Pinçon, C., Belkacem, K., Goupil, M. J., & Marques, J. P. 2017, A&A, 605, A31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  42. Reimers, D. 1975, Mem. Soc. Roy. Sci. Liège, 8, 369 [Google Scholar]
  43. Rogers, T. M. 2015, ApJ, 815, L30 [NASA ADS] [CrossRef] [Google Scholar]
  44. Rogers, T. M., & McElwaine, J. N. 2017, ApJ, 848, L1 [NASA ADS] [CrossRef] [Google Scholar]
  45. Rogers, F. J., Swenson, F. J., & Iglesias, C. A. 1996, ApJ, 456, 902 [NASA ADS] [CrossRef] [Google Scholar]
  46. Rüdiger, G., Hollerbach, R., Gellert, M., & Schultz, M. 2007, Astron. Nachr., 328, 1158 [NASA ADS] [CrossRef] [Google Scholar]
  47. Spada, F., Gellert, M., Arlt, R., & Deheuvels, S. 2016, A&A, 589, A23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  48. Spruit, H. C. 1999, A&A, 349, 189 [NASA ADS] [Google Scholar]
  49. Spruit, H. C. 2002, A&A, 381, 923 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  50. Suijs, M. P. L., Langer, N., Poelarends, A.-J., et al. 2008, A&A, 481, L87 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  51. Sweigart, A. V. 1999, in New Views of the Magellanic Clouds, eds. Y. H. Chu, N. Suntzeff, J. Hesser, & D. Bohlender, IAU Symp., 190, 370 [NASA ADS] [Google Scholar]
  52. Tayar, J., & Pinsonneault, M. H. 2013, ApJ, 775, L1 [NASA ADS] [CrossRef] [Google Scholar]
  53. Tayar, J., & Pinsonneault, M. H. 2018, ApJ, 868, 150 [NASA ADS] [CrossRef] [Google Scholar]
  54. Tayar, J., Ceillier, T., García-Hernández, D. A., et al. 2015, ApJ, 807, 82 [NASA ADS] [CrossRef] [Google Scholar]
  55. Wood, P. R., & Faulkner, D. J. 1986, ApJ, 307, 659 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A: Evolution of rotation from ZAMS to core helium burning

thumbnail Fig. A.1.

Hertzsprung-Russell diagram of 2.5 M models, one without a νadd of 107 (dashed line) and one with a νadd of 107 (solid line). Neither model includes the TS dynamo.

thumbnail Fig. A.2.

Time evolution of surface gravity g. Timescale is t* = tWD-t), with tWD being the age of the star at the end of the calculations. The offset in time comes from a slightly longer white dwarf phase for the 107 cm2 s−1 model compared to the nTS model.

thumbnail Fig. A.3.

Coupling made visible: shown here is the evolution of core (solid line) and envelope (dashed) rotation rates from the ZAMS to the start of the AGB phase. Differences between the two models become visible at the start of the hydrogen shell burning phase, where the model without TS dynamo and νadd shows that the core and envelope rotation rates move apart, while the model including a νadd of 107 cm2 s−1 shows the rotation rates are coupled.

thumbnail Fig. A.4.

Angular velocity and angular momentum profiles of the nTS model for four moments as described in the label.

Before comparing our MESA models to the data points, we introduce the models by discussing their rotational evolution up to the core helium burning phase. The Hertzsprung-Russell diagram (HRD) of two models, one without a νadd (labelled nTS) and one with a νadd (labelled 107 cm2 s−1) is shown in Fig. A.1. The two models do not include the TS dynamo. This figure shows that the two models are comparable. The same is true for the evolution of the surface gravity g shown in Fig. A.2, where log g is shown versus log10(t*) ≃ log(tWD-t). In this figure, the horizontal segments of the lines are the core hydrogen (MS) and helium (Core He b) burning phases. The hydrogen shell burning phase takes place in a short amount of time at log10(t*/yr) ≃ 10.160, the hydrogen/helium shell burning phase after the core helium burning phase at log10(t*/yr) ≃ 10.154. This paper focusses on the core helium burning phase, which starts at log10(g cm−1 s−2) ≃ 1.8 and a log10(t*/yr) ≃ 10.160. Then, in a relatively short amount of time, log10(g cm−1 s−2)≃2.9 is reached. From there, during the remaining core helium burning phase the log g evolves with a constant slope until log10(g cm−1 s−2)≃2.4 is reached. This loop is visible in all following log10 g vs Ω figures, where the lower halve of the curves is the long-lasting phase.

thumbnail Fig. A.5.

Same as Fig. A.4, but for the model that includes a νadd of 107 cm2 s−1.

Figure A.3 shows the time evolution of the angular velocity of core Ωc (solid lines) and envelope Ωe (dashed lines) from the start of the main sequence to the start of the AGB phase. During the core burning phases, the rotation rates of core and envelope are close to constant in both models, with the model including νadd showing a near solid body rotation trend during the main sequence. The nTS model, however, shows large differences between core and envelope rotation rates during the shell burning phases. These phases are characterised with core contraction and envelope expansion (also known as the mirror principle2, see Kippenhahn et al. 2013), resulting in a steeply increasing core rotation rate and steeply decreasing envelope rotation rate.

The model including νadd shows different trends during the shell burning phases. The coupling provided by νadd allows for transport of angular momentum even when the core is contracting. As a result, the core rotation rate follows the trends of the envelope rotation rate and decreases during the shell burning phases (orange lines in Fig. A.3). This trend is as observed by Aerts et al. (2017), who compare a compilation of rotation rates of main sequence stars to the rotation rates of more evolved stars by Mosser et al. (2012). They find that there must be a drop in core rotation before or during the end of hydrogen and the start of helium core burning phases.

The details of the angular velocity Ω and corresponding angular momentum j profiles from core to surface are given in Figs. A.4 and A.5. Both figures show this profile at four moments in the evolution: the start and end of the main sequence and the start and end of the core helium burning phase. The solid body start of the models is visible in both figures, and from there the differences appear. As mentioned before, the angular velocity of the core and envelope in the model without νadd (left panel of Fig. A.4) evolve separately and oppositely due to the mirror principle. This effect is already visible at the end of the main sequence, and results in a difference between core and envelope rotation rate of several orders of magnitude at the end of the core helium burning phase. In the right panel the j profiles are shown. A decrease in j in a region during a certain phase indicates transport of angular momentum. A sharp feature is usually the outer edge of a convective zone, which creates a barrier for transport of angular momentum. The general lack of transport of angular momentum in the nTS model is visible in the j profiles of Fig. A.4, because they largely overlap.

When an additional viscosity of νadd = 107 cm2 s−1 is added, the differences between core and envelope angular velocity are smaller than in the nTS model (left panel of Fig. A.5). The whole star is close to solid body rotation up to the end of the core helium burning phase, as also shown in Fig. A.3. In this model a large amount of angular momentum is transported out of the core between the end of the main sequence and the start of the core helium burning phase (right panel of Fig. A.5). This efficient transport is also able to overcome the edge of convective regions, resulting in a lack of sharp features in the j-profiles. The transport continues during the core helium burning phase, creating a short moment at the end of the core helium burning phase when the convective envelope rotates at a higher angular velocity than the rest of the star.

Appendix B: Model uncertainties

thumbnail Fig. B.1.

Model uncertainties: the first three models presented here have been calculated with the best fit parameters (50 km s−1, 107 cm2 s−1), while the initial mass is varied. The fourth model includes only the ES and SSI instabitity. The models labelled “2.5 M” and “ES+SSI” overlap.

thumbnail Fig. B.2.

Code comparison: the models presented here are calculated with GENEC to show our conclusions are independent from evolutionary code.

In Sect. 4 we were not able to mass the data points at the highest surface gravities corresponding to KIC 5184199 and KIC 4659821. Here, we show that this is a consequence of setting the initial mass to 2.5 M. When matching the initial mass to the masses listed in Table 1, we can indeed match the highest surface gravities, as shown in Fig. B.1. For all models in this comparison, we use νadd = 107 cm2 s−1. The model with the lowest initial mass (2.2 M) reaches the higher surface gravities of the two data points earlier unreached. These two data points correspond to the observations of stars with initial masses of 2.18 ± 0.23 and 2.21 ± 0.18 M, indeed matching the lower initial mass of 2.2 M. When comparing the model with the highest initial mass (2.9 M) to the data points, we find that the star with the highest mass, KIC 7581399, of 2.90 ± 0.34 M, has a log10(g cm−1 s−2) = 2.843 ± 0.013 and is located on the 2.5 M model. This might imply that the actual mass of KIC 7581399 is located near the lower end of the error margin.

The implementation of rotation in MESA allows for the inclusion and exclusion of individual rotationally induced instabilities. The dynamical instabilities (DSI and SH) are not part of the GENEC models as published by Eggenberger et al. (2012,2017). Here we investigate their effects on the transport of angular momentum in the MESA models studied in this paper. To test this, we calculated an extra model with an initial mass of 2.5 M and νadd = 107 cm2 s−1 with only the ES and SSI included, and added this model to Fig. 3 with the label “ES+SSI”. The overlap of this model and the 2.5 M model, which also includes the dynamical instabilities, shows that the SH and DSI do not contribute to the transport of angular momentum. Edelmann et al. (2017) have already shown issues with the 1D implementation of the DSI in stellar evolutionary codes, and therefore being able to exclude this instability in studies on angular momentum transport reduces the uncertainties of our results. They also confirm that the settings of the GENEC models are satisfactory.

In Sect. 1 we mention that we ran the best fit models of Eggenberger et al. (2012,2017), and found νadd values that were similar to theirs to explain the observations. Here, we show GENEC models (see Eggenberger et al. 2008, for a description of this code and their implementation of rotation) calculated to match the MESA models of this study. Three GENEC models are shown in Fig. B.2, with their νadd and initial mass, as labelled. The same trends can be identified in these models as in the MESA models of earlier sections: when the initial mass is reduced, the data points at high surface gravities can be reached. A νadd of 107 cm2 s−1 provides a better fit than 5 × 106 cm2 s−1. Therefore, our conclusions are code independent.

Appendix C: Rotation near the core

As mentioned in Sect. 2.2, the numbers in the core rotation rates column in Table 1 are actually “near core” rotation rates. Their location is 0.1–1% of the normalised radius away from the most central point, see Fig. 5 in Deheuvels et al. (2015). In this region the obtained rotation rate is constant despite the noise in this figure. In Fig. C.1 we show a similar figure for the nTS and 107 cm2 s−1 models, where the rotation rate at the start and end of the core helium burning phase is shown. We see that the model including the extra νadd shows a constant trend in the region of interest at both times, as needed for the comparison to the data of Deheuvels et al. (2015). However, the nTS model shows a strong decrease in this region, providing another argument against these standard models.

thumbnail Fig. C.1.

For comparison with Fig. 5 of Deheuvels et al. (2015). The region of interest is between r/R of 10−3 and 10−2.

All Tables

Table 1.

Properties of the seven KIC stars from Deheuvels et al. (2015).

Table 2.

Summary of all published values for νadd to date.

All Figures

thumbnail Fig. 1.

Effect of inclusion TS dynamo. Rotation rates of the core (solid line) and envelope (dashed line) of the models with the TS dynamo (wTS) and without (nTS). The initial rotation rates of the models are included in the legend. Data points are from Deheuvels et al. (2015).

In the text
thumbnail Fig. 2.

Determining the νadd. The ratio of core to surface rotation rate as a function of surface gravity for three models calculated with an initial rotational velocity of 50 km s−1, while the νadd is varied. The data points are from Deheuvels et al. (2015).

In the text
thumbnail Fig. 3.

Determining the νadd. Core and surface (solid and dashed line, respectively) rotation rates as a function of surface gravity for three models in Fig. 2. The data points are from Deheuvels et al. (2015).

In the text
thumbnail Fig. 4.

Dependence on inclusion time. The models presented here have been calculated with the best fit parameters (50 km s−1, 107 cm2 s−1), apart from the model labelled “2.5 × 106 cm2 s−1”, while varying the moment of including the νadd. The labels reflect the phase when the νadd is included. For the 2.5 × 106 cm2 s−1 model, the moment of inclusion is at the start of the core helium burning phase.

In the text
thumbnail Fig. 5.

WD periods as a function of mass. The data points are from Kawaler (2015) and Hermes et al. (2017). The black diamonds are the DAVs, the black crosses are other pulsating white dwarfs. All coloured symbols are our predicted WD periods: the sphere and hexagon are the models without νadd; the triangles are the models with different values of νadd; the square is the model that excludes νadd at the end of core helium burning; and the star, cross, and plus signs are the models that exclude νadd at different times during the core helium burning phase.

In the text
thumbnail Fig. A.1.

Hertzsprung-Russell diagram of 2.5 M models, one without a νadd of 107 (dashed line) and one with a νadd of 107 (solid line). Neither model includes the TS dynamo.

In the text
thumbnail Fig. A.2.

Time evolution of surface gravity g. Timescale is t* = tWD-t), with tWD being the age of the star at the end of the calculations. The offset in time comes from a slightly longer white dwarf phase for the 107 cm2 s−1 model compared to the nTS model.

In the text
thumbnail Fig. A.3.

Coupling made visible: shown here is the evolution of core (solid line) and envelope (dashed) rotation rates from the ZAMS to the start of the AGB phase. Differences between the two models become visible at the start of the hydrogen shell burning phase, where the model without TS dynamo and νadd shows that the core and envelope rotation rates move apart, while the model including a νadd of 107 cm2 s−1 shows the rotation rates are coupled.

In the text
thumbnail Fig. A.4.

Angular velocity and angular momentum profiles of the nTS model for four moments as described in the label.

In the text
thumbnail Fig. A.5.

Same as Fig. A.4, but for the model that includes a νadd of 107 cm2 s−1.

In the text
thumbnail Fig. B.1.

Model uncertainties: the first three models presented here have been calculated with the best fit parameters (50 km s−1, 107 cm2 s−1), while the initial mass is varied. The fourth model includes only the ES and SSI instabitity. The models labelled “2.5 M” and “ES+SSI” overlap.

In the text
thumbnail Fig. B.2.

Code comparison: the models presented here are calculated with GENEC to show our conclusions are independent from evolutionary code.

In the text
thumbnail Fig. C.1.

For comparison with Fig. 5 of Deheuvels et al. (2015). The region of interest is between r/R of 10−3 and 10−2.

In the text

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.