text:
# Random versus deterministic exponents in a rich family of diffeomorphisms

## 1 Introduction

Numerical experiments with area-preserving surface diffeomorphisms often produce the following dynamical picture: elliptical islands floating in ergodic seas. A reasonable guess is that these ergodic seas typically have positive measure, and further, that the Lyapunov exponents on these seas are on average nonzero. An example of tiny elliptical islands in the context of differential equations can be found in \[NSV\], where all rough numeric tests are in favor of ergodic behavior.

In this paper, we add to the pile of experimental evidence in favor of this conjecture.111Due to space limitations, this version of the paper does not contain figures. A version with the figures is available at either http$`://`$www.math.northwestern.edu$`/_\stackrel{~}{}`$wilkinso$`/`$papers ($`8\frac{1}{2}\times 11`$ format) or http$`://`$www.maia.ub.es$`/`$dsg$`/`$2002$`/`$index.html (A4 format).  We also discuss a possible theoretical approach to finding positive Lyapunov exponents in certain families of area-preserving diffeomorphisms of the sphere $`S^2`$. The possibility of such an approach was discussed in \[BuPuShWi\]. The families we consider are not obtained from a specific set of equations, but from the following construction. Let $`f:S^2S^2`$ be an area-preserving diffeomorphism of the round sphere, and let $`SO(3)`$ be the isometry group of $`S^2`$. Let

$$=\{gfgSO(3)\}$$

be the left $`SO(3)`$-coset of $`f`$ in $`\text{Diff }(S^2)`$, and let $`\nu `$ be the push-forward of Haar measure on $`SO(3)`$ to $``$. Provided that $`f`$ is not itself an isometry, the family $``$ has nonzero random Lyapunov exponents with respect to $`\nu `$ (see Proposition 2.2 below). The question this paper addresses is whether these random exponents can somehow be connected to the Lyapunov exponents of individual members of $``$, at least on $`\nu `$-average.

To test whether there might be such a connection, we chose $`f`$ to be a twist map, all of whose Lyapunov exponents are zero. The resulting family $``$ has similarities to the standard family on the $`2`$-torus. The dynamics of the individual elements of $``$ and how they depend on parameters is an interesting topic, but we only study here some key properties in the case of small $`\epsilon `$. We mainly focus on two quantities, the random exponent $`R(\nu )`$ and the average exponent $`\mathrm{\Lambda }(\nu )`$, which we now define.

Let $`\mu `$ be Lebesgue measure on $`S^2`$ normalized to be a probability measure. Suppose for now that $`\nu `$ is an arbitrary Borel probability measure supported on a subset $``$ of $`\text{Diff }_\mu (S^2)`$, the space of $`\mu `$-preserving diffeomorphisms of $`S^2`$. For $`f`$, the largest Lyapunov of $`f`$ at $`xS^2`$ is found by computing the limit:

$$\underset{n\mathrm{}}{lim}\frac{1}{n}\mathrm{log}T_xf^n=\lambda _1(x,f),$$
(1)
which exists for $`\mu `$-almost every $`x`$ by the subadditive ergodic theorem. We define the average exponent of $`f`$ to be

$$\lambda (f)=_{S^2}\lambda _1(f,x)𝑑\mu (x),$$
(2)
and the average exponent of $`\nu `$ to be

$$\mathrm{\Lambda }(\nu )=_{\text{Diff }_\mu ^r(S^2)}\lambda (f)𝑑\nu (f).$$
(3)
Rather than iterate a single diffeomorphism $`f`$, we might choose instead a sequence of diffeomorphisms $`\{f_1,f_2,\mathrm{}\}`$ and form their composition:

$$f^{(n)}:=f_nf_{n1}\mathrm{}f_1.$$

If the sequence is chosen to be independent and identically distributed with respect to $`\nu `$, then almost surely the limit

$$\underset{n\mathrm{}}{lim}\frac{1}{n}\mathrm{log}T_xf^{(n)}=:R(x,(f_i)_1^{\mathrm{}},\nu )$$
(4)
will exist, for $`\mu `$-almost every $`x`$. (This too follows from the subadditive ergodic theorem, applied in the appropriate context). Further, the integral of $`R(x,(f_i)_1^{\mathrm{}},\nu )`$ with respect to $`\mu `$ is almost surely independent of the sequence $`(f_i)_1^{\mathrm{}}`$. We define the random exponent of $`\nu `$ to be this integral:

$$R(\nu )=_{S^2}R(x,(f_i)_1^{\mathrm{}},\nu )𝑑\mu (x).$$
(5)
(see Kifer for an introduction to the subject of random diffeomorphisms and their exponents. We also give a self-contained introduction in section 2). The random exponent $`R(\nu )`$ is usually positive, unless $`\nu `$ is fairly degenerate \[Ca\].

The quantity $`\mathrm{\Lambda }(\nu )`$ is mysterious from a computational perspective, but useful from a dynamical one. The quantity $`R(\nu )`$ is relatively easy to estimate and is often positive.

Our goal is to understand in general if there is a notion of richness for a probability measure $`\nu `$ on the volume preserving diffeomorphisms of a closed manifold $`M`$ such that the positivity of $`R(\nu )`$ implies the positivity of $`\mathrm{\Lambda }(\nu )`$. Here we are investigating whether the $`SO(3)`$ invariance of the measure $`\nu `$ on $`\text{Diff }_\mu (S^2)`$ might provide such a notion of richness. In \[BuPuShWi\] we asked if even more might hold, that we might bound $`\mathrm{\Lambda }(\nu )`$ from below in terms of $`R(\nu ).`$

###### Question 1.1

Is there a positive constant $`C`$ — perhaps $`1`$ — such that $`\mathrm{\Lambda }(\nu )CR(\nu )`$?

Some motivation for Question 1.1 can be found in similar question for the iterates of linear maps (see \[DeSh\] where an affirmative answer to the analogue of Question 1.1 is proven with $`C=1`$). In Section 3, we describe a theoretical framework in which to address Question 1.1 and related questions. We discuss the linear case in Section 4.

Returning to the specific family of diffeomorphisms mentioned earlier, we now describe the experiment in more detail.

For $`\epsilon >0`$, we define a one-parameter family of twist maps $`f_\epsilon `$ as follows. Express $`S^2`$ as the sphere of radius $`1/2`$ centered at $`(0,0)`$ in $`𝐑\times 𝐂`$, so that the coordinates $`(r,z)S^2`$ satisfy the equation

$$|r|^2+|z|^2=1/4.$$

In these coordinates define a twist map $`f_\epsilon :S^2S^2`$, for $`\epsilon >0`$, by

$$f_\epsilon (r,z)=(r,\text{exp}(2\pi i(r+1/2)\epsilon )z).$$

Let $`_\epsilon `$ be the orbit $`SO(3)f_\epsilon `$. Let $`\nu `$ be the push-forward of Haar measure on $`SO(3)`$. We denote the resulting random and average Lyapunov exponents by $`R(\epsilon )`$ and $`\mathrm{\Lambda }(\epsilon )`$, respectively.

The numerical results are described in Section 6. To summarize these results, it appears that the inequality $`\mathrm{\Lambda }(\epsilon )R(\epsilon )`$ is satisfied for large $`\epsilon `$, and it is definitely not satisfied for small $`\epsilon `$. We now have rigorous results to confirm some of these observations. The strongest of these results is proved in Section 7: for $`\epsilon `$ close to $`0`$, there is no $`C>0`$ satisfying the inequality in Question 1.1: in fact, we show in Corollary 7.7 that for small $`\epsilon `$, $`\mathrm{\Lambda }(\epsilon )`$ is less than $`\epsilon ^3R(\epsilon )^{3/2}`$. The numerics support an upper bound on $`\mathrm{\Lambda }(\epsilon )`$ that is exponentially small, and we show in Theorem 7.6 that on most of $`_\epsilon `$ this is indeed the case. On the other hand, $`\mathrm{\Lambda }(\epsilon )`$ does appear to be positive for positive $`\epsilon `$ in all of the range where we can meaningfully compute. Section 7 also contains a study of the bifurcation structure of fixed points inside the family $`_\epsilon `$. We include this analysis because it sheds light on where, in both $`_\epsilon `$ and $`S^2`$, new elliptic periodic points and their surrounding islands are produced. Homoclinic bifurcations also give rise to horseshoes, which are associated (at least heuristically) with positive measure sets with nonzero exponents.

For the case of large $`\epsilon `$, we show in Section 8, that an inequality like that in Question 1.1 is satisfied when a small amount of noise is introduced. More precisely, we first prove in Section 3 some results about the quantities $`R(\epsilon )`$ and $`\mathrm{\Lambda }(\epsilon )`$ and a third quantity $`R(\epsilon ,\delta )`$, which measures the exponents of the “in-between” process in which each element of $`_\epsilon `$ has added noise in a $`\delta `$-ball inside $`_\epsilon `$ (see Section 3 for details). In particular, we prove that any element $`h`$ of a family $`=SO(3)f`$ described above will have average $`\delta `$-diffused exponents $`R(h,\nu ,\delta )`$ that are positive, unless $`f`$ is an isometry. In addition, there exists a stationary measure $`m_{h,\delta }`$ for such a process on the projective bundle $`PS^2`$ that is absolutely continuous with respect to Liouville measure and projects to Lebesgue measure $`\mu `$ on $`S^2`$; this measure is unique among stationary measures with these properties and is the unique fixed point of a “simple” linear operator. The integrated measure

$$m_\delta =_{}m_{h,\delta }𝑑\nu (h)$$

determines $`R(\nu ,\delta )`$, which is the average of the $`\delta `$-diffused exponents $`R(h,\nu ,\delta )`$ over $`h`$. Whenever $`m_\delta `$ is equal to Lebesgue measure $`m`$, we have the equality: $`R(\nu ,\delta )=\mathrm{\Lambda }(\nu )`$. For the family $`_\epsilon `$, denote by $`m_{\epsilon ,\delta }`$ this integrated measure. In Section 8, we prove (Theorem 8.1) that for any $`\delta >0`$, $`lim_\epsilon \mathrm{}m_{\epsilon ,\delta }=m`$. Using this result, we prove that if enough noise is introduced, then the inequality in Question 1.1 is satisfied as $`\epsilon \mathrm{}`$; in particular, we show in Corollary 8.2 that $`R(\epsilon ,\delta )R(\epsilon )`$ tends to $`0`$ as $`\delta 0`$ and $`\epsilon \mathrm{}`$ sufficiently quickly, for instance $`\epsilon >\delta ^{25}`$. Results of a similar nature for the standard family were obtained by Carleson and Spencer \[CaSp\] and are described in Section 8 below.

We suspect that a further study of these measures $`m_\delta `$ would be interesting. Even for an $`SO(n)`$\- or $`SU(n)`$-invariant family of matrices, the properties of the analogous “in-between” measures $`m_\delta `$ are, for the most part, unknown. In Section 4, we discuss what is known about these measures. For $`SO(2)`$, we prove that $`m_\delta `$ is Lebesgue measure for all $`\delta >0`$, and for $`SU(n)`$, $`m_\delta `$ is not Lebesgue measure if $`\delta `$ is sufficiently small.

We thank Roy Adler, Victor Klepstyn, Yuri Kifer, Marco Martens, Lai-Sang Young and especially Tolya Katok for conversations about this work. The computing facilities of the Dynamical Systems Group of the University of Barcelona have been widely used. We thank the supporting institutions. Michael Shub was partly supported by NSF Grant $`\mathrm{\#}`$DMS-9988809. Carles Simó was partly supported by grants DGICYT BFM2000-805, CIRIT 2001SGR-70 and INTAS00-221. Amie Wilkinson was partly supported by NSF Grant $`\mathrm{\#}`$DMS-0100314.

## 2 Background on random transformations and exponents

In this section we introduce some notation and gather together some facts and propositions. What we have to say in this section and the next is standard and can be found for example in \[Ki1\], \[Ki2\], \[Ki3\], \[Li-Qi\], \[Ca\], \[GoMa\] in most cases in greater generality. We have outlined proofs here in order to be reasonably self contained.

If $`\text{Diff }(M)`$, and $`\nu `$ is a probability measure on $``$, then $`(,\nu )`$ generates a random process given by selecting an independent, $`\nu `$-distributed sequence $`(h_i)_1^{\mathrm{}}`$ and forming the compositions:

$$h^{(n)}=h_nh_{n1}\mathrm{}h_1.$$

To study all possible outcomes of this experiment, we introduce the following auxiliary spaces and transformations: the shift space, $`^{\mathrm{}}:=\mathrm{\Pi }_{j=1}^{\mathrm{}}`$, the one-sided shift $`\sigma :^{\mathrm{}}`$ given by:

$$\sigma (h_1,h_2,\mathrm{})=(h_2,h_3,\mathrm{}),$$

and the skew product $`\tau :^{\mathrm{}}\times M`$ given by:

$$\tau ((h_i)_1^{\mathrm{}},x)=(\sigma ((h_i)_1^{\mathrm{}}),h_1(x)).$$

Then $`\sigma `$ has a natural invariant measure $`\nu ^{\mathrm{}}`$, the product measure induced by $`\nu `$, but a priori $`\tau `$ has no preferred invariant measure.

###### Definition 2.1

Let $`\nu `$ be a probability measure on $`\text{Diff }(M)`$. A measure $`\mu `$ on $`M`$ is stationary for the random process given by $`(,\nu )`$ if any of the following equivalent conditions is satisfied:

1. $`\tau _{}(\nu ^{\mathrm{}}\times \mu )=\nu ^{\mathrm{}}\times \mu `$
2. $`ev_{}(\nu \times \mu )=\mu `$, where $`ev:\times MM`$ is the evaluation map:

$$ev(h,x)=h(x)$$
3. $`\mu \nu =\mu `$, where $``$ is the convolution operator defined by:

$$\mu \nu (A)=_{}\mu (h^1(A))𝑑\nu (h),$$

for every $`\mu `$-measurable subset $`AM`$.

Stationary measures always exist \[Ki1\] and are the random analogue of invariant measures in the nonrandom setting. Part of the focus of this paper is to find natural stationary measures in the case where $`M=T_1S^2`$ and $``$ and $`\nu `$ are derived from Haar measure on $`SO(3)`$.

Given an injective linear map $`A:VW`$ between normed vector spaces we denote by $`A_{\mathrm{}}`$ the induced map from the unit sphere in $`V`$ to the unit sphere in $`W`$, which is defined by $`v\frac{A(v)}{A(v)}.`$ We use the same notation for the induced map on the projective space $`PV`$. We denote the tangent bundle of $`S^2`$ by $`TS^2`$, the unit tangent bundle by $`T_1S^2`$, and the projective bundle by $`PS^2`$. We let $`m`$ denote the normalized Liouville measure on $`PS^2`$, so $`m`$ is a probability measure which pushes forward under projection to $`S^2`$ to $`\mu `$. The fibers of $`TS^2`$ and $`PS^2`$ over a point $`zS^2`$ are denoted by $`T_zS^2`$ and $`P_zS^2`$. For any manifolds $`M,N`$ and differentiable map $`F:MN`$ the derivative of $`F`$ at $`xM`$ is denoted by $`T_xF`$; for $`vT_xM`$ we will usually write “$`TFv`$” instead of $`T_xF(v)`$. Finally, we denote by $`F_{\mathrm{}}:T_1MT_1N`$ the map that covers $`F`$ and is $`(T_xF)_{\mathrm{}}`$ on the fiber over $`xM`$. Since the tangent map to $`gSO(3)`$ preserves unit tangent vectors, we will write “$`g`$” for $`g_{\mathrm{}}`$.

Now let $`f\text{Diff }_\mu (S^2)`$, let $`=\{gfgSO(3)\}`$, and let $`\nu `$ be the push-forward to $``$ of Haar measure on $`SO(3)`$. Let $`m`$ be normalized Liouville measure on $`PS^2`$. Associated to $``$ we then have the set

$$_{\mathrm{}}=\{h_{\mathrm{}}|h\}=\{gf_{\mathrm{}}gSO(3)\},$$

and the measure $`\nu _{\mathrm{}}`$, the push-forward to $`\text{Diff }(PS^2)`$ of Haar measure on $`SO(3)`$. Let $`\sigma :^{\mathrm{}}`$, $`\tau :^{\mathrm{}}\times S^2`$, $`\sigma _{\mathrm{}}:_{\mathrm{}}^{}{}_{}{}^{\mathrm{}}`$, and $`\tau _{\mathrm{}}:_{\mathrm{}}^{}{}_{}{}^{\mathrm{}}\times PS^2`$ be the associated auxiliary transformations to the random processes generated by $`(,\nu )`$ and $`(_{\mathrm{}},\nu _{\mathrm{}})`$ respectively.

###### Lemma 2.1

The measures $`\mu `$ and $`m`$ are stationary for $`\nu `$ and $`\nu _{\mathrm{}}`$ respectively.

The transformations $`\tau `$, $`\sigma `$ and $`\sigma _{\mathrm{}}`$ are ergodic with respect to $`\nu ^{\mathrm{}}\times \mu `$, $`\nu ^{\mathrm{}}`$, and $`\nu _{\mathrm{}}^{\mathrm{}}`$, respectively.

Proof:  It is straightforward to check that they are stationary.

Ergodicity is not much harder to check. QED

Now we can compute $`R(\nu )`$ more explicitly:

###### Proposition 2.2

Let $`f\text{Diff }_\mu (S^2)`$, let $`=\{gfgSO(3)\}`$, and let $`\nu `$ be the push-forward to $``$ of Haar measure on $`SO(3)`$. Let $`m`$ be normalized Liouville measure on $`PS^2`$, the projective bundle of $`S^2`$. Then

$$R(\nu )=_{PS^2}\mathrm{log}Tfvdm(v).$$

Moreover, $`R(\nu )>0`$, unless $`f`$ is an isometry.

Proof:  We first apply Birkhoff’s Ergodic Theorem to the measure-preserving transformation $`\tau _{\mathrm{}}:_{\mathrm{}}^{\mathrm{}}\times PS^2`$ and the function

$$\psi ((h_i\mathrm{})_1^{\mathrm{}},v)=\mathrm{log}Tfv$$

to obtain that

$`\underset{p\mathrm{}}{lim}{\displaystyle \frac{1}{p}}\mathrm{log}Th_p\mathrm{}Th_1(v)`$ $`=`$ $`\underset{p\mathrm{}}{lim}{\displaystyle \frac{1}{p}}{\displaystyle \underset{j=1}{\overset{p}{}}}\mathrm{log}Tf(h_j\mathrm{}\mathrm{}h_1\mathrm{}(v))`$
$`=`$ $`\underset{p\mathrm{}}{lim}{\displaystyle \frac{1}{p}}{\displaystyle \underset{j=0}{\overset{p1}{}}}\psi (\tau _{\mathrm{}}^j((h_i\mathrm{})_1^{\mathrm{}},v))`$
$`=:`$ $`L((h_i\mathrm{})_1^{\mathrm{}},v)`$
exists a.e. in $`_{\mathrm{}}^{\mathrm{}}\times PS^2`$. The integral of this limit $`L`$ is

$`{\displaystyle L((h_i\mathrm{})_1^{\mathrm{}},v)d(\nu _{\mathrm{}}^{\mathrm{}}\times m)}`$ $`=`$ $`{\displaystyle _{_{\mathrm{}}^{\mathrm{}}\times PS^2}}\psi d(\nu _{\mathrm{}}^{\mathrm{}}\times m)`$
$`=`$ $`{\displaystyle _{PS^2}}\mathrm{log}Tfvdm(v).`$
Next, we apply Oseledec’s theorem to the map $`\tau :^{\mathrm{}}\times S^2`$ and the cocycle $`((h_i)_1^{\mathrm{}},x)T_xh_1`$. We obtain that for almost all $`(h_i)_1^{\mathrm{}}`$, almost all $`xS^2`$, and for almost all $`vT_{1,x}S^2`$, the limit

$$K((h_i)_1^{\mathrm{}},x):=\underset{p\mathrm{}}{lim}\frac{1}{p}\mathrm{log}Th_p\mathrm{}Th_1(v)$$

exists, is independent of $`v`$, and has $`\nu ^{\mathrm{}}\times \mu `$-integral equal to $`R(\nu )`$. The function:

$`K((h_i)_1^{\mathrm{}})`$ $`:=`$ $`{\displaystyle _{S^2}}K((h_i)_1^{\mathrm{}},x)𝑑\mu (x)`$
$`=`$ $`{\displaystyle _{xS^2}}{\displaystyle _{uT_{1,x}S^2}}\underset{p\mathrm{}}{lim}{\displaystyle \frac{1}{p}}\mathrm{log}Th_p\mathrm{}T_xh_1(u)dud\mu (x)`$
$`=`$ $`{\displaystyle _{PS^2}}\underset{p\mathrm{}}{lim}{\displaystyle \frac{1}{p}}\mathrm{log}Th_p\mathrm{}Th_1(v)dm(v)`$
is $`\sigma `$-invariant and has $`\nu ^{\mathrm{}}`$-integral equal to $`R(\nu )`$; ergodicity of $`\sigma `$ implies that it is a.e. constant and therefore equal to $`R(\nu )`$. We conclude that

