{"id":"cb736a59-bef0-481c-a9bd-cfe73688e576","arxiv_id":"2508.20558","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Dyonic black holes with non-monotonic metric functions support multiple geometrically distinct but topologically identical periodic orbits for the same rational label, plus bound orbits with energy above unity.","lead":"Timelike particles around dyonic black holes can share the same periodic-orbit label while following several different radial paths when the metric function is non-monotonic. The paper maps these multi-branch orbits and shows that bound orbits with energy above unity can also exist in such spacetimes.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central multi-branch claim rests on an undocumented numerical evaluation of Eq.","rationale":"I read the paper in good faith: it applies the standard Levin periodic-orbit framework to a known family of dyonic black holes, and the figures and tables are internally consistent with the claimed double-well phenomenology. The E>1 'bound' orbits are physically meaningful as classically confined states in a double-barrier potential, so I do not treat that as a semantic defect. The concern is not about interpretation but about whether the numerical solutions correspond to the actual spacetime. The metric function Eq. (3) is highly nontrivial exactly where the new structure appears, and the paper's silence on hypergeometric evaluation, quadrature, and root-finding accuracy makes the central numerical result unverifiable as written. Independent re-evaluation with arbitrary precision and a second analytic-continuation route is the check that would settle the issue. This is essentially the same load-bearing assumption the Pith reader identified, so the CONDITIONAL verdict remains appropriate; I see no basis to strengthen or weaken the verdict further, and no criticism of the authors beyond the need for documented numerical validation is implied.","tokens_in":48254,"tokens_out":18449,"duration_ms":186995,"concrete_test":"Independently recompute f(r) from Eq. (3) at α1=1.025755, M=1, α2=2.76, p=0.15, q=1.05 for r in [r_h, 1] using two independent methods: (i) arbitrary-precision evaluation of 2F1(1/4,1;5/4;z) (e.g., mpmath at 50 digits), and (ii) the Euler transformation 2F1(1/4,1;5/4;-y)=(1+y)^(-1/4) 2F1(1/4,1/4;5/4;y/(1+y)); verify agreement below 10^-12. Then recompute the MBO roots and the Table II entries for ϵ=0.1, q=(1,2,0), with root-finding tolerance below 10^-10. If the E1/E2 splitting persists and both orbits close with the same (z,w,v), the central claim survives; if the inner barrier or the E1/E2 pair disappears, the claim fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim is an existence claim: for the non-monotonic slice M=1, α2=2.76, p=0.15, q=1.05 and α1∈[1.025755,1.07317), the effective potential develops a double-barrier that supports several distinct periodic orbits for the same rational q, as well as confined orbits with E>1. The least secure link is the numerical evaluation of the metric function Eq. (3) in the region where the inner barrier lives. The hypergeometric argument is z=-4p²α2/(r⁴α1), which is strongly negative in the relevant near-horizon region: at r≈0.2, |z|≈155, and at r≈0.15, |z|≈490, far outside the radius of convergence |z|<1 of the defining series. The non-monotonicity of f(r), and therefore the double-well structure of Veff, depends on a correct analytic continuation of 2F1 on this sheet. The paper states no algorithm, tolerance, or validation for this continuation. The claimed branch splitting is also numerically delicate: in Table II, for ϵ=0.1 and q=(1,2,0), E1=0.917706569 and E2=0.918695150 differ by only about 10^-3, and these two nearby roots are the direct evidence for 'multiple branches'. If the hypergeometric evaluation is wrong at the 10^-4 level near the inner barrier, the double-barrier itself—and hence the multi-branch periodic orbits and the E>1 bound orbits—could be artifacts of the evaluation scheme rather than properties of the spacetime. The same fragility extends to the root finding for Eqs. (10)–(11) and the improper integral Eq. (15), for which no error estimates or convergence checks are supplied.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies bound periodic orbits of massive test particles in a static, spherically symmetric dyonic black hole arising from quasi-topological electromagnetism, with metric function given by Eq. (3). For fixed M=1, α2=2.76, p=0.15, q=1.05 and varying α1, the authors classify horizon structures and identify a non-monotonic regime 1.025755≤α1<1.07317 in which f(r) is non-monotonic outside the horizon and the effective potential develops a double-barrier structure. Within this regime they report three marginally bound orbit (MBO) branches, two of which (LMBO2 and LMBO3) are used to construct double-well effective potentials, and they present numerical evidence for multiple energy branches of periodic orbits with the same rational number q=(z,w,v), including bound orbits with E>1. For the monotonic regimes (3- and 4-horizon, and the second 1-/2-horizon range), they recover a single MBO and a single periodic orbit branch per q, consistent with previous results. The paper combines the standard Levin et al. periodic-orbit classification with a specific modified-gravity black hole family and provides extensive tables and trajectory figures.","tokens_in":48603,"tokens_out":8000,"duration_ms":78206,"significance":"If the numerical claims are robust, the paper reports a genuinely new orbital phenomenon: a non-monotonic metric function outside the horizon produces a double-barrier effective potential, multiple MBOs, multiple periodic-orbit branches sharing the same rational label q, and E>1 orbits confined by an outer barrier. These features are absent in Schwarzschild, Kerr, and Reissner-Nordström backgrounds and could be relevant for strong-field tests of modified gravity and for EMRI waveform modeling. The paper is clearly structured, uses standard geodesic and periodic-orbit formalism, and provides a large set of explicit trajectories and parameter tables. Its main weaknesses are numerical: the hypergeometric evaluation in Eq. (3) is undocumented, and the branch-splitting evidence rests on energy differences as small as a few times 10^-5 without error estimates or convergence criteria.","major_comments":[{"comment":"The central non-monotonicity of f(r) relies on numerical evaluation of 2F1(1/4,1;5/4; -4p^2 α2/(r^4 α1)) for large negative arguments. For the representative parameters, at r≈0.2 the argument is approximately -150 and at r≈0.15 it is approximately -480, far outside the |z|<1 radius of convergence of the defining series; the chosen sheet of the analytic continuation is therefore crucial. The paper states neither the evaluation algorithm nor any validation against an independent method, so the double-barrier structure and all subsequent claims rest on an undocumented numerical step. Please specify the hypergeometric evaluation method, provide convergence or cross-check tests (for example, against a power series in a different variable or known limits), and quantify the resulting uncertainty in f and its derivatives near the inner barrier.","section":"Section II A, Eq. (3)"},{"comment":"The multi-branch claim is supported by energy and angular momentum roots that are very close together; for example, in Table II for (1,2,0) at ϵ=0.1, E1=0.917706569 and E2=0.918695150 differ by about 10^-3, while for (2,3,1) the E1 and E2 values differ by about 3×10^-5. The paper gives no root-finding tolerances, no quadrature error estimates for Eq. (15), and no convergence criteria for the numerical integration of the geodesic equations. Please add error bars or tolerance-based statements demonstrating that these roots are distinct beyond numerical error; otherwise the branch splitting could be an artifact of the numerical scheme.","section":"Section IV A, Tables I-VI"},{"comment":"The analysis finds three MBO branches but excludes the first branch (LMBO1<LISCO) with the explanation that the corresponding potential shapes are 'not associated with normal closed bound orbits.' Since the abstract and conclusions emphasize 'multiple marginally bound orbits,' this exclusion needs a quantitative demonstration (for example, a scan for periodic solutions in that angular momentum range) rather than a post hoc statement. If no periodic orbits exist there, that should be shown; if they do exist, they should be included in the analysis.","section":"Section IV A, Fig. 3 and concluding summary"},{"comment":"Because Veff→1 from below at spatial infinity, an orbit with E>1 is not bound in the usual asymptotic sense; it is confined by the outer potential barrier. The paper should state this distinction explicitly and demonstrate for the plotted E>1 trajectories that the radial turning points and the barrier prevent escape. Currently the phrase 'bound orbits with E>1' is used without this caveat, and the figures alone do not show that the orbit cannot eventually cross the outer barrier over longer integration times.","section":"Section IV A, E>1 bound orbits"}],"minor_comments":[{"comment":"The stated condition 1≤v≤z−1 excludes the single-leaf cases (1,2,0) and (1,3,0) used throughout; the later text correctly uses v∈[0,z−1], so Eq. (17) should be corrected.","section":"Eq. (17)"},{"comment":"The value of α1 is inconsistent: Tables I and II use 1.025775, while the text and Table III use 1.025755, and Fig. 8(a) labels the single-horizon case as 1.012577 rather than 1.025755.","section":"Tables I, II, and Fig. 8"},{"comment":"There are numerous typographical errors, including 'V ertex number', 'emains identical', 'MBOLMBO2', 'LMB02', and 'rNBO3' in Fig. 3; these should be corrected before publication.","section":"Throughout"},{"comment":"The inner local minimum of the effective potential is described as 'a marginally stable circular orbit,' but a local minimum is a stable circular orbit; marginal stability corresponds to an inflection point where ∂rrVeff=0. The terminology should be fixed.","section":"Section IV A, Fig. 5 discussion"},{"comment":"Several figure cells contain 'No picture' placeholders instead of trajectories (for example, Fig. 9, (2,3,1) at ϵ=0.3); these should either be filled with the corresponding orbit or removed, since a placeholder is not an informative scientific result.","section":"Figures 9 and 15"},{"comment":"The caption states 'Lα1=1.8', which is presumably a typo for 'Lα1=1.08'; in addition, the two-horizon cases α1=1.5 and α1=10 in Table VIII are described as 'four-horizon black holes', inconsistent with the text of Section IV C.","section":"Table VII caption"}],"recommendation":"major_revision","confidential_remarks":"The manuscript fits the scope of the journal and the central physical idea is interesting. My main concern is numerical reproducibility: the hypergeometric evaluation and the very close energy roots underlying the multi-branch claim need explicit validation and error estimates. I do not see a circularity problem: the metric comes from Ref. [39] and the earlier works are used for context. If the authors add the requested numerical documentation, the paper would be publishable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The interesting thing here is real: for a slice of the dyonic black hole parameter space where f(r) is non-monotonic outside the horizon, the effective potential develops a double-barrier, and the paper shows several periodic orbits sharing the same rational label q, plus bound orbits with E > 1. I don't know another background where that has been shown. The paper does this with the standard Levin framework, no fitted parameters, and it carefully contrasts the non-monotonic regime with the monotonic single-branch cases. The table-and-trajectory evidence is substantial, and the contrast with the 3-/4-horizon monotonic cases is a good check on the logic.\n\nThe soft spot is exactly where the reader put it. The metric function is a hypergeometric 2F1 whose argument is strongly negative near the inner barrier, far outside the radius of convergence of the defining series. The paper states no analytic continuation, no evaluation scheme, no tolerances. Meanwhile the branch splitting is numerically delicate: in Table II, E1 and E2 for q=(1,2,0) differ by about 10^-3. If the hypergeometric evaluation is off at the 10^-4 level in that region, the double-barrier itself, and with it the multi-branch orbits and E>1 bound orbits, could be artifacts. That is a load-bearing gap, not a cosmetic one. I also agree the paper overclaims generality: 'unique signatures' is supported by one parameter slice, M=1, α2=2.76, p=0.15, q=1.05. The first MBO branch is excluded post hoc because LMBO1 < LISCO, which is a reasonable presentation choice, but it should be stated more plainly as a restriction. The paper is also sloppy in places: typos in table captions, mislabeled columns in Figure 9, and the claim of a 'new' winding number n = z(1+q) is really just the standard azimuthal count. These are minor individually, but they add friction.\n\nI don't see circularity. The metric comes from an externally published solution, the geodesic equations are standard, and the 'predictions' are direct numerical solutions rather than fits. The authors are not using self-citations to prop up the main result. If the numerics hold up, this is a solid subfield contribution.\n\nMy recommendation: send it to a serious referee. The phenomenon is novel enough to warrant the time, but I would make the request explicit — ask for the code or data, an error analysis, and a demonstration that the hypergeometric continuation is correct at the relevant |z| before accepting. I'd want to see that before citing it myself.","headline":"Genuinely new periodic-orbit phenomena in a non-monotonic dyonic black hole, but the central existence claim needs the numerics substantiated before I'd trust it.","tokens_in":49161,"tokens_out":1604,"would_cite":false,"duration_ms":19725,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Dyonic black holes with a non-monotonic metric function admit several closed bound orbits with the same rational label, including bound orbits with energy greater than one.","keywords":["periodic orbits","dyonic black holes","quasi-topological electromagnetism","double-barrier effective potential","marginally bound orbits","high-energy bound orbits","rational orbit labels","zoom-whirl-vertex classification"],"falsifier":"Recompute the two listed energies for $q=(1,2,0)$ at $\\alpha_1=1.025755$ and $L=2.778089089$ using an independent high-precision integrator for Eq. (15); if the two solutions collapse to one, or if a computed trajectory fails to close with winding number $n=z(1+q)=3$, the multi-branch claim is refuted.","tokens_in":48006,"feed_emoji":"🕳️","tokens_out":9284,"duration_ms":82539,"temperature":0.7,"pith_summary":"This paper argues that the shape of the metric function $f(r)$ outside a black hole horizon controls the bound-orbit menu, not just the number of horizons. In dyonic black holes from quasi-topological electromagnetism, once $f(r)$ becomes non-monotonic, the effective radial potential develops twin barriers and wells, so several marginally bound orbits exist and the same rational label $q=w+v/z$ can describe two or even three geometrically different periodic orbits. The paper also constructs closed bound orbits with energy $E>1$ when the outer barrier exceeds one, and finds up to three coexisting branches when the barrier peak reaches $E=1$. If these numerical findings hold, periodic-orbit tables and inspiral waveform searches in such spacetimes must allow several energy-angular-momentum solutions for each rational label.","feed_headline":"Three orbits can share one rational label in dyonic black holes","feed_subtitle":"Non-monotonic gravity creates double-well potentials and can hold particles with energy above one.","key_machinery":"The load-bearing object is the metric function $f(r)=1-2M/r+\\alpha_1 p^2/r^2+(q^2/(\\alpha_1 r^2))\\,{}_2F_1(1/4,1;5/4;-4p^2\\alpha_2/(r^4\\alpha_1))$, whose $\\alpha_1$-dependent non-monotonicity creates the double-barrier effective potential. The analysis then hangs on three conditions: marginally bound orbits at $V_{\\rm eff}=1$ and $\\partial_r V_{\\rm eff}=0$, ISCO at $V_{\\rm eff}=E^2$ with $\\partial_r V_{\\rm eff}=\\partial_{rr}V_{\\rm eff}=0$, and the frequency-ratio integral $\\Delta\\phi_r=2\\int L/(r^2\\sqrt{E^2-V_{\\rm eff}})\\,dr$ whose value, through $q=\\Delta\\phi_r/(2\\pi)-1$, assigns the $(z,w,v)$ rational label and winding number $n=z(1+q)$ to each closed orbit. This machinery converts the question of how many periodic orbits exist for a given $q$ into a root-finding problem for the energy and angular momentum that make the integral take a fixed rational value; a double-well potential makes that problem genuinely multi-valued.","core_discovery":"The central claim is that in the parameter window $1.025755\\le \\alpha_1<1.07317$ of the dyonic black hole family (with $M=1$, $\\alpha_2=2.76$, $p=0.15$, $q=1.05$), the metric function $f(r)$ is non-monotonic outside the event horizon and produces an effective potential $V_{\\rm eff}=f(r)(1+L^2/r^2)$ with two maxima and two minima. Because of that landscape, the equations defining marginally bound orbits have three solutions, and the frequency-ratio integral that assigns each periodic orbit its rational label $q=\\Delta\\phi_r/(2\\pi)-1$ becomes multi-valued in energy and angular momentum. The paper exhibits explicit closed trajectories: for example, at $\\alpha_1=1.025755$ with $\\epsilon=0.1$, the same $q=(1,2,0)$ appears at $E\\approx0.91771$ and $E\\approx0.91870$, and for some labels a third branch appears at higher $E$. All branches have the same winding number $n=z(1+q)$ and are topologically equivalent in the $(z,w,v)$ classification, yet they differ in radial extent and eccentricity; the innermost branch even becomes more circular as $E$ or $L$ increases while outer branches become more eccentric. The same analysis shows that bound orbits with $E>1$ exist whenever the outer barrier exceeds unity, a feature absent for monotonic $f(r)$.","pith_inferences":["Any static spherically symmetric metric whose $f(r)$ develops two local minima outside the horizon should show the same double-barrier orbit structure; testing a different theory, such as a compact object with a nested shell, would indicate whether the mechanism is generic.","The paper fixes $\\alpha_2=2.76$, $p=0.15$, and $q=1.05$; mapping the second and third branch boundaries across those parameters would show how robust the branches are and where they pinch off.","If real, $E>1$ bound orbits would be an unusual gravitational-wave source: particles could orbit in a well behind an outer barrier without escaping, and the paper does not yet compute such waveforms.","The reported inversion, where inner branches become more circular as energy rises while outer branches grow more eccentric, could serve as a geometric fingerprint for non-monotonic metric functions, provided an independent calculation reproduces it."],"forward_implications":["In the non-monotonic window, the standard one-periodic-orbit-per-rational-label tables are incomplete: the same $q$ can correspond to two or three closed orbits, so searches that invert $q$ to $(E,L)$ must allow multiple solutions.","Bound motion is not limited to $E<1$ in these spacetimes; particles with energy above the asymptotic value can still be trapped by the outer barrier, which changes how 'bound' should be defined near such objects.","All branches with the same $q$ share the same winding number and $(z,w,v)$ classification, so interpretations that rely only on orbit topology will not distinguish them; the radial profile matters.","Monotonic $f(r)$ cases, including the three- and four-horizon regimes, reduce to the known single-branch behavior, so the multi-branch phenomenon is tied to the shape of $f(r)$, not to the number of horizons alone."],"supporting_citations":[{"why":"Supplies the dyonic black hole solution and the hypergeometric metric function $f(r)$ whose non-monotonic behavior drives the whole analysis.","marker":"[39]"},{"why":"Supplies the rational-number labeling of periodic orbits, the $(z,w,v)$ zoom-whirl-vertex classification, and the winding-number definition used to establish topological equivalence.","marker":"[17]"},{"why":"Supplies the topological pairing of stable and unstable timelike circular orbits that underpins the double-well effective potential interpretation.","marker":"[46]"},{"why":"Documents multiple photon spheres in these spacetimes, providing prior evidence that non-monotonic $f(r)$ creates multiple barriers.","marker":"[40]"}],"fun_headline_variants":["Three distinct orbits share one rational label in dyonic holes","Double-well potential gives three orbit families per rational number","Non-monotonic metric reopens orbit branches for dyonic black holes","Super-energy bound orbits appear when potential barrier peaks","Same winding number, different radial paths in dyonic holes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument assumes the numerical evaluation of the metric function and of the orbit conditions is accurate enough to resolve energy differences as small as about one part in a thousand; if it is not, the claimed separate branches could merge or disappear.","fun_headline_variants_meta":{"raw":{"variants":["Three distinct orbits share one rational label in dyonic holes","Double-well potential gives three orbit families per rational number","Non-monotonic metric reopens orbit branches for dyonic black holes","Super-energy bound orbits appear when potential barrier peaks","Same winding number, different radial paths in dyonic holes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001176,"raw_usage":{"total_tokens":4976,"prompt_tokens":1177,"completion_tokens":3799,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":793,"completion_tokens_details":{"reasoning_tokens":3718}},"tokens_in":793,"tokens_out":3799,"duration_ms":23992,"temperature":1.0,"reasoning_tokens":3718,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T16:45:42.564541+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the two listed energies for $q=(1,2,0)$ at $\\alpha_1=1.025755$ and $L=2.778089089$ using an independent high-precision integrator for Eq. (15); if the two solutions collapse to one, or if a computed trajectory fails to close with winding number $n=z(1+q)=3$, the multi-branch claim is refuted.","supporting_citations":[],"review_version":2}