Small-vorticity phototactic orientation response (source code)

= Small-vorticity phototactic orientation response
{title2=$\langle\mathbf p\rangle=K_0\mathbf e_x+\varepsilon K_1\mathbf e_z+O(\varepsilon^2)$}

For light along $\mathbf e_x$, <vorticity> along $\mathbf e_y$ and $\varepsilon=\omega/(2D_r)$, spherical coordinates $\mathbf p=(\cos\theta,\sin\theta\cos\phi,\sin\theta\sin\phi)$ give a first orientation-density correction proportional to $\sin\phi\,g(\theta)$. The linearized alignment-diffusion operator preserves azimuthal harmonics and the <rotation> of the axisymmetric equilibrium forces only this sine harmonic. Angular integration makes the first correction to the mean purely vertical. Uniqueness follows from the weighted identity $\mathcal A(f_0h)=\nabla_p\cdot(f_0\nabla_ph)$: <integration by parts> shows the only homogeneous smooth solution is a multiple of $f_0$, removed by <normalization>.