Issue 
A&A
Volume 525, January 2011



Article Number  A5  
Number of page(s)  9  
Section  The Sun  
DOI  https://doi.org/10.1051/00046361/201015073  
Published online  26 November 2010 
Magnetic helicity fluxes in interface and flux transport dynamos
^{1}
NORDITA, AlbaNova University Center,
Roslagstullsbacken 23,
10691
Stockholm,
Sweden
email: piyalic@nordita.org
^{2}
Department of Astronomy, AlbaNova University Center, Stockholm
University, 10691
Stockholm,
Sweden
Received:
28
May
2010
Accepted:
15
September
2010
Context. Dynamos in the Sun and other bodies tend to produce magnetic fields that possess magnetic helicity of opposite sign at large and small scales, respectively. The buildup of magnetic helicity at small scales provides an important saturation mechanism.
Aims. In order to understand the nature of the solar dynamo we need to understand the details of the saturation mechanism in spherical geometry. In particular, we aim to understand the effects of magnetic helicity fluxes from turbulence and meridional circulation.
Methods. We consider a model with only radial shear confined to a thin layer (tachocline) at the bottom of the convection zone. The kinetic α owing to helical turbulence is assumed to be localized in a region above the convection zone. The dynamical quenching formalism is used to describe the buildup of mean magnetic helicity in the model, which results in a magnetic α effect that feeds back on the kinetic α effect. In some cases we compare these results with those obtained from a model with a simple algebraic α quenching formula.
Results. In agreement with earlier findings, the magnetic α effect has the opposite sign compared with the kinetic α effect and leads to a catastrophic decrease of the saturation field strength proportional to the inverse magnetic Reynolds number. At high latitudes this quenching effect can lead to secondary dynamo waves that propagate poleward because of the opposite sign of α. These secondary dynamo waves are driven by smallscale magnetic helicity instead of the smallscale kinetic helicity. Magnetic helicity fluxes both from turbulent mixing and from meridional circulation alleviate catastrophic quenching. Interestingly, supercritical diffusive helicity fluxes also give rise to secondary dynamo waves and grand minimalike episodes.
Key words: Sun: dynamo / Sun: activity / magnetic fields / magnetohydrodynamics (MHD)
© ESO, 2010
1. Introduction
The solar dynamo models developed so far that agree with solar magnetic field observations usually solve the αΩ mean field dynamo equations. The turbulent α effect first proposed by Parker (1955) is believed to be generated by helical turbulence in the convection zone of the Sun. Because α is generated by quadratic correlations of the smallscale turbulence, we need a closure in order to complete the set of mean field equations, e.g., the first order smoothing approximation (FOSA), and express the mean electromotive force in terms of the mean magnetic fields. This turbulent α encounters a critical problem when the energy of the mean field becomes comparable to the equipartition energy of the turbulence in the convection zone, and therefore it becomes increasingly difficult for the helical turbulence to twist rising blobs of magnetic field. The solar dynamo modellers have traditionally used what is referred to as algebraic alpha quenching to mimic this phenomena. This involves replacing α by , an expression introduced by Jepps (1975), or by , where α_{0} is the unquenched value and R_{m} is the magnetic Reynolds number, is the mean magnetic field and B_{eq} is the equipartition magnetic field. The latter expression has been discussed since the early work of Vainshtein & Cattaneo (1992). The R_{m} in the denominator is included because the smallscale fluctuating magnetic field reaches equipartition long before the mean magnetic field does. This has been supported by several numerical experiments to determine the saturation behaviour of α (e.g. Cattaneo & Hughes 1996; Ossendrijver et al. 2002). Given the large magnetic Reynolds numbers of astronomical objects, this phenomenon is referred to as catastrophic quenching.
After the discovery of the layer of strong radial shear (called the tachocline by Spiegel & Zahn 1992) at the bottom of the solar convection zone, Parker (1993) proposed a new class of solar dynamo models called the interface dynamo. In these models the shear is confined to a narrow overshoot layer immediately beneath the convection zone, which is also the region of α effect. The dynamo wave propagates in a direction given by the ParkerYoshimura rule at the interface between the two layers defined by a steep gradient in the turbulent diffusivity. The toroidal field produced by the stretching effect of the shear is much stronger than the poloidal field and remains confined in the overshoot layer, away from the region where the α effect operates. Note that the interface dynamo model may have serious problems when solarlike rotation with positive latitudinal shear is included (Markiel & Thomas 1999). Similarly, in the BabcockLeighton class of flux transport models (Choudhuri et al. 1995; Durney 1995) the toroidal and the poloidal fields are produced in two different layers. Unlike in the interface dynamo models, the coupling between the two layers is mediated both by diffusion and the conveyor belt mechanism of the meridional circulation.
It has been proposed that in interface and BabcockLeighton type dynamos, the α effect is not catastrophically quenched at high R_{m} because the strength of the toroidal field is very weak in the region of finite turbulent α (e.g. Tobias 1996; Charbonneau 2005). However, according to our knowledge, not much has been done to study the variation of the strength of the saturation magnetic field with the magnetic Reynolds number for these classes of αΩ dynamos. Zhang et al. (2006) made an attempt to reproduce the surface observations of current helicity in the Sun with a 2D mean field dynamo model in spherical coordinates coupled with the dynamical quenching equation. In a separate paper (Chatterjee et al. 2010) we demonstrated that interface dynamo models are also subject to catastrophic quenching.
Magnetic helicity has been identified as a key player in any dynamo process. Because the net magnetic helicity is conserved, any process can only create large and smallscale magnetic helicities of opposite signs. This was shown in direct numerical simulations of a twisted magnetic flux tube by Blackman & Brandenburg (2003), in good agreement with observations of sigmoid loops in the solar corona. However, the smallscale magnetic helicity backreacts on the helical turbulence and quenches the dynamo (Blackman & Field 2000; Kleeorin et al. 2000).
It has now been shown that this mechanism reduces the saturation strength of the magnetic field (B_{sat}) with increasing magnetic Reynolds number (R_{m}). Nevertheless this constraint may be lifted if the system is able to rid itself of smallscale helicity through at least one of several ways such as open boundaries, advective, diffusive and sheardriven fluxes (Shukurov et al. 2006; Zhang et al. 2006; Sur et al. 2007; Käpylä et al. 2008; Brandenburg et al. 2009; Guerrero et al. 2010). Even though the helicity constraint in direct numerical simulations (DNS) of dynamos with strong shear have been clearly identified, the results can be matched with mean field models with a weaker algebraic quenching than α^{2} dynamos (Brandenburg et al. 2001). It is possible to include this process in meanfield dynamo models through an equation describing the evolution of the smallscale current helicity. We shall refer to this equation as the dynamical quenching mechanism.
In this paper we perform a series of calculations with mean field αΩ models in spherical geometry along with a dynamical equation for the evolution of α for magnetic Reynolds numbers in the range 1 ≤ R_{m} ≤ 2 × 10^{5}. An important feature of the calculation is that the region of strong narrow shear is separated from the region of helical turbulence. In addition to providing detailed results not mentioned in Chatterjee et al. (2010), this paper is also aimed at studying somewhat more complicated models like flux transport (FT) models, which include meridional circulation (MC). The role of diffusive helicity fluxes modelled into the dynamical quenching equation by using a Fickian diffusion term is also discussed for various models. Helicity fluxes across an equator can indeed be modelled by this diffusion term, as was shown by Mitra et al. (2010). In Sects. 2.1 and 2.2 we discuss the features of the two classes of αΩ models used. The formulation of dynamical α quenching is given in Sect. 2.3. The results for 1 < R_{m} < 2 × 10^{5}, with or without smallscale helicity fluxes, are presented in Sects. 3.1–3.3, while Sect. 3.4 presents the results for the flux transport dynamo model. Finally, we draw the conclusions of this study in Sect. 4.
2. αΩ dynamo models with dynamical α quenching
2.1. Interface dynamo
We solve the induction equation in a spherical shell assuming axisymmetry. Our dynamo model consists of the induction equations for the mean poloidal potential A_{φ}(r,θ) and the mean toroidal field B_{φ}(r,θ) (see Eqs. (2) and (3) of Guerrero et al. 2010, GCB from now on). Axisymmetry demands that for all variables ∂/∂φ = 0. The turbulent coefficients α and η_{t} in these equations are estimated as follows: from the mixing length theory we know that the turbulent diffusivity in the convection zone is given by (cf. Sur et al. 2008) (1)where u_{rms} is the rms velocity of the turbulent eddies, k_{f} is the wavenumber of the energycarrying eddies, corresponding to the inverse pressure scale height near the base of the convection zone. In our models, we have taken k_{f} = 10k_{1}, where k_{1} is the longest wave that can fit the domain latitudinally. Below the convection zone the magnetic diffusivity has the molecular value, η_{r}. Both regions are smoothly matched by the step function (2)where r_{e} = 0.73 R_{⊙}, and d_{e} = 0.025 R_{⊙}. Thus the magnetic diffusivity profile is given by (3)The turbulent kinetic α effect considered in this model has the form (4)where g_{α} is a nondimensional coefficient equal to 1 or R_{m} depending on the assumed form of algebraic quenching in the models and is the value of the magnetic field at equipartition between the magnetic and the turbulent kinetic energies. The amplitude of the kinetic α effect is estimated by considering the first order smoothing approximation (FOSA) as being α_{0} = τϵ_{f}ω_{rms}u_{rms}/3, where ω_{rms} is the rms vorticity of the turbulence and τ ~ (k_{f}u_{rms})^{1} is the eddy correlation time scale. The prefactor ϵ_{f}, usually of the order of 0.1 or less, is used because (u·ω)_{rms} < u_{rms}ω_{rms}. The case ϵ_{f} = 1 means that the flow is maximally helical. These approximations give us an estimate of α_{0} in terms of eddy diffusivity η_{t} and forcing scale k_{f} as We would consider ϵ_{f} instead of α_{0} as a free parameter in the model apart from η_{t}. Even though the helical turbulence pervades almost the entire convection zone, here we take r_{a} = 0.77 R_{⊙} and d_{a} = 0.015 R_{⊙}, so that we can have a large spatial separation between the shear and turbulent layers. Consequently we consider a differential rotation profile like that in the high latitude tachocline of the Sun given by (5)where Ω_{0} = 14 nHz, r_{w} = 0.68 R_{⊙} and d_{w} = 0.015 R_{⊙}. The radial profiles of η_{t}, α and ∂Ω/∂r are plotted as a function of fractional radius r/R_{⊙} in Fig. 1. In that figure it is possible to see that the region of strong radial shear is separated from the region of helical turbulence and the diffusivity has a strong gradient at a radius lying between these two source layers. The reason of this is to decrease the cycle period T_{cyl} of the oscillatory dynamos to a reasonably small fraction of the diffusion time t_{diff}.
Fig. 1 Profiles of radial shear ∂Ω/∂r, α and η as a function of fractional solar radius. 
2.2. Flux transport dynamo
Apart from the interface dynamo, we also run numerical experiments for models that include meridional circulation and consider the socalled BabcockLeighton (BL) α effect. In these models the dynamo is thought to be operating in two separated regions, the toroidal component of the magnetic field is produced at the base of the convection zone and the poloidal field is produced by the decay of tilted bipolar active regions on the solar surface. The MC plays an important role in the dynamics of the system by advecting the magnetic flux and connecting both source regions. For this reason, these models are often referred to as fluxtransport (FT) dynamo models. In the literature FT models have been studied extensively by several authors (Dikpati & Charbonneau 1999; Chatterjee et al. 2004; Guerrero & Dal Pino 2008, and references therein). These models have now reached a stage where they are able to reproduce the butterfly diagram as well as the correct phase between the polar fields and the toroidal fields.
We recall that the α_{K} is now not due to the helical turbulence in the bulk of the convection zone, but due to a phenomenological BL α where the poloidal field is produced from the toroidal field by decay of tilted bipolar active regions. For the models in this section, the BL α is assumed to be concentrated only in the upper 0.05% of the convection zone and becomes maximum at ± 40° latitude and goes to zero at poles. The expression for α_{K} is (6)where d = 0.015 R_{⊙}.
For the meridional circulation we use the analytical expressions given by van Ballegooijen & Choudhuri (1988), where the radial and latitudinal components of the velocity field are given by (7)(8)where ζ = R_{⊙}/r−1, r_{b} = 0.71R_{⊙}, ζ_{b} = R_{⊙}/r_{b}−1, and . In the equations above, the flow is poleward at the surface with a maximum amplitude of u_{0} = 20 m s^{1}, and the equatorward return flow occurs at r_{b} with an amplitude of around 3 m s^{1}. Unlike in flux transport dynamo models, the meridional circulation does not reverse the direction of propagation of the dynamo wave in interface dynamo models as long as the meridional circulation is confined within the convection zone (Petrovay & Kerekes 2004). The shear is still radial and given by Eq. (5) with r_{w} = 0.7 R_{⊙}. Finally, the turbulent diffusivity has the same profile as in Eq. (3), but with η_{t} = 2 × 10^{11} cm s^{1} and r_{e} = 0.7 R_{⊙}.
2.3. Dynamical α quenching
Pouquet et al. (1976) first showed that the turbulent α effect is modified due to the generation of smallscale helicity in the way given by Eq. (9) below. The second term is sometimes referred to as the magnetic αeffect, (9)where ω, u, j, b denote the fluctuating component of the vorticity, velocity, current, and magnetic field in the plasma. Also, b = ∇ × a. It is possible to write an equation for the evolution of the magnetic part of α or α_{M} from the equation for the evolution of the smallscale magnetic helicity density using the relation (10)As the equation for is gaugedependent, it makes sense only to write an equation for the volumeaveraged quantity in order to avoid dependence on a specific gauge (Blackman & Brandenburg 2002). Our dynamo equations are independent of any gauge because we solve for the magnetic potential component A_{φ} with an axisymmetric constraint. It is important for us that the equation for α_{M} is also gaugeindependent. Subramanian & Brandenburg (2006) used the Gauss linking formula for the expression for h_{f} and wrote an equation independent of the gauge for the magnetic helicity density under the assumption that the correlation length for all fluctuating variables remains small compared to the system size at all times. Using Eq. (10) we write the same equation in terms of α_{M}, (11)where and are the mean field EMF and the mean magnetic field. We recall that in the mean field approach, we do not solve evolution equations for the fluctuating components a and b. Instead we model the effect of the small scales by writing an evolution equation for h_{f} similar to Eq. (11). Because the evolution of the magnetic field components (A_{φ}, B_{φ}) in the large scales is explicitly solved, we do not need a corresponding equation for the largescale helicity (). Moreover, and are in principle gaugedependent. However, if there is scale separation, can be written as a density of linkages and is therefore gaugeindependent (Subramanian & Brandenburg 2006), while is not. This was demonstrated in Hubbard & Brandenburg (2010).
The flux term, , consists of individual components, e.g., advection owing to the mean flow, VishniacCho fluxes (Vishniac & Cho 2001), effects of mean shear, diffusive fluxes, etc. When we make an ansatz for the flux of normalized smallscale helicity, , the effect on the largescale helicity flux is also taken into account by the nonlinear coupling between meanfield induction equations and Eq. (11). If smallscale magnetic helicity is removed by stellar winds in Eq. (11), a corresponding advection term in the mean field induction equation also changes the largescale helicity accordingly.
For the interface models we put F_{α} = 0 unless mentioned otherwise. For the flux transport dynamo (Sect. 2.2), the normalized helicity flux F_{α} in Eq. (11) is given by (12)where κ = κ_{0}η(r) is the diffusion coefficient for α_{M}. Let us define the diffusion time in the model as . The decay time in Eq. (11) is then . Note that we use g_{α} = 0 in Eq. (4) whenever we employ the dynamical quenching equation, because dynamical quenching is usually more important.
Our computational domain is defined to be the region confined by 0 ≤ θ ≤ π and 0.55 R_{⊙} ≤ r ≤ R_{⊙}. Unless otherwise stated, the boundary conditions for A_{φ} are given by a potential field condition at the surface (Dikpati & Choudhuri 1994) and A_{φ} = 0 at the poles. We also performed some calculations with the vertical field condition at the top boundary, which means that B_{θ} = B_{φ} = 0. At the bottom we use the perfect conductor boundary condition of Jouve et al. (2008) with A_{φ} = ∂(rB_{φ})/∂r = 0. But a more realistic perfect conductor boundary condition in our opinion would be ∂(rB_{θ})/∂r = ∂(rB_{φ})/∂r = 0. Also B_{φ} = 0 on all other boundaries. The equation for α_{M} is an initial value problem for F_{α} = 0. For finite fluxes we also set α_{M} = 0 at all boundaries. We checked that the results are not very sensitive to the different bottom boundary conditions given above, mainly because the bottom boundary is far removed from the dynamo region.
3. Results
3.1. Nonlinear interface dynamo without helicity fluxes
In order to study the R_{m} dependence of the saturation of the magnetic field in the interface dynamo of Sect. 2.1, we keep all the dynamo parameters the same for all runs and change η_{r} from 2 × 10^{5} cm^{2} s^{1} to 2 × 10^{10} cm^{2} s^{1} while keeping η_{t} fixed at 4 × 10^{10} cm^{2} s^{1}.
Fig. 2 Critical α in terms of a fraction of η_{t}k_{f} as a function of magnetic Reynolds number R_{m} for the interface dynamo model of Fig. 1. 
To be able to correctly compare the dynamo models for different R_{m}, it is important to first calculate the critical value of α_{0}, denoted by α_{c} for each model. This is done by solving the mean field dynamo equations with an algebraic quenching with g_{α} = 1. The values of α_{c} as a function of R_{m} are shown in Fig. 2. From this figure we observe that the model is most efficient for R_{m} ~ 20, because for this value the critical dynamo number is the lowest. A similar variation of α_{c} with the ratio η_{t}/η_{r} was obtained analytically for interface dynamos by MacGregor & Charbonneau (1997, see their Fig. 5A). All the simulations from now on are performed by setting α_{0} = 2α_{c}, corresponding to the R_{m} for each model. The cycle period of these dynamo models (T_{cyl}) slightly increases with R_{m} and is found to be of the order of 2t_{diff}, for the investigated range of R_{m}. Figure 3 shows the evolution of the volumeaveraged magnetic energy as a function of time for an algebraic form of quenching with g_{α} = 1. Note that in Fig. 3, the slopes in the kinematic phase are almost similar for all R_{m} within the errors in the numerical determination of α_{c}.
Fig. 3 Time evolution of the volumeaveraged magnetic energy in the domain scaled with the equipartition energy for the case with algebraic quenching and g_{α} = 1 for models with R_{m} = 1 (diamond), R_{m} = 20 (solid), R_{m} = 200 (dashed), R_{m} = 2 × 10^{3} and R_{m} = 2 × 10^{5} (triangles). 
Fig. 4 Same as Fig. 3 but with dynamical quenching and g_{α} = 0 for the R_{m} indicated in the figure. 
Fig. 5 Saturated value of the volumeaveraged magnetic energy scaled with the equipartition energy as a function of R_{m} for models with dynamical α quenching (triangles+solid) and algebraic quenching with g_{α} = 1 (squares + dashed) and with g_{α} = R_{m} (cross + dasheddotted). Note that all runs included in this figure used α_{0} = 2α_{c}. 
Now, we set g_{α} = 0 and consider a dynamically quenched α by solving Eq. (11) along with the linear induction equations for the mean fields A_{φ} and B_{φ}. The time evolution of the volumeaveraged magnetic energy for these systems for a range of magnetic Reynolds numbers is shown in Fig. 4. The strong R_{m} dependence, which is reminiscent of catastrophic quenching in all astrophysical dynamos, can be easily discerned from this figure. For high values of R_{m}, the saturation phase is clearly different from that obtained for algebraic quenching (see Fig. 3).
In this context it is important to recall an important difference between dynamos in periodic and open domains such as those considered here. In a periodic domain, the dynamo is expected to reach and sustain mean fields of the order of the inverse square root of the scale separation ratio (Blackman & Brandenburg 2002). However, in an open domain the final saturation field strength is always of the order of (Brandenburg & Dobler 2001; Brandenburg & Subramanian 2005), although the mean field may reach an early equipartition strength peak just at the time when the smallscale field reaches saturation (Fig. 3). The physical reality of this peak remains somewhat mysterious, although direct simulations also sometimes show this peak (Hubbard & Brandenburg 2010). In Fig. 5 we summarize the results of models with a simple algebraic quenching (Eq. (4)) for both g_{α} = 1 and g_{α} = R_{m}, as well as models with dynamical quenching (Eq. (11)). Both algebraically quenched models with g_{α} = R_{m} and those with dynamical α quenching exhibit a monotonic decrease of the saturation magnetic energy with R_{m}.
In our models we see that the region of strong toroidal field, B_{φ}, is different from that of the poloidal field produced by the α effect as shown in Fig. 6. We verified that the same dependence of the B_{sat} on R_{m} is reproduced in models where the two source regions are further spatially separated by setting r_{a} = 0.87 R_{⊙} instead of 0.77 R_{⊙}.
Fig. 6 Radial profiles of A_{φ} and B_{φ} at two different latitudes (λ) in the saturated phase for R_{m} = 2 × 10^{3}. 
The nature of the saturation curves of the magnetic energy (Fig. 4) is strongly governed by the ratio of t_{α} and T_{cyl}. For low R_{m}, t_{α} ≤ T_{cyl}, so that the time evolution of is flat in the saturation regime – similar to the algebraic quenching case (see Fig. 3). On the other hand, for R_{m} = 2 × 10^{3}, the decay time t_{α} ≫ T_{cyl} and so the system is underdamped, i.e., there are amplitude modulations of the magnetic energy before it settles to a final saturation value (see dashdotted line in Fig. 4). For R_{m} = 2 × 10^{5}, the underdamped oscillation has a long period (t_{α}) so that the model has to be run for more than 500 t_{diff} before the dynamo field starts becoming “strong” again. Because of the long computational time involved in this exercise we have not continued the calculation beyond 60 t_{diff}. Hence, the determination of saturation magnetic energy may be inaccurate.
Fig. 7 Timelatitude plot of α_{M}(0.72 R_{⊙},θ) for a) R_{m} = 20 and b) R_{m} = 200. 
A similar pattern is observed in the timelatitude plot of α_{M} as shown in Fig. 7. In this figure we note that for R_{m} = 20 (upper panel), α_{M} shows strong oscillations (see the modulation in the colour of the filled contours) since t_{α} ≪ T_{cyl}. For R_{m} = 200 (lower panel), t_{α} ~ T_{cyl}, so that the amplitude of the oscillation is weak because the α_{M} decays at the same rate at which it is produced. From the same figure it may be concluded that the smallscale current helicity, α_{M}, is predominantly negative (positive) in the Northern (Southern) hemisphere. For high values of R_{m} (≳2000) the models start showing changes in parity for t > 40t_{diff}. Nevertheless, the magnetic energy and the dynamo period, T_{cyl}, remain fairly constant even while the system fluctuates between symmetric and antisymmetric parity at irregular time intervals (see Fig. 8). This parity oscillation is absent in the corresponding models with algebraic quenching.
Fig. 8 a) Evolution of parity (purely dipolar = −1 and purely quadrupolar = +1) for an interface dynamo model with R_{m} = 2 × 10^{3} and no fluxes. b) A small part in the timelatitude diagram of B_{φ} indicated by dotted lines in a) where the parity is changing from quadrupolar to dipolar. 
Fig. 9 Timelatitude plot of the toroidal field a) and c) and α_{M}b) and d) for interface dynamo models with α = 4α_{c} and R_{m} = 20. 
3.2. Secondary dynamo waves with F_{α} = 0
An interesting result emerges when we increase the value of the kinetic α effect (α_{0} = 4α_{c} instead of 2α_{c} as in Sect. 3.1) in the interface dynamo model with dynamical α quenching and R_{m} = 20. In this case we observe that in addition to the primary dynamo wave travelling equatorward, a polewardpropagating secondary dynamo wave appears in the butterfly diagram for B_{φ} at r = 0.72 R_{⊙} (see fingerlike projections at high wave number in Fig. 9a). A weak signature of this secondary wave can also be seen in the butterfly diagram at 0.8 R_{⊙} in Fig. 9b.
The existence of these secondary dynamo waves can be explained by studying the source term in Eq. (11), i.e., . Because of a low value of R_{m}, the toroidal field generated in the shear layer diffuses into the overshoot layer, and by virtue of the term , generates a region of α_{M} that is positive (negative) in the northern (southern) hemisphere. Because α_{K} → 0 for r ≤ 0.73 R_{⊙}, the total α effect there has the sign of α_{M}. The new dynamo wave therefore travels poleward according to the ParkerYoshimura rule.
Vishniac & Cho (2001) also presented a case of an αΩ dynamo driven by a supercritical helicity flux. Their mechanism requires a finite initial magnetic field unlike here, where the initial field is ~10^{6}B_{eq}. The difference compared to the case above is that the mean field dynamo is not driven by supercritical Vishniac & Cho fluxes, but by a local generation of smallscale magnetic helicity. Note that in this case the secondary dynamo wave is energetically powered by the kinematic part of the helical convection. Blackman & Field (2004) also found a magnetichelicity driven dynamo (MHDD) in addition to a kinetichelicity driven dynamo (KHDD) by solving a coupled set of equations for largescale helicity and the mean smallscale helicity. They highlight an important difference between KHDD and MHDD, namely that in the former, the α_{K} produces smallscale and largescale magnetic helicity of opposite signs, whereas in the latter magnetic helicity has same sign in both scales. Our calculations in this section agree with their observation in the sense that α_{M} has the same sign as . Evidence for dominance of magnetically generated α in stratified magnetorotational turbulence in accretion discs has also been found by Gressel (2010).
We have not observed any evidence of chaotic behaviour in the range of magnetic Reynolds number 20 ≤ R_{m} ≤ 2 × 10^{5} for supercritical α ≤ 4α_{c} in agreement with Covas et al. (1997). However, if the α effect is highly supercritical, the dynamical quenching formula for α_{M} is insufficient for dynamo saturation, and additional algebraic quenching terms are needed (Kleeorin & Rogachevskii 1999). We continue the discussion of secondary dynamo waves driven by diffusive magnetic helicity fluxes in the next section.
3.3. Supercritical diffusive magnetic helicity fluxes
Recently, Brandenburg et al. (2009) showed that catastrophic quenching in onedimensional α^{2} dynamos can be alleviated by introducing a Fickian diffusive flux in Eq. (11) given by (13)There was an attempt to calculate the diffusion coefficient from direct numerical simulations and it was found that κ ~ 0.3η_{t} for R_{m} ~ 20 (Mitra et al. 2010). In the context of of αΩ dynamos GCB showed that even a very small diffusive flux for the mean smallscale helicity alleviates catastrophic quenching for a scale separation k_{f}/k_{1} ~ 10. Their Fig. 7 moreover shows that levels off at high R_{m} for finite diffusive fluxes even though the value of at which the curve levels off may vary with k_{f}.
Fig. 10 a) Time evolution of the volume averaged magnetic energy in the domain scaled with the equipartition energy for models with R_{m} = 2 × 10^{3} and κ_{0} = 10^{5} (∗ ) as well as with R_{m} = 2 × 10^{5} and κ_{0} = 10^{5} (solid line) for the interface dynamo model of Sect. 3.1 with dynamical α quenching. The saturation curve for zero fluxes for the model with R_{m} = 2 × 10^{5} has been shown by the dashed line; b) and c) show timelatitude diagrams for the toroidal field at the depths indicated for models with R_{m} = 2 × 10^{3} and κ_{0} = 10^{5}. 
Fig. 11 The timelatitude plot of toroidal field a) and c), and α_{M}b) and d) with α = 2α_{c} for R_{m} = 2 × 10^{5} and κ_{0} = 0.01. Meridional snapshots of e) sign(B_{φ})( B_{φ}  /B_{eq})^{1/2} and f) α_{m} × 10^{3} for the same case. 
For κ = 0, the saturation curves in Fig. 4 show that the B_{sat} either displays underdamped oscillations for R_{m} = 2 × 10^{3} or goes through very low values for R_{m} ~ 2 × 10^{5} and takes a long time to relax to a steady amplitude. In this section we introduce a diffusive flux with κ(r) = κ_{0}η(r). Figure 10 shows the results for models with two different values of R_{m} and κ. For R_{m} = 2 × 10^{3} and κ_{0} = 10^{5}, the saturation energy is comparable to the corresponding case with κ = 0 (see dashdotted line in Fig. 4). But now grand minimalike episodes with respect to the primary mode appear in the system, as may be seen in the oscillations in the volumeaveraged energy (represented by ∗ ’s in Fig. 10a) with a period ~5 times the period of the equatorwardpropagating mode long after saturation. Another interesting behaviour can be discerned from the butterfly diagram of the toroidal field (Fig. 10b, c). It appears that the dynamo here is governed by the competition between the equatorwardpropagating primary mode and the polewardpropagating secondary mode.
For R_{m} = 2 × 10^{5} and κ = 0.01, we observe underdamped oscillations with a final B_{sat} ~ 0.1B_{eq}. On studying the corresponding butterfly diagrams (Fig. 11a–d) we find a poleward propagating mode caused by the radial diffusion of the α_{M} into the stable layers. This was not possible for the model with the same R_{m} but κ = 0. Figures 11e, f show meridional snapshots of sign(B_{φ})( B_{φ}  /B_{eq})^{1/2} and α_{M} in order to get a clear idea of the distribution of magnetic fields.
The polewardpropagating mode is now driven by supercritical diffusive helicity fluxes, as opposed to supercritical Vishniac & Cho fluxes (see Brandenburg & Subramanian 2005, for examples of this behaviour). For the same model (R_{m} = 2 × 10^{5}) there exists a critical κ_{c} ~ 10^{5} such that the secondary dynamo wave does not appear if κ_{0} < κ_{c}. The volumeaveraged magnetic energy in this case decays eventually. Note that this threshold for κ is highly dependent on R_{m}. For instance a model with R_{m} = 2 × 10^{3} and κ_{0} = 10^{5} produces a dynamo with finite saturation magnetic energy and dynamo wave propagation governed by α_{M}, whereas for κ_{0} = 10^{4}, the dynamo shows a runaway growth.
Unlike in GCB, we used supercritical helicity fluxes in this section. Note that the critical value, κ_{c}, is much lower than that used in GCB where the maximum value of κ_{0} used is ~10. Interface dynamo models are a particular case of αΩ models where a strong toroidal field is produced below the convection zone. Even a low flux of α_{M} into the stable layers can produce very large magnetic fields and power a secondary dynamo. Direct numerical simulations of α^{2} dynamos have established that a largescale magnetic field is easily excited on the scale of the system i.e., for a high k_{f}/k_{1} ratio (Archontis et al. 2003). However the length scale of the magnetic field in Figs.10b, c and 11a, c, e is comparable to , which suggests that the degree of scale separation may have become insufficient to write the electromotive force as a simple multiplication, as is done in the expression . Then it may become necessary to write the electromotive force as a convolution, which essentially corresponds to a lowpass filter (see, e.g., Brandenburg et al. 2008). However, we have not pursued this aspect any further.
3.4. Flux transport BabcockLeighton dynamo
In this section we perform numerical experiments of the flux transport dynamo model described in Sect. 2.2. As in Sect. 3.1, we first find the critical α_{BL} required to have a selfexcited dynamo. In this case α_{c} = 5.1 m s^{1} for R_{m} = 2 × 10^{3}. We pursue the rest of the calculations with α_{BL} = 6.0 m s^{1} and including dynamical α quenching. This slightly supercritical choice of α_{BL} avoids the production of very large α_{M} leading to secondary dynamos as discussed in Sect. 3.2. We should emphasize that Eq. (9) represents a first order correction to α and should be treated with caution while in supercritical regimes.
Fig. 12 Time evolution of the volume averaged magnetic energy for the flux transport dynamo model of Sect. 3.2 for R_{m} = 2 × 10^{3} with κ_{0} = 0.3 (dashed); R_{m} = 2 × 10^{5} with κ_{0} = 0.3 (solid); R_{m} = 2 × 10^{3} with κ = 0 (dasheddotted); R_{m} = 2 × 10^{5} with κ_{0} = 0 (diamond+dashed). 
At first we do not consider any smallscale helicity fluxes in the model, i.e., F_{α} = 0 in Eq. (11). The saturation curve for R_{m} = 2 × 10^{3} in Fig. 12 is now underdamped, whereas the dynamo fails to generate a finite B_{sat} for R_{m} = 2 × 10^{5} even though it initially has the same growth rate. On increasing α_{BL} = 10 ms^{1} from 6 ms^{1} the saturation curve for R_{m} = 2 × 10^{5} also displays underdamped behaviour. This indicates that for α_{BL} = 6 ms^{1}, the total α in the domain was simply becoming subcritical and the dynamo was not able to sustain itself through the saturation phase. In Fig. 13a, b, we show the meridional snapshots of the magnetic fields and α_{M} for a saturated state. In this figure, the colour bars indicate that α_{M} increases while B_{φ} decreases with R_{m} for the same value of α_{BL}.
Now we include an advective flux of α_{M} caused by the meridional circulation in Eq. (11). In order to keep the system numerically stable, we also require a diffusive flux in Eq. (11). It is clear from Fig. 12 that the underdamped behaviour in the model without fluxes is suppressed due to a diffusive flux of α_{M}. The production of α_{M} is now countered by a diffusive decay in a time scale much shorter than that given by . The same figure shows that the saturation value of the magnetic energy is now almost independent of R_{m}. With diffusive and advective fluxes owing to the meridional circulation in Eq. (11) the smallscale helicity is distributed through out the convection zone as shown in Fig. 14a, b. It is instructive to compare the magnitudes of α_{M} in Figs. 13 and 14.
The diffusive flux of smallscale helicity is therefore crucial for the operation of a successful mean field αΩ dynamo. An interesting observation is the distribution of α_{M} in the concentrated region at the lower part of the convection zone (see Fig. 13) in contrast to Fig. 14. Even though α_{BL} is a surface phenomena, considerable magnetic helicity is generated in the entire convection zone when the meridional circulation sinks the poloidal field lines at high latitudes and brings them near the tachocline where the toroidal fields are generated.
Fig. 13 Meridional crosssections showing the distribution of toroidal field and α_{M} for a BabcockLeighton dynamo without MC and diffusive helicity fluxes in Eq. (11) for a) R_{m} = 2 × 10^{3} and b) R_{m} = 2 × 10^{5}. The streamlines of the positive and negative poloidal field are shown by solid and dashed lines respectively. Note that the magnetic field has decayed to very small values for R_{m} = 2 × 10^{5}. 
Fig. 14 Meridional crosssections showing the distribution of toroidal field and α_{M} for a BabcockLeighton dynamo with MC and diffusive helicity fluxes for R_{m} = 2 × 10^{3} at two different epochs. The streamlines of the positive and negative poloidal field are shown by solid and dashed lines respectively. 
4. Conclusions
We have performed calculations for αΩ dynamos in a spherical shell for spatially segregated α and Ω source regions. The two classes of models we studied resemble Parker’s interface dynamo and the flux transport dynamo with a BabcockLeighton α.
In agreement with earlier work (Chatterjee et al. 2010), we found that it is not possible to escape the catastrophic α quenching by merely separating the regions of shear and αeffect. The saturation value of magnetic energy in the interface model decreases as for both dynamical quenching with g_{α} = 0 and the algebraic quenching with g_{α} = R_{m} (Fig. 5). Nevertheless, we found that a richer dynamical behaviour emerges for the cases with dynamical α effect, in terms of parity fluctuations (Fig. 8) and the appearance of “secondary” dynamos (Fig. 9), namely dynamo waves governed by the magnetic α effect. We observe this interesting behaviour in two different situations. Firstly, for low R_{m}, the toroidal field is weakly confined below the overshoot layer and therefore α_{M} is generated there by the term in the expansion of the forcing term in Eq. (11). Secondly, for high R_{m}, the secondary dynamos are only observed if there is a supercritical diffusive helicity flux where the α_{M} diffuses into the stable layers. However, because of the lack of scale separation between the mean field and the forcing scale of the helical turbulence we refrain from interpreting this in terms of the poleward migration seen in the Sun. We have to be cautious about using the dynamical quenching equation for dynamo numbers that are not very large compared to the critical dynamo number. For highly supercritical α_{K}, the behaviour of the system begins to be governed by α_{M}. According to our knowledge, there are no direct numerical simulation of secondary dynamos or magnetichelicity driven dynamos applied to stellar magnetic fields. This will be explored in a future paper.
We do not see any evidence for chaotic behaviour in the time series of magnetic energy because the dynamo period and the saturation energy remain fairly constant. This may not be the case for diffusive helicity fluxes, which introduce further complexity to the system. Yet the addition of a diffusive helicity flux relaxes the catastrophic dependence of the saturation magnetic energy (Figs. 10a, 12).
We would expect that the magnetic field affects all turbulent coefficients including α and η. But for this analysis we did not included an equation for the variation of η_{t}. This is justified for the simple two layer model with a lower η_{t} in the region of the production of strong toroidal fields and a higher η_{t} in the region of weaker poloidal fields. Note also that by quenching the diffusivity inversely with magnetic energy in a nonlinear dynamo model, Tobias (1996) was able to produce a bonafide interface model where the magnetic field was restricted to a thin layer at an interface between a layer of shear and cyclonic turbulence. None of the previous interface models have used the dynamical quenching equation.
In the flux transport dynamo model, the magnetic energy fails to reach significant saturation when both the meridional circulation and the diffusive helicity fluxes are artificially turned off in the helicity evolution equation. This is expected and is demonstrated in Fig. 12. It is interesting that the BabcockLeighton dynamos, where α is concentrated only in a narrow layer at the surface, also produce considerable amount of magnetic helicity inside the convection zone when dynamical quenching is employed (Figs. 13, 14).
It remains to explore the role of solar wind and coronal mass ejections, which might help in disposing of smallscale helicity from the Sun and thus alleviate catastrophic quenching. An investigation of the effects of Vishniac & Cho fluxes found them to be of secondary importance compared to diffusive
helicity fluxes for αΩ mean field dynamos (Guerrero et al. 2010). Unfortunately the direct numerical simulations have not yet reached the modest Reynolds numbers used in this paper (~10^{4}), which are still much lower than the astrophysical dynamos. In order to verify if the equation for dynamical quenching works for the αΩ dynamos in the same way as it does in α^{2} dynamos, we need to embark upon systematic comparisons between DNS with shear and convection and mean field modelling.
Acknowledgments
We are grateful to Eric Blackman for comments that helped to improve the paper considerably. This work was supported in part by the European Research Council under the AstroDyn Research Project No. 227952 and the Swedish Research Council Grant No. 62120074064.
References
 Archontis, V., Dorch, S. B. F., & Nordlund, Å. 2003, A&A, 397, 393 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Blackman, E. G., & Brandenburg, A. 2002, ApJ, 579, 359 [Google Scholar]
 Blackman, E. G., & Brandenburg, A. 2003, ApJ, 584, L99 [NASA ADS] [CrossRef] [Google Scholar]
 Blackman, E. G., & Field, G. B. 2000, MNRAS, 318, 724 [NASA ADS] [CrossRef] [Google Scholar]
 Blackman, E. G., & Field, G. B. 2004, Phys. Plasmas, 11, 3264 [NASA ADS] [CrossRef] [Google Scholar]
 Brandenburg, A., & Dobler, W. 2001, A&A, 369, 329 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Brandenburg, A., & Subramanian, K. 2005, AN, 326, 400 [Google Scholar]
 Brandenburg, A., Bigazzi, A., & Subramanian, K. 2001, MNRAS, 325, 685 [NASA ADS] [CrossRef] [Google Scholar]
 Brandenburg, A., Rädler, K.H., & Schrinner, M. 2008, A&A, 482, 739 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Brandenburg, A., Candelaresi, S., & Chatterjee, P. 2009, MNRAS, 398, 1414 [NASA ADS] [CrossRef] [Google Scholar]
 Cattaneo, F., & Hughes, D. W. 1996, PRE, 54, R4532 [Google Scholar]
 Charbonneau, P., 2005, Liv. Rev. Sol. Phys., 2, http://www.livingreviews.org/lrsp20052 [Google Scholar]
 Chatterjee, P., Nandy, D., & Choudhuri, A. R. 2004, A&A, 427, 1019 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Chatterjee, P., Brandenburg, A., & Guerrero, G. 2010, GAFD, in press, [arXiv:1005.5708] [Google Scholar]
 Choudhuri, A. R., Schüssler, M., & Dikpati, M. 1995, A&A, 303, L29 [NASA ADS] [Google Scholar]
 Covas, E., Tworkowski, A., Brandenburg, A., & Tavakol, R. 1997, A&A, 317, 610 [NASA ADS] [Google Scholar]
 Dikpati, M., & Charbonneau, P. 1999, ApJ, 518, 508 [NASA ADS] [CrossRef] [Google Scholar]
 Dikpati, M., & Choudhuri, A. R. 1994, A&A, 291, 975 [NASA ADS] [Google Scholar]
 Durney, B. R. 1995, Sol. Phys., 160, 213 [Google Scholar]
 Gressel, O. 2010, MNRAS, 405, 41 [NASA ADS] [CrossRef] [Google Scholar]
 Guerrero, G., & de Gouveia Dal Pino, E. M. 2008, A&A, 485, 267 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Guerrero, G., Chatterjee, P., & Brandenburg, A. 2010, MNRAS, in press, [arXiv:1005.4818] (GCB) [Google Scholar]
 Hubbard, A., & Brandenburg, A. 2010, GAFD, submitted, [arXiv:1004.4591] [Google Scholar]
 Jepps, S. A. 1975, JFM, 67, 625 [Google Scholar]
 Jouve, L., Brun, A. S., Arlt, R., et al. 2008, A&A, 483, 949 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Käpylä, P. J., Korpi, M. J., & Brandenburg, A. 2008, A&A, 491, 353 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kleeorin, N., & Rogachevskii, I. 1999, PRE, 59, 6724 [Google Scholar]
 Kleeorin, N., Moss, D., Rogachevskii, I., & Sokoloff, D. 2000, A&A, 361, L5 [NASA ADS] [Google Scholar]
 MacGregor, K. B., & Charbonneau, P. 1997, ApJ, 486, 484 [NASA ADS] [CrossRef] [Google Scholar]
 Markiel, J. A., & Thomas, J. H. 1999, ApJ, 523, 827 [NASA ADS] [CrossRef] [Google Scholar]
 Mitra, D., Candelaresi, S., Chatterjee, P., Tavakol, R., & Brandenburg, A. 2010, Astron. Nachr., 331, 130 [Google Scholar]
 Ossendrijver, M., Stix, M., Brandenburg, A., & Rüdiger, G. 2002, A&A, 394, 735 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Parker, E. N. 1955, ApJ, 122, 293 [NASA ADS] [CrossRef] [Google Scholar]
 Parker, E. N. 1993, ApJ, 408, 707 [Google Scholar]
 Petrovay, K., & Kerekes, A. 2004, MNRAS, 351, L59 [NASA ADS] [CrossRef] [Google Scholar]
 Pouquet, A., Frisch, U., & Léorat, J. 1976, JFM, 77, 321 [NASA ADS] [CrossRef] [Google Scholar]
 Shukurov, A., Sokoloff, D., Subramanian, K., & Brandenburg, A. 2006, A&A, 448, L33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Spiegel, E. A., & Zahn, J.P. 1992, A&A, 265, 106 [NASA ADS] [Google Scholar]
 Subramanian, K., & Brandenburg, A. 2006, ApJ, 648, L71 [Google Scholar]
 Sur, S., Shukurov, A., & Subramanian, K. 2007, MNRAS, 377, 874 [NASA ADS] [CrossRef] [Google Scholar]
 Sur, S., Brandenburg, A., & Subramanian, K. 2008, MNRAS, 385, L15 [NASA ADS] [CrossRef] [Google Scholar]
 Tobias, S. M. 1996, ApJ, 467, 870 [NASA ADS] [CrossRef] [Google Scholar]
 van Ballegooijen, A. A., & Choudhuri, A. R. 1988, ApJ, 333, 965 [NASA ADS] [CrossRef] [Google Scholar]
 Vainshtein, S. I., & Cattaneo, F. 1992, ApJ, 393, 165 [NASA ADS] [CrossRef] [Google Scholar]
 Vishniac, E. T., & Cho, J. 2001, ApJ, 550, 752 [NASA ADS] [CrossRef] [Google Scholar]
 Zhang, H., Sokoloff, D., Rogachevskii, I., et al.2006, MNRAS, 365, 276 [NASA ADS] [CrossRef] [Google Scholar]
