You're watching a sinkhole open in the middle of a river. The water is flowing — the whole river is moving — but locally, right at the sinkhole, material is collapsing inward. Now here's the question that makes this hard: where does the sinkhole end and the river begin? The boundary isn't sharp because the river's current affects the collapse and the collapse affects the flow. Every textbook sinkhole model assumes it sits in a still pond (Schwarzschild in vacuum). But real sinkholes form in moving water. This paper builds the mathematical model of a sinkhole in a flowing river — a collapsing over-dense region that smoothly, asymptotically matches onto an expanding Friedmann-Robertson-Walker cosmology. The horizons, the singularity structure, and the causal geometry all change when the river is moving. The committed claim: exact solutions within the Lemaître-Tolbert-Bondi (LTB) dust collapse family that represent collapsing over-dense regions extending asymptotically to closed, open, or flat FRW cosmological models, with detailed analysis of the central singularity strength, nakedness, and the existence of both apparent and event horizons. Category of result: exact analytical solutions in general relativity. This is not a first — cosmological black holes and LTB collapse have been studied for decades — but the specific construction of solutions that smoothly interpolate between gravitational collapse and all three FRW geometries (k = -1, 0, +1), with a systematic treatment of horizon structure and singularity properties, advances the state of the art. The ladder here is the Schwarzschild black hole — the standard textbook solution that describes a black hole in otherwise empty space. Every astrophysics student learns Schwarzschild first, and most intuitions about event horizons, singularities, and causal structure come from that model. But Schwarzschild assumes asymptotic flatness: space is empty and static far from the hole. The real universe is expanding, filled with matter, and has a cosmological history. The Oppenheimer-Snyder collapse model (1939) is the classical prior art for collapse in a cosmological context, but it's highly idealized — uniform density, sharp boundary. LTB models generalize this to inhomogeneous dust collapse, and papers by Joshi, Dwivedi, Singh, and others have explored naked singularities in LTB settings. Firouzjaee and Mansouri's contribution is constructing LTB solutions where the density profile smoothly transitions from a collapsing over-density to an expanding FRW background across all three curvature cases, then carefully analyzing what happens to the horizons and singularities in that setting. Architecturally, this is classical general relativity — exact solutions of Einstein's field equations for pressureless dust (the LTB metric), which is a spherically symmetric, inhomogeneous generalization of FRW. No numerical relativity, no perturbation theory, no simulations. The method is constructing specific functional forms for the LTB free functions (energy function and bang-time function) that produce the desired asymptotic behavior. The analytical tractability comes from the dust assumption (zero pressure), which makes the Einstein equations integrable. This is the same family as Oppenheimer-Snyder but with inhomogeneous initial data designed to smoothly match cosmological expansion at large radius. Integrity is clean for what it is — mathematical physics. The results are exact solutions with analytical proofs about singularity strength (using the Tipler and Królak definitions), horizon existence, and causal structure. There's no experimental validation possible for these specific solutions — they're idealized models of a physical situation (gravitational collapse in an expanding universe). The 6 figures across 9 pages likely show density profiles, apparent horizon trajectories, and Penrose-Carter diagrams. The internal consistency check is whether the solutions actually satisfy Einstein's equations and whether the asymptotic matching to FRW is genuinely smooth. The weakness of all LTB models is the dust assumption: real astrophysical collapse involves pressure, radiation, angular momentum, and non-spherical perturbations. The results are exact but for an idealized matter model. The milestone to track: the field needs cosmological black hole solutions with pressure. Dust collapse is analytically tractable but physically incomplete — pressure gradients, shock formation, and thermodynamic effects change the singularity and horizon structure qualitatively. If someone constructs exact or high-precision numerical solutions for a perfect fluid (not just dust) collapsing region that asymptotically matches FRW, with the same level of causal analysis this paper provides, that would be the next real unlock. The technical barrier is that adding pressure makes the Einstein equations non-integrable in general. Estimated timeline: this has been an open problem for decades, so 'next 5-10 years with numerical methods' is honest. The obvious experiment they didn't run: numerical evolution of the same initial data with pressure included, to see whether the naked singularity results survive beyond the dust idealization. My read: (a) this is a different skill set and computational infrastructure. Firouzjaee and Mansouri are doing exact-solution mathematical GR; numerical relativity is a different community with different tools (Cactus, Einstein Toolkit, SpEC). The paper establishes the analytical scaffolding; the numerical community can now test robustness. It's a natural division of labor, not an omission.