$`R(\nu )`$ $`=`$ $`{\displaystyle _{PS^2}}\mathrm{log}Tf(v)dm(v)`$
It remains to see that $`_{PS^2}\mathrm{log}Tfvdm(v)>0`$ if $`f`$ is not an isometry. For this we use the following elementary lemma.

###### Lemma 2.3

Let $`ASL(2,𝐑)`$. Then the Jacobian of $`A_{\mathrm{}}`$ with respect to Lebesgue measure on $`S^1`$ is given by:

$$\text{Jac}(A_{\mathrm{}})(v)=Av^2,$$

for $`vS^1`$.

Since $`f`$ preserves $`\mu `$, it follows from this lemma that the Jacobian of $`f_{\mathrm{}}`$ with respect to $`m`$ at $`vT_1S^2`$ is $`Tf(v)^2`$. Since $`Tf_{\mathrm{}}`$ is a diffeomorphism,

$`{\displaystyle _{PS^2}}Tf(v)^2𝑑m(v)`$ $`=`$ $`{\displaystyle _{PS^2}}\text{Jac}(f_{\mathrm{}})(v)`$
$`=`$ $`1.`$
By Jensen’s inequality,

$`{\displaystyle _{PS^2}}\mathrm{log}Tf(v)^2dm(v)`$ $``$ $`\mathrm{log}({\displaystyle _{PS^2}}Tf(v)^2𝑑m(v))`$
$`=`$ $`0`$
with inequality holding unless $`\mathrm{log}Tfv`$ is constant and equal to $`0`$. Rearranging the inequality, we see that, unless $`f`$ is an isometry, we must have $`_{PS^2}\mathrm{log}Tf(v)dm(v)>0`$. QED

## 3 A theoretical framework

### 3.1 Connecting $`R(\epsilon )`$ to $`\mathrm{\Lambda }(\epsilon )`$

In this section, we attempt to interpolate between $`R(\epsilon )`$ and $`\mathrm{\Lambda }(\epsilon )`$ via a third quantity, $`R(\epsilon ,\delta )`$, which we call the random $`\delta `$-diffused exponent. When $`\delta `$ is greater than or equal to the radius of $`SO(3)`$, $`R(\epsilon ,\delta )`$ is equal to $`R(\epsilon )`$; as $`\delta `$ approaches $`0`$, $`R(\epsilon ,\delta )`$ approaches (in limsup) a lower bound for $`\mathrm{\Lambda }(\epsilon )`$. Roughly speaking, $`R(\epsilon ,\delta )`$ is the exponent (averaged over $`_\epsilon `$) obtained by introducing random perturbations (viewed as noise) of order $`\delta `$ to each element of $`_\epsilon `$, staying within the family $`_\epsilon `$. In Lemma 3.2 we show that $`lim\; sup_{\delta 0}R(\epsilon ,\delta )\mathrm{\Lambda }(\epsilon )`$. On the other hand, we derive in Proposition 3.3 a formula for $`R(\epsilon ,\delta )`$:

$`R(\epsilon ,\delta )`$ $`=`$ $`{\displaystyle _{PS^2}}\mathrm{log}Tfvdm_{\epsilon ,\delta }(v).`$
The probability measure $`m_{\epsilon ,\delta }`$ in this formula has nice properties: it projects to Lebesgue measure $`\mu `$ on $`S^2`$, and is absolutely continuous with respect to $`m`$, with smooth density.

Now, recall (Proposition 2.2) that

$`R(\epsilon )`$ $`=`$ $`{\displaystyle _{PS^2}}\mathrm{log}Tfvdm(v).`$
If it were the case that $`m_{\epsilon ,\delta }m`$ as $`\delta 0`$, then it would follow that:

$`\mathrm{\Lambda }(\epsilon )`$ $``$ $`\underset{\delta 0}{lim\; sup}R(\epsilon ,\delta )`$
$`=`$ $`\underset{\delta 0}{lim\; sup}{\displaystyle _{PS^2}}\mathrm{log}Tfvdm_{\epsilon ,\delta }(v)`$
$`=`$ $`{\displaystyle _{PS^2}}\mathrm{log}Tfvdm(v)`$
$`=`$ $`R(\epsilon ).`$
Hence the properties of this measure $`m_{\epsilon ,\delta }`$ are potentially quite interesting with regard to Question 1.1.

Here we collect some properties of $`m_{\epsilon ,\delta }`$ and $`R(\epsilon ,\delta )`$. First of all, $`R(\epsilon ,\delta )`$ is always positive for $`\delta >0`$ (in fact, we prove in Corollary 3.4 that this is true not just on average, but for individual elements of $`_\epsilon `$). In other words, introducing noise (no matter how small) to an element $`h_\epsilon `$ invariably produces positive exponents.

The measure $`m_{\epsilon ,\delta }`$ has additional properties as well. We prove that we can write:

$$m_{\epsilon ,\delta }=_{SO(3)}m_{g,\epsilon ,\delta }𝑑g,$$

where, for each $`gSO(3)`$, $`m_{g,\epsilon ,\delta }`$ is the unique probability measure on $`PS^2`$ with the properties:

1. $`m_{g,\epsilon ,\delta }`$ is stationary for the $`\delta `$-diffused process about $`gf_\epsilon `$;
2. $`m_{g,\epsilon ,\delta }`$ is absolutely continuous with respect to $`m`$, with smooth density;
3. $`m_{g,\epsilon ,\delta }`$ projects to Lebesgue measure $`\mu `$ on $`S^2`$.

As part of the proof, we show that each measure $`m_{g,\epsilon ,\delta }`$ is the unique fixed point of a “simple” linear operator.

Finally, from the way $`m_{\epsilon ,\delta }`$ is constructed, it follows that $`m_{\epsilon ,\delta }`$ shares all of the symmetries of $`f_\epsilon `$. In particular, the density $`\phi _{\epsilon ,\delta }`$ is invariant under all rotations that fix the North pole. The further study of these measures $`m_{g,\epsilon ,\delta }`$ and $`m_{\epsilon ,\delta }`$ might be of independent interest. We discuss the linear version of $`m_{\epsilon ,\delta }`$ in Section 4. In Section 8, we examine the behavior of $`m_{\epsilon ,\delta }`$ as $`\epsilon \mathrm{}`$.

We now turn to the proofs of assertions 1.-3.. We first prove a standard semicontinuity result for random exponents.

###### Lemma 3.1

Let $`\{\gamma _i\}`$ be a sequence of probability measures on $`\text{Diff }(M)`$ that converges weakly to a probability measure $`\gamma `$. Suppose that $`\mu `$ is stationary for the random process generated by $`(,\gamma _i)`$, for every $`i`$, (for example, if $`\text{Diff }_\mu (M)`$). Then

$$\underset{\gamma _i\gamma }{lim\; sup}R(\gamma _i)R(\gamma ).$$

Proof:  Let $`a_n((h_i)_1^{\mathrm{}},x)=\mathrm{log}T_xh^{(n)}`$. Then $`a_n:^{\mathrm{}}\times M𝐑`$ is subadditive with respect to $`\tau `$. By the subadditive ergodic theorem it then follows that

$`R(\gamma )`$ $`=`$ $`{\displaystyle _{^{\mathrm{}}\times M}}\underset{n\mathrm{}}{lim}{\displaystyle \frac{1}{n}}a_nd(\gamma ^{\mathrm{}}\times \mu )`$
$`=`$ $`\underset{n\mathrm{}}{lim}{\displaystyle \frac{1}{n}}{\displaystyle _{^{\mathrm{}}\times M}}a_nd(\gamma ^{\mathrm{}}\times \mu )`$
$`=`$ $`\underset{n}{inf}{\displaystyle \frac{1}{n}}{\displaystyle _{^{\mathrm{}}\times M}}a_nd(\gamma ^{\mathrm{}}\times \mu ).`$
Now for any fixed $`n`$ we have

$`\underset{\gamma _i\gamma }{lim\; sup}R(\gamma _i)`$ $``$ $`\underset{\gamma _i\gamma }{lim\; sup}{\displaystyle \frac{1}{n}}{\displaystyle _{^{\mathrm{}}\times M}}a_nd(\mu \times \gamma _i^{\mathrm{}})`$
$`=`$ $`{\displaystyle \frac{1}{n}}{\displaystyle _{^{\mathrm{}}\times M}}a_nd(\mu \times \gamma ^{\mathrm{}}).`$
So

$`\underset{\gamma _i\gamma }{lim\; sup}R(\gamma _i)`$ $``$ $`\underset{n}{inf}{\displaystyle \frac{1}{n}}{\displaystyle _{^{\mathrm{}}\times M}}a_nd(\mu \times \gamma ^{\mathrm{}})`$
$`=`$ $`R(\gamma ).`$
QED

We will apply this lemma to the situation where $`\gamma _i`$ is supported on a small ball in $`_\epsilon `$ converging, as $`i\mathrm{}`$, to a Dirac measure supported on an element of $`_\epsilon `$.

Let $`\delta >0`$, and let $`U_\delta `$ be a symmetric $`\delta `$ ball around the identity in $`SO(3)`$. Give $`U_\delta `$ the restriction of Haar measure, normalized to be a probability measure and similarly for $`_{g,\epsilon ,\delta }:=U_\delta gf_\epsilon `$, for every $`gSO(3)`$ and $`\epsilon >0.`$ We denote this last measure by $`\nu _{g,\epsilon ,\delta }`$. Let $`R(g,\epsilon ,\delta )=R(\nu _{g,\epsilon ,\delta })`$.

###### Definition 3.1

The ($`\delta `$-)diffused random exponent is the average over $`SO(3)`$ of $`R(g,\epsilon ,\delta )`$:

$$R(\epsilon ,\delta )=_{gSO(3)}R(g,\epsilon ,\delta )𝑑\nu (g).$$

###### Lemma 3.2

$$\underset{\delta 0}{lim\; sup}_{gSO(3)}R(g,\epsilon ,\delta )𝑑g\mathrm{\Lambda }(\epsilon ).$$

Proof:  Note that $`lim_{\delta 0}\nu _{g,\epsilon ,\delta }=\delta _{gf_\epsilon }`$, Dirac measure supported on $`gf_\epsilon `$. By the previous lemma,

$`\underset{\delta 0}{lim\; sup}R(g,\epsilon ,\delta )`$ $`=`$ $`\underset{\delta 0}{lim\; sup}R(\nu _{g,\epsilon ,\delta })`$
$``$ $`R(\delta _{gf_\epsilon })`$
$`=`$ $`\lambda (gf_\epsilon )`$
for each $`gSO(3)`$, so the same is true for the integral. QED

Now let $`h_0:S^2S^2`$ be any $`\mu `$-preserving diffeomorphism, and let $`=U_\delta h_0`$. As in the previous section, define the space

$$_{\mathrm{}}=\{h_{\mathrm{}}h\},$$

and evaluation maps

$$ev:\times S^2S^2,ev_{\mathrm{}}:\times PS^2PS^2.$$

Let $`\nu _\delta `$ and $`\nu _\delta \mathrm{}`$ be the push-forwards of normalized Haar measure to $``$ and $`_{\mathrm{}}`$, respectively. Note that $`\mu `$ is stationary for the process generated by $`(,\nu _\delta )`$.

###### Proposition 3.3

If $`h_0`$ is not an isometry, then for fixed $`\delta >0`$ the random process on $`PS^2`$ generated by $`(_{\mathrm{}},\nu _\delta \mathrm{})`$ has a stationary measure $`m_\delta `$ that is absolutely continuous with smooth density, covers $`\mu `$, and is the unique such stationary measure.

Moreover,

$$R(\nu _\delta )=_{PS^2}\mathrm{log}Th_0vdm_\delta (v)>0$$

###### Corollary 3.4

For fixed $`\delta >0`$, $`\epsilon 0`$ and $`gSO(3)`$, the random process on $`PS^2`$ generated by $`(_{g,\epsilon ,\delta },\nu _{g,\epsilon ,\delta })`$ has a stationary measure $`m_{g,\epsilon ,\delta }`$ that is absolutely continuous with smooth density, covers $`\mu `$, and is the unique such stationary measure.

Moreover,

$$R(g,\epsilon ,\delta )=_{PS^2}\mathrm{log}Tf_\epsilon vdm_{g,\epsilon ,\delta }>0$$

Proof of Corollary 3.4: Apply Proposition 3.3 to the case where $`h_0=gf_\epsilon `$.

Proof of Proposition 3.3: Let $`h_0`$ and $`\delta >0`$ be given. We break the proof into steps.

Step 1: Construction of $`m_\delta `$.

Recall that the convolution of a probability measure $`\pi `$ on $`_{\mathrm{}}`$ and a probability measure $`m`$ on $`PS^2`$ is a probability measure $`\pi m`$ on $`PS^2`$ defined by

$$\pi m(E)=_{_{\mathrm{}}}m(h_{\mathrm{}}^1E)𝑑\pi (h_{\mathrm{}})$$

for every $`m`$-measurable $`EPS^2`$. That a measure $`m`$ is stationary for the measure $`\pi `$ is equivalent to the fact that $`\pi _{\mathrm{}}m=m`$. For $`k>1`$ we let $`\pi ^km=\pi (\pi ^{k1}m)`$. For any probability measure $`m`$ on $`PS^2`$ any weak limit of the Cesàro sums $`\frac{1}{n}_1^n\pi ^km`$ is a stationary measure for $`\pi `$. Beginning with a measure $`m`$ which pushes forward under projection to $`\mu `$ produces an invariant measure by this process with the same property. If we start with $`m`$ as Liouville measure on $`PS^2`$ and $`\pi =\nu _\delta \mathrm{}`$ we call this limiting measure $`m_\delta `$.

Step 2: $`m_\delta `$ is absolutely continuous, with smooth density

For any measurable set $`APS^2`$, we have:

$`m_\delta (A)`$ $`=`$ $`(\nu _\delta \mathrm{}\times m_\delta )ev_{\mathrm{}}^1(A)`$
$`=`$ $`(\nu _\delta \mathrm{}\times m_\delta )\{(h_{\mathrm{}},v)|h_{\mathrm{}}(v)A\}`$
$`=`$ $`{\displaystyle _{vPS^2}}\nu _\delta \mathrm{}\{h_{\mathrm{}}|h_{\mathrm{}}(v)A\}𝑑m_\delta (v)`$
Now if the Liouville measure $`m(A)`$ equals zero, then $`\nu _\delta \{h_{\mathrm{}}|h_{\mathrm{}}(v)A\}`$ must also be zero, for every $`vPS^2`$. Thus $`m_\delta (A)`$ is zero. It follows that $`m_\delta (A)`$ is absolutely continuous with respect to $`m`$. So there is a non-negative integrable function $`\phi _\delta `$ defined on $`PS^2`$ so that for any measurable $`APS^2`$,

$$m_\delta (A)=_A\phi _\delta (x)𝑑m(x).$$

In Lemma 3.6 we will prove that $`\phi _\delta `$ satisfies the following formula:

$`\phi _\delta (x)`$ $`=`$ $`{\displaystyle \frac{1}{m(B(x,\delta ))}}{\displaystyle _{h_0\mathrm{}^1B(x,\delta )}}\phi _\delta (y)𝑑m(y)`$
$`=`$ $`{\displaystyle \frac{1}{m(B(x,\delta ))}}{\displaystyle _{B(x,\delta )}}\phi _\delta (h_0\mathrm{}^1z)\text{Jac}(h_0\mathrm{}^1)(z)𝑑m(z).`$
It follows now fairly directly that $`\phi _\delta `$ is as smooth as $`h_0`$, since the average over a $`\delta `$-ball of an $`L^1`$ function is continuous, and of a $`C^k`$ function, is $`C^{k+1}`$.

Step 3: $`R(\nu _\delta )`$ satisfies the integral formula, and the exponents of $`m_\delta `$ are nonzero.

The argument that $`R(\nu _\delta )=_{PS^2}\mathrm{log}Th_0vdm_\delta `$ is now the same as in the proof of proposition 2.2, using Birkhoff’s and Oseledec’s theorems, where $`_{\mathrm{}}`$ is replaced by $`_{\mathrm{}}`$, $`\nu `$ by $`\nu _\delta `$ and $`m`$ by $`m_\delta `$.

Next we will prove that the largest exponent is positive and from that we will deduce uniqueness. As in Proposition 2.2, for any $`h=gh_0`$ the Jacobian of $`h_{\mathrm{}}`$ with respect to $`m_\delta `$ at the vector $`v`$ in $`T_1S^2`$ is

$`\rho (h,v)`$ $`=`$ $`{\displaystyle \frac{\phi _\delta (h_{\mathrm{}}v)}{Thv^2\phi _\delta (v)}}={\displaystyle \frac{\phi _\delta (gh_0\mathrm{}v)}{Th_0v^2\phi _\delta (v)}}`$ (6)
provided that $`\phi _\delta (v)0`$.

We claim that the function $`\rho (h,v)`$ cannot be $`\nu _\delta \times m_\delta `$ \- almost everywhere equal to $`1`$. Suppose for the sake of contradiction that $`\rho (h,v)=1`$ a.e. Since $`\phi _\delta `$ is continuous, if we fix $`m_\delta `$-a.e. $`v`$, then the function $`\rho (,v)`$ is continuous. Thus, for $`m_\delta `$-a.e. $`v`$, we must have $`\rho (h,v)=1`$ for every $`hU_\delta `$.

Next, notice in the expression (6) for $`\rho (h,v)`$, that the only term that depends on $`gU_\delta `$ is the numerator $`\phi _\delta (gh_0\mathrm{}v)`$. Rewriting this expression, we have, for almost every $`v`$ in the support of $`\phi _\delta `$,

$`\phi _\delta (gh_0\mathrm{}v)`$ $`=`$ $`Th_0v^2\phi _\delta (v),`$ (7)
for every $`gU_\delta `$. Since $`m_\delta `$ projects to $`\mu `$, and $`\phi _\delta `$ is continuous, we have that for $`\mu `$-a.e. $`xS^2`$, the set $`O=\{v\phi _\delta (v)>0\}`$ is an open, $`h_0\mathrm{}`$-invariant set in $`PS^2`$ that intersects almost every fiber. Equation (7) implies $`\phi _\delta `$ must be constant on connected components of $`O`$, since varying $`g`$ inside of $`U_\delta `$, the vector $`gh_0\mathrm{}v`$ covers an open neighborhood of $`h_0\mathrm{}v`$ in $`PS^2`$.

But, again by equation (7), on each such component of $`O`$, we must have that $`Th_0`$ is constant. Since $`O`$ intersects almost every fiber of $`PS^2`$, we obtain that for almost every $`xS^2`$, there exists a connected open set $`I_x`$ in the fiber of $`PS^2`$ on which $`T_xh_0`$ has constant norm. Since $`h_0`$ preserves area, we must have $`T_xh_0=1`$ for $`\mu `$-a.e. $`x`$, contradicting the assumption that $`h_0`$ is not an isometry.

So $`\rho (h,v)`$ is not a.e. equal to $`1`$, and by Jensen’s inequality, we have:

$`{\displaystyle _{\times PS^2}}\mathrm{log}\rho (h,v)𝑑\nu _\delta (h)𝑑m_\delta (v)`$ $`<`$ $`\mathrm{log}{\displaystyle _{\times PS^2}}\rho (h,v)𝑑\nu _\delta (h)𝑑m_\delta (v)`$
$`=`$ $`0.`$
But

$`R(\nu _\delta )`$ $`=`$ $`{\displaystyle _{PS^2}}\mathrm{log}Th_0vdm_\delta (v)`$
$`=`$ $`{\displaystyle \frac{1}{2}}{\displaystyle _{\times PS^2}}\mathrm{log}Thv^2d\nu _\delta (h)𝑑m_\delta (v)`$
$`=`$ $`{\displaystyle \frac{1}{2}}{\displaystyle _{\times PS^2}}\mathrm{log}\left(Thv^2{\displaystyle \frac{\phi _\delta (h_{\mathrm{}}(v))}{\phi _\delta (v)}}\right)𝑑\nu _\delta (h)𝑑m_\delta (v),`$
$`=`$ $`{\displaystyle \frac{1}{2}}{\displaystyle _{\times PS^2}}\mathrm{log}\rho (h,v)𝑑\nu _\delta (h)𝑑m_\delta (v)`$
$`>`$ $`0.`$
(Here we used the stationarity of $`m_\delta `$ to conclude that the integral of $`\mathrm{log}(\phi _\delta (h_{\mathrm{}}(v))/\phi _\delta (v))`$ is $`0`$).

Part 4: $`m_\delta (A)`$ is unique

It remains to prove that the measure $`m_\delta `$ is unique among absolutely continuous stationary measures which cover $`\mu `$, now that we know that the random Lyapunov exponents are not zero.

Let $`\gamma `$ be any such measure. Then, as for $`m_\delta `$, there is a nonnegative smooth function $`\psi `$ such that for any measurable $`APS^2`$,

$$\gamma (A)=_A\psi (v)𝑑m(v).$$