All Figures
Fig. 1 Profiles of radial shear ∂Ω/∂r, α and η as a function of fractional solar radius. 

In the text 
Fig. 2 Critical α in terms of a fraction of η_{t}k_{f} as a function of magnetic Reynolds number R_{m} for the interface dynamo model of Fig. 1. 

In the text 
Fig. 3 Time evolution of the volumeaveraged magnetic energy in the domain scaled with the equipartition energy for the case with algebraic quenching and g_{α} = 1 for models with R_{m} = 1 (diamond), R_{m} = 20 (solid), R_{m} = 200 (dashed), R_{m} = 2 × 10^{3} and R_{m} = 2 × 10^{5} (triangles). 

In the text 
Fig. 4 Same as Fig. 3 but with dynamical quenching and g_{α} = 0 for the R_{m} indicated in the figure. 

In the text 
Fig. 5 Saturated value of the volumeaveraged magnetic energy scaled with the equipartition energy as a function of R_{m} for models with dynamical α quenching (triangles+solid) and algebraic quenching with g_{α} = 1 (squares + dashed) and with g_{α} = R_{m} (cross + dasheddotted). Note that all runs included in this figure used α_{0} = 2α_{c}. 

In the text 
Fig. 6 Radial profiles of A_{φ} and B_{φ} at two different latitudes (λ) in the saturated phase for R_{m} = 2 × 10^{3}. 

