{"id":"f620fc6d-cb3a-4158-9532-e405c94b6681","arxiv_id":"2504.17863","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Generic scalar perturbations near integrable singularities have non-integrable energy density, and extended objects cannot traverse them safely, challenging their viability as black hole alternatives.","lead":"This paper argues that 'integrable singularities', proposed mild black hole centers, are probably not safe to cross. Small wiggles in matter fields pile up with infinite energy, and real objects would be torn apart or heated enormously.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The non-integrability conclusion for l>0 rests on the unproven assumption that the divergent Frobenius solution is generically excited by smooth initial data; the paper itself concedes this possibility is not ruled out, and its l=0 analysis shows how discarding a solution changes the result.","rationale":"The reader's weakest assumption is exactly the load-bearing point. The Frobenius analysis gives local solutions, but a linear wave equation selects a unique solution from initial data; the divergent-mode coefficient is a linear functional of that data. The paper's conclusion that rho is non-integrable for l>0 only follows if that functional is nonzero generically. The paper acknowledges the opposite possibility in footnote 10. The l=0 subsection illustrates the stakes: after writing the general solution (V.15) with a Sigma b_n r^n term, the text keeps only r + r ln r, thereby discarding the Psi ~ const solution that would give Phi ~ 1/r and non-integrable spherical modes. Thus the same selection issue that could destroy the l>0 conclusion is already visible in the paper's l=0 treatment. A numerical evolution of the actual wave equation in the interior is the cleanest way to settle the issue. If the dominant mode is generically excited, the paper's l>0 claim stands, and the l=0 sector may be even worse than stated; if it is not, the central non-integrability argument collapses. This concern does not change the reader's CONDITIONAL verdict, because the concern is explicit and the paper itself flags it.","tokens_in":17874,"tokens_out":25424,"duration_ms":257098,"concrete_test":"Numerically solve the linear massless scalar wave equation in the interior of the integrable metric (II.8) with C0=3, M=1, using ingoing double-null coordinates. Launch smooth compact-support data for the l=1 (and l=0,2) modes across the horizon, evolve toward r=0, and fit the asymptotic radial profile to Psi(r)=A_+ r^{R1}+A_- r^{R2} (or the corresponding log solutions at l=0) over a window where higher orders are subdominant. Repeat for several pulse widths and centers. If the dominant divergent coefficient A_+ is nonzero for any generic family, the non-integrability claim is supported; if it vanishes identically for all such data, the conclusion fails. This directly tests the paper's unproven generic-excitation assumption and also reveals whether l=0 modes have the dropped Psi ~ const component.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that l>0 scalar perturbations have non-integrable energy density is a statement about the solution selected by Cauchy data, but the calculation in Sec. V.A only describes the two Frobenius solutions at r=0. Because r=0 is a future focusing boundary rather than a place where boundary conditions can be imposed, the coefficient of the divergent r^{-beta} solution must be shown to be generically nonzero. The paper does not show this; footnote 10 explicitly says the authors 'have not ruled out' that A1,A2 vanish for all l,k. The issue is concrete: in the l=0 sector, the general solution in (V.15) contains a Sigma b_n r^n term, but the text then retains only Psi ~ r + r ln r, discarding the Psi ~ const Frobenius solution whose Phi ~ 1/r would make even spherical modes non-integrable. The same kind of solution selection could in principle suppress the l>0 divergent mode. Thus the non-integrability argument hinges on an excitation assumption that is plausible but unverified.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies spherically symmetric black holes with integrable singularities, characterized by a mass function m(r)=m1 r+O(r^2) near r=0. It argues, first, that such singularities are focusing points whose inclusion in an extension would violate global hyperbolicity; second, that scalar test-field perturbations on these backgrounds have a non-integrable energy density for l>0 modes, whereas the paper claims l=0 modes remain integrable; and third, that extended physical observers would face severe obstacles due to deviations from geodesic motion and the large energy encountered near the singularity. The overall conclusion is that integrable singularity models face serious theoretical and practical challenges as alternatives to regular black holes.","tokens_in":18131,"tokens_out":8705,"duration_ms":80573,"significance":"If the central perturbative claim were established, the paper would be a significant contribution to the current debate on integrable singularities, because it would show that the integrability condition imposed on the background source does not automatically extend to dynamical fields. The manuscript is transparent about its assumptions, uses standard Frobenius and geodesic-deviation machinery, and explicitly flags residual possibilities. However, the main conclusion is currently conditional: the l>0 non-integrability result depends on an unproven genericity assumption about the excitation of the divergent Frobenius solution, and the l=0 analysis contains a concrete error. The paper's significance is therefore real but provisional.","major_comments":[{"comment":"The l=0 Frobenius solution is not handled correctly. The second independent solution in (V.15) is Ψ2 = ln(r) r Σ a_n r^n + Σ b_n r^n, so it contains a constant term b_0. Hence the generic leading behavior is Φ0 = Ψ0/r ∼ A2 b_0/r plus subleading A2 ln r, not Φ0 ∼ ln r as stated in (V.17). The corresponding energy density is ρ ∼ r^{-4}, so ∫ρ r^2 dr diverges; the claim that spherical modes are integrable is therefore wrong. This error does not weaken the general thrust of the paper; if anything, it shows that the constant Frobenius mode would make even l=0 modes non-integrable, but the l=0 discussion as written must be corrected.","section":"Sec. V.A, Eqs. (V.15)-(V.17)"},{"comment":"The central claim that l>0 scalar perturbations have non-integrable energy density depends on the divergent Frobenius solution being the one selected by generic Cauchy data. The manuscript computes only the two local solutions near r=0 and does not compute the excitation coefficients from horizon-crossing initial data; footnote 10 explicitly states that the authors 'have not ruled out' that A1,A2 vanish for all l,k. Because r=0 is a future focusing point rather than a boundary where regularity conditions can be imposed, the paper needs an argument that smooth, generic initial data produce a nonzero coefficient for the r^{-β} mode. Without such an argument, Eqs. (V.24)-(V.25) describe a possible mode, not necessarily the physical field selected by evolution. This is load-bearing: fine-tuned suppression of the divergent mode would invalidate the main conclusion.","section":"Sec. V.A and footnote 10"}],"minor_comments":[{"comment":"The prefactor in (V.22) is incorrect. Using (V.21) with f=1-2m(r)/r and m≈m1 r, the contribution from the radial derivative is 1/2(∂rΦ)^2 and the u^r term contributes m1(∂rΦ)^2, giving ρ ≈ (1+α/2)(∂rΦ)^2, not (1+α)/2(∂rΦ)^2. The divergence rates are unaffected.","section":"Sec. V.A, Eq. (V.22)"},{"comment":"The numerical estimate m(a=10 m)≈8.6×10^27 kg corresponds to C0=1 in (II.8), but the paper does not state the value used. Since Fig. 1 uses C0=3, the same calculation with that value gives roughly 2.6×10^28 kg for M≈10^30 kg; please clarify the choice of C0.","section":"Sec. VI.B"},{"comment":"The formula for the indicial roots is typographically unclear: the term involving α√α is written without indicating the division that makes the argument dimensionless. Please rewrite the roots in an unambiguous form.","section":"Eqs. (V.14) and (A.1)"},{"comment":"The notion of a 'Cauchy point' is introduced informally. Since the argument about the loss of global hyperbolicity relies on this concept, it would be helpful to give a precise definition or a supporting reference.","section":"Sec. IV.B"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the journal's scope and the perturbative idea is interesting. The main risk is the excitation assumption: if a concrete mode-excitation calculation or a precise genericity statement is provided, the paper's central claim would be substantially stronger. The l=0 error is fixable and does not by itself warrant rejection. I would not recommend accept at this stage because the current version overstates the certainty of the main conclusion."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The main new thing here is a Frobenius analysis of scalar test fields on integrable-singularity backgrounds. For l>0 modes the paper shows the field goes like r^{-β} with β≥1/2, making the energy density ρ~r^{-2(β+1)} non-integrable at r=0. That calculation is basically right, and it is a genuine new application to this background class. The paper is also honest: footnote 10 concedes the divergent solution might not be excited for all l,k, though they think fine-tuning would be needed.\n\nThe soft spots are real. First, the l=0 sector is mishandled. The general Frobenius solution (V.15) includes a constant term in Ψ from the second solution; that gives Φ~1/r, not ln r. The paper discards it without argument. If that constant is generically present, even spherical perturbations have non-integrable energy density, which contradicts the Introduction's claim that spherical modes are consistent with integrability. Second, the prefactor in Eq (V.22) is off by 1/2 — the coefficient should be 1/2 + m1, not m1. This doesn't affect integrability but is sloppy. Third, the mass estimate in Sec VI.B uses a length-to-mass conversion that looks wrong; I get a few times 10^28 kg for their example, not 8.6×10^27 kg. These are fixable.\n\nThe larger issue is conceptual. The non-integrability conclusion for l>0 depends on the divergent Frobenius solution being generically excited by smooth initial data at horizon crossing. Because r=0 in the interior is a focusing point rather than a boundary, the authors argue you can't impose regularity and must keep both solutions. That's plausible, and their contrast with stellar cores is well made, but it is still an assumption, and the l=0 error shows how much the result changes when you keep or drop a solution.\n\nWhat the paper does well, besides the core calculation, is frame r=0 as a Cauchy point — a Cauchy horizon shrunk to a point — and the discussion of C2-inextendibility is clear. The MPD and extended-object arguments are more qualitative but reasonable.\n\nI'd send this to peer review. The l>0 result is new and likely correct modulo the excitation assumption; the l=0 error and prefactor mistake are correctable. Just read the l=0 section with care. It deserves a serious referee.","headline":"A genuinely new Frobenius analysis of scalar perturbations on integrable-singularity backgrounds, with a sound l>0 non-integrability result, but marred by an l=0 sector error, a prefactor slip, and a load-bearing excitation assumption that is acknowledged but unproven.","tokens_in":18626,"tokens_out":5942,"would_cite":true,"duration_ms":48329,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83C57","83C75"],"pacs":["04.70.-s","04.20.-q"],"model":"deepseek-v4-flash","headline":"The paper argues that integrable black-hole singularities lose their integrability once generic scalar-field perturbations are included.","keywords":["integrable singularities","scalar field perturbations","Frobenius method","energy density non-integrability","black hole interiors","traversable singularity","conjugate points","extended bodies"],"falsifier":"Take the metric $f(r)=1-2(m_1r+\\dots)/r$ with $m_1>1/2$, place smooth initial data for a massless scalar field on a spacelike surface inside the horizon, and evolve the mode equation numerically for $l=0$ and $l>0$. If for generic data the divergent coefficients, such as the $r^{-1/2}$ or complex-power modes, vanish for every $l>0$, the non-integrability claim collapses; if they are nonzero, the energy-density integral over a small ball should diverge as $\\int_0^\\epsilon \\rho r^2 dr\\sim \\epsilon^{1-2\\beta}$ for $\\beta>1/2$.","tokens_in":17717,"feed_emoji":"🕳️","tokens_out":10768,"duration_ms":99285,"temperature":0.7,"pith_summary":"Black holes whose central singularity is only integrable — the metric stays finite and tidal forces stay finite, but curvature diverges — look like a middle path between Schwarzschild and regular black holes. This paper argues that the middle path closes when perturbations are included. In a spherically symmetric background with mass function $m(r)=m_1 r+O(r^2)$, a massless scalar test field grows like $\\ln r$ for spherical modes but like $r^{-\\beta}$ with $\\beta\\ge 1/2$ for all $l>0$ modes. The energy density of those non-spherical modes scales as $\\rho\\sim r^{-2(\\beta+1)}$, so its integral up to $r=0$ is infinite: the integrability that the background was built to satisfy fails for generic perturbing fields. The paper adds that the singularity is a focusing point for every causal geodesic, and that extended bodies will not follow the finely tuned radial geodesics for which the singularity is weak, so traversability is doubtful.","feed_headline":"Generic waves make 'integrable' singularities unintegrable","feed_subtitle":"Non-spherical scalar modes carry a diverging energy density, so the background's built-in integrability fails under perturbation.","key_machinery":"The load-bearing object is the radial wave equation for a massless scalar on the fixed background, written near $r=0$ as $$\\Psi''+\\frac{P}{r}\\Psi'+\\frac{Q}{$r^{2}$}\\Psi=0,$$ with $P$ and $Q$ obtained from $m(r)=m_1r+m_2r^2+O(r^3)$. Since $r=0$ is a regular singular point, the method of power-series solutions gives indicial roots $R_{1,2}=\\frac12\\left(1\\pm\\sqrt{1-\\frac{4l(l+1)}{\\alpha}}\\right)$ with $\\alpha=2m_1-1>0$, and hence $\\Phi_0\\propto r^{R-1}$. These roots are what convert the background's mild divergence into a field-mode divergence $r^{-\\beta}$, $\\beta\\ge1/2$, and then into the non-integrable density $\\rho\\sim r^{-2(\\beta+1)}$. The same geometry gives the null expansion $\\Theta\\to-\\infty$, making $r=0$ a focusing point and, in any extension, a conjugate point to every point of the trapped region, meaning a point where neighboring geodesics refocus.","core_discovery":"The paper's central claim is that the integrability condition encoded in the background spacetime does not survive contact with test fields. For $l=0$ the scalar field diverges only logarithmically, $\\Phi_0\\propto\\ln r$, and the observer-measured energy density $\\rho\\sim r^{-2}$ is still integrable. For every $l>0$, the power-series analysis of the radial mode equation gives a leading term $\\Phi_0\\propto r^{-\\beta}$ with $\\beta\\ge 1/2$; when $l$ is large the roots become complex and the field oscillates inside an $r^{-1/2}$ envelope. Either way $\\rho\\sim r^{-2(\\beta+1)}$ and $\\int \\rho r^2 dr$ diverges at $r=0$. Because $r=0$ is a moment in time inside the horizon rather than a spatial boundary, the divergent mode cannot be discarded by the regularity boundary condition used in stellar interiors; it can be removed only by fine-tuning the initial data. The paper concludes that generic perturbations would backreact strongly enough to threaten the integrable-singularity background, and that the focusing geometry makes the singularity a universal conjugate point and an effective barrier for extended objects.","pith_inferences":["Beyond the paper, the same regular-singular-point structure should make any massless higher-spin field, or any non-spherical extension of the metric, produce a similar or stronger divergence, since the effect comes from the indicial exponents rather than from scalar-field self-coupling.","If the universal-conjugate-point picture is correct, any extension through $r=0$ would behave like a Cauchy horizon collapsed to a point; one might expect nonlinear instabilities analogous to mass inflation, although the paper does not prove this.","A concrete extension the paper leaves open: evolving the scalar field from smooth initial data inside the horizon would settle whether the divergent coefficients are generically nonzero for all $l$, turning the present claim into a quantitative statement."],"forward_implications":["Any generic massless scalar perturbation with angular number $l>0$ carries a non-integrable energy density toward $r=0$, so the test-field approximation breaks down before the singularity and backreaction must be considered.","The focusing property makes $r=0$ a conjugate point for every causal geodesic in the trapped region, so a spacetime that includes $r=0$ in the manifold cannot be globally hyperbolic.","The weakness of the singularity for radial geodesics does not extend to realistic bodies: spin–curvature coupling and internal stresses send parts of an extended object onto non-radial trajectories for which the singularity is strong.","A falling extended object encounters a finite but enormous mass-energy concentrated near the singularity; for a stellar-mass example the enclosed energy exceeds the rest mass of a 10-meter water sphere by about 21 orders of magnitude.","The paper concludes that, without a backreaction or isotropization mechanism that eliminates non-spherical modes, integrable singularities are not a robust alternative to regular black holes."],"supporting_citations":[{"why":"Supplies the integrable-singularity background $m(r)=m_1r+O(r^2)$ and the trapped-interior condition $m_1>1/2$.","marker":"[39]"},{"why":"Provides the Schwarzschild-interior test-field analysis whose power-series approach is adapted here.","marker":"[57]"},{"why":"Establishes the finite tidal forces and possible traversability of integrable singularities that the paper tests.","marker":"[29]"},{"why":"Introduces the quantum-core picture in which $\\epsilon\\sim r^{-2}$ is treated as normalizable, motivating the integrability condition.","marker":"[31]"},{"why":"Checks that the singularity is deformationally weak for radial geodesics, the property called into question.","marker":"[44]"},{"why":"Supplies the conjugate-point theorem used to show that extending through $r=0$ conflicts with global hyperbolicity.","marker":"[47]"},{"why":"Shows the mass-inflation null singularity is weak for all geodesics, the contrast case for extended-body crossing.","marker":"[65]"},{"why":"Defines deformational strength, the criterion used to distinguish strong from potentially traversable singularities.","marker":"[42]"}],"fun_headline_variants":["Integrable singularities break under generic perturbation","Non-spherical modes render integrable singularities unstable","Higher scalar modes ruin integrable singularity viability","Tidal forces and divergent modes challenge integrable singularities","Integrability fails when test fields hit the singularity"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument assumes that the divergent power-series mode is genuinely excited by smooth, generic initial data at horizon crossing, rather than being absent because its coefficients are fine-tuned to zero for all $l$ and $k$, and that the interior really has $m_1>1/2$ so the trapped region reaches an integrable singularity.","fun_headline_variants_meta":{"raw":{"variants":["Integrable singularities break under generic perturbation","Non-spherical modes render integrable singularities unstable","Higher scalar modes ruin integrable singularity viability","Tidal forces and divergent modes challenge integrable singularities","Integrability fails when test fields hit the singularity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000226,"raw_usage":{"total_tokens":1471,"prompt_tokens":954,"completion_tokens":517,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":570,"completion_tokens_details":{"reasoning_tokens":445}},"tokens_in":570,"tokens_out":517,"duration_ms":5509,"temperature":1.0,"reasoning_tokens":445,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T10:34:09.737958+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the metric $f(r)=1-2(m_1r+\\dots)/r$ with $m_1>1/2$, place smooth initial data for a massless scalar field on a spacelike surface inside the horizon, and evolve the mode equation numerically for $l=0$ and $l>0$. If for generic data the divergent coefficients, such as the $r^{-1/2}$ or complex-power modes, vanish for every $l>0$, the non-integrability claim collapses; if they are nonzero, the energy-density integral over a small ball should diverge as $\\int_0^\\epsilon \\rho r^2 dr\\sim \\epsilon^{1-2\\beta}$ for $\\beta>1/2$.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the integrable-singularity background $m(r)=m_1r+O(r^2)$ and the trapped-interior condition $m_1>1/2$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Schwarzschild-interior test-field analysis whose power-series approach is adapted here."},{"cited_title":"Sbierski, Journal of Differential Geometry 108, 319 (2018)","cited_arxiv_id":null,"evidence_quote":"Establishes the finite tidal forces and possible traversability of integrable singularities that the paper tests."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Shows the mass-inflation null singularity is weak for all geodesics, the contrast case for extended-body crossing."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines deformational strength, the criterion used to distinguish strong from potentially traversable singularities."}],"review_version":1}