Let $`\gamma _x`$ be the disintegration of $`\gamma `$ on the fiber $`T_{1,x}S^2`$. The density of $`\gamma _x`$ with respect to Lebesgue measure on the fiber $`T_{1,x}S^2`$ is the restriction of $`\psi `$ to the fiber.

We need the notion of a natural extension of a non-invertible transformation. Let $`(\mathrm{\Omega },𝒜,p)`$ be a probability space with $`p`$-preserving transformation $`T:\mathrm{\Omega }\mathrm{\Omega }`$. Let

$$\widehat{\mathrm{\Omega }}=\{(\mathrm{},\omega _1,\omega _0)\mathrm{\Pi }_{\mathrm{}}^1\mathrm{\Omega }T(\omega _i)=\omega _{i+1},i1\}.$$

Let $`\widehat{T}:\widehat{\mathrm{\Omega }}`$ be the map:

$$\widehat{T}(\mathrm{},\omega _1,\omega _0)=(\mathrm{},\omega _1,\omega _0,T(\omega _0)),$$

and let $`\widehat{\pi }_0:\widehat{\mathrm{\Omega }}\mathrm{\Omega }`$ be the projection onto the first factor:

$$\widehat{\pi }_0(\mathrm{},\omega _1,\omega _0)=\omega _0.$$

Let $`\widehat{𝒜}`$ be the smallest $`\sigma `$-algebra on $`\widehat{\mathrm{\Omega }}`$ so that $`\widehat{\pi }_0`$ and $`\widehat{T}`$ are both measurable. On $`(\widehat{\mathrm{\Omega }},\widehat{𝒜})`$, there is a unique probability measure $`\widehat{p}`$, invariant under $`\widehat{T}`$, that pushes forward under $`\widehat{\pi }`$ to $`p`$. The measure-preserving system $`(\widehat{\mathrm{\Omega }},\widehat{𝒜},\widehat{p})\widehat{T}`$ is called the natural extension of $`(\mathrm{\Omega },𝒜,p)T`$. The natural extension $`\widehat{T}`$ is invertible, and ergodic if $`T`$ is ergodic.

Let $`\tau `$ and $`\tau _{\mathrm{}}`$ be the associated auxiliary maps to the processes generated by $`(,\nu _\delta )`$ and $`(_{\mathrm{}},\nu _\delta \mathrm{})`$, respectively. The natural extension of $`\tau _{\mathrm{}}`$ with respect to the measure $`\nu _\mathrm{}\delta ^{\mathrm{}}\times \gamma `$ fibers over the natural extension of $`\tau `$ with respect to $`\nu _\delta ^{\mathrm{}}\times \mu `$; the fiber over $`((h_i)_{\mathrm{}}^{\mathrm{}},x)`$ is $`T_{1,x}S^2`$. For $`\widehat{\nu }^{\mathrm{}}`$-almost every $`𝐡=(h_i)_{\mathrm{}}^{\mathrm{}}`$ and almost every fiber $`T_{1,x}S^2`$ there is a measure $`\widehat{\gamma }_{𝐡,x}`$ which is the disintegration of $`\widehat{\nu ^{\mathrm{}}\times \gamma }`$ along the fiber. Note that the extension $`\widehat{\nu ^{\mathrm{}}\times \gamma }`$ is determined by this system of measures, and therefore so is $`\gamma `$.

By (\[Li-Qi\], Proposition 1.1, p.131), if $`(𝐡=(h_i)_{\mathrm{}}^{\mathrm{}},x)_{\mathrm{}}^{\mathrm{}}\times S^2=\widehat{^{\mathrm{}}\times S^2}`$, then $`\widehat{\gamma }_x`$ is the limit of the push-forwards:

$`\widehat{\gamma }_{𝐡,x}`$ $`=`$ $`\underset{n\mathrm{}}{lim}(h_n\mathrm{}\mathrm{}h_1\mathrm{})_{}(\gamma _{h_n^1\mathrm{}h_1^1(x)})`$ (8)
###### Lemma 3.5

The limit (8) does not depend on the initial choice of absolutely continuous stationary measure $`\gamma `$ covering $`\mu `$.

Proof:  The average exponents of $`\tau `$ are nonzero, so the average exponents of $`\widehat{\tau }`$ are also nonzero. It is easy to see that $`\tau `$ is ergodic, and thus, so is $`\widehat{\tau }`$ and the exponents of $`\widehat{\tau }`$ are in fact nonzero $`\mu `$-a.e.

Let $`u(𝐡,x)PS^2`$ be the unstable Lyapunov direction for $`\widehat{\tau }`$ over $`(𝐡,x)`$. For every $`ϵ>0`$, there exists an $`n>0`$, and a set $`G_{\mathrm{}}^{\mathrm{}}\times S^2`$ such that

* $`\nu _{\mathrm{}}^{\mathrm{}}\mu (G)>1ϵ`$
* for every $`(𝐡,x)G`$, the $`\gamma _{𝐡,x}`$-measure of an $`ϵ`$-neighborhood of $`u(𝐡,x)`$ in $`T_{1,x}S^2`$ is at least $`1ϵ`$.

It follows that the limit in (8) is concentrated on the point $`u(𝐡,x)`$. QED

Thus the natural extension of $`\tau `$ and $`\gamma `$ is the same as the natural extension of $`\tau `$ and $`m_\delta `$, so $`\tau `$ and $`m_\delta `$ are themselves equal. QED

The next lemma completes the proof of Proposition 3.3.

###### Lemma 3.6

Let $`\gamma `$ be an absolutely continuous measure with respect to $`m`$ on $`PS^2`$ that is stationary for $`\nu _\delta `$. Then the density function $`\psi `$ defining $`\gamma `$ satisfies

$$\psi (x)=\frac{1}{m(B(x,\delta ))}_{h_0\mathrm{}^1B(x,\delta )}\psi (y)𝑑m(y).$$

Proof:  For any measurable set $`APS^2`$, we have:

$`{\displaystyle _A}\psi (y)𝑑m(y)`$ $`=`$ $`\gamma (A)`$
$`=`$ $`\nu _\delta \times \gamma (ev_{\mathrm{}}^1(A))`$
$`=`$ $`\nu _\delta \times \gamma \{(h,y)_{\mathrm{}}\times PS^2h(y)A\}`$
$`=`$ $`{\displaystyle _{gU_\delta }}\gamma \{h_0\mathrm{}^1g^1A\}𝑑\nu _\delta (g)`$
$`=`$ $`{\displaystyle _{gU_\delta }}{\displaystyle _{h_0\mathrm{}^1g^1A}}\psi (y)𝑑m(y)𝑑\nu _\delta (g)`$
$`=`$ $`{\displaystyle _{gU_\delta }}{\displaystyle _A}\psi (h_0\mathrm{}^1g^1(x))J(h_0\mathrm{}^1g^1(x))^1𝑑m(x)𝑑\nu _\delta (g)`$
$`=`$ $`{\displaystyle _A}{\displaystyle _{gU_\delta }}\psi (h_0\mathrm{}^1g^1(x))J(h_0\mathrm{}^1g^1(x))^1𝑑\nu _\delta (g)𝑑m(x)`$
So,