In the text 
Fig. 7 Timelatitude plot of α_{M}(0.72 R_{⊙},θ) for a) R_{m} = 20 and b) R_{m} = 200. 

In the text 
Fig. 8 a) Evolution of parity (purely dipolar = −1 and purely quadrupolar = +1) for an interface dynamo model with R_{m} = 2 × 10^{3} and no fluxes. b) A small part in the timelatitude diagram of B_{φ} indicated by dotted lines in a) where the parity is changing from quadrupolar to dipolar. 

In the text 
Fig. 9 Timelatitude plot of the toroidal field a) and c) and α_{M}b) and d) for interface dynamo models with α = 4α_{c} and R_{m} = 20. 

In the text 
Fig. 10 a) Time evolution of the volume averaged magnetic energy in the domain scaled with the equipartition energy for models with R_{m} = 2 × 10^{3} and κ_{0} = 10^{5} (∗ ) as well as with R_{m} = 2 × 10^{5} and κ_{0} = 10^{5} (solid line) for the interface dynamo model of Sect. 3.1 with dynamical α quenching. The saturation curve for zero fluxes for the model with R_{m} = 2 × 10^{5} has been shown by the dashed line; b) and c) show timelatitude diagrams for the toroidal field at the depths indicated for models with R_{m} = 2 × 10^{3} and κ_{0} = 10^{5}. 

