# Multiple forms of intermittency in PDE dynamo models

## I Introduction

Intermittency has been observed in a variety of real settings as well in a vast number of numerical models. A great deal of effort has therefore gone into understanding these modes of behaviour in the context of deterministic dynamical systems theory. These studies have demonstrated the existence of a number of different types of intermittency (such as Pomeau–Manneville , Crisis , On-off  intermittencies), each with their own associated signatures and scalings. Many of these forms of intermittency have in turn been concretely shown to be present in experiments and numerical studies of dynamical systems in a variety of settings (see  and references therein).

An important potential domain of applicability of such behaviour arises in understanding the mechanisms underlying the intermediate time scale variability in the Sun  \- the occurrence of the so called Maunder or grand minima \- during which solar activity (as deduced from the sunspot numbers) was greatly diminished . This behaviour is also confirmed by evidence coming from the analysis of proxy data . There is also some evidence for similar types of variability in solar-type stars .

The idea that some type of dynamical intermittency may under pin the grand minima type variability in the sunspot record (the intermittency hypothesis ) goes back at least to the late 1970’s . This idea has been the subject of intense study over the recent years and has involved the employment of various classes of dynamo models, including ordinary differential equations (ODE) (e.g.  as well as partial differential equations (PDE) models (e.g. ). In addition to the phenomenological evidence for the presence of intermittent-type behaviours in dynamo models , concrete evidence has recently been found for the presence of particular types of intermittency in both ODE dynamo models  as well as a recently discovered generalisation of on-off intermittency, referred to as in-out intermittency , in PDE models .

Here we wish to report concrete evidence for the occurrence of two other types of intermittency, namely the crisis–induced and Pomeau–Manneville Type-I intermittencies, in PDE mean–field dynamo models. The organisation of the paper is as follows. In Sec. II we briefly introduce the model studied here. Sec. III summarises our evidence demonstrating the presence of these types of intermittencies in this model and finally in Sec. IV we draw our conclusions.

## II Model

Ideally one would wish to employ the full 3-dimensional dynamo models with a minimum number of approximations and simplifying assumptions. Despite a number of important recent attempts , the difficulty of dealing with small scale turbulence makes a detailed and extensive self consistent study of such fully turbulent regimes in stars still computationally impractical (see e.g. .

In view of this an alternative approach in studies of stellar dynamos has been to employ mean–field models . We should mention that there is an ongoing debate regarding the nature and realistic value of such models . Nevertheless, 3-D turbulence simulations do seem to produce magnetic fields whose global properties (such as field parity and time dependence) are similar to those expected from corresponding mean–field dynamo models . In this way mean–field dynamo models seem to reproduce certain features of the more complicated models and allow the study of certain global properties of magnetic fields in the Sun and solar-type stars (see for example ).

The standard mean–field dynamo equation is given by

$$\frac{𝐁}{t}=\times \left(𝐮\times 𝐁+\alpha 𝐁\eta _t\times 𝐁\right),$$
(1)   
where $`𝐁`$ and $`𝐮`$ are the mean magnetic field and mean velocity respectively and the turbulent magnetic diffusivity $`\eta _t`$ and the coefficient $`\alpha `$ arise from the correlation of small scale turbulent velocities and magnetic fields ($`\alpha `$ effect) . We consider the usual algebraic form of $`\alpha `$–quenching namely

$$\alpha =\frac{\alpha _0\mathrm{cos}\theta }{1+|𝐁|^2},$$
(2)   
where $`\alpha _0=\mathrm{constant}`$ and $`\theta `$ is the co-latitude.

We solve Eq. 1 in an axisymmetric configuration and in the following, as is customary , we shall discuss the behaviour of the solutions by monitoring the total magnetic energy, $`E=\frac{1}{2\mu _0}𝐁^2𝑑V`$, where $`\mu _0`$ the induction constant, and the integral is taken over the dynamo region. We split $`E`$ into two parts, $`E=E_A+E_S`$, where $`E_A`$ and $`E_S`$ are respectively the energies of the antisymmetric and symmetric parts of the field with respect to the equator. The overall parity $`P`$ is given by $`P=[E_SE_A]/E`$, so $`P=1`$ denotes an antisymmetric (dipole-like) pure parity solution and $`P=+1`$ a symmetric (quadrupole-like) pure parity solution.

For the numerical results reported in the following section, we used a modified version of the axisymmetric dynamo code of Brandenburg et al. (1989)  employed recently in . These models are constructed from a complete sphere of radius $`R`$ by removing an inner concentric sphere of radius $`r_0`$ and a conical section of semi-angle $`\theta _0`$ about the rotation axis, from both the north and south polar regions (see  for details of the model and the relevant parameters). To test the robustness of the code we verified that no qualitative changes were produced by employing a finer grid and different temporal step length (we used a grid size of $`41\times \mathrm{\hspace{0.17em}81}`$ mesh points and a step length of $`10^4R^2/\eta _t`$ in the results presented in this paper). For the following results we use $`C_\mathrm{\Omega }=10^4`$, which give the magnitude of the differential rotation and $`\theta _0=45^{}`$. The magnitude of the $`\alpha `$-effect is given by the dynamo parameter $`C_\alpha `$.

In the next section we show in turn concrete evidence for the occurrence of crisis–induced and Pomeau–Manneville Type-I intermittencies.

## III Results

### A Crisis–induced Intermittency

As far as their detailed underlying mechanisms and temporal signatures are concerned, crises come in three varieties . Here we shall be concerned with only one of these types, referred to as “attractor merging crisis”, whereby as a system parameter is varied, two or more chaotic attractors merge to form a single attractor. There are both experimental and numerical evidence for this type of intermittency (see for example  and references therein). In particular, this type of behaviour has been discovered in a 6-dimensional truncation of mean–field dynamo models .

Fig. 1 shows the plots of the energy and parity for the above model as a function of time, calculated with $`r_0=0.2`$ and $`C_\alpha =25.202`$ which show a bimodal behaviour, switching intermittently between two different chaotic states.

To determine the nature of this behaviour more precisely, we have plotted in Fig. 2 the return maps for the PDE models (1), showing the attractors before and after the merging. As can be seen the resulting merged attractor is, as expected, larger than the superposition of the two pre-existing attractors.

These results can be taken as indications for the presence of crisis–induced intermittency in this model. To substantiate this further, we recall that another important signature of this type of intermittency is the way $`\tau `$, the average time between switches, scales with the system parameter, in this case, $`C_\alpha `$. According to Grebogi et al. , for a large class of dynamical systems this relation takes the form

$$\tau \left|C_\alpha C_\alpha ^{}\right|^\gamma ,$$
(3)   
where the real constant $`\gamma `$ is the critical exponent characteristic of the system under consideration and $`C_\alpha ^{}`$ is the critical value of $`C_\alpha `$ at which the two chaotic attractors merge.

The model under study here is a PDE system which is formally infinite dimensional. Such PDE models are numerically costly to integrate over long enough intervals of time (sometimes in excess of 5000 time units) necessary in order to obtain the scaling of the type (3). Furthermore, the demonstration of such scaling requires a precise determination of the critical value $`C_\alpha ^{}`$ which is difficult since as one approaches this value $`\tau `$ diverges and the integration time becomes prohibitive. Despite these difficulties, we have succeeded to obtain strong evidence for the presence of such a scaling as depicted in Fig. 3, with the corresponding $`\gamma =1.08\pm 0.05`$. Grebogi et al.  conjecture that there may be a general tendency for $`\gamma `$ to be larger for higher–dimensional attractors. We do have a value of $`\gamma `$ higher than the previous one found for a related six dimensional ODE dynamo model  but much lower than the value range suggested by Grebogi et al. Therefore, the conjectured range may need modification for large high–dimensional systems.

There is also evidence for an enlargement of the final attractor after merging, as shown by the larger amplitudes of variation in the parity, in the sense that the parity gets closer to $`1`$ after the merging, as depicted in Fig. 2. This helped us to numerically arrive at a better estimate for the critical value $`C_\alpha ^{}`$.

These indicators, taken together, amount to strong evidence for the presence of crisis–induced intermittency for this model.

### B Pomeau–Manneville Type-I Intermittency

This type of intermittency, which is brought about through a tangent bifurcation, results in the system switching back and forth between a “ghost” periodic orbit and sudden bursts of chaotic behaviour . There are both experimental and numerical evidence for this type of intermittency (see for example  and references therein). In particular this type of behaviour has been discovered in a 12-dimensional truncation of mean–field dynamo model .

To demonstrate the presence of this type of intermittency in the above PDE dynamo model, we have plotted in Fig. 4 the energy and parity as a function of time for the parameter values $`r_0=0.7`$ and $`C_\alpha =28.0`$, which clearly demonstrates switches between nearly periodic behaviour and sudden bursts. We note that interestingly the energy in this case shows strong modulation which could be of interest in accounting for the occurrence of grand type minima in sunspot activity.

Another signature of this type of intermittency is provided by the specific characteristics of its corresponding power spectrum. By employing finite dimensional maps , it has been shown that the corresponding spectra have a broad-band feature whose shape obeys approximately the inverse-power law $`1/f`$ for $`f>f_s`$, where $`f_s`$ is the saturation frequency. Below this frequency there is a flat plateau induced by noise that causes arbitrarily long laminar phases to become finite.

As further evidence for this type of intermittency in the model (1), we have plotted in Fig. 5 the power spectrum at $`C_\alpha =28.0`$, obtained by averaging over 16 different initial conditions corresponding to different initial parities. As can be seen, the power spectrum shows both the flat plateau and the $`1/f`$ power law scaling.

Taken together, these indicators amount to strong evidence for the presence of Pomeau–Manneville Type-I intermittency for this model.

## IV Conclusion

We have obtained concrete evidence, in terms of phase space signatures, spectra and scalings to demonstrate the presence of crisis–induced and the Pomeau–Manneville Type-I intermittencies in axisymmetric mean–field PDE dynamo models. Despite the rather idealised nature of these models, this is of potential importance since it shows the occurrence of two more types of intermittency (in addition to in–out intermittency recently discovered ) in these models which may in turn be taken as an indication that more than one type of intermittency may occur in solar and stellar dynamos. This suggests that any observational programme for identifying the mechanisms underlying grand minima type variability needs to take into account the possibility that multiple intermittency mechanisms may be operative in different stars of the similar type, or even in the same star over different epochs. This would also be of importance in the interpretation of proxy data. In this way a more appropriate hypothesis regarding such variability would be that of multiple–intermittency hypothesis.   

We would like to thank Axel Brandenburg for providing us with the original code and Andrew Tworkowski for helpful discussions. EC is supported by grant BD/5708/95 – PRAXIS XXI, JNICT. RT benefited from PPARC UK Grant No. L39094.