$$\psi (x)=_{gU_\delta }\psi (h_0\mathrm{}^1g^1(x)J(h_0\mathrm{}^1g^1(x))^1d\nu _\delta (g).$$

Since the integrand only depends on the point $`xPS^2`$ we push the measure $`\nu _\delta `$ forward to $`PS^2`$, use that $`J(g)^1=1`$ and that Haar measure on $`SO(3)`$ pushes forward to Liouville measure on $`PS^2`$ to obtain

$$\psi (x)=\frac{1}{m(B(x,\delta ))}_{yB(x,\delta )}\psi (h_0\mathrm{}^1(y))J(h_0\mathrm{}^1(y))𝑑m(y).$$

Finally, changing variables one more time, gives:

$$\psi (x)=\frac{1}{m(B(x,\delta ))}_{h_0\mathrm{}^1B(x,\delta )}\psi (y)𝑑m(y).$$

 QED

This completes the proof of Proposition 3.3.

Returning to discussion of the family $`_\epsilon `$, we have verified that properties 1.-3. hold for the measures $`m_{g,\epsilon ,\delta }`$. Now let

$$m_{\epsilon ,\delta }=_{gSO(3)}m_{g,\epsilon ,\delta }𝑑g.$$

Since

$$_{SO(3)}R(g,\epsilon ,\delta )𝑑g=_{PS^2}\mathrm{log}Tfvdm_{g,\epsilon ,\delta }(v)𝑑g,$$

it follows from Corollary 3.4 and Lemma 3.2 that:

###### Proposition 3.7

$$\mathrm{\Lambda }(\epsilon )\underset{\delta 0}{lim\; sup}_{PS^2}\mathrm{log}Tf_\epsilon (x)vdm_{\epsilon ,\delta }$$

### 3.2 What is this measure $`m_{ϵ,\delta }`$?

For the family of twist maps $`_\epsilon `$ under consideration, we prove in Section 7 that the answer to Question 1.1 is “no”, and for small $`\epsilon `$, we have $`\mathrm{\Lambda }(\epsilon )<R(\epsilon )`$. It then follows from Proposition 3.7 that for small $`\epsilon `$, Lebesgue measure $`m`$ is not a weak limit of $`m_{\epsilon ,\delta }`$ as $`\delta 0`$. At the opposite extreme, we show in Section 8 that as $`\epsilon `$ tends to infinity, the measures $`m_{\epsilon ,\delta }`$ do approach Lebesgue measure, for $`\delta >0`$ fixed. We hope that a future experiment will reveal more precisely how these measures behave in $`\delta `$, for moderate values of $`\epsilon `$. As will be seen in section 6, the tiny differences between $`\mathrm{\Lambda }(\epsilon )`$ and $`R(\epsilon )`$ for moderate and large values of $`\epsilon `$, are expected to give rise to numerical difficulties in estimating the behavior with respect to $`\delta .`$

In Section 4, we show that in $`SO(2)`$-invariant families of $`2\times 2`$ matrices, if $`m_\delta `$ is the analogous “in-between” measure, averaged over the family, then $`m_\delta `$ is Lebesgue measure on $`S^1`$, for all $`\delta >0`$. On the other hand, we also show that for unitarily-invariant families, as $`\delta 0`$, the $`m_\delta `$ do not approach the natural unitarily invariant measure on the appropriate Grassmannian manifold, but instead they limit on an even “better” measure, in the sense that this measure forces the (strict) inequality $`\mathrm{\Lambda }>R`$. We describe the construction of $`m_\delta `$ for matrices in Section 4.

## 4 The linear case

Question 1.1 was originally motivated by a result of Dedieu-Shub about random and deterministic exponents for families of matrices. In this section, we describe these results and apply the framework of the previous section to the matrix setting.

Let $`L_i`$ be a sequence of linear maps mapping finite dimensional normed vector spaces $`V_i`$ to $`V_{i+1}`$ for $`i𝐍`$. Let $`vV_0\backslash \{0\}`$. If the limit $`lim\frac{1}{k}\mathrm{log}L_{k1}\mathrm{}L_0(v)`$ exists it is called a Lyapunov exponent of the sequence. It is easy to see that if two vectors have the same exponent then so does every vector in the space spanned by them. It follows that there are at most $`\text{dim}(V_0)`$ exponents. We denote them $`\lambda _j`$ where $`jk\text{dim}(V_0)`$. We order the $`\lambda _i`$ so that $`\lambda _i\lambda _{i+1}`$.

Given a probability measure $`\mu `$ on $`GL(n,𝐂)`$, the space of invertible $`n\times n`$ complex matrices, we may form infinite sequences of elements chosen at random from $`\mu `$ by taking the product measure on $`GL(n,𝐂)^𝐍`$. Thus we may also talk about the Lyapunov exponents of sequences or almost all sequences in $`GL(n,𝐂)^𝐍`$.

For measures $`\mu `$ on $`GL(n,𝐂)`$ satisfying a mild integrability condition, we have, by Oseledec’s Theorem, $`n`$ Lyapunov exponents $`r_1r_2\mathrm{}r_n\mathrm{}`$ such that, for almost every sequence $`\mathrm{}g_k\mathrm{}g_1GL(n,𝐂)`$, the limit

$$lim\frac{1}{k}\mathrm{log}g_k\mathrm{}g_1v$$

exists for every $`v𝐂^n\{0\}`$ and equals one of the $`r_i`$, $`i=1\mathrm{}n`$, see Gol’dsheid and Margulis \[GoMa\] or Ruelle \[Ru\] or Oseledec \[Os\]. We may call the numbers $`r_1,\mathrm{},r_n`$ random Lyapunov exponents or even just random exponents. If the measure is concentrated on a point $`A`$, these numbers:

$$\lambda _i(A)=lim\frac{1}{n}\mathrm{log}A^nv,i=1\mathrm{}n,$$

are $`\mathrm{log}|e_1|`$, …, $`\mathrm{log}|e_n|`$, where $`e_i(A)=e_i`$, $`i=1\mathrm{}n`$, are the eigenvalues of $`A`$ written with multiplicity and $`|e_1||e_2|\mathrm{}|e_n|.`$

The integrability condition for Oseledec’s Theorem is

$$gGL(n,𝐂)\mathrm{log}^+(g)\text{is }\mu \text{integrable}$$

where for a real valued function $`f`$, $`f^+=\mathrm{max}[0,f].`$ Here we will assume more so that all our integrals are defined and finite, namely:

$$()gGL(n,𝐂)\mathrm{log}^+(g)\text{and }\mathrm{log}^+(g^1)\text{are }\mu \text{integrable}.$$

In this matrix setting, there are rigorous lower bounds for the average exponents (= logarithms of moduli of eigenvalues) of unitarily-invariant families in $`GL(n,𝐂)`$. In \[DeSh\], the following bound is proved:

###### Theorem 4.1

\[DeSh\] If $`\mu `$ is a unitarily invariant measure on $`GL(n,𝐂)`$ satisfying $`()`$ then, for $`k=1,\mathrm{},n,`$

$$_{AGL(n,𝐂)}\underset{i=1}{\overset{k}{}}\mathrm{log}|\lambda _i(A)|d\mu (A)\underset{i=1}{\overset{k}{}}r_i.$$

By unitary invariance we mean $`\mu (U(X))=\mu (X)`$ for all unitary transformations $`UU(n,𝐂)`$ and all $`\mu `$-measurable $`XGL(n,𝐂)`$.

We can rephrase a special case of this theorem in a form similar to Question 1.1. Fix $`AGL(n,𝐂)`$. As above, let $`\nu `$ be normalized Haar measure on $`U(n,𝐂)`$, and also denote by $`\nu `$ the push-forward of $`\nu `$ to the coset $`U(n,𝐂)AGL(n,𝐂)`$. Let $`R(A)=r_1(A)`$ be the largest random exponent of $`\nu `$, and let

$$\mathrm{\Lambda }(A)=_{BU(n,𝐂)A}\mathrm{log}|e_1(B)|d\nu (B).$$

Then we have:

###### Corollary 4.2

\[DeSh\] For $`n2`$, and for any $`AGL(n,𝐂)`$,

$$\mathrm{\Lambda }(A)R(A).$$

Equality holds if and only if $`AU(n,𝐂)`$.

Thus non-zero Lyapunov exponents for the family, i.e. non-zero random exponents, implies that at least some of the individual linear maps have non-zero exponents, i.e eigenvalues of modulus not equal to 1. Hence the question that we posed for diffeomorphisms has a positive answer for sufficiently rich (i.e. unitarily-invariant) families of matrices.

###### Remark 4.3

Theorem 4.1 is not true for general measures on $`GL(n,𝐂)`$ or $`GL(n,𝐑)`$ even for $`n=2`$. Consider

$$A_1=\left(\begin{array}{cc}1& 0\\ 1& 1\end{array}\right),A_2=\left(\begin{array}{cc}1& 1\\ 0& 1\end{array}\right),$$

and give probability $`1/2`$ to each. The left hand integral is zero but as is easily seen the right hand sum is positive. So, in this case the inequality goes the other way. We do not know a characterization of measures which make Theorem 4.1 valid.

We expect similar results for orthogonally invariant probability measures on $`GL(n,𝐑)`$ but the only case in which such a result has been proved is in dimension 2, where we have:

###### Theorem 4.4

\[DeSh\] Let $`\mu `$ be a probability measure on $`GL(2,𝐑)`$ satisfying

$$gGL(2,𝐑)\mathrm{log}^+(g)\text{and}\mathrm{log}^+(g^1)\text{are }\mu \text{integrable}.$$

a. If $`\mu `$ is a $`SO(2,𝐑)`$ invariant measure on $`GL^+(2,𝐑)`$ then,

$$_{AGL^+(2,𝐑)}\mathrm{log}|\lambda _1(A)|d\mu (A)=_{AGL^+(2,𝐑)}_{xS^1}\mathrm{log}Axdxd\mu (A).$$

b. If $`\mu `$ is a $`SO(2,𝐑)`$ invariant measure on $`GL^{}(2,𝐑)`$, whose support is not contained in $`O(2,𝐑)`$ i.e. in the set of scalar multiples of orthogonal matrices, then

$$_{AGL^{}(2,𝐑)}\mathrm{log}|\lambda _1(A)|d\mu (A)>_{AGL^{}(2,𝐑)}_{xS^1}\mathrm{log}Axdxd\mu (A).$$

Here $`GL^+(2,𝐑)`$ (resp. $`GL^{}(2,𝐑)`$) is the set of invertible matrices with positive (resp. negative) determinant. To rephrase Theorem 3 a., fix $`AGL^+(2,𝐑)`$, and let $`\nu `$ be the push forward of Haar measure on $`SO(2,𝐑)`$ to the coset $`SO(2,𝐑)A`$. Let $`R(A)`$ be the largest random exponent for the process induced by $`\nu `$, and let $`\mathrm{\Lambda }(A)=_{BSO(2,𝐑)A}\mathrm{log}|e_1(A)|d\nu (A)`$. The we have:

###### Corollary 4.5

\[DeSh, AvBo\] For any $`AGL^+(2,𝐑)`$,

$$\mathrm{\Lambda }(A)=R(A).$$

We give an alternate proof of Corollary 4.5 in the following subsection.

### 4.1 $`m_\delta `$ for matrices

Let $`AGL(n,𝐂)`$ or $`GL(n,𝐑)`$ and $`\mu `$ be the Haar measure on $`U(n,𝐂)`$ or $`SO(n,𝐑)`$, respectively, normalized to be a probability measure. Let $`G`$ denote $`GL(n,𝐂)`$ or $`GL(n,𝐑)`$. As we did for families of diffeomorphisms in Section 3, we now interpolate between random products and deterministic powers of matrices by changing $`\mu `$. Let $`\delta >0`$ and $`G_\delta `$ be the $`\delta `$-neighborhood of the identity in $`G`$. For $`gG`$, $`G_\delta gA`$ is a neighborhood of $`gA`$ in $`GA`$. We normalize Haar measure restricted to $`G_\delta `$ and push it forward to $`G_\delta gA`$. Let us call this measure $`\mu _{\delta ,g}.`$ Let $`r_1(\delta ,g)`$ be the largest random exponent for this measure. At the end of this subsection, we prove:

###### Proposition 4.6

$`lim_{\delta 0}r_1(\delta ,g)=\mathrm{log}|e_1(gA)|.`$

Let $`G_{n,k}(𝐂)`$ denote the Grassmannian manifold of $`k`$ dimensional vector subspaces in $`𝐂^n`$, and let $`m`$ be the natural unitarily invariant probability measure on $`G_{n,k}(𝐂)`$. Any $`n\times n`$ complex matrix acts on the homogeneous space $`G_{n,k}`$ by left-multiplication. For $`AGL(n,𝐂)`$ and $`PG_{n,k}(𝐂)`$, denote by $`A|P`$ the restriction of $`A`$ to the subspace $`P`$. Now let $`m_{\delta ,g}`$ be the stationary measure on the Grassmannian $`G_{n,1}(=𝐂P^{n1})`$ induced by $`\mu _{\delta ,g}.`$

###### Proposition 4.7

$`r_1(\delta ,g)=_{PG_{n,1}}\mathrm{log}(A|P)dm_{\delta ,g}`$

Let $`m_\delta =_{gG}\nu _{\delta ,g}𝑑\mu .`$ It follows that:

###### Proposition 4.8

$`_{gG}\lambda _1(gA)d\mu =lim_{\delta 0}_{PG_{n,1}}\mathrm{log}(A|P)dm_\delta .`$

Now

$$r_1=_{PG_{n,1}}\mathrm{log}(A|P)dm(P).$$

So a comparison of $`_{gG}\lambda _1(gA)𝑑\mu `$ and $`r_1`$ can be achieved via an understanding of the relationship between $`m`$ and $`m_\delta .`$ We have two results in this direction.

First, recall that inequality in Corollary 4.2 is strict when $`n2`$ unless $`A`$ is an isometry. So for $`G=SU(n)`$ the measures $`m_\delta `$ favor the expanding directions of $`A`$ as $`\delta 0`$. By Proposition 4.8, we obtain the immediate corollary:

###### Corollary 4.9

For $`G=SU(n)`$, $`n2`$,

$$\underset{\delta 0}{lim}m_\delta m,$$

unless $`A`$ is an isometry.

Experimentally the same seems to hold for $`SO(n)`$ when $`n>2`$, but we have not checked this very carefully.

By contrast, the equality in Corollary 4.5 is consistent with $`lim_{\delta 0}m_\delta =m`$. In fact, more is true:

###### Theorem 4.10

For $`G=SO(2)`$,

$$m_\delta =m,$$

for all $`\delta >0.`$

Combined with Proposition 4.8, Theorem 4.10 gives another proof of Corollary 4.5.

Proof of Proposition 4.6: Following the proof of Lemma 3.2, one obtains that

$$\mathrm{log}|e_1(gA)|\underset{\delta 0}{lim}r_1(\delta ,g).$$

If $`A`$ is replaced by $`cA`$, for $`c𝐂\{0\}`$, then both sides of the equality change by $`\mathrm{log}|c|`$. So we may assume that $`|detA|=1`$, and it will be enough to prove that $`\mathrm{log}|e_1(gA)|lim_{\delta 0}r_1(\delta ,g)`$ under the hypothesis that $`|e_1(gA)|>1`$.

Let $`E𝐂^n`$ be the generalized eigenspace of the eigenvalues of $`gA`$ whose modulus equals $`|e_1(gA)|`$. Then given $`\epsilon >0`$, there is a metric on $`𝐂^n`$, a closed cone $`K𝐂^n`$ containing $`E`$ in its interior, and a $`\delta >0`$ such that, for any $`BGL(n,𝐂)`$ in the $`\delta `$-neighborhood of $`gA`$, we have:

1. $`B(K)K`$, and
2. $`Bv(|e_1(gA)|\epsilon )v`$, for all $`vK`$.

It follows that, for any $`\epsilon >0`$, there is a $`\delta >0`$ such that,

$$B_n\mathrm{}B_1(|e_1(gA)|\epsilon )^n,$$

for all sequences $`B_1,\mathrm{},B_n`$ in $`U_\delta gA`$. Hence

$$r_1(\delta ,g)\mathrm{log}(|e_1(gA)|\epsilon ),$$

and so

$$\underset{\delta 0}{lim}r_1(\delta ,g)\mathrm{log}|e_1(ga)|.$$

 QED

We now turn to the proof of Theorem 4.10.

Proof of Theorem 4.10: By an argument presented in the proof of Proposition 4.6, we may assume that $`detA=1`$. The projective action of $`SL(2,𝐑)`$ on $`𝐑P^1`$ is conjugate to the standard action on the circle $`S^1=\{z𝐂||z|=1\}`$ by linear fractional transformations. The conjugacy sends the rotation by $`\theta `$ to rotation by $`2\theta `$. Let $`f:S^1S^1`$ be the linear fractional transformation induced by $`A`$, and let $`=\{\alpha f|\alpha S^1\}`$. Let

$$U_{\alpha ,\delta }=\{\beta \alpha f|\mathrm{arg}(\beta )(\delta ,\delta )\},$$

and let $`\nu _{\alpha ,\delta }`$ be normalized Lebesgue measure on $`U_{\alpha ,\delta }`$, pushed forward from $`(\delta ,\delta )`$. Denote by $`m_{\alpha ,\delta }`$ the stationary measure on $`S^1`$ induced by $`\nu _{\alpha ,\delta }`$. We will show that:

$$m_\delta =_{S^1}m_{\alpha ,\delta }𝑑\alpha $$

is Lebesgue measure on $`S^1`$.

The same argument as in the proof of Lemma 3.6 shows that for every $`\alpha S^1`$,

$$dm_{\alpha ,\delta }(z)=\phi _{\alpha ,\delta }(z)dz,$$

where

$$\phi _{\alpha ,\delta }(z)=\frac{1}{2\delta }_{y(\alpha f)^1B(z,\delta )}\phi _{\alpha ,\delta }(y)𝑑y,$$

and $`B(z,\delta )=\{\beta z|\mathrm{arg}(\beta )(\delta ,\delta )\}`$. Setting

$$k_\delta (z)=\frac{1}{2\delta }1_{B(1,\delta )}(z),$$

we have that

$$\phi _{\alpha ,\delta }(z)=_{yS^1}k_\delta (\alpha \overline{z}f(y))\phi _{\alpha ,\delta }(y)𝑑y,$$

where we use $`\overline{z}`$ to denote the multiplicative inverse of $`zS^1`$. Note that $`_{S^1}k_\delta (z)𝑑z=1`$.

Consider the following, more general setting. Let $`k`$ be a non-negative function on $`S^1`$ such that $`k=1`$ and $`k^2<\mathrm{}`$ (all integrals are with respect to the normalized Lebesgue measure on $`S^1`$). For example, $`k=k_\delta `$. Let $`f:S^1S^1`$ be a fractional linear transformation, so that we can write, for $`|z|=1`$:

$$f(z)=\underset{n0}{}c_nz^n.$$

For $`\alpha S^1`$ , define the operator $`L_{\alpha ,f}`$ on real functions on $`S^1`$ by:

$$L_{\alpha ,f}\phi (z)=k(\alpha \overline{z}f(y))\phi (y)𝑑y.$$

The operator $`L_{\alpha ,f}`$ is a positive operator, $`L_{\alpha ,f}\phi =\phi `$. There exists a unique function $`\phi _\alpha `$ satisfying $`L_{\alpha ,f}\phi _\alpha =\phi _\alpha `$ and $`\phi _\alpha =1`$. The function $`\phi _\alpha `$ is positive and continuous, upper and positive lower bounds for $`\phi _\alpha `$ can be chosen uniformly in $`\alpha `$. In the case $`k=k_\delta `$, we have $`\phi _\alpha =\phi _{\alpha ,\delta }`$.

###### Proposition 4.11

We have, for all $`zS^1`$, $`\phi _\alpha (z)𝑑\alpha =1`$.

The proposition follows directly from the following two claims:

Claim 1. For all $`m0`$, all $`zS^1`$, $`L_{\alpha ,f}^m1(z)𝑑\alpha =1`$.

Claim 2. The sequence $`\frac{1}{n}_{m=1}^nL_{\alpha ,f}^m1(z)`$ converges to $`\phi _\alpha (z)`$ in $`L^1(d\alpha ,dz)`$ as $`n\mathrm{}`$.

Claim 2 follows from the ergodic theorem for Markov (i.e. $`L1=1`$) operators applied to the operator $`\psi (\alpha ,z){\displaystyle \frac{1}{\phi _\alpha (z)}}L_{\alpha ,f}(\phi _\alpha \psi (\alpha ,.))(z)`$ and the initial function $`{\displaystyle \frac{1}{\phi _\alpha (z)}}`$.

In order to prove Claim 1, we compute, for a function $`\phi L^2`$, $`\phi (z)=_n\gamma _nz^n`$, the Fourier coefficients $`\gamma _n^{}`$ of the function $`L_{\alpha ,f}\phi `$. We find, for $`n0`$:

$`\gamma _n^{}`$ $`=`$ $`{\displaystyle \overline{z}^nL_{\alpha ,f}\phi (z)𝑑z}`$
$`=`$ $`{\displaystyle \overline{z}^nk(\alpha \overline{z}f(y))\phi (y)𝑑y𝑑z}`$
$`=`$ $`\alpha ^n\widehat{k}(n){\displaystyle (f(y))^n\phi (y)𝑑y}`$
$`=`$ $`\alpha ^n\widehat{k}(n){\displaystyle \underset{k0}{}c_k^{(n)}y^k\phi (y)dy}`$
$`=`$ $`\alpha ^n\widehat{k}(n){\displaystyle \underset{k0}{}}c_k^{(n)}\gamma _k`$
and, analogously:

$$\gamma _n^{}=\alpha ^n\widehat{k}(n)\underset{k0}{}c_k^{(n)}\gamma _k,$$

where we wrote $`(f(z))^n=_{k0}c_k^{(n)}z^k,(f(z))^0=1`$, and $`\widehat{k}`$ is the Fourier transform of $`k`$.

Iterating these formulas, we obtain for a function $`\phi L^2`$, $`\phi (z)=_n\gamma _nz^n`$, the Fourier coefficients $`\gamma _n^{(m)}`$ of the function $`L_{\alpha ,f}^m\phi `$:

$$\gamma _n^{(m)}=\alpha ^n\widehat{k}(n)\underset{n_1,\mathrm{},n_m0}{}(\mathrm{\Pi }_{s=1}^{m1}\alpha ^{n_s}\widehat{k}(n_s))c_{n_1}^{(n)}c_{n_2}^{(n_1)}\mathrm{}c_{n_m}^{(n_{m1})}\gamma _{n_m}$$

for $`n0`$, and

$$\gamma _n^{(m)}=\alpha ^n\widehat{k}(n)\underset{n_1,\mathrm{},n_m0}{}(\mathrm{\Pi }_{s=1}^{m1}\alpha ^{n_s}\widehat{k}(n_s))c_{n_1}^{(n)}c_{n_2}^{(n_1)}\mathrm{}c_{n_m}^{(n_{m1})}\gamma _{n_m}$$

for $`n0`$.

To get the Fourier coefficients of the bounded continuous function $`L_{\alpha ,f}^m1(z)𝑑\alpha `$, we integrate in $`\alpha `$ the Fourier coefficients of the bounded continuous functions $`L_{\alpha ,f}^m1(z)`$. In the above sum, all terms vanish, except the ones with $`n+n_1+\mathrm{}+n_{m1}=0`$. Since all $`n_i`$ have the same sign, the one nonzero integral corresponds to $`n=n_1=\mathrm{}=n_{m1}=0`$. Claim 1 follows. QED

###### Remark 4.12

Theorem 4.10 holds even without randomization. Suppose that $`A`$ has determinant equal to $`1`$ and let $`O`$ vary over $`SO(2,𝐑).`$ Then for almost every $`O`$ the eigenvalues of $`OA`$ are either complex with irrational argument or real and there is one eigenvalue of modulus bigger than one. In the first case Cesaro sums of the push forward of Lebesgue measure by $`OA_{\mathrm{}}`$ converge to the unique invariant measure of $`OA_{\mathrm{}}`$. In the second case to the Dirac measure supported on the expanding eigenspace. Call these measures $`m_{OA}`$, then $`m_{OA}𝑑O`$ is Lebesgue measure. The proof is the same, but easier.

## 5 Experimental Method

### 5.1 Haar Measure on $`SO(3)`$

It is clear that an element of $`SO(3)`$ is determined by its axis and angle of rotation. Here we describe how to pick axis and angle uniformly with respect to Haar measure on $`SO(3)`$.

Let $`S^2`$ be the usual two sphere with measure $`\mu `$. Let $`S^1`$ be the usual unit circle of angles from $`0`$ to $`2\pi `$ given the probability measure with density function $`(1\mathrm{cos}\theta )/(2\pi )`$. Let $`S^2\times S^1`$ be the product space with the product measure which we denote by $`\gamma `$. There is a natural map $`P:S^2\times S^1SO(3)`$ which maps a vector $`x`$ and an angle $`\theta `$ to the orthogonal transformation which fixes $`x`$ and rotates by angle $`\theta `$ around $`x`$ according to the right hand rule. The map $`P`$ sends $`(x,0)`$ to the identity in $`SO(3)`$ for all $`xS^2`$ and is two to one when $`\theta 0`$, $`(x,\theta )`$ and $`(x,\theta )`$ map to the same point. $`P`$ maps $`\gamma `$ to the Haar measure on $`SO(3)`$. If we identify $`(x,\theta )(x,\theta )`$ we obtain $`S^2\times S^1/`$ which is a circle bundle over real projective 2-space. $`P`$ induces a map $`P^{}:S^2\times S^1/SO(3)`$.

###### Proposition 5.1

$`P^{}`$ is one-one off of the zero section of $`S^2\times S^1/`$ and gives a measurable isomorphism between $`(S^2\times S^1/,\gamma )`$ and $`(SO(3),\text{Haar})`$.

Proof:  That $`P^{}`$ is one-one off the zero section of $`S^2\times S^1/`$ is easily verified.

Fix the standard product metric on $`S^2\times S^1`$, normalized so that each factor has total volume 1. Normalized Haar measure on $`SO(3)`$ is Riemannian volume with respect to the bi-invariant metric we now describe. The Lie algebra of $`SO(3)`$ is the space of anti-symmetric matrices $`so(3)`$; on this algebra, we put the inner product:

$$=\frac{1}{2c^2}tr(AB^t),$$

where $`c=2\pi ^{2/3}`$. An orthonormal basis for $`so(3)`$ is $`\{X,Y,Z\}`$, where

$$X=\left(\begin{array}{ccc}0& 0& 0\\ 0& 0& c\\ 0& c& 0\end{array}\right),Y=\left(\begin{array}{ccc}0& 0& c\\ 0& 0& 0\\ c& 0& 0\end{array}\right),$$

and

$$Z=\left(\begin{array}{ccc}0& c& 0\\ c& 0& 0\\ 0& 0& 0\end{array}\right).$$

(Note $`[X,Y]=cZ`$, etc). In the bi-invariant metric induced by this inner product, $`SO(3)`$ has constant sectional curvatures, see \[DoC\]:

$$\lambda =\frac{c^2}{4}=\frac{[X,Y]^2}{4}=\frac{[X,Z]^2}{4}=\frac{[Y,Z]^2}{4}.$$

The diameter of $`SO(3)`$ in this metric is $`r=\pi /c`$ and the total volume is

$$\frac{\pi }{\lambda }\left(2r\frac{\mathrm{sin}(2\sqrt{\lambda }r)}{\sqrt{\lambda }}\right)=1.$$

Let $`\rho (x,\theta )d\mu (x)d\theta `$ be the pullback of the volume form on $`SO(3)`$ to $`S^2\times S^1`$ under $`P`$. To prove that $`P^{}`$ is an isomorphism, it suffices to show that $`\rho `$ is the density of $`\gamma `$, that is, to show that:

$$\rho (x,\theta )=(1\mathrm{cos}\theta )/\pi ,$$

for all $`(x,\theta )S^2\times S^1`$.

For any $`x,yS^2`$, if $`N=A\times \mathrm{\Theta }S^2\times S^1`$ is a product neighborhood of $`(x,\theta )`$ of Lebesgue measure $`\delta `$, then there exists a $`gSO(3)`$ such that $`\widehat{N}=gA\times \mathrm{\Theta }`$ is a neighborhood of $`(y,\theta )`$ of Lebesgue measure $`\delta `$. From the definition of $`P`$, it follows that $`P(\widehat{N})=gP(N)g^1`$. Since $`\rho (x,\theta )d\mu (x)d\theta `$ is the pullback of an $`SO(3)`$-invariant form, we obtain that for any $`x,yS^2`$, $`\theta S^1`$,

$$\rho (x,\theta )=\rho (y,\theta )=:\rho (\theta ).$$

Finally, we compute $`\rho (\theta )`$. Since the geodesics of $`SO(3)`$ through $`I`$ are precisely the one-parameter subgroups, the image under $`P^{}`$ of the curve $`t(p,t)`$ is a geodesic through $`I`$ of speed $`1/c`$. It follows that for any $`\theta (0,\pi )`$, $`P^{}`$ sends $`S^2\times (0,\theta )/`$ diffeomorphically onto $`B_{c^1\theta }(I)\{0\}`$, the punctured ball of radius $`c^1\theta `$ about the identity in $`SO(3)`$. The volume of such a ball is

$`{\displaystyle _0^\theta }\rho (t)𝑑t`$ $`=`$ $`\text{vol}(B_{c^1\theta }(I))`$
$`=`$ $`{\displaystyle \frac{\pi }{\lambda }}\left(2c^1\theta {\displaystyle \frac{\mathrm{sin}(2\sqrt{\lambda }c^1\theta }{\sqrt{\lambda }}}\right)`$
$`=`$ $`{\displaystyle \frac{1}{\pi }}(\theta \mathrm{sin}\theta ).`$
By the Fundamental Theorem of Calculus,

$$\rho (\theta )=\frac{1}{\pi }(1\mathrm{cos}(\theta )),$$
(9)
which completes the proof. QED

We call $`P`$ or more appropriately $`P^{}`$ polar coordinates on $`SO(3)`$.

### 5.2 Computing Random Exponents

Recall that $`R(\epsilon )=_{PS^2}\mathrm{log}Tf_\epsilon vdm(v).`$ Using the results of \[AvBo\] we can reduce the right hand integral to a one variable integral which we can then evaluate numerically very accurately. This is how the random Lyapunov exponents are computed.

###### Proposition 5.2

$$R(\epsilon )=_0^{\frac{1}{2}}\mathrm{log}(1+(2\pi \epsilon x(1x))^2)𝑑x.$$

Proof. Consider the “inverse Archimedean projection”

$$\mathrm{\Psi }:(\theta ,x)(g(x)\mathrm{cos}\theta ,g(x)\mathrm{sin}\theta ,x\frac{1}{2}),$$

where $`g(x)=\sqrt{x(1x)}`$. This map sends the cylinder $`C=S^1\times [0,1]`$ onto the sphere $`S^2`$ and is area-preserving: the pullback $`\mathrm{\Psi }^{}d\mu `$ is a multiple of the Lebesgue volume form $`d\theta dx`$ on $`C`$. The Riemannian metric on $`S^2`$ pulls back to the metric:

$$<v,w>_{(\theta ,x)}=v^tB(x)^2w,$$

where

$$B(x)=[D\mathrm{\Psi }(x,\theta )^tD\mathrm{\Psi }(x,\theta )]^{1/2}=\left(\begin{array}{cc}g(x)& 0\\ 0& (2g(x))^1\end{array}\right).$$

Setting $`\stackrel{~}{f}_\epsilon =\mathrm{\Psi }^1f_\epsilon \mathrm{\Psi }`$, we have that:

$$\stackrel{~}{f}_\epsilon (\theta ,x)=(\theta +2\pi \epsilon x,x).$$

We next compute:

$`R(\epsilon )`$ $`=`$ $`{\displaystyle _{PS^2}}\mathrm{log}Tf_\epsilon vdm(v).`$
$`=`$ $`{\displaystyle _{T_1C}}\mathrm{log}T\stackrel{~}{f}_\epsilon v\mathrm{\Psi }_{\mathrm{}}^{}dm(v).`$
In the second equation, the unit tangent bundle $`T_1C`$ and the quantity $`T\stackrel{~}{f}_\epsilon v`$ are defined with respect to the $`\mathrm{\Psi }`$-pullback Riemannian metric on $`C`$.

If $`vTC`$ is a unit vector with respect to the pullback metric, then $`u=Bv`$ is a unit vector with respect to the Euclidean metric, and $`T\stackrel{~}{f}_\epsilon v=BT\stackrel{~}{f}_\epsilon B^1u_{\text{Eucl.}}`$. Hence we can write:

$`R(\epsilon )`$ $`=`$ $`{\displaystyle \frac{1}{2\pi }}{\displaystyle _{(x,\theta )C}}{\displaystyle _{uS^1}}\mathrm{log}BT_{(x,\theta )}\stackrel{~}{f}_\epsilon B^1u_{\text{Eucl.}}dxd\theta du.`$
$`=`$ $`{\displaystyle \frac{1}{2\pi }}{\displaystyle _{(x,\theta )C}}{\displaystyle _{uS^1}}\mathrm{log}\left(\begin{array}{cc}1& 4\pi \epsilon x(1x)\\ 0& 1\end{array}\right)u_{\text{Eucl.}}dxd\theta du.`$
For $`ASL(2,𝐑)`$, \[AvBo\] show that

$$_{uS^1}\mathrm{log}Au_{\text{Eucl.}}du=\mathrm{log}((A_{\text{Eucl.}}+A_{\text{Eucl.}}^1)/2).$$

Applying this to $`A=\left(\begin{array}{cc}1& \alpha \\ 0& 1\end{array}\right)`$, we obtain that

$$_{uS^1}\mathrm{log}Au_{\text{Eucl.}}du=\frac{1}{2}\mathrm{log}(1+\alpha ^2/4).$$

It follows that:

$`R(\epsilon )`$ $`=`$ $`{\displaystyle _0^1}{\displaystyle _{uS^1}}\mathrm{log}\left(\begin{array}{cc}1& 4\pi \epsilon x(1x)\\ 0& 1\end{array}\right)u_{\text{Eucl.}}dudx`$
$`=`$ $`{\displaystyle \frac{1}{2}}{\displaystyle _0^1}\mathrm{log}(1+(4\pi \epsilon x(1x))^2/4)𝑑x`$
$`=`$ $`{\displaystyle _0^{1/2}}\mathrm{log}(1+(2\pi \epsilon x(1x))^2)𝑑x.`$
## 6 Experimental Results

In this section we describe several numerical experiments and their results. First we obtain experimental values of $`R(\epsilon )`$ for random maps as introduced in Section 2. We use a method similar to those we use later to compute estimates for $`\mathrm{\Lambda }(\epsilon )`$. Proposition 5.2 allows to us to check the accuracy of the computed estimate of $`R(\epsilon )`$ against the precise value of $`R(\epsilon )`$ given there. Then we pass to the computation of $`\mathrm{\Lambda }(\epsilon ).`$ Three different approaches to the computation of $`\mathrm{\Lambda }(\epsilon )`$ are presented and discussed which will allow us to obtain accurate enough values to draw conclusions. Finally, a sample of numerical estimates of $`\mathrm{\Lambda }(\epsilon )`$ is shown.

### 6.1 The case of random maps

To obtain experimental values of $`R(\epsilon )`$ for different $`\epsilon `$ we proceed as in Proposition 2.2 . A random point $`x`$ in $`S^2`$ and a random vector $`\xi `$ in $`T_{1,x}S^2`$ is chosen. A random sequence $`g_i`$ in $`SO(3)`$ is selected and the derivative of the maps $`g_if_\epsilon `$ are applied to the tangent vector $`\xi `$. The rate of increase of the logarithm of $`T_xf^{(n)}(\xi )`$ is described, where $`f^{(n)}=g_nf_\epsilon \mathrm{}g_1f_\epsilon .`$ The results are the same with probability 1. For brevity, we refer loosely to the use of formula (4).

Let us describe the selection of random elements:

* An initial point can be described in polar coordinates by a longitude $`\lambda _x`$ and a latitude $`\beta _x`$. The value of $`\lambda _x`$ is chosen at random in $`[0,2\pi ]`$ with uniform probability. Concerning $`\beta _x`$, a random value $`z[1,1]`$ is selected with uniform probability and then we let $`\beta _x=\mathrm{sin}^1z.`$ This gives the uniform probability for $`xS^2.`$
* A tangent vector $`\xi T_{1,x}S^2`$ is generated by choosing an angle $`\psi [0,2\pi ]`$ with uniform probability and letting $`\xi `$ make an angle $`\psi `$ with the unit tangent vector to the latitude through $`x`$ taken in the positive sense. We call this last vector the horizontal vector at $`x`$. It is not defined at the poles, but as the poles have measure zero this is irrelevant in the current context.
* A random rotation $`gSO(3)`$ is determined by an axis and an angle of rotation. As described in subsection 5.1 one can take these as elements in $`S^2\times S^1.`$ It is not necessary to carry out the identification described there. The axis is selected just as the point $`xS^2`$ was above. Let $`\theta `$ be the rotation angle. To select it according to formula (9) pick a random value $`z[0,2\pi ]`$ with uniform probability and solve the equation $`z=\theta \mathrm{sin}(\theta )`$ for $`\theta .`$ The equation is nothing other than the well known Kepler equation with eccentricity equal to 1. There are efficient solvers for it.

Then, given initial values of $`(x,\xi )`$ one can apply formula (4) to approximate $`R(x,\epsilon )`$ by using a finite number of iterates, $`N.`$ In turn, to approximate $`R(\epsilon )`$ the integral in formula (5) can be computed using a sample of size $`M`$ in $`PS^2.`$ Let $`R_{N,M}(\epsilon )`$ be the value obtained.

This value is compared to the one given by Proposition 5.2 which has been computed using a Simpson method with iterative mesh refinement (the values of $`R(\epsilon )`$ are shown in figure XXX). The following limit approximations are straightforward to derive:

$$\begin{array}{ccccc}\hfill R(\epsilon )& =& \frac{\pi ^2}{15}\epsilon ^2\frac{2\pi ^4}{315}\epsilon ^4+O(\epsilon ^6)\hfill & \text{for}\hfill & \epsilon 0,\hfill \\ & & & & \\ \hfill R(\epsilon )& =& \mathrm{log}(2\pi \epsilon )2+\frac{1}{2\epsilon }+O(\frac{\mathrm{log}(\epsilon )}{\epsilon ^2})\hfill & \text{for}\hfill & \epsilon \mathrm{}.\hfill \end{array}$$
(12)
Skipping the $`O`$ terms one has relative errors less than $`0.01`$ if $`\epsilon <0.30`$ and $`\epsilon >3.19`$, respectively.

Let

$$\mathrm{\Delta }_{N,M}(\epsilon )=R_{N,M}(\epsilon )R(\epsilon ).$$
(13)
We have verified experimentally that $`\mathrm{\Delta }_{N,M}(\epsilon )`$ has essentially zero average and a standard deviation like

$$\sigma _{N,M}(\epsilon )\frac{\kappa (\epsilon )}{\sqrt{MN}},$$

provided $`N`$ is large enough.

Tests have been done for several choices of $`N,M`$ and $`\epsilon .`$ Using $`M=10^3`$ and $`N=10^k,k=2,\mathrm{},6,`$ the values of $`\kappa (\epsilon )`$ have been estimated for $`\epsilon `$ ranging from $`10^1`$ to $`10^3.`$ There are no significant differences from $`k=4`$ on. Figure xxx displays, for different values of $`\epsilon `$, the interval $`[3\kappa (\epsilon ),3\kappa (\epsilon )]`$ and the results of single runs (i.e., taking $`M=1`$) for $`N=10^k,k=4,\mathrm{},9.`$ More concretely, the plotted values are the deviations $`\mathrm{\Delta }_{N,1}(\epsilon )`$ given by formula (13) multiplied by $`\sqrt{N}.`$

These results indicate that $`R_{N,M}(\epsilon )`$ and $`R(\epsilon )`$ agree to order $`10^5`$ taking $`NM=10^{10}.`$ Further checks have been done for larger values of $`\epsilon `$ (up to $`\epsilon =10^6`$) which show no significant variation of $`\kappa (\epsilon )`$ between $`\epsilon =10^3`$ and $`\epsilon =10^6.`$

### 6.2 Computing the Lyapunov exponent in the deterministic case

In principle one can follow a similar scheme to compute $`\mathrm{\Lambda }(\epsilon )`$. That is, using (1) with a finite number of iterates, $`N`$, an approximation of $`\lambda _1(x,gf_\epsilon )`$ is obtained. Then, (2) is computed using a Montecarlo method sampling $`xS^2`$ as described with samples of size $`M_p.`$ Finally an estimate of $`\mathrm{\Lambda }(\epsilon )`$ is computed by applying again a Montecarlo method to (3), sampling $`gSO(3)`$ as explained above and using samples of size $`M_r.`$ In any case the samples are taken according to the appropriate measures. The total number of iterates of the maps and their differential is, hence, $`M_rM_pN.`$ Let us denote by $`\mathrm{\Lambda }_{M_r,M_p,N}(\epsilon )`$ a value obtained in this way.

Results of this approach are shown in figure xxx for different values of $`\epsilon .`$ They require some explanation. For a fixed $`gSO(3)`$ the values of $`\lambda _1(x,gf_\epsilon )`$ are estimated for $`M_p`$ random values of $`x.`$ The standard deviation of the values of $`\lambda _1(x,gf_\epsilon )`$ is then computed. This value, $`\sigma _g,`$ depends of the choice of $`g.`$ Let $`\sigma _{S^2}`$ be the average value of $`\sigma _g`$ when a full sample of $`gSO(3)`$ is considered. On the other hand all the $`M_rM_p`$ determinations of $`\lambda _1(x,gf_\epsilon )`$ can be used to estimate a global standard deviation, $`\sigma _{\text{total}}.`$

It is clear that $`\sigma _{S^2}`$ measures the average dispersion of the maximal Lyapunov exponent when different points are taken in the phase space. The dispersion depends on the concrete rotation $`g`$ taken. Typically, for the $`g`$ such that relatively small values of the average $`\lambda (gf_\epsilon )`$ of $`\lambda _1(x,gf_\epsilon )`$ on $`S^2`$ are obtained, it is seen that $`\sigma _g`$ is larger. It should also be mentioned that the errors in the determination of $`\lambda _1(x,gf_\epsilon )`$, due to the finiteness of the number of iterates $`N`$, also contribute to this dispersion.

On the other hand there is also a dispersion in $`\lambda (gf_\epsilon )`$ when different $`g`$ are taken. The standard deviation $`\sigma _{\text{total}}`$ measures the cumulative effect of both dispersions.

Now let us make several comments on the observed behavior, based on computations carried out with quite different values of $`M_r,M_p`$ and $`N.`$

* The estimates of $`\mathrm{\Lambda }(\epsilon )`$ are close to $`R(\epsilon )`$ for $`\epsilon `$ large. For instance, for $`\epsilon >3`$ one has that $`|\mathrm{\Lambda }_{5000,5000,8192}(\epsilon )R(\epsilon )|`$, in the runs done, is below $`0.00325`$. The situation is worse for small $`\epsilon `$ because the difference reaches the value $`0.105`$ for $`\epsilon =0.55.`$ If we proceed to compute the relative error (r.e.) $`=(\mathrm{\Lambda }_{5000,5000,8192}(\epsilon )R(\epsilon ))/R(\epsilon )`$ the agreement is even worse. For $`\epsilon =1`$ one has r.e. $`0.14`$ and r.e. $`<0.9`$ if $`\epsilon <0.42.`$ In these comparisons one should take into account that the discrepancies also include the errors done in the estimates of $`\mathrm{\Lambda }(\epsilon )`$ (see Sections 6.3 and 6.4). But the present results already indicate that the differences for $`\epsilon `$ small are not only due to statistical errors. To make this more evident some additional computations have used a total number of iterates (for some selected $`\epsilon `$) with $`M_rM_pN`$ largely exceeding $`10^{12}.`$

Furthermore it seems also clear that for $`\epsilon <0.3`$ there is a “saturation” in the behavior of the estimates of $`\mathrm{\Lambda }_{M_r,M_p,N}(\epsilon )`$ and of the standard deviations. Indeed, it can be seen that the trend on the figure changes completely (this is also the purpose to use logarithmic scales). Systematic errors occur which completely invalidate the statistical results. To make this more evident some values, for $`\epsilon `$ in the range $`[0.2,0.3]`$, computed also in a probabilistic way but with a different estimator (see Section 6.3) are also shown as dots in the lower left part. For these computations $`M_r=14400,M_p=16384`$ and $`N=16384`$ have been used.

The maximal value $`\lambda _{1,\text{max}}(\epsilon )`$ of $`\lambda _1(x,gf_\epsilon )`$ for $`xS^2`$ and $`gSO(3)`$ is larger than $`R(\epsilon ).`$ This requires samples with $`M_r,M_p`$ large to be detected if $`\epsilon `$ is small. It will be clear from Section 7 and the upper formula in (12) that the quotient $`q(\epsilon )=\lambda _{1,\text{max}}(\epsilon )/R(\epsilon )`$ is unbounded if $`\epsilon 0.`$ If large values of $`\epsilon `$ are considered, it is observed that $`q(\epsilon )`$ tends slowly to 1 when $`\epsilon \mathrm{}.`$ In fact it follows from the lower formula in (12) and the analysis in Section 7 (which is partly valid for any $`\epsilon `$) that the difference $`\lambda _{1,\text{max}}(\epsilon )R(\epsilon )`$ is bounded by $`2\mathrm{log}(2)+O(\epsilon ^1).`$ To see differences close to the bound one has to use very large values of $`M_r,M_p`$.
* The estimates of $`\sigma _{S^2}`$ are mildly sensitive to the concrete values of $`M_r,M_p`$ and $`N`$, provided these values are not too small, and assuming $`0.3<\epsilon <3.`$ For $`\epsilon >30`$ a clear dependence with respect to $`N`$, of the form $`N^{1/2}`$, is seen. From $`\epsilon =3`$ to $`\epsilon =30`$ there is a gradual increase in the dependence with respect to $`N.`$ For $`\epsilon <0.3`$ a tendency towards a behavior of the form $`N^1`$, which increases when $`\epsilon `$ decreases, is clear.

If extremal values of $`\sigma _{S^2}`$ are considered when a sample of $`g`$ is taken, it is clear that the minimum must be zero. But there is a significant difference for $`\epsilon <2,`$ because the minimum is already close to zero for samples of moderate size, while for larger $`\epsilon `$ the minimum, which is almost insensitive to $`\epsilon `$ goes to zero slowly when $`M_r`$ increases. On the other hand, the maximum of the observed values of $`\sigma _{S^2}`$ increases until $`\epsilon 10.`$ It only stabilizes to a value with small dependence on $`\epsilon `$ for $`\epsilon 100.`$
* The values of $`\sigma _{\text{total}}`$ are much larger than $`\sigma _{S^2}`$ for $`\epsilon >1`$, while for $`\epsilon <0.3`$ they are essentially equal. In particular they have a small dependence with respect to $`M_r,M_p`$ and $`N`$, if these are not too small and assuming $`\epsilon >0.3.`$

An analysis of the reasons of the observed behavior is useful because it helps in three different aspects: a) to understand the different contributions to the errors in the estimates of $`\mathrm{\Lambda }(\epsilon )`$; b) to see the main differences between the cases of random and deterministic maps; c) to suggest alternative methods to obtain better estimates.

