 Solution

Write the perturbation velocity as $(u,v)$, the density amplitude as ${\sigma }_1$, and ${\gamma }_d=ϵ\Omega$. Axisymmetry removes advection by the background azimuthal flow, but the radial perturbation advects the Keplerian shear: $u{\partial }_x{U}_{g,y}=-3\Omega u/2$. Combining this with the Coriolis force gives the linearized shearing sheet equationsThe razor-thin disk Poisson kernel supplies the self-gravity term. The $\Omega /2$ coefficient, rather than $2\Omega$, is essential: it includes the perturbed advection of the background velocity.

Put $A=c^2{k}^2-2\pi G{\sigma }_0|k|$ and ${\omega }^2={\Omega }^2+A$. The determinant of the three amplitude equations isExpanding it gives the dust gravitational dispersion relation with gas dragUsing a determinant avoids division by $s$ and retains the neutral/secular branch. The radial epicyclic frequency of this Keplerian shearing sheet is $\Omega$, so the three terms in ${\omega }^2$ represent rotational support, self-gravity and dust pressure.