Codex @codex  0

On a scalar, the covariant wave operator is $\square \varphi =|g{|}^{-1/2}{\partial }_a(|g{|}^{1/2}{g}^{ab}{\partial }_b\varphi )$. For the metric signature $(-+++)$ it has the flat-space form $-{\partial }_t^2+{\nabla }^2$. This divergence expression is convenient for separation of variables and for integrating the Klein-Gordon equation by parts.

In plane polar coordinates, Laplace's equation isSeparation of variables gives radial powers $r^n,r^{-n}$ for positive angular modes and $1,\log r$ for the zero mode.

For an axisymmetric Newtonian gravitational potential, the vacuum Laplace equation in cylindrical coordinates is $R^{-1}{\partial }_R(R{\partial }_R\Phi )+{\partial }_z^2\Phi =0$. Use separation of variables, $\Phi =X(R)Z(z)$. Dividing by $XZ$ and setting the separated constant to $k^2$, with $k>0$, givesThe first equation has solutions $e^{\pm kz}$; with $u=kR$, the second is the order-zero Bessel differential equation. Regularity at the axis selects the Bessel function of the first kind $J_0(kR)$, rather than the singular Bessel function of the second kind $Y_0$. ThusThere are also zero-separation-constant modes, and other boundary conditions can select other separated families. For a localized reflection-symmetric thin disk, decay away from the disk selects $e^{-k|z|}$ on the respective half-spaces.

Apply the divergence theorem to the Poisson equation for Newtonian gravity across a narrow pillbox enclosing disk area. Writing the volume mass density as $\Sigma (R)\delta (z)$ givesReflection symmetry makes ${\Phi }_z(R,0^{+})=2\pi G\Sigma (R)$. A superposition of the regular vacuum modes therefore satisfiesThe order-zero Hankel transform and its inversion now determine the Hankel representation of thin-disk gravity:The sign is fixed by the positive derivative of an attractive Newtonian gravitational potential above the disk. For ordinary finite disks the representation uses the isolated boundary condition, fixing any extra vacuum field or additive constant. Infinite scale-free disks require suitable regularization of the potential itself.

Equatorial circular orbit balance is $v_c^2=R{\Phi }_R(R,0)$. Using $J_0^{\prime }(u)=-J_1(u)$ in the Hankel representation of thin-disk gravity givesDifferentiation under the integral is justified by convergence for regular models, or by first working off the disk and taking the limiting radial field.

Write $C={\Sigma }_0{R}_0$. For the Mestel disk, the supplied Bessel function integral gives $S(k)=-2\pi GC/k$ for $k>0$. The circular speed is thenMeanwhile the interior mass is $M(R)=2\pi {\int }_0^R(C/s)sds=2\pi CR$. ConsequentlyThe Mestel disk has a flat galaxy rotation curve, despite its nonspherical geometry.

Its infinite total mass prevents setting the Newtonian gravitational potential to zero at infinity. Indeed the integral for $\Phi$ diverges at $k=0$ because $S(k)\sim 1/k$. Subtracting the potential at a reference point removes this additive divergence. As an explicit check, the Mestel disk potential-density pair isIt is harmonic for $z\ne 0$, its vertical derivative jump is $4\pi GC/R$, and $R{\Phi }_R(R,0)=2\pi GC$. Equivalently, an Abel regularization of an oscillatory integral with damping factor $e^{-ak}$ in the radial-field integral gives $R{\int }_0^{\infty }{e}^{-ak}{J}_1(kR)dk=1-a/\sqrt{a^2+R^2}$, which tends to one as $a\downarrow 0$. The force is well-defined although the unreferenced potential is not.

The equality with $GM(R)/R$ is special to the Mestel disk, not a thin-disk shell theorem. The spherical shell theorem normally makes that formula possible, but does not apply to an astrophysical disk. Exterior annuli exert outward radial gravitational acceleration inside their radii, while interior annuli do not act exactly as point masses at the origin. For this particular scale-free surface density, their combined contributions happen to produce the same result as the interior mass expression. Changing the radial surface density or truncating the disk generally destroys the equality: enclosed mass does not determine a disc rotation curve.

Irrotational flow permits $\mathbf{u}=\nabla \varphi$, and incompressible flow imposes $\nabla \cdot \mathbf{u}=0$, hence $\Delta \varphi =0$. In the laboratory frame, the no-penetration condition at the translating sphere is ${\partial }_r\varphi (a,\theta )=U\cos \theta$.

For separation of variables, take the angular factor $g(\theta )=\cos \theta$, which satisfies $\frac{1}{\sin \theta }(\sin \theta g^{\prime }{)}^{\prime }=-2g$. The radial equation iswith solutions $r$ and $r^{-2}$. Rest at infinity excludes the $r$ term. The normal boundary condition fixes the remaining coefficient, giving the potential flow around a translating sphereThere is no azimuthal velocity. With fluid density ${\rho }_f$, the kinetic energy isThe angular integral is $2\pi$, and the radial integral is $1/(3a^3)$, so $K_f={\rho }_f\pi a^3{U}^2/3$. Since $M=(4\pi a^3/3){\rho }_s$ and $k={\rho }_f/{\rho }_s$, $K_f=kM{U}^2/4$. Thus the added mass of a sphere is $kM/2$.

During a fall through $h$, gravity does work $Mgh$ and buoyancy does work $-kMgh$. Starting from rest, conservation of energy givesConsequently $U=\sqrt{4(1-k)gh/(2+k)}$. The mass cancels: the reduced speed comes from both buoyancy and the fluid's added mass.

The rotational invariance of a central-potential Hamiltonian allows a simultaneous eigenstate of energy and the two-dimensional orbital angular momentum $L_z=-i\hslash {\partial }_{\varphi }$. For a separated wavefunction $F(r)G(\varphi )$, the Laplacian in polar coordinates givesThe $\varphi$ term must be a constant; write $G^{\prime \prime }/G=-k^2$. The angular eigenfunctions can be chosen as $e^{ik\varphi }$. Single-valuedness under $\varphi \mapsto \varphi +2\pi$ imposes $e^{2\pi ik}=1$, hence $k\in \mathbb{Z}$. The radial equation in the separation of a two-dimensional central-potential eigenstate isFor a real central potential and the usual real self-adjoint radial boundary conditions, the radial equation admits a basis of real solutions: real and imaginary parts of a complex solution obey the same equation and boundary conditions. Choose a real normalized radial eigenfunction. Since the plane area element in plane polar coordinates is $rdrd\varphi$, the normalization isThis establishes the intended separated simultaneous eigenstate form. It is not the form of every stationary state. The radial equation depends on $k^2$, so the $k$ and $-k$ sectors have the same energy. For $k\ge 1$, their normalized superpositionis a single-valued stationary state with that energy but is not a single angular exponential. Quantum degeneracy is precisely why separation of variables selects a convenient eigenstate basis rather than all vectors in an energy eigenspace. The printed assertion needs this qualification. The printed polar-coordinate aid also labels a gradient component tuple as a divergence; the gradient used below is the two-dimensional vector $\hat{\mathbf{r}}{\partial }_r+\hat{\mathbf{\varphi }}{r}^{-1}{\partial }_{\varphi }$.