* In the deterministic case the dynamics on $`S^2`$ for a given rotation $`g`$ is relevant. This dynamics is “destroyed” (or “smoothed”, “averaged”) in the case of random maps. Hence, the initial point is irrelevant and the estimates improve in a probabilistic way, depending on the total number of iterates $`MN`$, in the random case.
* Changing $`g`$ in the deterministic case produces dramatic changes in the dynamics and, hence, on $`\lambda (gf_\epsilon ).`$ This is specially clear for $`g`$ close to the identity (axis close to the pole or small rotation) or for rotations of angle very close to $`\pi `$ around an axis of small latitude. The variability of $`\lambda (gf_\epsilon )`$ with respect to $`g`$ is a major source of dispersion in the results, specially for large $`\epsilon .`$ A standard deviation $`\sigma _{\text{total}}`$ around $`0.4`$, mainly due to the variation of $`\lambda (gf_\epsilon ),`$ would require sampling with $`M_r`$ of the order of $`10^8`$, at least, to have accurate results.
* For a fixed $`g`$, changing $`xS^2`$ has a very mild effect for large $`\epsilon .`$ Despite the possible existence of tiny islands (see Section 7) the dynamics ”looks ” ergodic. A similar behavior has been observed for standard–like maps in \[GLS\] and in the case of volume preserving flows it is seen in \[NSV\], where an analysis of the places where the islands should be expected is carried out before finding them explicitly.

Figure xxx displays a sample of orbits in $`S^2`$ for fixed $`g`$ and two different values of $`\epsilon `$. For $`\epsilon =0.3`$ the dynamics is mainly dominated by an integrable behavior, with many invariant curves and small chaotic seas, the largest one seen in the front part. This is persistent with respect to changes in $`g.`$ The system is even more integrable (that is, invariant curves fill up a larger area) for most of the rotations $`gSO(3).`$ For $`\epsilon =2`$ (which is not so large!) only minor islands subsist, and they can even disappear for different $`g.`$ For values like $`\epsilon =10`$ it is hard to see any island unless $`g`$ is selected on a set of small measure. In the random case one observes a uniform distribution of iterates in $`S^2`$ and the same is essentially true for large $`\epsilon `$ in the deterministic case.
* On the other hand, for fixed $`g`$ and small $`\epsilon `$ the value of $`\lambda _1(x,gf_\epsilon )`$ depends strongly on $`x`$. But the behavior is typically rather sharp. Either one obtains a moderate value of the order of $`\epsilon ^{1.5}`$ or it is zero. The smaller the value of $`\epsilon ,`$ the larger the measure of the $`x`$ with exponent zero, of course. The average value can be very small and despite the standard deviation $`\sigma _{S^2}`$ also being small, large samples with respect to $`x`$ have to be taken if small relative errors are desired.
* The worst point concerning accuracy, especially for small $`\epsilon `$, are the errors in the computation of $`\lambda _1(x,gf_\epsilon ).`$ Indeed, for an integrable motion (e.g., $`x`$ in an invariant curve) $`\lambda _1`$ is zero, but the estimates $`\frac{1}{n}\mathrm{log}T_xf^n`$ are, generically, of the order of $`\frac{\mathrm{log}(n)}{n}.`$ This implies that the convergence to zero is slow. For large $`\epsilon `$ this effect is relatively not so dramatic, but oscillations in the behavior of the quotients and different trends can be expected.

### 6.3 Improved procedures

Several alternative procedures have been used in previous computations to determine the maximal Lyapunov exponent of a given map, averaged on the phase space:

* If the dynamics has a uniform hyperbolicity but with superimposed strong periodic or quasiperiodic oscillations, the following strategy has been used in \[BS\]. It consists in detecting, by an iterative procedure, an upper envelope of the plot of the quotients $`\frac{1}{n}\mathrm{log}T_xf^n`$ as a function of $`n.`$ Then, and after skipping a transient regime, one fits a function of the form $`\alpha +\beta /n`$ to the envelope. The value of $`\alpha `$ is a good estimator for $`\lambda _1(x,f)`$ and, as the system in \[BS\] is a skew product with linear action on the fibers, the value of $`x`$ is irrelevant.
* If the values of $`\lambda _1(x,f)`$ depend strongly on $`x`$, it is possible to divide the phase space in pixels of a given size (in general, $`d`$-dimensional pixels) and start the computations at a point in each pixel. However, if the number of pixels is large and the system depends on several additional parameters, the method can be prohibitive from a computational point of view. Then, together with each initial pixel one considers all the pixels visited by the orbit. The estimated Lyapunov exponent is assigned to all of them. One requires each pixel to be visited a minimal number of times (in case of need one takes several initial points in the pixel) and an averaged Lyapunov exponent is assigned to the pixel. Later on this is averaged over the full phase space.

This method has been used in \[SS\] to study the classical Hill’s problem and how the degree of chaos behaves with respect to the energy.

One should also take into account the stickiness of invariant curves. An initial point in a chaotic sea can remain close to an island for a large number of iterates. Hence, it is a good strategy to take a larger number of initial points even if one has to decrease the number of iterates for each one, provided this number is not too small. Furthermore, the local slope of $`\mathrm{log}T_xf^n`$ can have quite different trends if the number of iterates is large. Statistically this is not a problem because the interesting magnitude is the average behavior.
* The previous method still suffers from slow convergence of the quotients to $`\lambda _1(x,f)`$. An alternative method has been used in the context of flows with applications to galactic potentials in \[CiSi\] and later on extended to discrete transformations in \[CGS\], where references to other applications can be found. It is mainly intended to discriminate between regular and chaotic motion (that is, to decide if one can accept $`\lambda _1(x,f)=0`$ or not), but it also supplies an estimate of the Lyapunov exponent.

Given a map $`f`$, an initial point $`x`$ on a manifold $``$ and a random vector $`\xi T_{1,x}`$, let $`\xi _0=\xi `$ and define $`\xi _k=(T_{f^{k1}(x)}f)(\xi _{k1}).`$ For fixed integers $`m`$ and $`n`$ and $`N>0`$ let

$$Y_{m,n}(N)=N^n\underset{k=1}{\overset{N}{}}\mathrm{log}\left(\frac{\xi _k}{\xi _{k1}}\right)k^m,\overline{Y}_{m,n}(N)=\underset{k=1}{\overset{N}{}}Y_{m,n}(k).$$

Then, for a chaotic orbit the estimator

$$\widehat{Y}_{m,n}(N)=(m+1)(m+n+2)\frac{\overline{Y}_{m,n}(N)}{N^{n+m+2}}$$
(14)
tends to $`\lambda _1(x,f)`$, while for a generic regular orbit it behaves like $`\frac{(m+1)(m+n+2)}{m(m+n+1)}\frac{1}{N}.`$ The basic idea is to average the exponential rate of increase of the length of $`\xi _k`$ so that the transience has small relevance and to smooth out the irregularities of the quotients. Hence, it is a measure of the mean exponential growth of nearby orbits (MEGNO) and depends on the couple $`(m,n).`$ Suitable values (according to numerical experience) are $`m=2,n=0,`$ and then it is denoted as MEGNO20.

For $`(m,n)=(2,0)`$ a slightly better estimator for $`\lambda _1(x,f)`$ is obtained by using $`12\frac{\overline{Y}_{2,0}(N)}{N^4+4N^3+5N^2}.`$ Furthermore, when this method is used with these $`(m,n)`$, one can check for a behavior of the form $`\frac{2}{N}`$ to decide $`\lambda _1(x,f)0.`$ Typically $`\widehat{Y}_{2,0}(N)\frac{2}{N}=O(N^2)`$ for regular orbits.

Given a maximal number of iterates $`N_{\text{max}}`$ to be used in the estimates, an additional question is whether it can be better to use another value $`N<N_{\text{max}}`$ as a better choice to estimate $`\lambda _1(x,f)`$. In \[CGS\] a “right stop” criterion is introduced. It is specially relevant if the orbit is close to be regular, to prevent an overestimate of $`\lambda _1(x,f)`$, but it has not been used in the present computations.

Figure xxx illustrates the different behavior of MEGNO20 and the quotients in formula (1) in several cases. Details on the parameters used for the plots are as follows.

On the upper row, left plot, for $`\epsilon =0.3`$ a rotation of angle $`\theta `$ and axis of latitude $`\beta `$ with $`\frac{\beta }{2\pi }=\frac{3}{25}g,\frac{\theta }{2\pi }=\frac{18}{25}g`$ where $`g=(\sqrt{5}1)/2`$ has been selected. Random initial conditions are chosen. After a transient of 512 iterates, estimates of $`\lambda _1(x,f)`$ are produced and plotted for the next 2048 iterates. Solid lines (the lower ones) correspond to the estimates using MEGNO20, while discontinuous lines are produced by the classical formula. The middle part displays a similar plot for different initial conditions and the right part is similar to the middle one but for $`\epsilon =10.`$

