Oscillatory magnetic response in first-order smoothing ($a=\eta k^2,d=a^2+{\Omega }^2$)

For a velocity Fourier mode with wavenumber $k$, define $a=\eta k^2$, $d=a^2+{\Omega }^2$ and $\mathbf{C}=i({\mathbf{B}}_0\cdot \mathbf{k})\hat{\mathbf{u}}$. A cosine-forced mode has periodic response coefficients $(\mathbf{p},\mathbf{q})=(a\mathbf{C},\Omega \mathbf{C})/d$, while a sine-forced mode has $(-\Omega \mathbf{C},a\mathbf{C})/d$. Homogeneous transients decay as $e^{-at}$. The periodic response is a particular solution, not arbitrary initial data.