Solution (source code)

= Solution

At stationary onset, $\mu=0$. The lowest stress-free vertical mode is $\widehat W=\sin\pi z$, for which $D^2=-\pi^2$ and $\widehat D^2=-(\pi^2+k^2)$. Substitution gives
$$
\boxed{
Ra(k)=\frac{(\pi^2+k^2)^3+\pi^2(\lambda/\sigma)^2}{k^2}}.
$$
The rotation term is positive, so rotation raises the critical Rayleigh number and is stabilizing.