{"id":"a326a98d-bb72-479f-8f8f-521123fa2f5c","arxiv_id":"2412.06122","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Chaos grows with energy, magnetic field, and the MOG parameter, and shrinks with spin and angular momentum, for a charged particle in a magnetized Kerr-MOG black hole.","lead":"This paper builds explicit symplectic integrators for a charged particle moving around a Kerr-MOG black hole with a magnetic field, then maps how five parameters shift the motion between regular and chaotic. It matters because it offers a numeric tool and phase diagrams for a modified-gravity black hole, extending an established integrator program to the MOG spacetime.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Chaos-trend claims are not robust: FLI thresholds vary per figure, no step-size convergence is shown, and the two-parameter scan contains an L sign-flip that undermines the 'a and L dominate' conclusion.","rationale":"The explicit splitting construction follows the established time-transformation recipe, and the per-piece flows in Eq. (48) are internally consistent with the decomposed Hamiltonians in Eq. (45); I do not see a fatal derivation error. The comparison of S2, S4, and PRK64 is plausible as a numerical demonstration, but it is only one parameter set and does not support the physical conclusions by itself. The load-bearing part of the paper is the phase-space classification. That classification is fragile because FLI thresholds are chosen per figure with no universal criterion, only one step size (h=1) is used with no convergence check, and the two-parameter dominance claims are qualitative visual assessments rather than a quantitative sensitivity analysis. The reported L sign-flip in Fig. 7(b) is particularly damaging because it directly contradicts the single-parameter monotonic claim and the abstract's assertion that increasing L contracts chaos. These issues do not prove the trends are wrong, but they mean the central claims are not established at the claimed level of generality. The reader's conditional verdict already captures this risk; my additional observation about the two-parameter scan strengthens the same concern without changing the verdict.","tokens_in":20718,"tokens_out":5985,"duration_ms":62406,"concrete_test":"Recompute the two-parameter FLI maps of Figs. 6-8 with a fixed, pre-registered threshold (e.g., FLI=10) and with at least two smaller step sizes, h=0.5 and h=0.1, for w=10^7, and also for w=10^8 on a 20x20 subsample of the (beta,L) and (alpha,L) panels. Additionally, isolate L's effect at r=110 by varying L alone at the fixed beta values used in Fig. 7(b) (beta=3e-4 and 7e-4) while holding E, a, and alpha constant; if the fraction of chaotic orbits rises with L in either slice, the monotone 'L shrinks chaos' conclusion and the 'L dominates' claim fail in that regime.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central physical conclusions are classification statements, not integrator statements: the abstract claims that the chaotic region grows with E, beta, and alpha and shrinks with a and L, and that a and L dominate in two-parameter scans. These conclusions rest on FLI maps binarized with thresholds that change per figure (8, 10, 12, 20, 30, 50) with no calibration to reference orbits and no fixed rule. A different threshold would shift the red/cyan boundaries; for example, the E-dependence at r=11 already shows alternating chaotic and ordered windows between E=0.9933 and E=0.9966, so the claimed monotonicity is not uniform. The simulations use a single step size h=1 in virtual time for w=10^7, with no convergence study reported; because FLI growth rates are numerical quantities, a smaller step size could reclassify marginal orbits and change the claimed trends. The two-parameter analysis is purely visual, with all third parameters fixed at arbitrary values, and no quantitative measure is used to compare parameter influence. This matters concretely: in Fig. 7(b), increasing L at r=110 is reported to promote chaos, contrary to the single-parameter conclusion in Fig. 5(b); the explanation that beta dominates is asserted, not derived. Thus the concluding claim that a and L play the major role is not supported by the evidence shown.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper constructs explicit symplectic integrators for a charged test particle around a magnetized Kerr-MOG black hole. The authors apply a time transformation to the Hamiltonian, split it into five integrable components, and build three explicit integrators (S2, S4, and PRK64), reporting that PRK64 has the smallest Hamiltonian error. They then use Poincar\\'e sections and Fast Lyapunov Indicator (FLI) maps to classify orbits as regular or chaotic for various parameter values. The central physical claims are that the chaotic region expands with increasing energy E, magnetic field parameter β, and MOG parameter α, while it contracts with increasing angular momentum L and spin a, and that in simultaneous two-parameter scans a and L play the dominant role.","tokens_in":20968,"tokens_out":7594,"duration_ms":72956,"significance":"If the numerical evidence is sound, the paper would provide a useful explicit symplectic integrator for this modified-gravity spacetime and a first systematic phase-diagram study of its chaotic dynamics. The Hamiltonian splitting in Eqs. (43)-(48) is a nontrivial construction, and the explicit analytic solutions for the five sub-Hamiltonians are a genuine technical contribution. The energy-error comparison of S2, S4, and PRK64 at fixed step size is a useful benchmark. However, the paper does not ship code or machine-checked proofs, and its main qualitative conclusions rest on FLI threshold choices that are not justified or calibrated; the numerical robustness of the chaos trends is therefore not yet established.","major_comments":[{"comment":"The FLI thresholds are chosen separately for each figure (values 8, 10, 12, 20, 30, and 50) with no calibration to reference orbits and no fixed, universal criterion. Since the labels 'regular' and 'chaotic' are obtained by thresholding the FLI maps, and the paper's central conclusions are qualitative statements about the size and location of chaotic regions, a different threshold could reclassify marginal orbits and change the reported trends. Please provide a calibration procedure using reference regular and chaotic orbits, and add a robustness check showing that the conclusions are stable over a range of thresholds.","section":"§IV.A, §IV.B, Figs. 3 and 6–8"},{"comment":"All simulations use a single step size h=1 and an integration time of w=10^7, and no step-size convergence study or comparison with an independent integrator is reported. The FLI is a numerical quantity, so a smaller step size or a different integration method could reclassify weakly chaotic or sticky orbits. Please add convergence tests (for example h=0.5, 0.25 or energy-error scaling) and, where possible, an independent adaptive-integrator baseline for a representative subset of orbits.","section":"§III.B and §IV"},{"comment":"The conclusion that 'a and L play a major role' in the simultaneous variation case is not supported by the evidence shown. Fig. 7(b) reports that increasing L promotes chaos at r=110, which is opposite to the single-parameter conclusion of Fig. 5(b); the explanation that β dominates is asserted rather than derived. In addition, the two-parameter scans are compared only by visual inspection, with no quantitative measure of parameter influence, and the third parameters are fixed at nominal values without sensitivity analysis. Please quantify the sensitivity (for example, the fraction of chaotic grid points as a function of each parameter) and reconcile the L behavior before drawing the parameter-dominance conclusion.","section":"§IV.B and §V"}],"minor_comments":[{"comment":"With the definition FLI = log10(d(w)/d(0)), a regular orbit has d(w) growing at most linearly and hence FLI growing logarithmically with time, while a chaotic orbit has d(w) growing exponentially and FLI growing linearly. The text says the opposite; please correct this description.","section":"§IV, Fig. 1(d)"},{"comment":"The statement that symplectic algorithms 'rigorously preserve energy' is an overstatement; symplectic integrators preserve the symplectic form and have bounded energy error, not exact energy conservation. Please rephrase.","section":"Introduction and §III.A"},{"comment":"The thresholds quoted in the text are inconsistent: the values E≤0.993 and E≥0.9969 do not match the stated transition values E=0.9933 and E=0.9966. Please make these values consistent.","section":"§IV.A, Fig. 3(a)"},{"comment":"The notation 'P RK64' appears with a space; use 'PRK64' consistently.","section":"Throughout"},{"comment":"Given the numerical nature of the claims, depositing the code used to produce the FLI maps and Poincar\\'e sections would improve reproducibility; the current statement only promises availability on reasonable request.","section":"Data and Code Availability"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the integrable splitting appears to be a reasonable technical contribution. The main risk is numerical robustness: the qualitative chaos trends and the parameter-dominance claim are built on thresholded FLI maps without calibration or convergence studies. I see no evidence of misconduct or citation problems; the revision should focus on quantitative support for the physical conclusions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Punchline: useful integrator, shaky chaos claims. The genuinely new contribution is the five-term decomposition of the magnetized Kerr-MOG Hamiltonian and the explicit analytic solutions for each term, following Wu et al.'s time-transformation framework. The PRK64 integrator is shown to preserve energy better than S2 and S4 over 10^7 steps, which is a solid numerical result. The single-parameter scans are internally consistent: Poincaré sections and FLI both show chaos growing with E, β, α and shrinking with a, L. So the qualitative picture is probably right.\n\nThe soft spots are where the stress-test note lands. The FLI thresholds are chosen per figure (8, 10, 12, 20, 30, 50) with no calibration to reference orbits; that makes the binarization of the two-parameter maps arbitrary boundary placement. There is no step-size convergence test; h=1 is used throughout, and for chaotic systems a smaller step can change classifications of marginal orbits. The text's description of FLI growth in Fig. 1(d) is reversed: regular orbits give logarithmic growth of FLI, chaotic orbits give linear growth. And the two-parameter analysis is purely visual, with no quantitative measure of chaotic fraction, so the claim that 'a and L play a major role' is not supported by the evidence: in the (β,L) panel of Fig. 7, increasing L promotes chaos and β dominates, contradicting the single-parameter L trend. The authors do note this exception, but the conclusion still overstates.\n\nWho benefits: researchers working on explicit symplectic integrators for black hole spacetimes and on chaos in MOG/modified-gravity backgrounds. The integrator part is citable and correct as far as I can tell. The chaos classification needs more work before the qualitative trends are treated as robust. I would not desk reject this; it deserves a serious referee. The referee should ask for a fixed FLI threshold or a calibration procedure, a step-size convergence test for at least one representative orbit, and a more careful wording of the two-parameter conclusions.","headline":"Useful incremental symplectic integrator for Kerr-MOG; chaos-trend claims rest on per-figure FLI thresholds and a single step size, and the 'a and L dominate' summary overreaches.","tokens_in":21530,"tokens_out":4832,"would_cite":false,"duration_ms":46679,"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":"This paper claims that in a magnetized Kerr-MOG black hole, charged-particle orbits become more chaotic as energy, magnetic field, or MOG parameter rise, while spin and angular momentum suppress chaos, and it maps this with new explicit…","keywords":["chaotic dynamics","Kerr-MOG black hole","explicit symplectic integrator","Hamiltonian splitting","Fast Lyapunov Indicator","charged particle motion","magnetic field","Poincaré sections"],"falsifier":"Recompute the FLI maps for the $(\\beta,E)$ and $(\\alpha,a)$ planes with step sizes $h=0.1$ and $h=0.5$, and label orbits with a threshold-free measure such as the 0-1 test or by comparing Lyapunov exponents over two different time windows; if any claimed transition, such as chaos onset near $E\\approx 0.9933$ at $r=11$ or the ordered region at $a>0.6$, moves or disappears, the qualitative claims are not robust.","tokens_in":20501,"feed_emoji":"🕳️","tokens_out":4510,"duration_ms":41627,"temperature":0.7,"pith_summary":"The paper aims to show that the long-term dynamics of a charged test particle around a Kerr-MOG black hole in a uniform magnetic field can be studied reliably with explicit symplectic integrators, and that the resulting regular or chaotic character of orbits has a systematic dependence on five parameters. It constructs three integrators by splitting a time-transformed Hamiltonian into five integrable pieces and reports that the PRK64 version keeps the Hamiltonian error smallest over $10^{7}$ steps. Using Poincaré sections and the Fast Lyapunov Indicator, it claims that the chaotic region expands as the particle energy $E$, the magnetic field strength $\\beta$, or the MOG parameter $\\alpha$ increases, and contracts as the black hole spin $a$ or the angular momentum $L$ increases. When two parameters act together, the paper concludes that $a$ and $L$ dominate the response while the MOG parameter has relatively little influence. A sympathetic reader would care because this supplies a numerical tool and a first phase diagram for chaos in a modified-gravity black hole that is otherwise nonintegrable.","feed_headline":"Chaos grows with energy and magnetic field near Kerr-MOG black holes","feed_subtitle":"Explicit symplectic integrators map how five parameters push charged orbits from regular to chaotic.","key_machinery":"The central object is the five-term splitting of the time-transformed Hamiltonian $H = H_1+H_2+H_3+H_4+H_5$, where $H_1 = (\\Sigma/r^2)(H_p+p_0)$, $H_2 = [a^2+\\alpha(1+\\alpha)]/(2r^2)\\,p_r^2$, $H_3 = p_r^2/2$, $H_4 = -(1+\\alpha)/r\\,p_r^2$, and $H_5 = p_\\theta^2/(2r^2)$. Each subsystem has an explicit analytic solution in the new time $w$, so composing their flows with symmetric coefficients yields the second-order S2, fourth-order S4, and optimized PRK64 integrators. The splitting is what converts a nonintegrable Hamiltonian into a sequence of exactly solvable kicks, and the PRK64 composition is what keeps the energy error small over long integrations.","core_discovery":"The central claim is that the magnetized Kerr-MOG Hamiltonian, after the time transformation $d\\tau = (\\Sigma/r^2)\\,dw$, splits into five analytically solvable parts, so explicit symplectic integrators S2, S4, and PRK64 can be constructed; PRK64 preserves the Hamiltonian to roughly two orders of magnitude better than S4 over $10^7$ steps. With that integrator, the paper states that FLI and Poincaré sections show order-to-chaos transitions driven by $E$, $\\beta$, and $\\alpha$, and chaos-to-order transitions driven by $a$ and $L$. The main conclusions are qualitative: the chaotic area increases as $E$, $\\beta$, or $\\alpha$ increases, but $a$ and $L$ act in the opposite direction, and in two-parameter scans $a$ and $L$ play the major role.","pith_inferences":["Editorial: applying a threshold-independent chaos indicator such as the 0-1 test or finite-time Lyapunov exponents to the same orbits could test whether the monotonic trends survive without the ad hoc FLI cutoffs of 8, 10, 12, 20, 30, and 50.","Editorial: if the trends hold, one could try to connect the phase diagram to astrophysical observables such as variability of accretion flows or hotspots, where spin would dominate the chaotic signature and $\\alpha$ would be hard to constrain.","Editorial: the same five-term splitting with PRK64 coefficients could be reused for other modified-gravity or external-field backgrounds, provided their Hamiltonians reduce to the same split structure.","Editorial: a step-size convergence study with $h<1$ would clarify whether the reported thresholds, such as $E\\approx 0.9933$ at $r=11$, are physical boundaries or numerical artifacts of the single step size $h=1$."],"forward_implications":["Long integrations of charged-particle motion around Kerr-MOG black holes can be performed explicitly and symplectically, avoiding the pseudo-chaos that energy drift causes in ordinary integrators.","The claimed phase diagram predicts where regular and chaotic orbits sit as functions of $E$, $L$, $\\beta$, $a$, and $\\alpha$, so later studies can select parameter regions that isolate a single dynamical mechanism.","Two-parameter scans imply that spin and angular momentum, rather than the MOG parameter, control whether simultaneous parameter changes suppress chaos.","Because $\\alpha$ affects chaos only mildly, orbital chaos is likely to be a weak probe of MOG compared with changes in spin and magnetic field strength.","The same splitting-and-composition strategy should carry over to other stationary axisymmetric spacetimes whose Hamiltonians admit a similar time transformation."],"supporting_citations":[{"why":"Supplies the Kerr-MOG black hole solution that the whole study is built on.","marker":"[21]"},{"why":"Introduces the splitting of a black-hole Hamiltonian into integrable sub-Hamiltonians, the method adapted here.","marker":"[44]"},{"why":"Provides the time-transformation technique that makes the Kerr-MOG Hamiltonian splittable into five integrable parts.","marker":"[47]"},{"why":"Gives the charged-particle Hamiltonian and parameterization for Kerr-MOG in an external magnetic field that the paper starts from.","marker":"[55]"},{"why":"Provides Wald's electromagnetic four-potential used to model the asymptotically uniform magnetic field.","marker":"[56]"},{"why":"Supplies Yoshida's symmetric composition method used to build the fourth-order integrator S4.","marker":"[59]"},{"why":"Provides the optimized PRK64 coefficients that give the most accurate integrator in the paper.","marker":"[60]"},{"why":"Introduces the Fast Lyapunov Indicator as the chaos-detection method used to classify orbits.","marker":"[62]"},{"why":"Supplies the two-nearby-trajectory computation of Lyapunov indicators in curved spacetime that underlies the FLI values.","marker":"[63]"}],"fun_headline_variants":["Chaos grows with energy and B-field, shrinks with spin and L","Energy and β drive chaos; spin and L tame it in Kerr-MOG","Five parameters control chaos around Kerr-MOG black holes","Symplectic integrators map chaos drivers near Kerr-MOG","Charged particle chaos: energy and field increase, spin and L decrease"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The conclusions rest on classifying each orbit by comparing its Fast Lyapunov Indicator to a threshold chosen separately for each figure (8, 10, 12, 20, 30, or 50) after a single run with step size $h=1$ and integration time $w=10^7$; if those thresholds or that step size mislabel even some orbits, the claimed monotonic trends in chaos could shift.","fun_headline_variants_meta":{"raw":{"variants":["Chaos grows with energy and B-field, shrinks with spin and L","Energy and β drive chaos; spin and L tame it in Kerr-MOG","Five parameters control chaos around Kerr-MOG black holes","Symplectic integrators map chaos drivers near Kerr-MOG","Charged particle chaos: energy and field increase, spin and L decrease"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000904,"raw_usage":{"total_tokens":3927,"prompt_tokens":1022,"completion_tokens":2905,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":638,"completion_tokens_details":{"reasoning_tokens":2812}},"tokens_in":638,"tokens_out":2905,"duration_ms":19949,"temperature":1.0,"reasoning_tokens":2812,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T19:59:59.728089+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the FLI maps for the $(\\beta,E)$ and $(\\alpha,a)$ planes with step sizes $h=0.1$ and $h=0.5$, and label orbits with a threshold-free measure such as the 0-1 test or by comparing Lyapunov exponents over two different time windows; if any claimed transition, such as chaos onset near $E\\approx 0.9933$ at $r=11$ or the ordered region at $a>0.6$, moves or disappears, the qualitative claims are not robust.","supporting_citations":[{"cited_title":"Induction of chaotic fluctuations in particle dynamics in a uniformly accelerated frame","cited_arxiv_id":"1904.11760","evidence_quote":"Supplies the Kerr-MOG black hole solution that the whole study is built on."},{"cited_title":"Construction of explicit symplectic integrators in general relativity. III. Reissner-Nordstrom-(anti)-de Sitter black holes","cited_arxiv_id":"2103.12272","evidence_quote":"Provides the time-transformation technique that makes the Kerr-MOG Hamiltonian splittable into five integrable parts."},{"cited_title":"Charged particle motion and acceleration around Kerr-MOG black hole","cited_arxiv_id":"2311.16936","evidence_quote":"Provides Wald's electromagnetic four-potential used to model the asymptotically uniform magnetic field."},{"cited_title":"On the Structure of Symplectic Mappings. The Fast Lya- punov Indicator: a Very Sensitive Tool,","cited_arxiv_id":null,"evidence_quote":"Supplies the two-nearby-trajectory computation of Lyapunov indicators in curved spacetime that underlies the FLI values."}],"review_version":1}