{"id":"4426240b-7908-491c-816b-18fdbdd06ccf","arxiv_id":"2501.03178","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A pseudo-Newtonian simulation of accretion onto a Reissner-Nordstrom naked singularity shows the disk thickening at the maximum of the orbital frequency and pushing matter over the top into a torus near the zero-gravity radius.","lead":"We ran the first simulations of a thin gas disk around a highly charged point mass described by general relativity as a naked singularity. The gas piles up into a rotating ring near a special radius, offering a possible observational signature if such objects exist.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Pseudo-Newtonian potential unvalidated in the thick-disk regime; the over-the-top accretion and torus may be approximation artifacts.","rationale":"The reader's weakest assumption—that the pseudo-Newtonian potential faithfully captures the relevant RN gravity for thin-disk accretion, including vertical structure and behavior near the zero-gravity sphere—is the most load-bearing concern. The central mechanism is a 3D thick-disk effect (bulge formation and high-latitude flow) occurring in a regime where the approximation has not been validated. The paper's own statement in Section 5 that the calculations may be superseded by full GR is an admitted limitation, not a resolved one. A quantitative comparison to exact GR equilibria or, better, a full-GR hydrodynamical rerun, would settle whether the over-the-top accretion and the torus are physical. I considered whether the isotropic α-viscosity implementation or the lack of a convergence study could be more decisive, but those concern the numerical realism of a standard disk model; the pseudo-Newtonian potential is the foundational modeling assumption on which the central claim rests. Since the reader already flagged this and the verdict was CONDITIONAL, my read does not change the verdict.","tokens_in":10056,"tokens_out":23915,"duration_ms":236520,"concrete_test":"Run the same q=1.5, α=5e-3, thin-disk setup in a general-relativistic hydrodynamics code on the RN spacetime (e.g., the code used by Kluźniak & Krajewski 2024) and compare the meridional mass-accretion-rate profile dMdot/dθ and the density map with Figs. 2 and 4. If the inward flux at latitudes θ<80° or θ>100° inside r<4/3 r0 does not persist, or if the torus radial extent differs from the pseudo-Newtonian result by more than the local grid scale, the central claim is not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that the disk thickens at the zero-torque circle and accretes over the top at high latitudes, forming a torus near the zero-gravity sphere—is computed entirely with Newtonian hydrodynamics in PLUTO under the pseudo-Newtonian potential V(r) = -M/r + Q^2/(2r^2) (Eq. 5). That potential is constructed to reproduce the radial Keplerian frequency Ω_RN(r), and because it is spherically symmetric it also reproduces the midplane vertical epicyclic frequency κ_z = Ω_RN for any RN-like circular orbit. However, the mechanism operates in the thick-disk regime (h/r ~ 0.5 at the inner edge, as the paper states in Section 4), where the off-midplane shape of the effective potential and the relativistic corrections to the fluid momentum equation (lapse, connection coefficients, energy normalization) enter. These are not captured by Eq. 5. The paper's own limitation in Section 5 ('until such time as our pseudo-Newtonian calculations are superseded by ones in full GR') acknowledges the vulnerability. The comparison to the GR fluid equilibria of Mishra et al. (2024a) is only qualitative ('reminiscent'), so the potential is unvalidated precisely in the region where the central mechanism operates. If the effective potential for a thick torus in the RN metric differs substantially from U(r) in Eq. 6—for example, in the radial extent of the bulge or in the location of the inner edge of the torus—then both the high-latitude accretion path and the quasi-toroidal ring could be numerical artifacts of the approximation rather than features of accretion onto a naked singularity.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents the first pseudo-Newtonian simulations of a thin accretion disk around a Reissner-Nordström naked singularity. The gravity is modeled by the potential V(r) = -M/r + Q^2/(2r^2), which reproduces the RN Keplerian orbital frequency and its maximum at r = 4r0/3. Using the PLUTO code in axisymmetric 2D Newtonian hydrodynamics with α-viscosity, the authors follow a disk initially at r > 10M for several values of q > sqrt(9/8), focusing on q = 1.5. They find that the disk thickens near the radius of maximum angular frequency, that accretion continues inward at high latitudes through this bulge, and that matter accumulates in a quasi-toroidal structure near the zero-gravity radius r0. They suggest that such a rotating ring could be an observational signature of a naked singularity and that, for charges close to unity, the inner edge of the structure would lie well inside the Schwarzschild ISCO.","tokens_in":10372,"tokens_out":7851,"duration_ms":80418,"significance":"If the central mechanism is correct, the paper resolves a conceptual obstruction: a thin disk with outward angular-momentum transport cannot cross the radius where the viscous torque vanishes, yet the simulations show accretion continuing over the top of a pressure-supported bulge. This is an interesting and nontrivial scenario with a testable prediction (a compact rotating ring). The authors are appropriately modest in Section 5 about the pseudo-Newtonian approximation, and Fig. 4 directly addresses the angular-momentum transport claim by comparing fluxes in different polar-angle wedges. However, the load-bearing assumption that Eq. (5) captures the relevant vertical and thick-disk dynamics is not independently validated, and no grid-convergence study is reported. The paper is therefore a useful first step rather than a definitive statement about RN accretion; with added validation it could be of real interest to the accretion and alternative-gravity communities.","major_comments":[{"comment":"The pseudo-Newtonian potential in Eq. (5) is constructed specifically to reproduce the RN Keplerian frequency, so the simulated equatorial Ω(r) matching Eq. (4) in Fig. 3 is a consistency check rather than an independent validation. The central mechanism operates in the thick-disk regime, where the paper states h/r ~ 1/2 at the inner edge, and there the off-midplane shape of the effective potential, pressure gradients, and relativistic terms absent from Newtonian hydrodynamics can all alter the result. The authors' own caveat in Section 5 that the calculations may be superseded by full GR acknowledges this vulnerability. I ask for a quantitative test of Eq. (5) in this regime: for example, compare the simulated torus density and velocity field with the full-GR fluid equilibria of Mishra et al. (2024a) beyond the qualitative remark that the structure is 'reminiscent,' or repeat one run with a formulation that includes the RN lapse and connection terms. Without such a test, the over-the-top accretion and the torus could be artifacts of the pseudo-Newtonian approximation.","section":"Section 2, Eq. (5); Section 4, Fig. 3"},{"comment":"No grid-convergence study is reported. All runs use a single computational grid with R × θ = 217 × 200 cells, and the central conclusions (the bulge, over-the-top accretion, and the torus) are drawn from that one resolution. The claim that the disk 'negotiates' the zero-torque radius by thickening should be demonstrated to persist with higher resolution, especially in the polar-angle direction where the bulge is resolved. I also ask for a time-convergence or steady-state assessment: the averages are taken over t ∈ [19000, 21000] tg, but the text says the inner disk 'builds up' during the simulation, so the structure may still be evolving rather than quasi-steady.","section":"Section 3"},{"comment":"The low-density, non-rotating background fluid is present throughout the grid, and the authors explicitly show that it accretes radially from both sides into a spherical shell around r0 (Fig. 2, right panel, and Fig. 3, left panel). Because this shell forms at the same radius where the quasi-toroidal accretion structure accumulates, the paper should quantify whether the inner torus and the high-latitude fluxes in Fig. 4 are contaminated by the background. A direct test would be to recompute Mdot and Jdot after masking the background (for example, by a density or angular-momentum threshold), or to run at least one case with a much lower background density. The statement in Section 4 that the shell is 'not related to the accretion disk at all (at least at the time interval used for computing the average)' is not backed by a quantitative separation of the two components.","section":"Section 4 and Fig. 2"}],"minor_comments":[{"comment":"There is a typo in the text: 'Fiq. 3' should be 'Fig. 3'.","section":"Section 4"},{"comment":"The text refers to a 'wedge of azimuthal angle range [68°,112°],' but in spherical coordinates θ is the polar/co-latitudinal angle; the terminology should be made consistent with the Fig. 4 caption, which uses 'co-latitudinal intervals.'","section":"Section 4 and Fig. 4"},{"comment":"The sentence 'in the co-latitudinal direction we set a uniform grid with two different resolutions' is confusing; the grid has uniform spacing within each of three zones but different spacings between zones. Please rephrase to describe the three-zone angular grid more clearly.","section":"Section 3"},{"comment":"It would be helpful to state explicitly that the spherical potential also reproduces the midplane vertical epicyclic frequency (because in a spherical Newtonian potential the vertical frequency equals the orbital frequency), while noting that the off-midplane effective potential and the missing lapse/connection terms are the main sources of uncertainty. This would clarify the scope of the approximation.","section":"Section 2, Eq. (5)"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses a genuinely interesting question and the numerical setup is reasonable for a first exploration. My main concern is that the pseudo-Newtonian potential is validated only for the radial Keplerian frequency, while the claimed mechanism depends on vertical/thick-disk behavior at h/r ~ 0.5. A quantitative comparison with full-GR equilibria or a GR simulation of one case would substantially strengthen the case. A grid-convergence study and a background-contamination test are also necessary before the ring/over-the-top scenario can be considered robust. I would support publication after these additions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear [Colleague],\n\nThe one thing to know: this is a genuine first—the first numerical simulation of a thin accretion disk around an RN naked singularity—and the central result is a plausible mechanism for how the disk negotiates the vanishing viscous torque: it thickens into a bulge at the maximum of the orbital frequency and accretes over the top at high latitudes. The flux decomposition in Fig. 4 does support that claim, and the paper is honest about its main weakness, which is the pseudo-Newtonian potential itself. The potential V(r) = -M/r + Q^2/(2r^2) is a simple construction that exactly reproduces the RN Keplerian frequency, and the simulated Omega(r) matches that prediction well. That is not circular in a damaging sense for the equatorial profile, since the simulation is doing real hydrodynamics; the circularity concern is limited to the fact that the potential was built to produce that frequency.\n\nThe soft spot is real and the authors name it: the mechanism operates in a thick-disk regime (h/r ~ 0.5) where the off-midplane shape of the potential and the relativistic corrections to the momentum equation matter, and those are not captured by the spherical pseudo-Newtonian potential. The comparison with the GR equilibrium figures is only qualitative. So the over-the-top accretion and the torus could shift quantitatively, or even disappear, in full GR. That is a limitation, not a refutation: the paper is a proof-of-concept in a simplified potential, and the authors say so in Section 5.\n\nWhat is missing on the reproducibility side is more annoying: no code or parameter files, no grid convergence study, and only one resolution. For a numerical paper, that is a legitimate referee request. I would not desk-reject; I would send it to a referee with instructions to ask for the convergence study and the code release, and to probe the potential's behavior off the midplane.\n\nWho this is for: people working on accretion disk theory, compact-object observables, and the naked singularity debate. They will get a concrete, usable potential and a new mechanistic picture even if the quantitative details need GR confirmation.\n\nMy recommendation: accept for peer review, with substantive revisions expected. It deserves serious referee time. I'd bring it to reading group and I'd cite it if I worked in this area.","headline":"First thin-disk simulation around an RN naked singularity gives a plausible over-the-top accretion mechanism, but the pseudo-Newtonian potential needs GR validation and the numerics need a convergence study.","tokens_in":10919,"tokens_out":2420,"would_cite":true,"duration_ms":20264,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A thin accretion disk can keep feeding a Reissner-Nordström naked singularity past the radius where viscous torque vanishes, by bulging over the top and forming a rotating ring near the zero-gravity sphere.","keywords":["accretion disks","naked singularity","Reissner-Nordström spacetime","pseudo-Newtonian potential","numerical simulation","zero-gravity radius","toroidal structure","thin disk"],"falsifier":"A concrete check is to compare the vertical and epicyclic frequencies of the pseudo-Newtonian potential against the exact RN values at the bulge radii; a significant mismatch would mean the bulge and ring are artifacts of the potential. A full general-relativistic hydrodynamic simulation of the same thin disk, run to the same times, would settle whether the over-the-top accretion and the toroidal structure survive in the true RN spacetime.","tokens_in":9824,"feed_emoji":"🌀","tokens_out":8233,"duration_ms":58922,"temperature":0.7,"pith_summary":"The paper reports the first numerical simulations of a thin accretion disk around a Reissner-Nordström naked singularity, using a pseudo-Newtonian potential that exactly reproduces the RN Keplerian orbital frequency. The central difficulty is that the orbital frequency peaks at $r = 4r_0/3$ (where $r_0 = Q^2/M$ is the zero-gravity radius), so the viscous torque, which drives outward angular momentum transport, vanishes there. The simulations find that the disk does not stall: it thickens into a bulge at the zero-torque circle, and fluid is pushed over the top at higher latitudes, allowing accretion to continue inward. Material eventually accumulates in a rotating toroidal structure near the zero-gravity sphere. If real, such a ring could be an observational signature distinguishing a naked singularity from a black hole.","feed_headline":"Thin disk climbs over a zero-torque ridge to feed a naked singularity","feed_subtitle":"First simulations show how matter climbs over the stalled orbit and forms a ring around the singularity.","key_machinery":"The central object is the pseudo-Newtonian potential $V(r) = -M/r + Q^2/(2r^2)$, which reproduces in Newtonian mechanics the exact radial dependence of the RN Keplerian orbital frequency $\\Omega_{RN}(r) = \\sqrt{1 - r_0/r}\\,\\sqrt{M/r^3}$ for all $r > r_0$, including the zero at $r_0 = Q^2/M$ and the maximum at $4r_0/3$. The effective potential $U(r) = V(r) + \\ell^2/(2r^2)$ yields the same circular-orbit relation, so the simulated disk follows the test-particle frequency. Because the viscous stress is proportional to $d\\Omega/dr$, the sign change of $d\\Omega/dr$ is what creates the zero-torque obstacle, and the disk's vertical response to the resulting pressure build-up is the mechanism that carries the over-the-top accretion. The equations are evolved with Newtonian hydrodynamics and an $\\alpha$-viscosity prescription.","core_discovery":"The central discovery is that a thin disk around an RN naked singularity negotiates the vanishing viscous torque by changing its shape. At the radius of maximum orbital frequency, $r = 4r_0/3$, the disk thickens, and accretion proceeds across that radius at high latitudes, well above and below the equatorial plane. Beyond the maximum, the flow continues down to the vicinity of the zero-gravity sphere, where matter gathers into a quasi-toroidal, rotating structure. The simulated angular velocity profile matches the test-particle Keplerian frequency, its maximum sits close to $4r_0/3$, and the inner torus resembles fluid figures of equilibrium derived analytically in full general relativity.","pith_inferences":["The same 'bulge and go over the top' mechanism should appear in any axisymmetric potential whose orbital frequency has an interior maximum, so it may generalize to other exotic compact objects or modified gravity models.","The paper does not compute radiation from the torus; computing its luminosity and spectrum would give a concrete prediction testable with current interferometers.","If the vertical gravity of the pseudo-Newtonian potential differs from RN, the height of the bulge and the ring location could shift; a full-GR run is the direct test.","Because the simulation uses a non-radiative, purely hydrodynamic disk, the role of magnetic stresses (MRI) in the over-the-top crossing remains an open question."],"forward_implications":["Accretion onto a naked singularity can continue even where the viscous torque vanishes, through high-latitude flow over a thickened disk bulge.","For charge-to-mass ratios $q \\gtrsim 1$, the inner edge of the accretion structure lies inside the Schwarzschild ISCO, and the maximum orbital frequency is much higher than the Schwarzschild ISCO frequency.","The rotating ring near the zero-gravity sphere could serve as an observational signature that distinguishes an RN naked singularity from a black hole in horizon-scale images.","The simulated torus matches the shape of fluid figures of equilibrium obtained in full GR, suggesting the pseudo-Newtonian model captures the equilibrium geometry of the inner flow."],"supporting_citations":[{"why":"Derives the RN Keplerian orbital frequency and its maximum at $4r_0/3$, which the pseudo-Newtonian potential is designed to reproduce.","marker":"Pugliese et al. 2011"},{"why":"Identifies the zero-gravity radius $r_0$ and the levitating atmosphere, used to interpret the shell and the inner ring.","marker":"Vieira & Kluźniak 2023"},{"why":"Provides the analytic fluid figures of equilibrium for RN naked singularities that the simulated toroidal structure is compared against.","marker":"Mishra et al. 2024a"},{"why":"Earlier analytic equilibria for uniform angular momentum that motivate the shape of the inner structure.","marker":"Mishra & Kluźniak 2023"},{"why":"The hydrodynamic code used to evolve the disk in the simulations.","marker":"Mignone et al. 2007"},{"why":"The $\\alpha$-viscosity prescription that transports angular momentum outward in the disk.","marker":"Shakura & Sunyaev 1973"}],"fun_headline_variants":["Disk thickens at stalled orbit to feed naked singularity","Naked singularity fed by disk that flares over zero-torque radius","Simulation shows toroidal ring around charged naked singularity","Disk overcomes zero-torque ridge by thickening, forms ring","First simulation: disk flares to feed ring around naked singularity"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the pseudo-Newtonian potential $V(r) = -M/r + Q^2/(2r^2)$ with Newtonian hydrodynamics correctly captures the gravity of the RN naked singularity for a thin disk, in particular the vertical structure of the disk and the behavior of the flow near the zero-gravity sphere.","fun_headline_variants_meta":{"raw":{"variants":["Disk thickens at stalled orbit to feed naked singularity","Naked singularity fed by disk that flares over zero-torque radius","Simulation shows toroidal ring around charged naked singularity","Disk overcomes zero-torque ridge by thickening, forms ring","First simulation: disk flares to feed ring around naked singularity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000702,"raw_usage":{"total_tokens":3160,"prompt_tokens":931,"completion_tokens":2229,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":547,"completion_tokens_details":{"reasoning_tokens":2145}},"tokens_in":547,"tokens_out":2229,"duration_ms":87453,"temperature":1.0,"reasoning_tokens":2145,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T21:52:38.104400+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A concrete check is to compare the vertical and epicyclic frequencies of the pseudo-Newtonian potential against the exact RN values at the bulge radii; a significant mismatch would mean the bulge and ring are artifacts of the potential. A full general-relativistic hydrodynamic simulation of the same thin disk, run to the same times, would settle whether the over-the-top accretion and the toroidal structure survive in the true RN spacetime.","supporting_citations":[],"review_version":1}