Gravitational Collapse for Polytropic Gaseous Stars: Self-similar Solutions

The paper rigorously proves the existence of global, real-analytic Yahil-type self-similar solutions for 1<γ<4/3 with a single sonic point, strict monotonicity, and the velocity bound −2/3 y < u < 0 (Theorem 1.3), via a careful sonic-point Taylor analysis, LPH-branch selection, invariant-region arguments to the right, and a shooting construction to the left that identifies a critical sonic time ȳ*; see the formulation of (1.13), boundary conditions (1.14)–(1.15), and methodology outline (Steps 1–4) in the paper. By contrast, the model’s Phase-2 solution has multiple flaws: (i) it asserts local analytic existence at y=0 from standard analytic ODE theory despite the explicit 1/y singularity in ω′, which the paper treats delicately (the paper only establishes C1 regularity at y=0 and remarks that smoothness can be obtained with extra work); (ii) it gives an incorrect far-field decay exponent for 2−γ−ω, claiming α=γ/(2−γ), whereas the paper proves the sharp asymptotic y^{1/(2−γ)}(2−γ−ω)→k̄1>0; and (iii) its transonic-passage/shooting argument (via dτ=dy/D and an intermediate-value argument in ρ0) is not justified at the level of detail required by the paper, which instead uses a precise sonic Taylor expansion on the LPH branch, invariant structures, and a carefully defined fundamental set Y leading to ȳ*. Therefore the paper’s argument is correct and complete for the stated result, while the model’s solution is partially incorrect/incomplete.