pith. machine review for the scientific record. sign in

arxiv: 2605.14065 · v1 · submitted 2026-05-13 · ⚛️ physics.plasm-ph · astro-ph.SR

Recognition: no theorem link

Boris and Exponential Integrators in the Theory of Particles Interacting with Magnetic Turbulence

Andreas Shalchi
Authors on Pith no claims yet

Pith reviewed 2026-05-15 02:25 UTC · model grok-4.3

classification ⚛️ physics.plasm-ph astro-ph.SR
keywords charged particlesmagnetic turbulenceexponential integratorsRodrigues schemeBoris integratorNewton-Lorentz equationtest particle simulationsplasma physics
0
0 comments X

The pith

Rodrigues scheme from exponential integrators matches Boris accuracy for particle motion in magnetic turbulence without extra cost

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper derives both the Rodrigues scheme and the Boris integrator from exponential integrators for solving the Newton-Lorentz equation governing charged particle motion in magnetic fields. For test-particle simulations in turbulent magnetic fields where the field is generated anew at every time step, the two methods produce nearly identical results in practice. Although the Rodrigues approach uses the exact matrix exponential without approximation, this does not lead to measurable improvements in the tested scenarios. The key insight is that the Rodrigues integrator serves as a strong alternative because it achieves this performance without requiring longer computing times.

Core claim

The central claim is that exponential integrators yield the Rodrigues scheme, which evaluates the matrix exponential exactly, and this provides comparable accuracy and speed to the Boris integrator in simulations of particles interacting with magnetic turbulence. When the turbulent field is recreated at each step, the practical differences vanish, allowing the Rodrigues method to be used as an efficient alternative for exploring particle transport in magnetized plasmas.

What carries the argument

The Rodrigues formula for the exact matrix exponential in the velocity update of the Newton-Lorentz equation, which rotates the particle velocity precisely under the magnetic force.

Load-bearing premise

The assumption that recreating the magnetic field at each time step equalizes the practical performance of the Rodrigues and Boris schemes despite the exact exponential in the former.

What would settle it

A simulation in a static uniform magnetic field over thousands of gyroperiods, measuring if the perpendicular velocity magnitude remains exactly constant in Rodrigues while showing drift in Boris.

Figures

Figures reproduced from arXiv: 2605.14065 by Andreas Shalchi.

