{"id":"47d09a7b-9e7f-48a9-b7d5-2c3f2a866c15","arxiv_id":"1908.01886","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Timelike geodesics in the rotating Hayward regular black hole admit a Carter-like fourth constant of motion, and out-of-plane orbits differ from Kerr mainly at low spin and large length parameter g.","lead":"This paper studies the paths of massive particles around a rotating version of the Hayward regular black hole and finds a fourth conserved quantity, analogous to the Carter constant in Kerr spacetime. It numerically compares orbits with Kerr and reports the largest differences for slowly rotating holes with a large length parameter.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Hamilton equation (4.7) has the wrong sign for the E L_z cross term; the out-of-plane numerical orbits in Figs. 15-16 are therefore not reproducible from the printed equations.","rationale":"The reader's weakest_assumption points to the physical validity of the Bambi-Modesto metric and the unproved separation ansatz. I checked the separation: substituting the inverse metric (2.5) into the Hamilton-Jacobi equation with the ansatz (3.4) separates exactly because M(r) appears only in r-dependent coefficients proportional to 1/Δ, while all θ-dependence is independent of M(r). Thus the existence of the Carter-like constant Q is not a weak point. The physical-source caveat (Refs. [16-18] and Sec. I) is acknowledged by the authors and would affect the interpretation of g as a physical parameter, but not the mathematical geodesic structure of metric (2.3). The most concrete load-bearing concern is instead the sign error in Eq. (4.7): it is internal, it affects the numerical out-of-plane results that are presented as new, and it cannot be checked against code because none is shipped. This supports the reader's CONDITIONAL verdict, so I recommend no change to the verdict, but for a more specific reason than the one the reader emphasized.","tokens_in":13602,"tokens_out":20122,"duration_ms":191844,"concrete_test":"Recompute ∂H/∂θ with p_t = -E and p_φ = L_z; the coefficient of E L_z should be -2g^{tφ}_{,θ}, not +2g^{tφ}_{,θ}. With the corrected sign in Eq. (4.7), integrate the first-order system (4.4), (4.7)-(4.11) using the initial data stated for Figs. 15 (a = 0.9, v0 = 0.96, r0 = 5, Q = -0.2) and 16 (a = 0.6, v0 = 0.52, r0 = 6, Q > 0) for each displayed value of g. Compare the θ-range and orbit shapes with the published figures, and check conservation of H and p_φ to machine precision. If the corrected trajectories differ qualitatively from those in the paper, the numerical out-of-plane comparison is unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing weakness is an internal sign inconsistency in the Hamilton system used for the numerical out-of-plane trajectories. From H = 1/2(g^{αβ}p_αp_β) with p_t = -E and p_φ = L_z, the cross term is -2g^{tφ} E L_z, so ∂H/∂θ contains -2g^{tφ}_{,θ} E L_z. Equation (4.7) instead writes +2g^{tφ}_{,θ} E L_z inside the bracket, and all g's in Sec. IV are inverse metric components. Since the θ-equation is what makes the motion leave the equatorial plane, the trajectories shown in Figs. 15 and 16, and the claim that departures from Kerr are most visible for slowly rotating holes with large g, are not reproducible from the manuscript as written. This is more specific than the physical-source caveat: the Hamilton-Jacobi separation itself is sound, so the fourth constant is not in doubt; the unsupported part is the numerical orbit comparison. No code or full initial data are shipped, so the printed sign is the only specification available to a reader.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies timelike geodesics in the rotating Hayward regular black hole of Bambi and Modesto. It analyzes the horizon and ergosphere structure, claims the existence of a fourth Carter-like constant of motion via Hamilton-Jacobi separation, classifies radial and polar motion using effective potentials, and presents numerical trajectories in and out of the equatorial plane, comparing them with the Kerr spacetime. The central assertion is that departures from Kerr are most pronounced for slowly rotating holes with large values of the characteristic length parameter g.","tokens_in":13818,"tokens_out":31936,"duration_ms":268438,"significance":"If the derivation is correct, the paper provides a useful extension of Carter's integrability to a family of regular rotating metrics and one of the first systematic comparisons of out-of-plane geodesics in such spacetimes with Kerr, a topic relevant for gravitational-wave dephasing estimates. The Hamilton-Jacobi machinery is standard, the effective-potential analysis is pedagogically useful, and no parameters are fitted to data, so the fourth-constant claim is not circular. However, several inconsistencies in the printed equations currently prevent the results from being reproduced as written, and at least one typo appears to affect every subsequent formula.","major_comments":[{"comment":"The mass function is printed as M(r)=m r^2/(r^3+g^3). With this definition the g=0 limit gives M(r)=m/r and Δ=r^2+a^2-2m, not the Kerr Δ=r^2+a^2-2mr, contradicting the statement that the metric reduces to Kerr for g=0; the a=0 limit also differs from the Hayward f(r) in Eq. (2.2). The correct Bambi-Modesto mass function is M(r)=m r^3/(r^3+g^3) (with the usual denominator r^3+2mℓ^2 in Hayward's notation). Because M(r) enters Δ, R, the effective potentials, and all numerical integrations, this is a load-bearing error and every subsequent equation must be re-examined after the correction.","section":"II, Eq. (2.4)"},{"comment":"The text states Q=k-(aE-L_z)^2, but the printed R(r) and Θ(θ) are not the separation results for that relation: with Q=k-(aE-L_z)^2 the angular equation acquires an extra -2aE L_z term and the radial equation differs by -2aE L_z Δ. The pair (3.7)-(3.8) is instead consistent with Q=k-a^2E^2-L_z^2. Since the separation step is not shown, the reader cannot resolve the discrepancy. Please present the separation explicitly and use a consistent definition of the Carter constant.","section":"III.1, Eqs. (3.5)-(3.8)"},{"comment":"In the Hamiltonian equations, p_t=-E and p_φ=L_z, so the cross term in H is -2g^{tφ}E L_z; hence ∂H/∂θ contains -2g^{tφ}_{,θ}E L_z and ˙pθ should have -2g^{tφ}_{,θ}E L_z inside the bracket after the minus sign. Equation (4.7) prints +2g^{tφ}_{,θ}E L_z. Since the θ-momentum equation is the term that makes the motion leave the equatorial plane, the out-of-plane trajectories in Figs. 15 and 16 are not reproducible from the printed system. Correct the sign and regenerate the affected figures, or supply the code and initial data.","section":"IV, Eq. (4.7)"},{"comment":"The quantity written as f'(1) in the Q<0 case is not f'(1); the expression (Q+L_z^2-a^2(E^2-1))^2+4a^2(E^2-1)Q is the discriminant of the quadratic B v^2+A v+Q in v=u^2. Moreover, the stated inequality '≤0' appears to have the wrong sign for the existence of two positive roots (for a downward-opening parabola one needs a nonnegative discriminant). The classification of Q<0 polar motion is therefore not correctly derived.","section":"III.B, Eq. (3.23) and following"}],"minor_comments":[{"comment":"The separation ansatz (3.4) and the passage from Eq. (3.5) to Eqs. (3.6)-(3.8) are skipped; please show the algebra so the claimed fourth constant can be verified directly.","section":"III.1"},{"comment":"The Q=0 case is described in a confusing sentence: 'There is a solution in which θ is constant at the equatorial plane u=0 if L_z<a^2(E^2-1).' This condition is not derived and should be clarified or removed.","section":"III.B"},{"comment":"The caption reads 'Potential function, V+ for a = with g = 0.2 and Lz = −2.5'; the value of a is missing and should be supplied.","section":"Fig. 8 caption"},{"comment":"The parameter is ℓ in Eq. (2.2) and g elsewhere; the paper should state the identification used (for example g^3=2mℓ^2) or use a single symbol consistently.","section":"II, Eq. (2.2) vs (2.4)"},{"comment":"The abstract and Section V claim that departures from Kerr are most visible for slowly rotating holes with large g, but no quantitative measure is given; a quantitative comparison (e.g., precession-angle difference or trajectory separation) would strengthen this claim.","section":"V"},{"comment":"The paper acknowledges in Refs. [16-18] that the matter source is only approximate; this caveat should be stated explicitly where the metric is introduced, since it delimits the physical interpretation of the subsequent geodesic results.","section":"II, after Eq. (2.3)"}],"recommendation":"major_revision","confidential_remarks":"The error in Eq. (2.4) looks like a typographical slip (r^2 instead of r^3) but it is pervasive and must be corrected before the paper can be evaluated. I recommend asking the authors to verify that the numerical code uses the corrected mass function and to state the corrected definition explicitly. Given that no code or complete initial data are provided, the numerical section is currently not independently checkable; requiring code or at least a table of initial data would help. The physical caveat about the matter source [16-18] is not disqualifying for a geodesic-structure paper, but the authors should frame the metric as effective/phenomenological rather than as an exact solution of a known theory."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe genuinely new thing here is the out-of-equatorial-plane geodesic analysis for the Bambi-Modesto rotating Hayward metric: a fourth constant of motion analogous to Carter's, an effective-potential classification for timelike particles, and a numerical catalog of non-equatorial orbits. The Hamilton-Jacobi separation is standard—this metric is Kerr with M(r)—and the resulting R and Θ forms are almost certainly correct. The paper also gives a reasonable treatment of horizons and ergospheres. Credit where due: the analytic core is likely sound and the prose is clear.\n\nNow the soft spots. The stress-test sign error in Eq. (4.7) is real. With H = 1/2 g^{αβ}p_αp_β and p_t = -E, p_φ = L_z, the cross term in H is -g^{tφ}EL_z, so ∂H/∂θ contributes +g^{tφ}_{,θ}EL_z to ˙pθ. The printed equation has the opposite sign, and since this is the θ-equation that moves the orbit off the equator, Figs. 15 and 16—and the claim that departures from Kerr are most visible for slow rotation and large g—are not reproducible from the manuscript as written. No code or full initial data are shipped, so the printed equation is the only specification. This is a load-bearing flaw for the numerical comparison, though not for the existence of the fourth constant.\n\nThere is also a smaller inconsistency in the Q < 0 classification: Eq. (3.23) gives f'(1) = 2(2B + A), while the bullet says f'(1) equals a different expression that is not even dimensionally consistent. Probably a typo, but it needs fixing. The separation step itself is skipped—acceptable in a short paper, but a referee should ask to see it.\n\nOne broader caveat the authors themselves acknowledge: the Bambi-Modesto metric is only an approximate matter solution (their refs. [16–18]). That does not invalidate a geodesic study of the geometry, but it weakens any astrophysical constraint on g.\n\nBottom line: the analytic result is worth taking seriously, the numerical comparison is not currently supported. This deserves a serious referee—it is a legitimate extension with a fixable flaw. I would send it out, asking for corrected equations, a shown separation, and either code or enough initial data to reproduce the orbits. The likely readers are people working on regular black hole shadows, accretion, or EMRI dephasing. I would not cite the numerical comparison as it stands, though I might cite the fourth-constant result once verified.\n\nBest,\n[You]","headline":"The Carter-type constant for rotating Hayward is probably solid, but the printed sign in the θ-evolution equation undercuts the out-of-plane numerical comparison.","tokens_in":14313,"tokens_out":3526,"would_cite":false,"duration_ms":35227,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["04.70.Dy","04.62.+v","11.10.Kk"],"model":"deepseek-v4-flash","headline":"The geodesic motion around the rotating Hayward regular black hole is integrable through a Carter-like fourth constant, and its largest departures from Kerr occur for slowly rotating holes with a large regularization length.","keywords":["rotating regular black holes","Hayward metric","Carter constant","geodesic integrability","effective potential","timelike orbits","Kerr comparison","gravitational-wave emission"],"falsifier":"Take a generic timelike geodesic off the equatorial plane in metric (2.3), integrate the full second-order equations, and evaluate the separated expression $k=Q+(aE-L_z)^2$ along the orbit; if it drifts by more than numerical error, the separation ansatz is false. A complementary check is to insert the metric into the field equations with an exact nonlinear-electrodynamics stress-energy tensor and see whether $M(r)=m r^3/(r^3+g^3)$ remains a solution without extra terms.","tokens_in":13401,"feed_emoji":"🕳️","tokens_out":7925,"duration_ms":73325,"temperature":0.7,"pith_summary":"This paper establishes that timelike geodesics in the rotating Hayward regular black hole spacetime are completely integrable: in addition to energy, axial angular momentum, and rest-mass normalization, there is a fourth constant of motion analogous to the Kerr Carter constant. The authors derive the separated radial and polar equations from a Hamilton-Jacobi ansatz, classify the possible polar motion by the sign of this constant, and simulate orbits both in and out of the equatorial plane. The payoff is a concrete way to compare this regular black hole with Kerr: differences in trajectories appear mainly for slowly rotating holes with large values of the regularization length $g$, where particles that would fall into a Kerr black hole can instead be trapped in bound orbits. If correct, this gives observational routes, such as accretion dynamics and gravitational-wave emission, to constrain $g$.","feed_headline":"Rotating Hayward black holes inherit Kerr's fourth constant","feed_subtitle":"Off-equatorial orbits stay integrable, and departures from Kerr grow at low spin and large regularization length.","key_machinery":"The central object is the Hamilton-Jacobi separation ansatz $S=\\tfrac12\\tau-Et+L_z\\varphi+S_r(r)+S_\\theta(\\theta)$ applied to the rotating Hayward metric (2.3). This ansatz produces the Carter-like constant $Q$ and reduces the geodesic problem to first-order quadratures in $r$ and $\\theta$. Two effective-potential tools carry the analysis: the radial potential $V_\\pm$ obtained by squaring the radial equation, and the polar quartic $f(u)=Q+Au^2+Bu^4$ in $u=\\cos\\theta$, whose roots decide whether an orbit crosses the equatorial plane. For numerical trajectories the paper switches to the Hamiltonian form of the geodesic equations to avoid the sign ambiguities of square roots at turning points.","core_discovery":"For the rotating Hayward metric proposed by Bambi and Modesto, which has the Kerr form with a radial mass function $M(r)=m r^3/(r^3+g^3)$, the Hamilton-Jacobi equation separates in Boyer-Lindquist-like coordinates. The paper claims there exists a fourth integral $Q$, defined through the separation constant by $Q=k-(aE-L_z)^2$, that governs out-of-equatorial-plane motion just as Carter's constant does for Kerr. The radial and polar motions are governed by the squared potentials $R(r)$ and $\\Theta(\\theta)$, and the polar equation reduces to a quartic in $u=\\cos\\theta$ whose turning-point structure classifies the orbits according to the sign of $Q$. Numerical integration of the first-order equations shows that, compared with Kerr, the largest trajectory differences occur for slowly spinning black holes and large $g$; the precession rate increases with $g$, and bound-versus-plunge outcomes can switch purely by changing $g$.","pith_inferences":["A direct test the paper leaves implicit is to integrate the full second-order geodesic equations for a generic off-equatorial orbit and verify numerically that the separated combination $k=Q+(aE-L_z)^2$ remains constant; this would confirm the fourth constant without relying on the separation ansatz.","If the Bambi-Modesto metric is treated as an effective spacetime, the same $Q$-based machinery applies to black-hole shadow calculations, where off-equatorial photon orbits would probe $g$ through the shadow's shape and size.","Because the paper notes the matter source is only an approximation to the Hayward magnetic-monopole solution, the integrability result may not survive in an exact rotating solution of nonlinear electrodynamics; the fourth constant could be specific to the chosen mass function.","The bound-versus-plunge switching induced by $g$ suggests that accretion-disk inner edges and quasi-periodic oscillations in low-spin systems could encode $g$, a testable extension the paper does not pursue."],"forward_implications":["Geodesic motion in this rotating regular black hole spacetime is completely integrable, so every timelike orbit can be labeled by the four constants $(E,L_z,Q,\\mu)$ in the same way as in Kerr.","Out-of-equatorial-plane orbits are classified by the sign of $Q$: $Q>0$ oscillates through the equatorial plane, $Q<0$ stays on one side, and $Q=0$ confines motion to the equator or to a constant polar angle.","For a fixed spin, increasing $g$ shifts the horizons and reshapes the effective potential so that a particle with fixed energy and angular momentum can switch from plunging into the hole to following a bound precessing orbit.","Measurable differences from Kerr, such as faster precession, shifted polar turning angles, and altered bound-plunge boundaries, are largest for slowly rotating black holes with large $g$, making $g$ the natural target for observational constraints.","Because these trajectory differences would modulate gravitational-wave emission from a small companion, the comparison provides a template for testing regular black holes against Kerr with future observations."],"supporting_citations":[{"why":"Supplies the rotating Hayward metric (2.3) that is the spacetime under study.","marker":"[14]"},{"why":"Carter's separation method that the paper adapts to obtain the fourth constant of motion.","marker":"[30]"},{"why":"Introduces the original Hayward regular black hole whose rotating version is analyzed.","marker":"[11]"},{"why":"Earlier study of equatorial-plane motion and effective potentials for the rotating Hayward metric that this work extends.","marker":"[29]"},{"why":"Gives the regularity analysis used to justify treating the rotating metric as a regular black hole.","marker":"[25]"},{"why":"Shows the matter content of the rotating metric is only an approximation to the Hayward magnetic-monopole source, limiting physical interpretation.","marker":"[16–18]"}],"fun_headline_variants":["Hayward black holes keep a Carter-like fourth constant","Rotating Hayward inherits Kerr's integrability","Off-plane orbits stay integrable in Hayward spacetime","At low spin, Hayward orbits wander from Kerr","Regular black hole gets fourth integral of motion"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The analysis assumes the Bambi-Modesto rotating metric (2.3) is a physically valid regular black hole spacetime; the paper itself notes that its matter source is only an approximation to the Hayward magnetic-monopole solution, so the orbit results inherit that approximation.","fun_headline_variants_meta":{"raw":{"variants":["Hayward black holes keep a Carter-like fourth constant","Rotating Hayward inherits Kerr's integrability","Off-plane orbits stay integrable in Hayward spacetime","At low spin, Hayward orbits wander from Kerr","Regular black hole gets fourth integral of motion"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000505,"raw_usage":{"total_tokens":2408,"prompt_tokens":834,"completion_tokens":1574,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":450,"completion_tokens_details":{"reasoning_tokens":1500}},"tokens_in":450,"tokens_out":1574,"duration_ms":12142,"temperature":1.0,"reasoning_tokens":1500,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:02:25.878722+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a generic timelike geodesic off the equatorial plane in metric (2.3), integrate the full second-order equations, and evaluate the separated expression $k=Q+(aE-L_z)^2$ along the orbit; if it drifts by more than numerical error, the separation ansatz is false. A complementary check is to insert the metric into the field equations with an exact nonlinear-electrodynamics stress-energy tensor and see whether $M(r)=m r^3/(r^3+g^3)$ remains a solution without extra terms.","supporting_citations":[{"cited_title":"Ay´ on-Beato and A","cited_arxiv_id":null,"evidence_quote":"Supplies the rotating Hayward metric (2.3) that is the spacetime under study."},{"cited_title":"On regular rotating black holes","cited_arxiv_id":"1611.03654","evidence_quote":"Carter's separation method that the paper adapts to obtain the fourth constant of motion."},{"cited_title":"Dymnikova","cited_arxiv_id":null,"evidence_quote":"Introduces the original Hayward regular black hole whose rotating version is analyzed."}],"review_version":1}