{"id":"7c43a015-bd8c-4734-8a0b-0c6c1b963a53","arxiv_id":"2507.11616","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A semi-implicit integrator for covariant guiding center equations tracks charged particles in GRMHD black hole accretion backgrounds over macroscopic timescales.","lead":"This paper presents a new numerical method for following charged particles in the plasma around black holes without resolving every fast spin, by tracking the drift of the particle's guiding center on a magnetized fluid simulation. If it works, researchers can study how high-energy particles are accelerated and transported in black hole accretion flows over realistically long timescales.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Fig. 3 divergence diagnosis rests on magnetic-moment equation (30), whose Appendix C derivation fixes a placeholder term by a consistency condition and whose covariant form is ill-defined; if wrong, the divergence is an unexplained mismatch.","rationale":"The reader's weakest assumption identifies exactly the same load-bearing point: the t≈15 M divergence in Fig. 3 is explained using eq. (30), and if that equation is wrong, the divergence is not shown to be a field-interpolation artifact. I agree with that assessment. The concern is load-bearing because Fig. 3 is one of only two GRMHD accuracy demonstrations, and it is the only one that shows a clear failure of trajectory matching in a turbulent region. The paper's response to that failure is the magnetic-moment diagnostic, so the validity of that diagnostic determines whether the accuracy claim survives. The Appendix C derivation is genuinely circular: ξ is a placeholder for electric-field derivative terms, and it is fixed by insisting that the total expression vanishes when Faraday's law holds. That ensures consistency with the conservation law in the Maxwell case but does not derive the coefficient for the non-Maxwellian terms that the interpolation is supposed to produce. The covariant form in eq. (30) compounds the problem, since ζ^λ = (U_∥, b) is not a well-defined four-vector contraction. This is a specific, testable weakness rather than a vague demand for more tests. A flat-spacetime test with an analytically prescribed Maxwell violation would settle whether eq. (26)/(30) correctly predicts the measured magnetic-moment evolution. I do not think this concern requires changing the reader's CONDITIONAL verdict: the semi-implicit construction is clearly described, the tricubic interpolation is a reasonable choice, and the dipole scaling test in §4.1 is a good idea. But the paper's central accuracy claim currently rests on an under-derived diagnostic equation, so conditional acceptance remains the right posture.","tokens_in":16869,"tokens_out":21138,"duration_ms":258933,"concrete_test":"In flat spacetime with no gravity, choose an analytic electromagnetic field that violates the homogeneous Maxwell equations in a controlled way: for example, B = (x, 0, 0) with E = 0 (nonzero ∇·B), or a uniform, time-dependent B(t) with E = 0 (nonzero ∂tB). Integrate the full Lorentz force with a high-accuracy RK4 scheme and compute dμ/dt from the instantaneous definition of the magnetic moment, eq. (7). Independently evaluate the right-hand side of eq. (26) (and, separately, eq. (30)) along the same trajectory. If the predicted rate disagrees with the measured rate beyond truncation error, the Fig. 3 divergence explanation in §3.2 is unsupported and the evolving-mu trajectory is not a valid diagnostic.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central accuracy claim is not fully established because the one GRMHD test that shows a clear trajectory mismatch is explained away using eq. (30), the magnetic-moment evolution equation. That equation is derived in Appendix C by introducing a placeholder term ξ (eq. C35), asserting that 'we will be able to solve for this term,' and then fixing it by requiring the magnetic moment to be invariant when Faraday's law holds (eqs. C37-C41). This is a consistency condition, not a derivation: it guarantees the final expression vanishes on Maxwell solutions, but it does not uniquely determine the response to the specific Maxwell violations introduced by tricubic interpolation unless the perturbation calculation is actually carried through. The covariant translation in eq. (30) is also suspect, since ζ^λ is defined as (U_∥, b), mixing a scalar parallel velocity with a unit spatial vector; this is not a proper four-vector contraction. Consequently, the paper has not demonstrated that the t≈15 M divergence in Fig. 3 is an interpolation/static-snapshot Maxwell artifact rather than an integrator or guiding-center error. The dashed green evolving-mu trajectory in Fig. 3 is therefore not a valid accuracy demonstration, leaving an unexplained counterexample to the claimed trajectory-matching performance.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a semi-implicit numerical integrator for the covariant guiding-center equations of Trent et al. (2024), designed to evolve charged-particle trajectories in GRMHD backgrounds over macroscopic timescales. The scheme treats the Lorentz-force term implicitly and the remaining terms explicitly, with a second-order Runge-Kutta-like half-step/full-step structure derived in Appendix A. The authors apply the method to a static Athena++ SANE accretion-flow snapshot around a spinning black hole, using tricubic interpolation of the background fields and a 3+1 reconstruction of the electromagnetic field tensor. They validate the method by comparing guiding-center trajectories with full Lorentz-force integrations for a small set of initial conditions, and they examine the impact of grid resolution in a flat-space magnetic-dipole setup, finding an error scaling of roughly ξρ0 + δr^3. The paper also derives an evolution equation for the magnetic moment and uses it to argue that a trajectory divergence seen at t≈15M is caused by interpolation-induced violations of the homogeneous Maxwell equations rather than by an integrator failure.","tokens_in":17111,"tokens_out":10330,"duration_ms":122871,"significance":"If the central claims hold, the algorithm would be a useful tool for hybrid kinetic-GRMHD simulations of black-hole accretion, because it would allow particle trajectories to be followed over macroscopic times while avoiding the gyroperiod timestep restriction. The paper has several genuine strengths: the semi-implicit integrator is derived analytically rather than merely asserted; the flat-space stability analysis is clear and shows why explicit schemes fail; the comparison against full Lorentz-force integration is an appropriate external benchmark and not circular; and the grid-resolution study in Fig. 6 provides a concrete, falsifiable error scaling. The paper also honestly identifies a real and often neglected difficulty, namely that interpolated electromagnetic fields from GRMHD snapshots do not satisfy the homogeneous Maxwell equations and can artificially evolve the magnetic moment. However, the validation rests on only a handful of GRMHD trajectories, the counterintuitive gyroradius dependence in Fig. 2 is left unexplained, and the magnetic-moment evolution equation used to diagnose the Fig. 3 divergence is not derived rigorously.","major_comments":[{"comment":"The magnetic-moment evolution equation used to diagnose the Fig. 3 divergence is not derived in a way that supports the conclusion. In Eq. (C35) a placeholder term ξ is introduced and later fixed by requiring that ⟨dμ/dt⟩ vanish when Faraday's law holds (Eqs. C37–C41). This is a consistency condition, not a derivation: it guarantees the final expression vanishes on Maxwell solutions, but it does not determine the response to the specific interpolation-induced Maxwell violations unless the perturbation calculation is actually carried through. Moreover, the covariant form in Eq. (30) is ill-defined: ζ^λ ≡ (U_∥, b) mixes a scalar parallel velocity with a spatial unit vector, so the contraction is not manifestly covariant, and the denominator sqrt(2F_{ιη}F^{ιη}) does not equal the sqrt(B²−E²) used in Eq. (26) under the stated definitions. Since the claim that the t≈15M divergence in Fig. 3 is an interpolation/static-snapshot artifact rests on Eq. (30), that claim is currently unsupported. Please replace this with a direct evaluation of ∇·B and ∂B/∂t+∇×E along the actual guiding-center trajectory (as is already done for the full trajectory in the bottom panels of Figs. 4 and 5), or provide a systematic derivation of Eq. (30).","section":"§3.2, Eq. (30), Appendix C"},{"comment":"The discrete update of the magnetic moment is never specified. The text says the scheme \"allows for the magnetic moment to be evolved\" and Fig. 3 presents a trajectory \"calculated for an evolving magnetic moment,\" but there is no equation showing how μ_{n+1/2} and μ_{n+1} are obtained from Eq. (30), nor how the evolving μ is coupled back into the semi-implicit position/velocity update in Eqs. (14)–(18). Without this information, the dashed green curve in Fig. 3 is not reproducible and the accuracy of the evolving-μ variant cannot be assessed. Please provide the full discrete update and state explicitly whether μ is evolved with the same timestep and whether it is held fixed or advanced during the half-step and full-step stages.","section":"§2.3 and Fig. 3"},{"comment":"The result that the trajectory with the larger initial gyroradius (ρ0 = 10^{-5} M) matches the full trajectory more closely than the smaller one (ρ0 = 10^{-7} M) is explicitly called counterintuitive but is never explained. This ordering is important because it bears directly on the convergence of the guiding-center approximation and on the error budget in Eq. (32), which predicts a smaller absolute error for smaller ρ0 when the δr^3 interpolation term is subdominant. If the δr^3 term dominates for these trajectories, that should be stated and demonstrated; if the explanation involves the static-snapshot Maxwell violation, it should be quantified along each trajectory. As written, the paper leaves a central accuracy result unexplained, and this should be resolved before the method's accuracy claims can be accepted.","section":"§4, Fig. 2"}],"minor_comments":[{"comment":"The abstract claims that the method demonstrates \"accuracy and efficiency,\" but no wall-clock times, step counts, or comparisons with a fully explicit integrator are reported. Please add at least a quantitative measure of the efficiency gain, or soften the claim.","section":"§4, Abstract"},{"comment":"In Eq. (C41), the second term is written first as −(μ/B)∂B/∂t and then as −(μ/B)∂B/∂t·b in the same displayed equation; this appears to be a typo and should be corrected for clarity.","section":"Eq. (C41)"},{"comment":"The statement that the tricubic interpolation \"works for non-Cartesian coordinate grids as well\" needs clarification, because the interpolated quantities are coordinate components of vector fields in Kerr-Schild spherical coordinates. The interpolation error near coordinate singularities and the transformation to Boyer-Lindquist coordinates should be discussed, especially since the guiding-center equation depends on field derivatives.","section":"§3.1"},{"comment":"In Fig. 6, the error floor is attributed to the guiding-center approximation, but the text does not state how many trajectories were used per grid resolution or what the statistical uncertainty in the plotted differences is. Please add error bars or repeat runs to support the claimed scaling.","section":"§4.1, Fig. 6"},{"comment":"The statement that a time-evolving GRMHD background would counteract the spurious curl of the electric field via ∂B/∂t is plausible but speculative; it should either be tested with a time-dependent simulation or clearly marked as a forward-looking remark rather than a conclusion of this work.","section":"§4"},{"comment":"The timestep formula in Eq. (19) uses ξ = 10^{-3}, but no convergence study with respect to ξ is shown for the GRMHD backgrounds. Please add a brief study or a statement of how this value was established.","section":"Eq. (19)"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the journal's scope and addresses a relevant problem in hybrid kinetic-GRMHD modeling. The main obstacle to acceptance is that the validation of the method relies on the magnetic-moment evolution equation (30), which is not rigorously derived and whose covariant form is questionable. The Fig. 3 divergence could be a genuine counterexample to the claimed trajectory matching, and the current manuscript does not rule out an integrator or guiding-center error. I would encourage the editor to request a revision in which the authors either provide a direct calculation of the interpolation-induced Maxwell violations along the guiding-center trajectories or remove the evolving-μ diagnostic from the central accuracy claim. The self-citation to Trent et al. (2024) is legitimate because the guiding-center equations are the foundation of the algorithm, but the paper should make sure it does not depend on unpublished details from that work."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Everything you need to know: this is a genuine new numerical method—a semi-implicit integrator for the covariant guiding center equations, coupled to GRMHD data with tricubic interpolation—and the core derivation is coherent. It deserves a serious referee. But the accuracy demonstrations are thin, and the one dramatic trajectory divergence is explained using an equation whose derivation in Appendix C is hand-wavy.\n\nWhat's new: the integrator (eqs. 14–18) treats the unstable linear term implicitly and the flat-space stability analysis is sound. The coupling to Athena++ (Appendix B) is careful: interpolate B and velocity first, then construct the force-free field tensor. The scaling epsilon ~ xi rho0 + delta-r^3 is practically useful. They also deserve credit for flagging interpolation-induced magnetic-moment artifacts and quantifying them, and for validating against full Lorentz-force trajectory integration, which is the right external benchmark. The self-citations to Trent et al. (2023, 2024) are appropriate since the covariant equations are their own prior result.\n\nThe soft spots. The tests are only a handful of trajectories. Fig. 2 shows the counterintuitive result that a larger gyroradius gives better agreement, and they never explain it. That is a red flag. The efficiency claim—timestep set by the drift scale rather than the gyroperiod—is never demonstrated against an explicit integrator; they assert explicit schemes fail but show no comparison. Most substantively, the Fig. 3 divergence at t ≈ 15 M is blamed on interpolation Maxwell violations, but eq. (30), which supports that diagnosis, comes from adding a placeholder term xi in Appendix C and fixing it via a consistency condition. The covariant form zeta^lambda = (U_parallel, b) mixes a scalar and a unit vector; it does not look like a proper four-vector contraction. The empirical match between the mu evolution from eq. (30) and from the full trajectory in Figs. 4–5 is encouraging, but the derivation needs to be tightened before I would fully trust the diagnosis. No code or public data is provided, which limits independent verification.\n\nWho this is for: anyone building hybrid kinetic/GRMHD tools for black hole accretion flows. The method is plausible and the limitations are honestly stated. I would send it to review, asking the authors to add an explicit-integrator comparison, clarify Fig. 2, and fix or re-derive eq. (30). If those land, it becomes a reference method.","headline":"A coherent semi-implicit integrator for covariant guiding center motion coupled to GRMHD, but thin trajectory tests and a hand-wavy magnetic-moment derivation keep it from being a finished community tool.","tokens_in":17640,"tokens_out":6220,"would_cite":true,"duration_ms":59716,"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":"A second-order semi-implicit integrator tracks charged-particle guiding centers in black-hole accretion flows without resolving gyration, and exposes interpolation-driven magnetic-moment artifacts.","keywords":["guiding center","drift-kinetic","GRMHD","semi-implicit integrator","tricubic interpolation","magnetic moment","black hole accretion","test-particle dynamics"],"falsifier":"Re-run the Fig. 3 case at the same grid resolution but reconstruct the field by interpolating a vector potential, so that $\\nabla \\cdot B = 0$ and Faraday's law hold by construction; if the full and guiding-center trajectories still separate at $t \\approx 15 M$, the divergence cannot be blamed on interpolation Maxwell violations, and the paper's diagnosis fails.","tokens_in":16642,"feed_emoji":"🕳️","tokens_out":11956,"duration_ms":130789,"temperature":0.7,"pith_summary":"The paper claims that the recently derived covariant guiding-center equations of motion for charged particles can be integrated over macroscopic times in realistic general-relativistic magnetohydrodynamics (GRMHD) black-hole accretion-flow backgrounds, without resolving the gyroperiod. The key move is a second-order semi-implicit integrator: the unstable linear Lorentz term is treated implicitly and inverted analytically, while the stable gravitational and mirror terms remain explicit, so the timestep can be set by the drift scale rather than the gyroperiod. Comparing against full Lorentz-force integrations on a GRMHD accretion-flow snapshot, the authors find agreement at the level of a combined error budget $\\epsilon \\simeq \\xi \\rho_0 + \\delta r^3$, where $\\rho_0$ is the gyroradius and $\\delta r$ the grid spacing. They further show that interpolation of the discretized fields violates the homogeneous Maxwell equations and artificially changes the magnetic moment, and they derive a covariant equation for that spurious growth. If this is right, the method supplies a practical bridge from fluid GRMHD simulations to kinetic transport and particle acceleration in black-hole accretion flows.","feed_headline":"Semi-implicit integrator follows guiding centers in black hole flows","feed_subtitle":"Charged-particle orbits in GRMHD flows advance with drift-scale steps; interpolation artifacts are exposed.","key_machinery":"The load-bearing object is the two-stage semi-implicit integrator in eqs. (14)-(18), which treats only the linear velocity term $(q/m)F^{\\alpha}_{\\beta} U^{\\beta}$ implicitly, averages it between old and new velocity, and solves the resulting linear system analytically, while keeping the Christoffel term, the $\\mu \\nabla^{\\alpha} \\omega$ drift/mirror term, and the position update explicit. The timestep is set by eq. (19), the time for the guiding center to drift into a region of appreciably different field, with a choice $\\xi = 10^{-3}$; the conserved norm $U^{\\alpha} U_{\\alpha} + 2 \\mu \\omega = -1$ is used to update the time component of the four-velocity. Around this core sit the tricubic interpolation of the background magnetic field and fluid velocity (with the electromagnetic tensor reconstructed afterward to enforce the force-free condition), and the covariant magnetic-moment evolution equation (30), whose derivation fixes a free term by demanding that $\\mu$ be constant when Faraday's law holds. That equation is what lets the paper attribute trajectory errors to interpolation-induced Maxwell violations rather than to the integrator.","core_discovery":"The central claim is that the covariant guiding-center acceleration equation $dU^{\\alpha}/d\\tau = -\\Gamma^{\\alpha}_{\\beta\\nu} U^{\\beta} U^{\\nu} + (q/m) F^{\\alpha}_{\\beta} U^{\\beta} - \\mu \\nabla^{\\alpha} \\omega$ can be evolved stably with a second-order semi-implicit scheme (eqs. 14-18) whose timestep is set by the drift scale, not the gyroperiod. In the GRMHD tests shown, guiding-center trajectories track full Lorentz-force integrations within an error budget $\\epsilon \\simeq \\xi \\rho_0 + \\delta r^3$, with $\\xi$ of order a few; the gyroradius $\\rho_0$ sets the floor and interpolation error falls as the cube of the grid spacing. The paper also establishes that interpolating the electromagnetic field from a discretized GRMHD snapshot introduces violations of $\\nabla \\cdot B = 0$ and of Faraday's law, that these violations drive a spurious evolution of the magnetic moment, and that the covariant equation (30) quantitatively accounts for the divergence seen between full and guiding-center trajectories at $t \\approx 15 M$ in one test case. The same artifact affects full-trajectory integrations on interpolated fields, so the diagnosis is not specific to the guiding-center approximation.","pith_inferences":["The authors leave implicit that eq. (30) can be turned into an online monitor: a hybrid simulation could flag or compensate for interpolation-induced pitch-angle diffusion as it integrates.","A corollary of the error scaling is a grid-design rule for kinetic post-processing: to keep interpolation error below the gyroradius floor one wants $\\delta r \\lesssim (\\xi \\rho_0)^{1/3}$, a constraint that may bind more tightly than the guiding-center validity conditions themselves.","If the diagnosis of the $t \\approx 15 M$ divergence is correct, then vector-potential-based interpolation, which preserves the homogeneous Maxwell equations by construction, should eliminate the divergence at fixed resolution; this is a natural test the paper does not run.","The same machinery, applied to time-dependent backgrounds, might reveal whether the spurious magnetic-moment growth seen here is largely a static-snapshot artifact, and whether prior hybrid particle-MHD studies of turbulent plasmas were biased by interpolation Maxwell violations."],"forward_implications":["Charged-particle trajectories in inner accretion flows can now be integrated for durations of tens of $M$ with steps set by the drift scale, making kinetic post-processing of GRMHD snapshots practical.","The scaling $\\epsilon \\simeq \\xi \\rho_0 + \\delta r^3$ means the guiding-center approximation itself sets the accuracy floor, and increasing grid resolution pays off only cubically once interpolation dominates.","Any simulation that interpolates electromagnetic fields from a discretized grid is exposed to the same artifact: magnetic-moment conservation is broken by interpolation-induced violations of the homogeneous Maxwell equations, regardless of whether a guiding-center reduction is used.","Using a single static GRMHD snapshot leaves $\\partial B/\\partial t = 0$, so the interpolated field has an unbalanced curl of the electric field; time-evolving backgrounds should largely remove this particular contribution to the spurious magnetic-moment growth.","With time-dependent backgrounds, the method provides a route to studying non-thermal electron acceleration and the confinement of hot spots in flares around Sgr A*."],"supporting_citations":[{"why":"derives the covariant guiding-center equations of motion, eq. (8), that the new integrator solves","marker":"Trent et al. 2024"},{"why":"provides the base magnetic-moment evolution result for non-uniform magnetic fields that Appendix C extends to time-dependent perpendicular electric fields","marker":"Brizard & Markowski 2022"},{"why":"supplies the tricubic interpolation scheme used to reconstruct GRMHD fields and their derivatives","marker":"Lekien & Marsden (2005)"},{"why":"provides the GRMHD code used to generate the black-hole accretion-flow background","marker":"Stone et al. 2020"},{"why":"supplies the GRMHD accretion-flow snapshot used as the background for the trajectory tests","marker":"Avara et al. 2025"},{"why":"provides the definitions of electric and magnetic fields in terms of the electromagnetic field tensor used to make the magnetic-moment equation covariant","marker":"Komissarov 2011"},{"why":"documents that magnetic moment is not conserved on discretized backgrounds, motivating the paper's diagnostic","marker":"Ripperda et al. 2018"},{"why":"supplies the drift solution and perturbation analysis used to identify the unstable linear term in the guiding-center equations","marker":"Northrop 1963"},{"why":"establishes the GRMHD formulation and field-construction conventions that Appendix B follows","marker":"Gammie et al. 2003"}],"fun_headline_variants":["Semi-implicit integrator for guiding centers in black hole flows","Drift-scale particle tracking in black hole accretion flows","Exposing interpolation artifacts in black hole kinetic tracking","Guiding center orbits at drift scale in black hole accretion"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the paper's derived equation (30) correctly describes how interpolation-induced violations of Maxwell's equations change a particle's magnetic moment; if that equation is wrong, the paper has not shown that the $t \\approx 15 M$ trajectory divergence is an interpolation artifact rather than a failure of the integrator.","fun_headline_variants_meta":{"raw":{"variants":["Semi-implicit integrator for guiding centers in black hole flows","Drift-scale particle tracking in black hole accretion flows","Exposing interpolation artifacts in black hole kinetic tracking","Guiding center orbits at drift scale in black hole accretion"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001292,"raw_usage":{"total_tokens":5289,"prompt_tokens":971,"completion_tokens":4318,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":587,"completion_tokens_details":{"reasoning_tokens":4252}},"tokens_in":587,"tokens_out":4318,"duration_ms":34990,"temperature":1.0,"reasoning_tokens":4252,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T17:06:17.435948+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-run the Fig. 3 case at the same grid resolution but reconstruct the field by interpolating a vector potential, so that $\\nabla \\cdot B = 0$ and Faraday's law hold by construction; if the full and guiding-center trajectories still separate at $t \\approx 15 M$, the divergence cannot be blamed on interpolation Maxwell violations, and the paper's diagnosis fails.","supporting_citations":[{"cited_title":"Covariant Guiding Center Equations for Charged Particle Motions in General Relativistic Spacetimes","cited_arxiv_id":"2404.01391","evidence_quote":"derives the covariant guiding-center equations of motion, eq. (8), that the new integrator solves"},{"cited_title":"J., & Markowski, D","cited_arxiv_id":null,"evidence_quote":"provides the base magnetic-moment evolution result for non-uniform magnetic fields that Appendix C extends to time-dependent perpendicular electric fields"},{"cited_title":"2005, International Journal for Numerical Methods in Engineering, 63, 455, https://doi.org/10.1002/nme.1296","cited_arxiv_id":null,"evidence_quote":"supplies the tricubic interpolation scheme used to reconstruct GRMHD fields and their derivatives"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the GRMHD accretion-flow snapshot used as the background for the trajectory tests"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the drift solution and perturbation analysis used to identify the unstable linear term in the guiding-center equations"}],"review_version":1}