Figure 1
Figure 1. Figure 1: — The phase defined via Eq. (36) for the case of a particle interacting with a constant magnetic field pointing in the z-direction. Shown is the ratio of the numerically obtained phase for a step size of ∆t = 0.1 and the exact analytical result for a Rodrigues-based integrator (blue line) and the Boris scheme (red line). as it was done above (see Eqs. (7)-(10)). As in Sect. 3 we, therefore, have employed t… view at source ↗
Figure 2
Figure 2. Figure 2: — The parallel diffusion coefficient obtained from the simulations by employing the Rodrigues integrator. We have performed the simulations for different values of the step size, namely ∆t = 1 (green line), ∆t = 0.1 (blue line), and ∆t = 0.01 (red line). Here we have considered the case of small rigidities R = 0.1. Time 0 500 1000 1500 2000 2500 3000 3500 4000 4500 5000 0 0.5 1 1.5 2 2.5 3 Dif f u sio n C … view at source ↗
Figure 3
Figure 3. Figure 3: — Caption is as in [PITH_FULL_IMAGE:figures/full_fig_p010_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: — The parallel diffusion coefficient obtained from the simulations by employing the Rodrigues integrator. We have performed the simulations for different values of the step size, namely ∆t = 1 (green line), ∆t = 0.1 (blue line), ∆t = 0.01 (red line), and ∆t = 0.001 (black line). Here we have considered the case of large rigidities R = 10. Time 500 1000 1500 2000 2500 3000 3500 4000 4500 5000 Dif f u sio n … view at source ↗
Figure 5
Figure 5. Figure 5: — Caption is as in [PITH_FULL_IMAGE:figures/full_fig_p011_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: — Caption is as in [PITH_FULL_IMAGE:figures/full_fig_p011_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: — Rigidity conservation for a Rodrigues integrator. Time 0 500 1000 1500 2000 2500 3000 3500 4000 4500 5000 Dif f u sio n C o e f ficie n t 0 0.05 0.1 0.15 0.2 0.25 0.3 [PITH_FULL_IMAGE:figures/full_fig_p012_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: — The parallel diffusion coefficient obtained from the simulations by employing the Boris integrator. We have performed the simulations for different values of the step size, namely ∆t = 1 (green line), ∆t = 0.1 (blue line), and ∆t = 0.01 (red line). Here we have considered the case of small rigidities R = 0.1. Note, red and blue lines are in perfect coincidence. Time 0 500 1000 1500 2000 2500 3000 3500 40… view at source ↗
Figure 9
Figure 9. Figure 9: — Caption is as in [PITH_FULL_IMAGE:figures/full_fig_p012_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: — The parallel diffusion coefficient obtained from the simulations by employing the Boris integrator. We have performed the simulations for different values of the step size, namely ∆t = 1 (green line), ∆t = 0.1 (blue line), ∆t = 0.01 (red line), and ∆t = 0.001 (black line). Here we have considered the case of large rigidities R = 10. Time 0 500 1000 1500 2000 2500 3000 3500 4000 4500 5000 Dif f u sio n C… view at source ↗
Figure 11
Figure 11. Figure 11: — Caption is as in [PITH_FULL_IMAGE:figures/full_fig_p013_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: — Caption is as in [PITH_FULL_IMAGE:figures/full_fig_p013_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: — Rigidity conservation for a Boris integrator. particles interacting with magnetic turbulence numerically. Support by the Natural Sciences and Engineering Research Council (NSERC) of Canada is acknowledged. APPENDIX THE MOORE-PENROSE INVERSE For the Rodrigues approach the updated position can also be obtained by using Z ∆t 0 dτ eMτ = M+ [PITH_FULL_IMAGE:figures/full_fig_p014_13.png] view at source ↗
read the original abstract

The interaction of electrically charged particles with magnetic fields is a fundamental problem in several areas of physics. An example is the motion of energetic particles through a magnetized plasma. The most accurate and reliable way to explore theoretically the interactions between particles and fields is via test-particle simulations. In such simulations one creates the turbulent magnetic field and solves the Newton-Lorentz equation numerically by employing an integration scheme. In the current article we discuss exponential integrators and derive systematically from this the Rodrigues scheme as well as the famous Boris integrator. For an approach where one creates the magnetic field anew at each time step, both integrators are overall comparable. In theory the Rodrigues approach should be more accurate due to the fact that the occurring matrix exponential is evaluated without further approximations. Practically, both methods provide very similar results. It is argued in the current article that a Rodrigues based integrator is a very strong alternative because for the specific problem discussed here, it does not require longer computing times.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Circularity Check

0 steps flagged

Derivation from general exponential-integrator framework recovers known schemes without circular reduction

full rationale

The paper begins with the general exponential-integrator framework and systematically derives both the Rodrigues scheme and the Boris integrator as special cases. This is a forward mathematical construction that does not reduce any claimed result to its own inputs by definition, fitting, or self-citation. No load-bearing steps invoke prior self-citations as uniqueness theorems, smuggle ansatzes, or rename empirical patterns as new predictions. The practical comparison of run times for the specific turbulence problem is presented as an empirical observation within the described simulation setup rather than a derived necessity. The central claim therefore remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

The abstract does not identify any free parameters, axioms, or invented entities; the work rests on standard mathematical properties of matrix exponentials and the Newton-Lorentz equation.

pith-pipeline@v0.9.0 · 5463 in / 1053 out tokens · 132540 ms · 2026-05-15T02:25:33.072196+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

51 extracted references · 51 canonical work pages

  1. [1]

    J.420, 294 (1994)

    Kallenrode, M.-B., Wibberenz, G.: Astrophys. J.420, 294 (1994)

    work page 1994

  2. [2]

    Bieber, J.W., Wanner, W., Matthaeus, W.H.: J. Geophys. Res. 101, 2511 (1996) 16

    work page 1996

  3. [3]

    Burden, R.L., Faires, J.D.: Numerical Analysis, 9th ed., Brooks/Cole, Cengage Learning, Boston (2011)

    work page 2011

  4. [4]

    Space Res.461208 (2010)

    Dosch, A., Shalchi A.: Adv. Space Res.461208 (2010)

    work page 2010

  5. [5]

    Effenberger, F., Walter, D., Fichtner, H., et al.: Space Sci. Rev. 221, 75 (2025)

    work page 2025

  6. [6]

    J.969, 1 (2024)

    Els, P., Engelbrecht, N.E.: Astrophys. J.969, 1 (2024)

    work page 2024

  7. [7]

    Space Res.66, 2722 (2020)

    Engelbrecht, E.N., Wolmarans, C.P.: Adv. Space Res.66, 2722 (2020)

    work page 2020

  8. [8]

    J.908, 167 (2021)

    Engelbrecht, E.N., Moloto, K.D.: Astrophys. J.908, 167 (2021)

    work page 2021

  9. [9]

    J.792, 133 (2014)

    Ferrand, G., Danos, R., Shalchi, A., Safi-Harb, S., Mendygral, P.: Astrophys. J.792, 133 (2014)

    work page 2014

  10. [10]

    Hairer, E., Lubich, C., Wanner, G.: Geometric numerical integration: Structure-preserving algorithms for ordinary differential equations (2nd ed.), Springer (2006)

    work page 2006

  11. [11]

    Space Sci.361, 308 (2016)

    Heusen, M., Shalchi, A.: Astrophys. Space Sci.361, 308 (2016)

    work page 2016

  12. [12]

    Thesis, University of Manitoba (2017)

    work page 2017

  13. [13]

    Hu, J., Li, G., Ao, X., Zank, G.P., Verkhoglyadova, O.: J. Geophys. Res.122, 938 (2017)

    work page 2017

  14. [14]

    Hussein, M., Tautz, R.C., Shalchi, A.: J. Geophys. Res.120, 4095 (2015)

    work page 2015

  15. [15]

    J.817, 136 (2016)

    Hussein, M., Shalchi, A.: Astrophys. J.817, 136 (2016)

    work page 2016

  16. [16]

    Ivascenko, A., Lange, S., Spanier, F., Vainio, R.: Astrophys. J. 833, 223 (2016)

    work page 2016

  17. [17]

    Kolmogorov, A.N.: Dokl. Akad. Nauk. SSSR30, 301 (1941)

    work page 1941

  18. [18]

    Li, G., Zank, G.P., Rice, W.K.M.: J. Geophys. Res.108, 1082 (2003)

    work page 2003

  19. [19]

    Li, G., Zank, G.P., Rice, W.K.M.: J. Geophys. Res.110, A06104 (2005)

    work page 2005

  20. [20]

    Space Res.49, 1067 (2012)

    Li, G., Shalchi, A., Ao, X., Zank, G.P., Verkhoglyadova, O.P.: Adv. Space Res.49, 1067 (2012)

    work page 2012

  21. [21]

    Springer-Verlag, New York, (1999)

    work page 1999

  22. [22]

    Matthaeus, W.H., Goldstein, M.L., Aaron, R.D.: J. Geophys. Res. 95, 20673 (1990)

    work page 1990

  23. [23]

    J.590, L53 (2003)

    Matthaeus, W.H., Qin, G., Bieber, J.W., Zank, G.P.: Astrophys. J.590, L53 (2003)

    work page 2003

  24. [24]

    J.894, 121 (2020)

    Moloto, K.D., Engelbrecht, N.E.: Astrophys. J.894, 121 (2020)

    work page 2020

  25. [25]

    Moore, E.H.: Bulletin of the American Mathematical Society26, 394 (1920)

    work page 1920

  26. [26]

    Space Res.69, 2330 (2022) Pad´ e, H.: Sur la r´ epresentation approch´ ee d’une fonction par des fractions rationelles (Thesis), Ann

    Ramokgaba, I.I., Ndiitwani, D.C.: Adv. Space Res.69, 2330 (2022) Pad´ e, H.: Sur la r´ epresentation approch´ ee d’une fonction par des fractions rationelles (Thesis), Ann. ´Ecole Nor. 3, 3-93 (1892)

    work page 2022

  27. [27]

    Penrose, R.: A generalized inverse for matrices, Proceedings of the Cambridge Philosophical Society, 51, 406 (1955)

    work page 1955

  28. [28]

    J.671, L177 (2007)

    Perri, S., Zimbardo, G.: Astrophys. J.671, L177 (2007)

    work page 2007

  29. [29]

    Plasma14, 012311 (2007)

    Pommois, P., Zimbardo, G., Veltri, P.: Phys. Plasma14, 012311 (2007)

    work page 2007

  30. [30]

    J.640, L103 (2006)

    Qin, G., Matthaeus, W.H., Bieber, J.W.: Astrophys. J.640, L103 (2006)

    work page 2006

  31. [31]

    J.860, 3 (2018)

    Qin, G., Kong, F.-J., Zhang, L.-H.: Astrophys. J.860, 3 (2018)

    work page 2018

  32. [32]

    Plasmas20, 084503 (2013)

    Qin, H., Zhang, S., Xiao, J., Liu, J., Sun, Y., Tang, W.M.: Phys. Plasmas20, 084503 (2013)

    work page 2013

  33. [33]

    G., et al.: Mon

    Reichherzer, P., Becker Tjus, J., Zweibel, E. G., et al.: Mon. Not. R. Astron. Soc.514, 2658 (2022)

    work page 2022

  34. [34]

    J.235, 21 (2018)

    Sironi, L., Lapenta, G., Keppens, R.: Astrophys. J.235, 21 (2018)

    work page 2018

  35. [35]

    J.854, 137 (2018)

    Shen, Z.-N., Qin, G.: Astrophys. J.854, 137 (2018)

    work page 2018

  36. [36]

    J.256, 18 (2021)

    Shen, Z., Qin, G., Zuo, P., Wei, F., Xu, X.: Astrophys. J.256, 18 (2021)

    work page 2021

  37. [37]

    362, Springer, Berlin (2009)

    Shalchi, A.: Nonlinear Cosmic Ray Diffusion Theories, Astrophysics and Space Science Library, vol. 362, Springer, Berlin (2009)

    work page 2009

  38. [38]

    J.720, L127 (2010)

    Shalchi, A.: Astrophys. J.720, L127 (2010)

    work page 2010

  39. [39]

    J.881, L27 (2019)

    Shalchi, A.: Astrophys. J.881, L27 (2019)

    work page 2019

  40. [40]

    Rev.216, 23 (2020)

    Shalchi, A.: Space Sci. Rev.216, 23 (2020)

    work page 2020

  41. [41]

    J.923, 209 (2021)

    Shalchi, A.: Astrophys. J.923, 209 (2021)

    work page 2021

  42. [42]

    J.932, 127 (2022)

    Snodin, A.P., Jitsuk, T., Ruffolo, D., Matthaeus, W.H.: Astrophys. J.932, 127 (2022)

    work page 2022

  43. [43]

    Tautz, R.C.: Comput. Phys. Commun.181, 71 (2010)

    work page 2010

  44. [44]

    Astrophys.555, A101 (2013)

    Tautz, R.C., Lerche, I., Kruse, F.: Astron. Astrophys.555, A101 (2013)

    work page 2013

  45. [45]

    Plasmas15, 056701 (2008)

    Vay, J.-L.: Phys. Plasmas15, 056701 (2008)

    work page 2008

  46. [46]

    J.651, 211 (2006)

    Webb, G.M., Zank, G.P., Kaghashvili, E.Kh., le Roux, J.A.: Astrophys. J.651, 211 (2006)

    work page 2006

  47. [47]

    Yoshida, H.: Phys. Lett. A150, 262 (1990)

    work page 1990

  48. [48]

    Zank, G.P., Rice, W.K.M., Wu, C.C.: J. Geophys. Res.105, 25079 (2000)

    work page 2000

  49. [49]

    Zank, G.P., Li, G., Florinski, V., Hu, Q., Lario, D., Smith, C.W.: J. Geophys. Res.111, A06108 (2006) Zank G.P., Nakanotani M., Webb, G.M.: Astrophys. J.88, L116 (2019)

    work page 2006

  50. [50]

    J.639, L91 (2006)

    Zimbardo, G., Pommois, P., Veltri, P.: Astrophys. J.639, L91 (2006)

    work page 2006

  51. [51]

    Space Res

    Zimbardo, G., Perri, S., Pommois, P., Veltri, P.: Adv. Space Res. 49, 1633 (2012) This paper was built using the Open Journal of Astrophysics LATEX template. The OJA is a journal which provides fast and easy peer review for new papers in theastro-phsection of the arXiv, making the reviewing process simpler for authors and referees alike. Learn more at htt...

    work page 2012