 Solution

Let $f$ be probability density per unit solid angle, with ${\int }_{S^2}fd\Omega =1$ and isotropic initial condition $f(\mathbf{p},0)=1/(4\pi )$. In still fluid, bottom-heavy spherical-cell orientation dynamics gives $\dot{\mathbf{p}}=B^{-1}[\mathbf{k}-(\mathbf{k}\cdot \mathbf{p})\mathbf{p}]$. Adding isotropic rotational diffusion gives the gyrotactic orientation Fokker-Planck equationThe derivatives are tangential to the unit sphere. Rotational symmetry about $\mathbf{k}$ is preserved by both the equation and the initial condition, so $f$ depends only on $\theta ,t$. Since $\dot{\theta }=-\sin \theta /B$, the polar probability flux is $j_{\theta }=-\sin \theta f/B-D_r{f}_{\theta }$. Expanding ${\partial }_{t}f+{\sin }^{-1}\theta {\partial }_{\theta }(\sin \theta j_{\theta })=0$ yieldsRegular solutions have no probability flux through the poles. The mean swimming velocity is ${\mathbf{V}}_c=V_s⟨\mathbf{p}⟩$, with only a vertical component by symmetry.