I've been trying to understand the asymptotic behavior of Ricci flow, and there are two facts which I am unable to square away. I'm interested in higher dimensional manifolds, but my question is easier to state for Riemann surfaces. I suspect that the solution for surfaces will also solve the case in higher dimensions, so I'll focus on the two-dimensional case.
Short version
Given a Riemann surface which is topologically a sphere, the normalized Ricci flow is known to converge exponentially to a round metric. However, when we analyze the evolution of the scalar curvature for a small deformation of the round metric, it appears that there is a non-trivial center which seems to prevent the metric from converging exponentially. As such, I'm trying to determine the correct asymptotics for the scalar curvature, and why my calculation is going wrong.
Long version
Fact 1: In [1], Richard Hamilton proved that given a Riemann surface of positive curvature, if one evolves the metric by normalized Ricci flow \begin{equation} \frac{\partial}{\partial t} g_{i j}=(r-R) g_{i j}, \end{equation} then the metric will converge exponentially quickly to a metric of constant positive curvature (i.e., a round sphere). Here, $R$ is the scalar curvature and $r$ is the average scalar curvature, which acts to normalize the flow.
"Fact 2": The round unit sphere is a fixed point for normalized Ricci flow, so we study asymptotics of the flow for small (volume-preserving) deformations of the round metric. Under Ricci flow, the scalar curvature $R$ evolves via the reaction-diffusion equation \begin{equation} \frac{\partial R}{\partial t}=\Delta R+R^{2}-r R, \end{equation} where $\Delta$ is the Laplace-Beltrami operator with respect to $g(t)$. Since we are considering the unit sphere, the Gauss-Bonnet theorem implies that $r = 2$, so we consider the function $\phi = R-2$ and have that \begin{equation} \frac{\partial \phi}{\partial t}=\Delta \phi+R \phi. \end{equation}
If we linearize this equation at the round metric, we find that \begin{equation} \frac{\partial \phi}{\partial t}=\Delta_{\mathbb{S}^2} \phi+ 2 \phi, \end{equation}
To understand the asymptotic behavior of $\phi$, we can decompose it as $\phi = \sum_{\ell, m} a^{\ell}_m Y^m_{\ell}$, where $Y^m_{\ell}$ are spherical harmonics. Doing so, the flow under the preceding equation is simply \begin{equation} \phi(z,t) = \sum_{\ell, m} a^{\ell}_m e^{(-\ell(\ell+1)+2)t} Y^m_{\ell} \end{equation}
Here, this is using the fact that the spectrum of $\Delta_{\mathbb{S}^2}$ is $\ell (\ell+1)$. However, this seems to introduce a problem; the reaction term $2 \phi$ exactly cancels the diffusion effect on the $Y^m_{1}$-components of $\phi$. In other words, it seems like if $\phi$ is a principle eigenfunction of the Laplace-Beltrami operator (on the round sphere), then it is unchanged by Ricci flow, which seems to preclude exponential convergence of the flow.
Reintroducing the non-linearity
One might suspect that the issue is that I've linearized the flow at the round metric. However, the flow converges exponentially quickly to a round metric and a metric's principle eigenvalue is Lipschitz in the metric (which can be seen using the Rayleigh quotient property). As such, for a small deformation of the round metric where $\phi$ is close to a principal eigenfunction, we have that $ \left| \frac{\partial \phi}{\partial t}-\Delta_{\mathbb{S}^2} \phi- 2 \phi \right | < C e^{-\delta t} \phi $ for some $C$ and $\delta>0$. In other words, it seems that we can construct $\phi$ so that $$ \left| \frac{\partial \phi}{\partial t} \right | < C e^{-\delta t} \phi,$$ which again seems to preclude exponential convergence of the metric.
What's wrong with Fact 2?
Clearly, something is wrong with the analysis in the second argument. There is something that I have missed or computed incorrectly. However, I can't find the flaw in the argument, and I'd really like to know where I am going astray here.
[1] Hamilton, Richard S., The Ricci flow on surfaces, Mathematics and general relativity, Proc. AMS-IMS-SIAM Jt. Summer Res. Conf., Santa Cruz/Calif. 1986, Contemp. Math. 71, 237-262 (1988). ZBL0663.53031.