On the lower row estimates of the largest Lyapunov exponent for random initial conditions are shown. On the left part $`\epsilon =0.3`$ and $`\theta `$ and $`\beta `$ as before have been used, and 512 random initial points are plotted. Solid (resp. discontinuous) lines correspond to MEGNO20 (resp. classical) estimators. Transient and number of iterates are as before. The left (middle) plot in the upper row corresponds to number 4 (43) of these points. The middle part is similar but for $`\epsilon =3.`$ For completeness estimates in the random case using finite $`n`$ in (1) are also shown on the right for $`\epsilon =0.3.`$ Each iterate uses a random element in $`SO(3)`$ with density $`\widehat{\rho }(\beta ,\theta )=\mathrm{cos}(\beta )(1\mathrm{cos}(\theta ))/(2\pi ).`$ The solid line shows the value of $`R(0.3)0.0547518.`$

See also \[S2\] for additional methods and applications.

Due to the good properties of the procedure, estimates of $`\lambda _1(x,gf_\epsilon )`$ have been computed using procedure 3) above. For the integrations in $`S^2`$ a Montecarlo method has been used. This is good enough for large $`\epsilon `$, because of the mild dependence of $`\lambda _1(x,gf_\epsilon )`$ with respect to $`x`$ (for, say, $`\epsilon >1`$) for most of the $`xS^2`$ and most of the $`gSO(3).`$ For $`\epsilon `$ small and especially if $`\epsilon <0.3,`$ a method such as the one presented in 2) would be suitable, but there are additional problems, due to the smallness of the Lyapunov exponents, to be discussed later. Furthermore it will turn out that it is relevant to compute $`\mathrm{\Lambda }(\epsilon )`$ with a small relative error for large $`\epsilon ,`$ to allow for a careful comparison with $`R(\epsilon ).`$ But for $`\epsilon `$ small it will be clear from the results, even those obtained with a moderate accuracy, that $`\mathrm{\Lambda }(\epsilon )`$ is far below $`R(\epsilon ).`$ In any case, it seems that numerical estimates of $`\mathrm{\Lambda }(\epsilon )`$ for $`\epsilon <0.3`$ with small relative error require an enormous computational effort.

Finally, for the integrations in $`SO(3)`$ and taking into account the large standard deviation found for moderate and large values of $`\epsilon `$, it has been found more convenient to use numerical quadrature formulas based on a grid of points. More concretely, a product Simpson method has been used with respect to the latitude $`\beta _g`$ of the axis of rotation and the rotated angle $`\theta _g.`$ The longitude of the axis is irrelevant: changes in this longitude are equivalent to changes in the longitude $`\lambda _x`$ of $`xS^2.`$ Using a grid with $`\theta _g[0,2\pi ],\beta _g[0,\pi /2]`$ requires that the estimate of $`\lambda (gf_\epsilon )`$ be multiplied by the factor $`\mathrm{cos}(\beta _g)(1\mathrm{cos}(\theta _g)).`$

Initial estimates for the results shown in the next Section use, for the elements in $`SO(3)`$ the following data: $`\lambda =0`$, $`(\theta ,\beta )`$ on a grid of $`N_g\times N_g`$ with $`N_g=200`$. Then a sample of $`N_p=1600`$ random initial points and tangent vectors in $`PS^1`$ and $`M=8000`$ iterates are used. The programs have been done in such a way that beyond the estimates for these values, also estimates using grids with $`N_g=100,\mathrm{\hspace{0.17em}50}`$ or using samples with $`N_p=800,\mathrm{\hspace{0.17em}400},\mathrm{\hspace{0.17em}200},\mathrm{\hspace{0.17em}100}`$ and doing a number of iterates equal to $`M=4000,\mathrm{\hspace{0.17em}2000},\mathrm{\hspace{0.17em}1000},\mathrm{\hspace{0.17em}500}`$ are computed. This allows for a check on the internal consistency of the results.

It turns out that the use of different grids in $`SO(3)`$ stabilizes quickly. Concerning the dependence with respect to $`N_p`$ and $`M`$, it is clearly seen that there is no need for very large values of $`N_p`$ except in the case of small $`\epsilon `$ and one is interested in having small relative error. The dependence with respect to $`M`$ is clearly of the form ctant/$`M.`$ Hence, extrapolations with respect to $`M`$ have been used. The initial estimates allow for a fine tuning of the most suitable values for the grid, $`N_p`$ and $`M`$. For instance, assuming that one can accept a total of $`2^{41}`$ iterates (for every value of $`\epsilon `$), for large $`\epsilon `$ a typical choice is $`N_g=2^8,N_p=2^9,M=2^{16},`$ while for small $`\epsilon `$ it is $`N_g=2^7,N_p=2^{10},M=2^{17}.`$ Even with this large $`M`$ the results start to be not very good if $`\epsilon <0.2.`$

Figure xxx shows 3D views of $`h(\theta _g,\beta _g,\epsilon )=\lambda (gf_\epsilon )\mathrm{cos}(\beta _g)(1\mathrm{cos}(\theta _g))\frac{\pi }{2}`$ as a function of $`(\theta _g,\beta _g)`$ for different values of $`\epsilon .`$ Level lines of these surfaces are displayed in figure xxx. The plots give a good evidence of the smooth behavior of $`h(\theta _g,\beta _g,\epsilon )`$ for moderate and large values of $`\epsilon `$, and how the behavior becomes wilder, with sharp changes for small $`\epsilon .`$ It is clear that the results are the same if $`(\theta ,\beta )`$ is replaced by $`(2\pi \theta ,\beta )`$. Furthermore, note that for integer $`\epsilon `$ the symmetries of $`gf_\epsilon `$ imply that the results should be the same if $`\theta `$ is changed to $`2\pi \theta `$, as clearly seen in the first three plots. For $`\epsilon =0.5`$ the lack of symmetry $`\theta 2\pi \theta `$ is clear, but for large non–integer $`\epsilon `$ this lack of symmetry is harder to detect.

### 6.4 A sample of results

The results of applying the methodology just described are shown in figure xxx. Typical values for the number of iterates, initial data and grid have been given before.

On the top plots general views of $`R(\epsilon )`$ and estimates of $`\mathrm{\Lambda }(\epsilon )`$, to be denoted as $`\mathrm{\Lambda }(\epsilon )_{\text{num}}`$, can be seen. In particular, for large and small $`\epsilon `$ it is easy to check the limit behavior of $`R(\epsilon )`$ predicted by (12). On top left only values $`\epsilon 10`$ are shown. On this scale no differences can be seen between $`R(\epsilon )`$ and $`\mathrm{\Lambda }(\epsilon )_{\text{num}}`$ for $`\epsilon >3.`$ For small $`\epsilon `$ the differences are clear and they are quite dramatic for $`\epsilon <0.4.`$

On the bottom left plot the tiny differences $`\mathrm{\Lambda }(\epsilon )_{\text{num}}R(\epsilon )`$ are displayed. It seems that they tend to 0 as $`\epsilon \mathrm{},`$ in agreement with the second possibility in Section 3. Finally the behavior of $`\mathrm{\Lambda }(\epsilon )_{\text{num}}`$ for $`\epsilon `$ small is seen in detail on the bottom right part. A logarithmic scale has been used to reveal that $`\mathrm{log}(\mathrm{\Lambda }(\epsilon )_{\text{num}})`$ is dominated by a function of the form $`c/\epsilon `$ for some $`c>0.`$ Due to the smallness of the estimates and to the fact that it would require an enormous effort to estimate $`\mathrm{\Lambda }(\epsilon )`$ for $`\epsilon `$ close to $`0.1`$ (unless other methods, from deterministic analysis, are used), it is not completely clear what the correct behavior is. Even with a reduced set of data a fit of the values obtained for $`0.16\epsilon 0.3`$ gives a result of the form

$$\mathrm{log}(\mathrm{\Lambda }(\epsilon )_{\text{num}})2.45\frac{3.16}{\epsilon }$$
(15)
which must be taken with caution, but seems to give the correct trend. This suggests that the inequality in Question 1.1 is not satisfied for small $`\epsilon `$, which will be confirmed theoretically in Section 7. This fact is not a surprise, because similar facts occur in generic analytic families of area-preserving diffeomorphisms. The smallness of the Lyapunov exponent is related to the area of the chaotic seas which in turn is related to the splitting of the separatrices of fixed and periodic point for maps close to the identity. See \[FS\], where general upper estimates can be found. In fact, as with many a priori exponentially small upper bounds, this result can also be obtained as a corollary of averaging theory for analytic systems, see \[N\] and also Section 7.2.

Finally it should be mentioned that some computations have been done for large values of $`\epsilon `$ (up to $`10^6`$). Due to the strong chaotic properties it is enough to take small values (say $`N_p=256`$) of the number of initial points in $`S^2`$. But the grid in the parameters $`(\theta _g,\beta _g)`$ has to contain more points. Typical values of $`N_g`$ to have a good determination of $`\mathrm{\Lambda }(\epsilon )_{\text{num}}`$ are $`2^9`$ and $`2^{10}`$. The results confirm what is seen in the left lower part of figure xxx, that is, $`\mathrm{\Lambda }(\epsilon )_{\text{num}}>R(\epsilon )`$ and the difference goes to zero slowly.

## 7 The case of small $`\epsilon `$

As it is clear that the greatest problems occur for small $`\epsilon `$, it is worth it to carry out a preliminary analysis of the dynamics in this case. The first item to be studied is the location and stability of fixed points. This can be carried out, with the same effort, for any $`\epsilon `$. Furthermore this allows us to see how bifurcations give rise to new elliptic fixed points with the corresponding creation of islands. Later on the global behavior of $`gf_\epsilon `$ on $`S^2`$ is discussed. In what follows it is assumed that $`\epsilon >0.`$

### 7.1 Fixed points and their stability

To look for fixed points of $`gf_\epsilon `$ it is enough to consider axes of rotation with zero longitude. Let $`(\beta ,\theta )`$ be the latitude of the axis and the angle of the rotation. As follows from Section 5.1, it is not restrictive to assume $`\beta [0,\pi /2].`$ Then a fixed point $`A`$ is mapped by $`f_\epsilon `$ to a point $`A^{}`$ which by $`g`$ returns to $`A`$. Let $`b`$ be the latitude of $`A`$. It is clear that $`A`$ and $`A^{}`$ must have symmetric longitudes, $`\delta (b)`$ and $`\delta (b)`$, respectively, where $`\delta (b)=\frac{\pi }{2}\epsilon (1+\mathrm{sin}(b)).`$ It is easy to derive the condition for the fixed points

$$\mathrm{sin}\left(\frac{\theta }{2}\right)\mathrm{sin}(\beta )\mathrm{cos}(b)\mathrm{cos}(\delta (b))\mathrm{sin}\left(\frac{\theta }{2}\right)\mathrm{cos}(\beta )\mathrm{sin}(b)+\mathrm{cos}\left(\frac{\theta }{2}\right)\mathrm{cos}(b)\mathrm{sin}(\delta (b))=0.$$
(16)
For $`\epsilon `$ small one has $`\mathrm{sin}(\delta (b))=O(\epsilon )`$ and $`\mathrm{cos}(\delta (b))=1O(\epsilon ^2).`$ For the analysis of this case it is more convenient to write (16) in the form

$`\mathrm{sin}\left({\displaystyle \frac{\theta }{2}}\right)\mathrm{sin}(\beta b)`$ $``$ $`\epsilon ^2{\displaystyle \frac{\pi ^2}{4}}\mathrm{sin}\left({\displaystyle \frac{\theta }{2}}\right)\mathrm{sin}(\beta )\mathrm{cos}(b)(1+\mathrm{sin}(b))^2`$ (17)
$`+`$ $`\epsilon {\displaystyle \frac{\pi }{2}}\mathrm{cos}\left({\displaystyle \frac{\theta }{2}}\right)\mathrm{cos}(b)(1+\mathrm{sin}(b))+O(\epsilon ^3)=0.`$
From (17) it is clear that if $`\mathrm{sin}\left(\frac{\theta }{2}\right)=O(1),`$ that is, $`\theta `$ is not too close to 0 from the positive or the negative side, then one has solutions for $`b`$ of the form

$$b=\beta +\gamma \epsilon +O(\epsilon ^2)\text{or}b=\pi +\beta +\gamma \epsilon +O(\epsilon ^2),$$

where $`\gamma `$ is independent of $`\epsilon .`$ The value for $`\gamma `$ is given by

$$\gamma =\frac{\pi }{2}\mathrm{cos}(\beta )(1\pm \mathrm{sin}(\beta ))/\mathrm{tan}\left(\frac{\theta }{2}\right),$$

where the $`+`$ sign is used in the first case and the $``$ sign in the second. In both cases $`b`$ is close to either $`\beta `$ or $`\beta +\pi `$ and an important thing is that there are exactly two fixed points for $`gf_\epsilon .`$

Otherwise one can write $`\theta =m\pi \epsilon ,m>0`$ or $`\theta =2\pi +m\pi \epsilon ,m<0.`$ The dominant terms in (17) become in both cases

$$m\mathrm{sin}(b\beta )\mathrm{cos}(b)(1+\mathrm{sin}(b))=0.$$
(18)
Equation (18) has to be seen as an equation for $`b`$ depending on $`\beta `$ and $`m`$ (which accounts for $`\theta `$). As it has zero average it should have at least two different zeros. To look for more solutions it is relevant to compute the lines (in $`(\beta ,m)`$) where double zeros occur. The angle $`b`$ can be used to parameterize these lines. They are represented as

$$\begin{array}{ccc}\hfill m^2& =& 2+\mathrm{sin}^3(b)\frac{3}{4}\mathrm{sin}^2(2b),\hfill \\ & & \\ \hfill \beta & =& b\mathrm{arg}\left[\mathrm{sin}(b)\mathrm{cos}(2b)\sqrt{1}\left(\mathrm{cos}(b)+\frac{1}{2}\mathrm{sin}(2b)\right)\right].\hfill \end{array}$$
(19)
As $`\beta [0,\pi /2],`$ inspection of (18) shows that no double zeros can occur in the case $`m>0.`$ Hence the value of $`m`$ is confined to $`[2,0].`$ The bounds on $`\beta `$ also imply that the parameter $`b`$ in (19) has to be in $`[\pi /2,3\pi /2].`$ It is elementary to discuss the behavior of $`(m,\beta )`$ as a function of $`b`$. It is better seen by looking at figure xxx. It can just be said that two curves of double zeros appear with $`\beta [\pi /4,\pi /2].`$ They meet at $`m=\sqrt{2},\beta =\pi /4,`$ where a triple zero appears. Between both curves there are exactly four zeros.

When additional powers of $`\epsilon `$ are included, a routine application of the Implicit Function Theorem permits us to conclude the same behavior for the full equation (16). It should be noted that $`\beta =\pi /2`$ has to be excluded from the previous analysis: in that case the axis of rotation is also the axis of the twist.

The preceding analysis can be summarized as

###### Proposition 7.1

For $`\epsilon `$ small enough and any $`gSO(3)`$ there are always at least two fixed points of $`gf_\epsilon .`$ Bifurcations to exactly four fixed point appear for any longitude of the rotation axis and for latitude of the axis and angle of rotation $`(\beta ,\theta =\pi (2m\epsilon +O(\epsilon ^2)))`$ along two lines described by formulae (19) when the parameter $`b`$ ranges in $`(\pi /2,\pi )`$ and $`(\pi ,3\pi /2)`$, respectively.

To discuss bifurcations of the fixed points for general values of $`\epsilon `$ is an elementary but cumbersome task. As an illustration the case of bifurcations appearing on $`\beta =0`$ is presented. Then (16) reduces to

$$\mathrm{cos}\left(\frac{\theta }{2}\right)\mathrm{cos}(b)\mathrm{sin}(\delta (b))\mathrm{sin}\left(\frac{\theta }{2}\right)\mathrm{sin}(b)=0$$
(20)
and the condition for a double root becomes

