Singularity formation and self-shrinkers for the parabolic Bernoulli problem

We study singularity formation for classical solutions of the parabolic Bernoulli free boundary problem and the associated self-similar profiles. Inspired by Ricci and mean curvature flows, we introduce a Type I assumption on the blow-up rate at a singular free boundary point, under which the parabolic rescalings are compact and every tangent flow is a non-trivial self-shrinking solution of the elliptic profile equation. The proof establishes convergence of the positivity sets and recovers the Bernoulli condition in the limit, despite the absence of a minimizing structure. In the subcritical regime, every tangent flow is a double plane. We next classify all radial self-shrinkers. Besides the known ball and exterior profiles, there is a unique annular profile. We obtain dimension-uniform bounds and sharp asymptotics for its inner and outer radii, including a quantitative outward bias of its midpoint. These delicate estimates yield a complete classification of the signs of the linear spectrum for the annular solution in every dimension. As a consequence, the ball is dynamically stable modulo ambient symmetries, whereas the annulus is dynamically unstable, with genuine unstable modes of angular degrees $0,1,2,3$. We finally prove, from the spectral nondegeneracy of the compact radial profiles, that whenever one tangent flow is a ball or an annular solution, the tangent flow is unique. One borderline spectral sign in dimensions $2\le n\le 4000$ is verified by a rigorous computer-assisted argument using interval arithmetic.