Past exam of the mathematics course of the University of Cambridge / 2019 / iii / Paper 317 / 4

For $P=K{\rho }^2$, hydrostatic equilibrium gives $2K\nabla \rho =-\nabla \Phi$. The Poisson equation for Newtonian gravity therefore gives the Helmholtz equationFor spherical symmetry, the regular polytrope of index one solution isThe singular $\cos (kr)/r$ solution is excluded by finite central mass density. Taking the first zero as the stellar surface yieldsThe mass integral isDividing by the volume $4\pi R^3/3$ gives

The same polytrope of index one interior equation is ${\nabla }^2\rho +(2\pi G/K)\rho =0$. A positive separated solution with zero density on all six cube faces is the cubic polytropic interiorIts Laplacian is $-3{\pi }^2\rho /L^2$, so the required side length isWith $\Phi =C-2K\rho$ and $P=K{\rho }^2$, it satisfies the interior hydrostatic equilibrium and Poisson equation for Newtonian gravity. Integrating each sine factor yieldsThis formal interior solution is not an isolated physical cubic star. The cubic polytrope fails isolated gravitational matching: at a vertex, the product of sines has $\nabla \rho =0$, so the interior potential predicts zero gravitational acceleration. At the vertex $(0,0,0)$, however, the field generated by its own positive mass isand each component is strictly positive. There can be no continuous matching to the isolated external field without additional forces or mass sources. Thus solving the interior density equation and imposing zero face values is insufficient. Fluid stars also have no rigid structure to maintain sharp cubic faces, and observed stellar shapes are approximately spherical or rotationally flattened rather than cubic.