$$\frac{2}{\pi \epsilon }\mathrm{tan}(\frac{\pi \epsilon }{2}(1+\mathrm{sin}(b))=\mathrm{sin}(b)\mathrm{cos}^2(b).$$
(21)
The degenerate cases $`|b|=\pi /2`$ must be excluded in (21). It is immediate that new double fixed points appear on $`S^2`$ with $`\beta =0`$ if and only if $`\epsilon `$ is a positive integer. The number of double fixed points with $`\beta =0`$ (and some $`\theta `$) increases with $`\epsilon .`$ Also from (20) it follows that new zeros appear near $`\theta =0,`$ one for $`\theta >0`$ and the other for $`\theta <0.`$ These zeros move towards $`\theta =\pi `$ without ever reaching it. So, it is a simple matter to state how many fixed points exist for $`\beta =0`$ (except at the bifurcation values of $`\epsilon `$): there are at most $`2(1+E(\epsilon ))`$, where $`E`$ denotes the integer part of $`\epsilon .`$ For a given non-integer $`\epsilon `$ there are always values of $`\theta `$ such that this number $`2(1+E(\epsilon ))`$ is the exact number of fixed points. This has an elementary dynamical interpretation: new fixed points emanate from the north pole of $`S^2`$ when the rotation number of $`f_\epsilon `$ at the north pole (defined by continuity) passes through 0 (mod 1).

To study the stability of the fixed points we recall that they are generically elliptic (eigenvalues $`\mu `$ in $`S^1\{\pm 1\}`$), hyperbolic (real positive eigenvalues) and hyperbolic with reflection (real negative eigenvalues). Let $`E,H,R`$ denote the number of fixed points of each type. Euler–Poincaré formula gives $`EH+R=2`$ (for simple fixed points). At the creation of new fixed points $`E`$ and $`H`$ increase by 1. When double eigenvalues are equal to $`1`$ then, generically, $`E`$ decreases by 1 and $`R`$ increases by 1.

An analytic discussion on the stability of the fixed points is elementary (at least for small $`\epsilon `$) but cumbersome. It is worth mentioning that, for any $`\epsilon `$ the maximum eigenvalue at a fixed point is achieved on $`b=0`$ and has the expression

$$\mu _{\text{max}}=\frac{\pi \epsilon }{2}+\left(1+\left(\frac{\pi \epsilon }{2}\right)^2\right)^{1/2}.$$
(22)
A sample of illustrations is shown in figure xxx having $`\theta /2\pi `$ as horizontal variable and $`\beta /\pi `$ as vertical one. A region containing $`i`$ (resp. $`j,k`$) fixed points of type $`E`$ (resp. $`H,R`$) is denoted as $`R^kH^jE^i.`$ On the top left plot and for $`\epsilon =0.1`$ the two curves on the upper part of the plot are the curves of double zeros given by (19). Only in the region bounded by them there are 4 fixed points; the code is $`H^1E^3.`$ The codes for the black, dark grey and light grey regions are, respectively, $`R^2,E^2`$ and $`R^1E^1.`$ On the top right plot, the value of $`\epsilon `$ is $`3.456789.`$ The region containing the point $`(0.5,0)`$ has exactly 2 fixed points while the regions which contact with this one through arcs have 4. The darker region has 6 and the small region near the upper right corner has 10. The regions around $`\theta =0`$ have 8 fixed points. The solid lines give the location of all bifurcations and changes of stability.

In the bottom left plot, computed for $`\epsilon 9.876543`$, all the lines of bifurcation or change of stability are plotted. The number of fixed points, NFP, in the major regions is shown. The typical transitions are as follow: Consider, for instance, a passage from NFP$`=16`$ to NFP$`=18`$ near $`\beta =0`$ with increasing $`\theta .`$ First a line of creation of an elliptic and hyperbolic points is found. One passes from a code $`R^9H^7`$ to $`R^9H^8E^1.`$ This is followed by a change of stability by passing to $`R^{10}H^8.`$ Later on, inside the region with NFP$`=20`$, the points of $`R`$ type become again of type $`E`$. So the code passes, in the different changes, from $`R^{11}H^9`$ to $`H^9E^{11}.`$

Finally, in the bottom right plot, for $`\epsilon =0.5`$, regions similar to the case $`\epsilon =0.1`$ can be seen, with a different configuration. The level lines of figure xxx are also displayed. It is checked that the highest levels correspond to domains where the map has exactly one elliptic and one hyperbolic fixed points. This fact is also present for smaller values of $`\epsilon .`$

### 7.2 The maps as perturbed twists

The object of interest is the global dynamics of $`gf_\epsilon `$ in $`S^2`$. To this end it is convenient to write these maps in a slightly different, but equivalent, way. For this study $`S^2`$ will be taken as the sphere of radius 1 centered at the origin. Instead of considering $`f_\epsilon `$ as a twist around the $`z`$-axis of angle $`\pi \epsilon (1+z)`$, it will be taken as a twist around the axis of zero longitude and latitude $`\beta `$ with angle of rotation around this axis of a point of coordinates $`(x,y,z)`$ equal to $`\pi \epsilon (1+x\mathrm{cos}(\beta )+z\mathrm{sin}(\beta )).`$ Then, the rotation $`g`$ is simply a rotation of angle $`\theta `$ around the $`z`$-axis, to be denoted by $`R_\theta ^{(z)}.`$ Up to the substitution of $`\pi /2\beta `$ for $`\beta `$, the relative positions of the axes in this formulation and in the previous one are equivalent.

It is instructive to first consider the case $`\theta =2\pi \frac{p}{q},`$ where $`p,q`$ are coprime integers. Let us introduce $`\delta =\pi (1+x\mathrm{cos}(\beta )+z\mathrm{sin}(\beta )).`$ Then

$$f_\epsilon \left(\begin{array}{c}x\\ y\\ z\end{array}\right)=\left(\begin{array}{c}x\\ y\\ z\end{array}\right)+\epsilon \delta \left(\begin{array}{c}y\mathrm{sin}(\beta )\\ x\mathrm{sin}(\beta )z\mathrm{cos}(\beta )\\ y\mathrm{cos}(\beta )\end{array}\right)+O(\epsilon ^2).$$
(23)
The next step is the computation of the map $`M_{q,\theta ,\beta ,\epsilon }:=(R_\theta ^{(z)}f_\epsilon )^q,`$ the parameter $`\beta `$ being the latitude of the axis of $`f_\epsilon .`$ It is clear that at order zero in $`\epsilon `$ one has $`M_{q,\theta ,\beta ,0}=\text{Id}.`$ An elementary computation using formula (23) for $`f_\epsilon `$ and the expression of $`\delta `$ as a function of $`x,z`$ and $`\beta ,`$ gives

$$M_{q,\theta ,\beta ,\epsilon }\left(\begin{array}{c}x\\ y\\ z\end{array}\right)=\left(\begin{array}{c}x\\ y\\ z\end{array}\right)+R_\gamma ^{(z)}\left(\begin{array}{c}x\\ y\\ z\end{array}\right)+O(\epsilon ^2),$$
(24)
where $`R_\gamma ^{(z)}`$ is now a rotation around the $`z`$-axis in each one of the horizontal planes with angle of rotation depending on $`z`$ as follows

$$\gamma =\pi \epsilon q\left(\mathrm{sin}(\beta )+zP_2(\mathrm{sin}(\beta ))\right),$$
(25)
where $`P_2`$ denotes the second Legendre polynomial ($`P_2(w)=\frac{3}{2}w^2\frac{1}{2}`$).

This result tells us that the rotation $`R_\theta ^{(z)}`$ averages the effect of the map $`f_\epsilon `$ in a good way. Let us remark that in (24) it is assumed $`p/q`$ fixed and $`\epsilon `$ sufficiently small. From (25) it follows that the angle $`\gamma `$ is still small provided that $`\epsilon q`$ is small. In the trivial case $`p=0,q=1`$ the rotation is the identity and then the twist can be be taken with $`\beta =\pi /2,`$ recovering in formula (25) the angle rotated in the twist.

To pass to the general case for $`\theta `$ one needs a preliminary lemma.

###### Lemma 7.2

Let $`\rho (0,1)`$ and $`N𝐍.`$ Let

$$S_{\rho ,N}=\underset{1qN,\mathrm{\hspace{0.17em}0}pq,(p,q)=1}{}[\frac{p\rho }{q},\frac{p+\rho }{q}].$$

Then, if $`N+1\rho ^1`$ one has $`[0,1]S_{\rho ,N}.`$

Proof:  Let $`\alpha [0,1].`$ If $`\alpha =r/s𝐐`$ with $`(r,s)=1`$ and $`sN,`$ there is nothing to prove. Otherwise consider the approximants to $`\alpha `$ given by the continued fraction algorithm. Assume that $`\frac{p_1}{q_1}`$ and $`\frac{p_2}{q_2}`$ are consecutive approximants with $`q_1N`$ and $`q_2>N.`$ Then

$$\left|\alpha \frac{p_1}{q_1}\right|\frac{1}{q_1q_2}\frac{1}{q_1(N+1)}\frac{\rho }{q_1}.$$

 QED

Hence $`\frac{\theta }{2\pi }`$ can be written as $`\frac{p}{q}+\frac{\mu }{q}`$, where $`|\mu |\rho ,q+1\rho ^1,`$ where $`\rho `$ is not specified for the moment. Hence one can represent the map $`gf_\epsilon `$ as something similar to the previous case, that is, a rotation whose angle is a rational multiple of $`2\pi `$, composed with a map close to the identity, by writing

$$gf_\epsilon =R_{2\pi p/q}^{(z)}R_{2\pi \mu /q}^{(z)}f_\epsilon .$$

By a direct computation one obtains expression (24) again with the following modifications:

* If the rational which approximates $`\frac{\theta }{2\pi }`$ is 0, then $`q=1`$ and there is no average, so that we keep the map $`gf_\epsilon ,`$
* The rotation is now $`\gamma =\pi \epsilon q\left(\mathrm{sin}(\beta )+zP_2(\mathrm{sin}(\beta ))\right)+2\pi \mu `$,
* The error terms are, uniformly in $`\theta `$, of the form $`O(\epsilon ^2q+\frac{\rho ^2}{q}).`$

If one takes $`\rho =\epsilon ^{2/3}`$ then the maps $`(gf_\epsilon )^q`$ are, in all cases, $`\epsilon ^{1/3}`$–close to the identity and the error terms are at most $`\epsilon ^{4/3}`$. Note that besides the choice $`\rho =\epsilon ^{2/3}`$ there are other possibilities, but $`\rho =\epsilon ^{2/3}`$ is good enough to prove Corollary 7.7. Finally, it is clear that (25), or the modification just mentioned adding $`2\pi \mu `$, is a twist except for $`\beta =\beta _{\text{crit}}`$ such that $`P_2(\mathrm{sin}(\beta _{\text{crit}}))=0`$ ($`\beta _{\text{crit}}=\mathrm{sin}^1(1/\sqrt{3})`$).

To summarize we state the following

###### Proposition 7.3

If $`\theta `$ is not $`\epsilon ^{2/3}`$–close to zero and $`\beta \beta _{\text{crit}}`$ the maps $`gf_\epsilon `$ for $`\epsilon `$ small enough, have a power which is $`\epsilon ^{1/3}`$–close to the identity. This power satisfies a twist condition of order at least $`\epsilon P_2(\mathrm{sin}(\beta ))`$.

###### Remark 7.4

Rotations $`g`$ with small $`\theta `$ are irrelevant, for the present purpose, due to the Haar measure in $`SO(3)`$, and a single exceptional case (non-twist) is also unimportant. This will be seen later in detail.

###### Remark 7.5

In the case of small $`\theta `$ it is still possible to show that $`gf_\epsilon `$ produces a twist effect on each meridian in $`S^2.`$ A problem which appears though is that the angle rotated by the different points can pass through an extremum, losing in this way the twist property. In fact this is not so important because the existence of invariant curves when the twist condition is lost at some point has been established in \[S1\]. But this refinement is not necessary in the present context.

###### Theorem 7.6

With the possible exclusion of an open set $``$ in $`SO(3)`$ of small measure, there exist $`\epsilon _0`$ such that for $`\epsilon <\epsilon _0`$ the maps $`gf_\epsilon `$ have a dynamics exponentially close to an integrable flow in $`S^2.`$

Proof:  The proof is divided into steps.

1. The maps $`M_{q,\theta ,\beta ,\epsilon }`$, being a power of $`gf_\epsilon `$, have the same dynamics as $`gf_\epsilon `$. In all cases (including $`q=1`$ and the exceptional value of $`\beta `$) they are $`\epsilon ^{1/3}`$–close to the identity. Hence there exists a suspension given by the flow of a $`1`$–periodic vector field in $`S^2`$ such that the time–$`1`$ map associated to this flow coincides with $`M_{q,\theta ,\beta ,\epsilon }`$. The vector field is “slow” (of the order of $`\epsilon ^{1/3}`$) and the dominant terms do not depend on time. See \[BRS\] for details and an explicit construction.

It is relevant to note that the vector field is analytic with respect to the phase space variables (that is, the points in $`S^2`$) while the dependence in $`\epsilon `$ is discontinuous in $`SO(3)`$ (moving $`\theta [0,2\pi ]`$ changes the value of $`q`$), but the relevant thing is that it is bounded in $`\epsilon `$. Furthermore the dependence with respect to $`t`$ can be made of class $`𝒞^r`$ for any $`r>0,`$ but continuity in $`t`$ is sufficient for what follows. Furthermore the vector field is Hamiltonian.
2. The next step is to “average” the vector field with respect to $`t`$. This is the content of Neishtadt’s theorem \[N\]. See \[S0\] for a detailed proof. As a consequence the vector field can be written as an autonomous part and a remainder which is exponentially small in the current small parameter; that is, the remainder is bounded by $`\text{exp}(c\epsilon ^{1/3})`$ for some $`c>0`$. Furthermore, the averaged vector field is still Hamiltonian (see \[SV\] for a sketch of the proof).
3. As the averaged system is a Hamiltonian in $`S^2`$, it is integrable and, hence, foliated by invariant curves except on the separatrices, which are a set of zero measure. Most of the invariant curves subsist as a consequence of Moser’s twist theorem. To this end one should have that the perturbation is small compared with the twist condition. Hence, it is enough to exclude a neighborhood of the critical latitude $`\beta _{\text{crit}}`$ which can be taken also exponentially small. Furthermore the set of points in $`S^2`$ not covered by invariant curves of the full system has a measure bounded by the square root of the perturbation, again exponentially small in $`\epsilon .`$

Summarizing, when arbitrary $`gSO(3)`$ are considered the dynamics in $`S^2`$ is ordered (the points lie on invariant curves) except for points in a subset of $`S^2`$ of exponentially small measure. Furthermore $``$ consists of a neighborhood of the identity of size $`O(\epsilon ^{2/3})`$ and a neighborhood of $`\beta _{\text{crit}}`$ which is exponentially small in $`\epsilon `$. QED

###### Corollary 7.7

For $`\epsilon `$ small enough $`\mathrm{\Lambda }(\epsilon )<A\epsilon ^3.`$

Proof:  It is sufficient to make remark 7.4 more explicit. The differentials of the maps $`f_\epsilon `$ increase the length of the vectors in $`TS^2`$ by a factor of the form $`1+O(\epsilon )`$ and composing with $`gSO(3)`$ produces no essential changes in the factor. Hence the values of $`\lambda _1(x,gf_\epsilon )`$ are bounded by $`C\epsilon ,`$ where $`C`$ is a positive constant. The contribution to $`\mathrm{\Lambda }(\epsilon )`$ of the $`g`$ to which Theorem 7.6 applies is bounded by $`C\epsilon `$ times an exponentially small amount. On the other hand the contribution of the excluded neighborhood of $`\beta _{\text{crit}}`$ is also exponentially small.

Therefore the main contribution to $`\mathrm{\Lambda }(\epsilon )`$ can only come from the neighborhood of the identity excluded in Theorem 7.6. But the Haar measure of this set is of the order of

$$_0^{\epsilon ^{2/3}}(1\mathrm{cos}(\theta ))𝑑\theta =O(\epsilon ^2).$$

This bound and the previous one on $`\lambda _1(x,gf_\epsilon )`$ give the result. QED

If we want to consider more “realistic” upper bounds it is possible to proceed along the ideas in remark 7.5. A further consideration is that the largest stochastic zones are typically associated to the splitting of separatrices of the hyperbolic fixed points. From \[FS\] it follows that the splitting can be bounded by

$$\text{exp}\left(\frac{c}{\mathrm{log}(\mu _{\text{max}})}\right),$$

where $`\mu _{\text{max}}`$ is the maximal eigenvalue at the fixed points and $`c>0.`$ From (22) one has that for $`\epsilon `$ small $`\mathrm{log}(\mu _{\text{max}})=\frac{\pi \epsilon }{2}+O(\epsilon ^2).`$ This “heuristic” prediction is in good agreement with the observed behavior for $`\epsilon `$ small.

## 8 The case of large $`\epsilon `$

Let us call the subbundle of $`PS^2`$ tangent to the invariant circles of $`f_\epsilon `$ the horizontal bundle and denote it by $`H`$. As $`\epsilon \mathrm{}`$, a large portion of $`PS^2`$ is sucked into a small neighborhood of $`H`$ under $`f_{\epsilon \mathrm{\#}}`$. The measure $`\delta _{g(H)}`$ on $`PS^2`$ that is atomic in each tangent space and supported on $`g(H)`$ looks more and more like an invariant measure for $`gf_{\epsilon ,\mathrm{\#}}`$. These measures integrate to give Lebesgue:

$$m=_{gSO(3)}\delta _{g(H)}𝑑\nu (g).$$

This yields a heuristic argument for why the inequality in Question 1.1 should hold when $`\epsilon =\mathrm{}`$. In this section, we make this argument rigorous by adding some $`\delta `$-noise in $`_\epsilon `$, and replacing invariant measures with stationary measures. We prove:

###### Theorem 8.1

Let $`m_{\epsilon ,\delta }`$ be defined as in Section 3, and let $`\phi _{\epsilon ,\delta }`$ be the density of $`m_{\epsilon ,\delta }`$:

$$dm_{\epsilon ,\delta }=\phi _{\epsilon ,\delta }dm.$$

There exists $`C>0`$ such that, for all $`\epsilon ,\delta >0`$,

$$\phi _{\epsilon ,\delta }1_1<C\delta ^{11}\epsilon ^{1/2}$$

where $`_1`$ is the $`L^1`$-norm with respect to Lebesgue measure $`m`$ on $`PS^2`$.

This has the corollary:

###### Corollary 8.2

There exists a $`C>0`$ such that for all $`\epsilon ,\delta >0`$,

$$|R(\epsilon ,\delta )R(\epsilon )|C\delta ^{11}\epsilon ^{1/2}\mathrm{log}\epsilon .$$

In particular, for all $`\delta >0`$,

$$\underset{\epsilon \mathrm{}}{lim}|R(\epsilon ,\delta )R(\epsilon )|=0,$$

where $`R(\epsilon ,\delta )`$ is the random diffused exponent defined in Section 3.

As we were finishing this paper, unpublished work of L. Carleson and T. Spencer came to our attention \[CaSp\]. For the standard map on the 2-torus:

$$g_\epsilon :(x,y)(2x+\epsilon \mathrm{sin}(2\pi x)y,x)$$

where $`\epsilon `$ measures the strength of the nonlinearity, they prove that by adding a noise of strength $`\text{exp}\epsilon ^2`$ to the element $`g_\epsilon `$, a Lyapunov exponent of order $`\mathrm{log}\epsilon `$ can be established.

Proof of Corollary 8.2. From the definitions,

$`R(\epsilon ,\delta )`$ $`=`$ $`{\displaystyle _{PS^2}}\mathrm{log}Tf_\epsilon vdm_{\epsilon ,\delta }(v)`$
$`=`$ $`{\displaystyle _{PS^2}}\mathrm{log}Tf_\epsilon v\phi _{\epsilon ,\delta }(v)𝑑m(v),`$
whereas $`R(\epsilon )`$ is the integral of $`\mathrm{log}Tf_\epsilon v`$ with respect to $`m`$. Hence,

$`|R(\epsilon ,\delta )R(\epsilon )|`$ $`=`$ $`|{\displaystyle _{PS^2}}\mathrm{log}Tf_\epsilon v(\phi _{\epsilon ,\delta }(v)1)dm(v)|`$
$``$ $`\mathrm{log}Tf_\epsilon _{\mathrm{}}\phi _{\epsilon ,\delta }1_1`$
$``$ $`C\delta ^{11}\epsilon ^{1/2}\mathrm{log}\epsilon ,`$
by Theorem 8.1. QED

Proof of Theorem 8.1.

$`SO(3)`$ acts transitively on $`T_1S^2`$ by isometries and with trivial stabilizer. From now on we identify points in $`T_1S^2`$ with elements of $`SO(3)`$, and use the group structure in writing our formulas. We will use $`x,y,z`$ to denote elements of $`PS^2`$, $`p,q`$ for elements of $`S^2`$, and $`(p,v),(q,w)`$ for elements of $`PS^2`$ (or $`T_1S^2`$).

Recall from Section 3 that

$$\phi _{\epsilon ,\delta }=_{PS^2}\phi _{\epsilon ,\delta ,g}𝑑\nu (g),$$

where $`\phi _{\epsilon ,\delta ,g}`$ is the fixed point of the operator $`L_{\epsilon ,\delta ,g}`$ defined by:

$`L_{\epsilon ,\delta ,g}\psi (x)`$ $`=`$ $`{\displaystyle \frac{1}{m(B(x,\delta ))}}{\displaystyle _{(gf_{\epsilon }^{}{}_{\mathrm{}}{}^{})^1B(x,\delta )}}\psi (y)𝑑m(y)`$
Setting

$$k_\delta (x)=\frac{1}{m(B(e,\delta ))}1_{B(e,\delta )}(x),$$

where $`e`$ is the identity element of $`SO(3)`$, we rewrite $`L_{\epsilon ,\delta ,g}`$ as:

$`L_{\epsilon ,\delta ,g}\psi (x)`$ $`=`$ $`{\displaystyle _{PS^2}}k_\delta (x^1gf_{\epsilon }^{}{}_{\mathrm{}}{}^{}(y))\psi (y)𝑑m(y).`$ (26)
Let $`L_{\epsilon ,\delta }=L_{\epsilon ,\delta ,e}`$. It is clear from (26) that $`L_{\epsilon ,\delta ,g}\psi (x)=L_{\epsilon ,\delta }\psi (g^1x)`$

Denote by $`K_{\epsilon ,\delta }:PS^2\times PS^2𝐑_+`$ the kernel of the operator $`L_{\epsilon ,\delta }`$, so that

$$K_{\epsilon ,\delta }(x,y)=k_\delta (x^1f_{\epsilon }^{}{}_{\mathrm{}}{}^{}(y)).$$

Let $`\pi :PS^2S^2`$ be the projection along tangent fibers. By averaging along fibers, we shall approximate $`K_{\epsilon ,\delta }`$ by a new kernel $`\widehat{K}_{\epsilon ,\delta }`$ that is constant along fibers of the second $`PS^2`$ \- factor. Define $`\widehat{K}_{\epsilon ,\delta }:PS^2\times PS^2𝐑_+`$ by

$$\widehat{K}_{\epsilon ,\delta }(x,y)=_{\pi ^1\pi y}K_{\epsilon ,\delta }(x,z)𝑑m_{\pi y}(z),$$

where, for $`pS^2`$, $`m_p`$ denotes the disintegration of $`m`$ along the fiber $`\pi ^1p`$. For $`gSO(3)`$, we obtain a new operator $`\widehat{L}_{\epsilon ,\delta ,g}`$ on $`L^{\mathrm{}}(PS^2)`$, given by:

$$\widehat{L}_{\epsilon ,\delta ,g}\varphi (x)=_{PS^2}\widehat{K}_{\epsilon ,\delta }(g^1x,y)\varphi (y)𝑑m(y).$$

Let $`\widehat{L}_{\epsilon ,\delta }=\widehat{L}_{\epsilon ,\delta ,e}`$.

The next lemma shows that the operators $`\widehat{L}_{\epsilon ,\delta ,g}`$ have a good averaging property when applied to densities of measures that project to Lebesgue measure $`\mu `$ on $`S^2`$.

###### Lemma 8.3

Let $`\widehat{K}:PS^2\times PS^2𝐑_+`$ be any $`L^1`$ function such that:

1. for all $`pS^2`$,

$$_{\pi ^1p}_{PS^2}\widehat{K}(x,y)𝑑m_p(x)𝑑m(y)=1;$$
2. if $`\pi (y)=\pi (z)`$, then for all $`x`$, $`\widehat{K}(x,y)=\widehat{K}(x,z)`$.

For $`\varphi L^{\mathrm{}}(PS^2)`$, and $`gSO(3)`$, define $`\widehat{L}_g\varphi L^{\mathrm{}}(PS^2)`$ by

$$\widehat{L}_g\varphi (x)=_{PS^2}\widehat{K}(g^1x,y)\varphi (y)𝑑m(y).$$

Then:

1. $`\widehat{\phi }_g=\widehat{L}_g\phi `$ is the unique fixed point of $`\widehat{L}_g`$, for any $`\phi L^{\mathrm{}}(PS^2)`$ that satisfies:

$$_{\pi ^1p}\phi (z)𝑑m_p(z)=1,$$

for all $`pS^2`$.
2. for all $`xPS^2`$, we have:

$$_{SO(3)}\widehat{\phi }_g(x)𝑑\nu (g)=1.$$

###### Remark 8.4

Lemma 8.3 applies to the operators $`\widehat{L}_{\epsilon ,\delta ,g}`$. An example of a function $`\phi `$ that satisfies the hypotheses of Lemma 8.3 is the density $`\phi _{\epsilon ,\delta ,g}`$, for any $`\epsilon ,\delta ,g`$. We do not use that $`\widehat{\phi }_g=\widehat{L}_g\widehat{\phi }_g`$ below.

Proof of Lemma 8.3: Since it is constant along fibers of the second factor, $`\widehat{K}`$ projects to a function on $`PS^2\times S^2`$, which we shall also call $`\widehat{K}`$.

For $`gSO(3)`$, define $`\widehat{\phi }_g`$ by:

$$\widehat{\phi }_g(x)=\widehat{L}_g1(x)=_{pS^2}\widehat{K}(g^1x,p)𝑑\mu (p).$$

We compute directly that, for any $`\phi `$ satisfying the hypotheses of a.,

$`\widehat{L}_g(\phi )(x)`$ $`=`$ $`{\displaystyle _{yPS^2}}\widehat{K}(g^1x,y)\phi (y)𝑑m(y)`$
$`=`$ $`{\displaystyle _{pS^2}}{\displaystyle _{z\pi ^1p}}\widehat{K}(g^1x,p)\phi (z)𝑑m_p(z)𝑑\mu (p)`$
$`=`$ $`{\displaystyle _{pS^2}}K(g^1x,p)𝑑\mu (p)`$
$`=`$ $`\widehat{\phi }_g(x).`$
To see that $`\widehat{L}_g\widehat{\phi }_g=\widehat{\phi }_g`$ and finish the proof of a it is now sufficient to verify that

$$_{\pi ^1p}\widehat{\phi }_g(z)𝑑m_p(z)=1.$$

But $`\widehat{\phi }_g(z)=\widehat{L}_g1(z)=L_g1(z)`$ since the function $`1`$ is constant and hence constant on fibers. Now $`L_g(1)`$ is the density function of a measure on $`PS^2`$ which covers Lebesgue measure on $`S^2`$ and hence its integral on fibers equals $`1`$. This proves a.

Integrating $`\widehat{\phi }_g(x)`$ with respect to $`g`$ we obtain:

$`{\displaystyle _{gSO(3)}}\widehat{\phi }_g(x)𝑑\nu (g)`$ $`=`$ $`{\displaystyle _{gSO(3)}}{\displaystyle _{pS^2}}\widehat{K}(g^1x,p)𝑑\mu (p)𝑑\nu (g)`$
$`=`$ $`{\displaystyle _{yPS^2}}{\displaystyle _{pS^2}}\widehat{K}(y,p)𝑑\mu (p)𝑑m(y)`$
$`=`$ $`1,`$
completing the proof of b. QED

Returning to the proof of Theorem 8.1, let $`\widehat{\phi }_{\epsilon ,\delta ,g}=\widehat{L}_{\epsilon ,\delta ,g}1=\widehat{L}_{\epsilon ,\delta ,g}\phi _{\epsilon ,\delta ,g}`$ be the unique fixed point of $`\widehat{L}_{\epsilon ,\delta ,g}`$ given by Lemma 8.3. We now have:

$`\phi _{\epsilon ,\delta }1_1`$ $`=`$ $`{\displaystyle _{gSO(3)}}(\phi _{\epsilon ,\delta ,g}1)𝑑\nu (g)_1`$
$`=`$ $`{\displaystyle _{gSO(3)}}(\phi _{\epsilon ,\delta ,g}\widehat{\phi }_{\epsilon ,\delta ,g})𝑑\nu (g)_1+{\displaystyle _{gSO(3)}}(\widehat{\phi }_{\epsilon ,\delta ,g}1)𝑑\nu (g)_1`$
$``$ $`{\displaystyle _{gSO(3)}}(L_{\epsilon ,\delta ,g}\phi _{\epsilon ,\delta ,g}\widehat{L}_{\epsilon ,\delta ,g}\phi _{\epsilon ,\delta ,g})𝑑\nu (g)_1`$
$`+{\displaystyle _{gSO(3)}}(\widehat{\phi }_{\epsilon ,\delta ,g}1)𝑑\nu (g)_1`$
$`=`$ $`{\displaystyle _{gSO(3)}}(L_{\epsilon ,\delta }\phi _{\epsilon ,\delta ,g}g^1\widehat{L}_{\epsilon ,\delta }\phi _{\epsilon ,\delta ,g}g^1)𝑑\nu (g)_1`$
$``$ $`{\displaystyle _{gSO(3)}}L_{\epsilon ,\delta }\phi _{\epsilon ,\delta ,g}g^1\widehat{L}_{\epsilon ,\delta }\phi _{\epsilon ,\delta ,g}g^1_1𝑑\nu (g)`$
where we used Lemma 8.3 to obtain the second to last inequality.

Propositions 8.5 and 8.6, which we state and prove below, now imply that

$`{\displaystyle _{gSO(3)}}L_{\epsilon ,\delta }\phi _{\epsilon ,\delta ,g}g^1\widehat{L}_{\epsilon ,\delta }\phi _{\epsilon ,\delta ,g}g^1_1𝑑\nu (g)`$ $``$ $`C_1\delta ^3\epsilon ^{1/2}{\displaystyle _{gSO(3)}}\phi _{\epsilon ,\delta ,g}_{\mathrm{}}𝑑\nu (g)`$
$``$ $`C_2\delta ^3\epsilon ^{1/2}{\displaystyle _{gSO(3)}}\delta ^8𝑑\nu (g)`$
$`=`$ $`C\delta ^{11}\epsilon ^{1/2},`$
completing the proof of the theorem. It remains to state and prove Propositions 8.5 and 8.6.

###### Proposition 8.5

There is a $`C>0`$ such that, if $`\epsilon \delta ^4`$, then for any $`\varphi L^{\mathrm{}}(PS^2)`$,

$$L_{\epsilon ,\delta }(\varphi )\widehat{L}_{\epsilon ,\delta }(\varphi )_1C\varphi _{\mathrm{}}\delta ^3\epsilon ^{1/2}.$$

###### Proposition 8.6

There exists $`C>0`$ such that, for all $`\epsilon >0`$ and $`\delta >0`$,

$$\phi _{ϵ,\delta ,g}_{\mathrm{}}C\delta ^8.$$

Proof of Proposition 8.5:

$`L_{\epsilon ,\delta }(\varphi )\widehat{L}_{\epsilon ,\delta }(\varphi )_1`$ $``$ $`\varphi _{\mathrm{}}K_{\epsilon ,\delta }\widehat{K}_{\epsilon ,\delta }_1`$
$`=`$ $`\varphi _{\mathrm{}}{\displaystyle |K_{\epsilon ,\delta }(x,y)\widehat{K}_{\epsilon ,\delta }(x,y)|𝑑m(x)𝑑m(y)}`$
$`=`$ $`\varphi _{\mathrm{}}{\displaystyle |K_{\epsilon ,\delta }(x,y)_{\pi ^1\pi y}K_{\epsilon ,\delta }(x,z)𝑑m_{\pi y}(z)|𝑑m(x)𝑑m(y)}.`$
Let $`x=(p,v),y=(q,w)`$ be elements of $`PS^2`$. For a fixed $`p,v,q`$, the map $`wK_{\epsilon ,\delta }((p,v),(q,w))`$ is a constant multiple $`c\delta ^3`$ of the characteristic function for $`\pi ^1qf_{\epsilon }^{}{}_{\mathrm{\#}}{}^{1}B((p,v),\delta )`$. Note that, for any measurable set $`B`$ in a probability space with measure $`m`$, the average value of the function $`|m(B)1_B|`$ is $`2m(B)(1m(B))`$. It follows that:

$`K_{\epsilon ,\delta }\widehat{K}_{\epsilon ,\delta }_1`$ $`=`$ $`2c\delta ^3{\displaystyle _{(p,v)PS^2}}{\displaystyle _{qf_\epsilon ^1B(p,\delta )}}|\beta (p,v,q)(1\beta (p,v,q))|𝑑m((p,v))𝑑\mu (q),`$
where

$$\beta (p,v,q)=m_q(f_{\epsilon }^{}{}_{\mathrm{\#}}{}^{1}B((p,v),\delta )).$$

We next show that for some $`k>0`$, and for $`\epsilon >\delta ^4`$, there is a set $`GPS^2`$, with $`m(G)1\epsilon ^{1/2}`$ such that, for all $`(p,v)G`$, there is a set $`G^{}=G^{}(p,v)f_\epsilon ^1B(p,\delta )`$ with $`\mu (G^{})\mu (B(p,\delta ))\epsilon ^{1/2}`$, such that, for $`qG^{}`$,

$`\beta (p,v,q)`$ $``$ $`k\epsilon ^{1/2}\text{or}\beta (p,v,q)1k\epsilon ^{1/2}.`$ (27)
This implies that

$`K_{\epsilon ,\delta }\widehat{K}_{\epsilon ,\delta }_1`$ $``$ $`2c\delta ^3\left(k\epsilon ^{1/2}\mu (B(p,\delta ))+2\epsilon ^{1/2}\right)`$
$``$ $`C\delta ^3\epsilon ^{1/2},`$
which implies the result.

Fix $`\delta <1/2`$, and assume that $`\epsilon >\delta ^4`$. For $`\alpha >0`$, denote by $`C_\alpha `$ the $`\alpha `$-neighborhood of the horizontal bundle $`HPS^2`$. In other words, $`C_\alpha `$ is the set of lines in $`PS^2\pi ^1\{NP,SP\}`$ at angle $`\alpha `$ with the latitudinal circles.

It is not difficult to see that if the distance from $`p`$ to the poles is greater than $`\epsilon ^{1/4}`$, then

$`m_{f_\epsilon ^1p}(f_{\epsilon \mathrm{\#}}^1(C_{\epsilon ^{1/2}}\pi ^1p))`$ $``$ $`1\epsilon ^{1/2}.`$ (28)
Let $`G_0`$ be the set of $`p`$ for which (28) holds, and let $`G=\pi ^1G_0`$. Then $`m(G)1\epsilon ^{1/2}`$.

Fix $`(p,v)G`$, and consider a point $`f_\epsilon qB(p,\delta )G_0`$. The intersection of $`B((p,v),\delta )`$ with the fiber $`\pi ^1f_\epsilon q`$ is an interval $`I`$. If the endpoints of $`I`$ are disjoint from the interval $`J=C_{\epsilon ^{1/2}}\pi ^1f_\epsilon q`$, then either $`IJ`$ or $`IJ=\mathrm{}`$. In the former case, the length of $`f_{\epsilon \mathrm{\#}}^1I`$ is greater than the length of $`f_{\epsilon \mathrm{\#}}^1J`$, which by (28) is greater than $`1\epsilon ^{1/2}`$. In the latter case, the length of $`f_{\epsilon \mathrm{\#}}^1I`$ is less than $`\epsilon ^{1/2}`$. Hence if we let $`G^{}=G^{}(p,v)`$ be the set of $`qf_\epsilon ^1B(p,\delta )G_0`$ satisfying:

$`B((p,v),\delta )C_{\epsilon ^{1/2}}\pi ^1f_\epsilon q`$ $`=`$ $`\mathrm{},`$ (29)
then (27) holds for all $`qG^{}`$. It remains to show that $`\mu (G^{})\mu (B(p,\delta ))\epsilon ^{1/2}`$.

For $`(p,v)PS^2`$, denote by $`S((p,v),\delta )`$ the geodesic sphere of radius $`\delta `$ centered at $`(p,v)`$, so $`S((p,v),\delta )=B((p,v),\delta )`$. We will use the following lemma here and later in the proof of Proposition 8.6.

###### Lemma 8.7

There exists $`C>0`$ such that, for all $`(p,v)\pi ^1p`$, and all $`\alpha >0`$,

$$\mu \left(\pi (S((p,v),\delta )C_\alpha )\right)C\delta \alpha .$$

Proof. The claim follows from the following facts:

1. On $`T_1S^2\pi ^1\{NP,SP\}`$, the subbundle $`H`$ (regarded as a submanifold) is uniformly transverse to the fibers of $`PS^2`$.
2. There exists a $`C>0`$ such that for all $`\delta `$ sufficiently small, and for all $`(p,v)PS^2`$, the intersection $`S((p,v),\delta )H`$ is contained in a smooth curve of length $`C\delta `$.

The verification of these facts is left as an exercise. QED

From Lemma 8.7 it follows that:

$`\mu (G^{})`$ $``$ $`\mu (B(p,\delta ))\mu (G_0)\mu (f_\epsilon ^1\pi (S((p,v),\delta )C_{\epsilon ^{1/2}})`$
$``$ $`\mu (B(p,\delta ))\epsilon ^{1/2}`$
This completes the proof of Proposition 8.5. QED

Proof of Proposition 8.6: We know that $`\phi _{ϵ,\delta ,g}`$ is a function in $`L^1`$ which satisfies, for all $`(p,v)PS^2`$,

$`\phi _{ϵ,\delta ,g}(p,v)`$ $`=`$ $`c\delta ^3{\displaystyle _{f_{\epsilon \mathrm{\#}}^1B(g^1(p,v),\delta )}}𝑑\phi _{ϵ,\delta ,g}(q,w)𝑑\mu (q)𝑑m_q(w)`$ (30)
and, for all $`pS^2`$,

$`{\displaystyle _{\pi ^1p}}\phi _{ϵ,\delta ,g}(p,v)𝑑m_p(v)`$ $`=`$ $`1.`$ (31)
Then, by (30), the function $`\phi _{ϵ,\delta ,g}`$ is continuous, and therefore has a maximum $`M_{ϵ,\delta ,g}`$ that we denote by $`M`$. The idea is that a Hölder constant for $`\phi _{ϵ,\delta ,g}`$ can be estimated in terms of $`M`$. Reporting in (31) gives a bound for $`M`$ which is independent of $`\epsilon ,g`$. Since we want to use (31) at the end, it suffices to consider the Hölder constant along the fibers. So, let $`(p,v),(p,v^{})\pi ^1p`$. We have:

$$|\phi _{ϵ,\delta ,g}(p,v)\phi _{ϵ,\delta ,g}(p,v^{})|cM\delta ^3m(f_{\epsilon \mathrm{\#}}^1B(g^1(p,v),\delta )\mathrm{\Delta }f_{\epsilon \mathrm{\#}}^1B(g^1(p,v^{}),\delta )),$$

where $`m`$ is Lebesgue measure on $`T_1S^2`$, and $`A\mathrm{\Delta }B`$ stands for the set of points which belong to only one of the subsets $`A`$ or $`B`$.

###### Lemma 8.8

There exists $`C>0`$ such that, for all $`v,v^{}T_{1,p}S^2`$,

$$m(f_{\epsilon \mathrm{\#}}^1B(g^1(p,v),\delta )\mathrm{\Delta }f_{\epsilon \mathrm{\#}}^1B(g^1(p,v^{}),\delta ))C\delta d(v,v^{})^{1/4}.$$

###### Remark 8.9

In Proposition 3.3., we prove that $`\phi _{ϵ,\delta ,g}`$ is as smooth as $`f_\epsilon `$. The point of the arguments here is to get a Hölder constant independent of $`\epsilon `$.

From this lemma, it then follows that $`M5(Cc)^4\delta ^8/2`$, since:

$`1`$ $`=`$ $`{\displaystyle _{T_pS^2}}\phi _{ϵ,\delta ,g}(p,v)𝑑v`$
$``$ $`M{\displaystyle _{\mathrm{}}^{\mathrm{}}}(1Cc\delta ^2|t|^{1/4})^+𝑑t`$
$`=`$ $`{\displaystyle \frac{2M}{5}}(Cc)^4\delta ^8.`$
We now prove Lemma 8.8.

We have the two balls $`B((p,v),\delta ))`$ and $`B((p,v^{}),\delta )`$. Let $`\alpha =\sqrt{d(v,v^{})}`$. We may assume that $`\alpha <<\delta `$. The set $`B((p,v),\delta ))\mathrm{\Delta }B((p,v^{}),\delta )`$ meets the fiber $`\pi ^1q`$ in a pair of intervals, each of length $`\alpha ^2<<\sqrt{\alpha }`$. If the endpoints of these intervals do not lie in $`C_\sqrt{\alpha }`$, then the entire intervals must be disjoint from $`C_\sqrt{\alpha }`$. In other words, if

$`\pi ^1q(S((p,v),\delta )S((p,v^{}),\delta ))C_\sqrt{\alpha }=\mathrm{},`$ (32)
then

$`\pi ^1qB((p,v),\delta ))\mathrm{\Delta }B((p,v^{}),\delta )C_\sqrt{\alpha }=\mathrm{}.`$ (33)
Let $`GS^2`$ be the set of $`q`$ satisfying (32). By Lemma 8.7, $`\mu (G)12C\delta \sqrt{\alpha }`$.

###### Claim 8.10

There exists a $`C>0`$ such that, for all $`\alpha 1`$, $`pS^2`$, and $`\epsilon 0`$, if $`(p,v)C_\sqrt{\alpha }`$, then

$$T_{(p,v)}f_{\epsilon ,\mathrm{\#}}^1|_{T_{1,p}S^2}C\alpha ^1.$$

Proof. Recall that $`T_{(p,v)}f_{\epsilon ,\mathrm{\#}}^1|_{\pi ^1p}=T_pf_\epsilon ^1v^2`$. With respect to the orthonormal basis of $`T_pS^2`$ of the form $`\{e_1(p),e_2(p)\}`$, where $`e_1(p)H`$ points in the direction of $`f_\epsilon `$-twist and $`e_2(p)`$ points toward the north pole $`NP`$, $`T_pf_\epsilon `$ takes the form:

$$T_pf_\epsilon ^1=\left(\begin{array}{cc}1& \beta \\ 0& 1\end{array}\right),$$

for some $`\beta 0`$. A direct computation shows that there exists a constant $`C>0`$ such that, for all $`\alpha ,\beta 0`$, if the angle between a unit vector $`v𝐑^2`$ and the $`x`$-axis is greater than $`\sqrt{\alpha }`$, then:

$$\left(\begin{array}{cc}1& \beta \\ 0& 1\end{array}\right)v^2C\alpha ^1$$

From this the claim follows. QED

Claim 8.10 and (33) imply that for $`qG`$, the derivative of $`f_{\epsilon \mathrm{\#}}`$ on
$`\pi ^1qB((p,v),\delta )\mathrm{\Delta }B((p,v^{}),\delta )`$ is bounded:

$$T_{(q,w)}f_{\epsilon \mathrm{\#}}|_{T_{1,q}S^2}\alpha ^1,$$

for all $`w`$ such that $`(q,w)B((p,v),\delta ))\mathrm{\Delta }B((p,v^{}),\delta )`$. But for $`qG`$,

$`m_{f_\epsilon ^1q}(f_{\epsilon \mathrm{\#}}^1(B((p,v),\delta ))\mathrm{\Delta }B((p,v^{}),\delta ))`$ $``$ $`\alpha ^1m_q(B((p,v),\delta ))\mathrm{\Delta }B((p,v^{}),\delta ))`$
$``$ $`\alpha ^1d(v,v^{})`$
$`=`$ $`d(v,v^{})^{1/2}.`$
But then

$`m(B((p,v),\delta ))\mathrm{\Delta }B((p,v^{}),\delta )`$ $``$ $`2C_1\delta \alpha ^{1/2}+C_2\delta ^2d(v,v^{})^{1/2}`$
$``$ $`C\delta d(v,v^{})^{1/4},`$
completing the proof of Proposition 8.6 and of Theorem 8.1. QED

## 9 Discussion

We have wondered \[BuPuShWi\] about the relationship of the random Lyapunov exponent of a measure on the space of volume preserving diffeomorphisms of a manifold to the mean of the Lyapunov exponents of the individual members. The point of the question we raised was to be able to conclude that in a rich enough family of diffeomorphisms there must be some with positive Lyapunov exponents, that is to say positive entropy. At question is what sort of notion of richness would make such a conclusion valid. We even proposed that much more might conceivably be true, a lower bound for the mean of the Lyapunov exponents in terms of the random exponents for orthogonally invariant measures on volume preserving diffeomorphisms of the sphere. The orthogonal invariance of the measure was to provide the necessary “richness”.

In the studied family strong numerical evidence has been found about the existence of such a lower bound when the values of the stretching parameter $`\epsilon `$ are not too small. In some sense strong stretching has an effect similar to randomization, but it depends in a clear way on the concrete map. More concretely

* Even moderate values of $`\epsilon `$ like $`\epsilon 10`$ are enough to have an average of the metric entropy larger than the one corresponding to the random map.
* There exist unbounded parameters $`\epsilon `$ for which islands are born. The range of existence of these islands is small, but only the islands associated to fixed points have been considered.
* For small $`\epsilon `$ the estimated average entropy seems positive and definitely to be much less than the one of the random map. The numerical evidence is in favor of the existence of exponentially small lower and upper bounds (in the present example, with an analytic family).

The problems in numerically estimating exponents and how to overcome them have been discussed. A partial analysis of the family of maps has been done for $`\epsilon `$ small. Even a rough estimate of an upper bound of the averaged entropy is enough to show that the this averaged entropy falls below any constant multiple of the entropy of the randomized system, if $`\epsilon `$ is sufficiently small.

Finally, the effect of a small randomization of fixed size $`\delta `$ of the individual elements of the family $`_\epsilon `$ is considered. Now the mean of the local random exponents of the family is indeed asymptotic to the random exponent of the entire family as $`\epsilon `$ tends to infinity; that is, $`R(\epsilon ,\delta )`$ and $`R(\epsilon )`$ are asymptotic.