In the text 
Fig. 11 The timelatitude plot of toroidal field a) and c), and α_{M}b) and d) with α = 2α_{c} for R_{m} = 2 × 10^{5} and κ_{0} = 0.01. Meridional snapshots of e) sign(B_{φ})( B_{φ}  /B_{eq})^{1/2} and f) α_{m} × 10^{3} for the same case. 

In the text 
Fig. 12 Time evolution of the volume averaged magnetic energy for the flux transport dynamo model of Sect. 3.2 for R_{m} = 2 × 10^{3} with κ_{0} = 0.3 (dashed); R_{m} = 2 × 10^{5} with κ_{0} = 0.3 (solid); R_{m} = 2 × 10^{3} with κ = 0 (dasheddotted); R_{m} = 2 × 10^{5} with κ_{0} = 0 (diamond+dashed). 

In the text 
Fig. 13 Meridional crosssections showing the distribution of toroidal field and α_{M} for a BabcockLeighton dynamo without MC and diffusive helicity fluxes in Eq. (11) for a) R_{m} = 2 × 10^{3} and b) R_{m} = 2 × 10^{5}. The streamlines of the positive and negative poloidal field are shown by solid and dashed lines respectively. Note that the magnetic field has decayed to very small values for R_{m} = 2 × 10^{5}. 

In the text 
Fig. 14 Meridional crosssections showing the distribution of toroidal field and α_{M} for a BabcockLeighton dynamo with MC and diffusive helicity fluxes for R_{m} = 2 × 10^{3} at two different epochs. The streamlines of the positive and negative poloidal field are shown by solid and dashed lines respectively. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.