pith. sign in

arxiv: 2509.25888 · v1 · submitted 2025-09-30 · 🌌 astro-ph.HE · math-ph· math.MP· nlin.CD

The Stochastic-Dissipative St\"{o}rmer Problem-Trajectories and Radiation Patterns

Pith reviewed 2026-05-18 12:36 UTC · model grok-4.3

classification 🌌 astro-ph.HE math-phmath.MPnlin.CD
keywords Störmer problemmagnetic dipole fieldLangevin equationelectromagnetic radiationpower spectral densitycharged particle motiondissipationstochastic forces
0
0 comments X

The pith

Friction and stochastic forces in the Störmer problem produce electromagnetic radiation spectra with peaks at specific frequency intervals.

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

The paper extends the classical nonrelativistic Störmer problem of charged particle motion in a magnetic dipole field by adding velocity-proportional friction and stochastic forces. These effects are incorporated through a generalized Langevin equation that modifies the standard Lorentz force. Numerical simulations track the resulting particle trajectories and compute the electromagnetic radiation emitted during the motion. The power spectral density of this radiation consistently displays peaks within certain frequency intervals for every Störmer-type model examined.

Core claim

When dissipation and stochastic forces are included, the motion of charged particles in a purely magnetic dipole field follows a generalized Langevin equation. Detailed numerical analysis of trajectories and emitted electromagnetic power shows that the power spectral density exhibits peaks in the radiation spectrum corresponding to specific frequency intervals across all considered models, yielding a full description of the spectrum and motion properties.

What carries the argument

The generalized Langevin equation that augments the Lorentz force with friction proportional to velocity and stochastic terms, followed by numerical computation of trajectories and power spectral densities of the emitted radiation.

If this is right

  • Particle trajectories in the dipole field deviate from classical paths once friction and stochastic forces are active.
  • The emitted electromagnetic power varies with the chosen dissipation coefficient and stochastic force amplitude.
  • Peaks in the power spectral density appear at particular frequency intervals and characterize the radiation for every model.
  • Numerical integration provides a complete mapping of the radiation spectrum and associated physical motion properties.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same framework could be used to model radiation signatures from particles trapped in planetary or stellar magnetic dipole fields.
  • Observable peaks in real astrophysical spectra might indicate the presence of dissipative processes not captured by ideal models.
  • Testing the same setup with relativistic corrections would show whether the frequency locations of the peaks shift at high speeds.

Load-bearing premise

Dissipation is modeled specifically as friction proportional to velocity, with stochastic forces added through a standard Langevin extension that omits relativistic or quantum corrections.

What would settle it

Perform a laboratory experiment with charged particles in a controlled magnetic dipole field, apply measured friction and random forces, and record the emitted radiation spectrum to check for the predicted peaks at specific frequencies.

Figures

Figures reproduced from arXiv: 2509.25888 by Gabriela Raluca Mocanu, Tiberiu Harko.

Figure 1
Figure 1. Figure 1: FIG. 1. Periodic motion, radiation power, and PSD in the Clas [PITH_FULL_IMAGE:figures/full_fig_p008_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Chaotic motion, radiation power and PSD in the Classi [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Periodic motion in the Classical Dissipative St¨orm [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Chaotic motion in the Classical Dissipative St¨orme [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Periodic motion in the Stochastic-Dissipative St¨o [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Periodic motion in the Stochastic-Dissipative St¨o [PITH_FULL_IMAGE:figures/full_fig_p010_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Periodic motion in the Stochastic-Dissipative St¨o [PITH_FULL_IMAGE:figures/full_fig_p010_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Chaotic motion in the Stochastic-Dissipative St¨or [PITH_FULL_IMAGE:figures/full_fig_p010_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Chaotic motion in the Stochastic-Dissipative St¨or [PITH_FULL_IMAGE:figures/full_fig_p011_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Chaotic motion in the Stochastic-Dissipative St¨o [PITH_FULL_IMAGE:figures/full_fig_p011_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11. Brownian Motion in the Dissipative-Stochastic St¨o [PITH_FULL_IMAGE:figures/full_fig_p013_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12. Stochastic-Dissipative St¨ormer Problem escape r [PITH_FULL_IMAGE:figures/full_fig_p013_12.png] view at source ↗
read the original abstract

We consider a generalization of the classical nonrelativistic St\"{o}rmer problem, describing the motion of charged particles in a purely magnetic dipole field, by taking into account the effects of the dissipation, assumed to be of friction type, proportional to the velocity of the particle, and of the presence of stochastic forces. In the presence of dissipative/stochastic effects, the motion of the particle in the magnetic dipole field can be described by a generalized Langevin type equation, which generalizes the standard Lorentz force equation. We perform a detailed numerical analysis of the dynamical behavior of the particles in a magnetic dipolar field in the presence of dissipative and stochastic forces, as well as of the electromagnetic radiation patterns emitted during the motion. The effects of the dissipation coefficient and of the stochastic force on the particle motion and on the emitted electromagnetic power are investigated, and thus a full description of the spectrum of the magnetic dipole type electromagnetic radiation and of the physical properties of the motion is also obtained. The power spectral density of the emitted electromagnetic power is also obtained for each case, and, for all considered St\"{o}rmer type models, it shows the presence of peaks in the radiation spectrum, corresponding to certain intervals of the frequency.

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.

Referee Report

2 major / 1 minor

Summary. The manuscript generalizes the classical nonrelativistic Störmer problem for charged particles in a magnetic dipole field by adding velocity-proportional friction (dissipation) and stochastic forces, yielding a generalized Langevin equation. It reports numerical trajectories, electromagnetic radiation patterns, and power spectral densities, claiming that PSD peaks appear at specific frequency intervals for all considered Störmer-type models.

Significance. If the numerical results hold after verification, the work would extend classical charged-particle dynamics to include dissipative and stochastic effects, potentially informing models of radiation in astrophysical magnetospheres. The explicit computation of spectra under varying dissipation and stochastic strengths provides a concrete parameter study.

major comments (2)
  1. [Numerical methods] The integration procedure for the generalized Langevin equation (integrator type, time step, ensemble size, and convergence tests) is not described. Because the central claim—that PSD peaks exist for all models—rests entirely on these trajectories, the absence of such checks leaves open the possibility that reported peaks arise from aliasing of stochastic noise rather than the underlying dynamics.
  2. [Results (PSD analysis)] PSD results are presented without error bars, comparison to the deterministic (zero-stochastic) baseline, or quantitative measures of peak significance. This weakens the assertion that the peaks are robust features corresponding to distinct frequency intervals.
minor comments (1)
  1. [Abstract and introduction] Notation for the dissipation coefficient and stochastic force amplitude should be introduced once with explicit symbols and then used consistently when discussing their effects on trajectories and spectra.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments, which will help improve the clarity and robustness of our results. We address each major comment below.

read point-by-point responses
  1. Referee: [Numerical methods] The integration procedure for the generalized Langevin equation (integrator type, time step, ensemble size, and convergence tests) is not described. Because the central claim—that PSD peaks exist for all models—rests entirely on these trajectories, the absence of such checks leaves open the possibility that reported peaks arise from aliasing of stochastic noise rather than the underlying dynamics.

    Authors: We agree that the description of the numerical integration procedure was incomplete in the original submission. In the revised manuscript we will add a dedicated subsection specifying the integrator employed for the stochastic differential equation (an adapted stochastic Runge-Kutta scheme), the fixed time step, the number of independent realizations used for ensemble averaging, and the convergence tests performed with respect to time step and ensemble size. These additions will allow readers to assess the reliability of the trajectories and to exclude numerical artifacts such as aliasing in the computed power spectral densities. revision: yes

  2. Referee: [Results (PSD analysis)] PSD results are presented without error bars, comparison to the deterministic (zero-stochastic) baseline, or quantitative measures of peak significance. This weakens the assertion that the peaks are robust features corresponding to distinct frequency intervals.

    Authors: We acknowledge that the PSD figures would be strengthened by the inclusion of statistical measures. In the revision we will augment the PSD plots with error bars obtained from the standard deviation across the ensemble, add direct overlays of the deterministic (zero-stochastic-force) spectra for comparison, and report quantitative indicators of peak significance such as peak-to-background ratios. These changes will better substantiate that the reported spectral peaks are robust features of the radiation emitted in all considered Störmer-type models. revision: yes

Circularity Check

0 steps flagged

No circularity: results follow from direct numerical integration of the stated generalized Langevin equation

full rationale

The paper defines a generalized Langevin equation that augments the Lorentz force with a velocity-proportional friction term and a stochastic force, then integrates the resulting ODE numerically to generate trajectories and compute the power spectral density of the emitted radiation. The reported peaks in the PSD for all considered Störmer-type models are direct numerical outputs of this integration; no parameters are fitted to a subset of data and then relabeled as predictions, no self-definitional loops exist between inputs and outputs, and no load-bearing self-citations or uniqueness theorems are invoked to force the result. The derivation chain remains self-contained as a straightforward numerical exploration of the explicitly stated dynamical system.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The model depends on choosing specific functional forms for dissipation and noise, plus numerical parameters whose values are varied but not derived from first principles.

free parameters (2)
  • dissipation coefficient
    Controls strength of velocity-proportional friction; varied across simulations to study effects.
  • stochastic force strength
    Amplitude of random term in the Langevin equation; scanned to explore noise impact.
axioms (2)
  • domain assumption Dissipation is friction-type and strictly proportional to velocity
    Invoked when generalizing the Lorentz force to a Langevin equation.
  • domain assumption Stochastic forces admit a standard additive Langevin representation
    Used to incorporate random effects into the particle equation of motion.

pith-pipeline@v0.9.0 · 5766 in / 1369 out tokens · 51582 ms · 2026-05-18T12:36:42.269423+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

65 extracted references · 65 canonical work pages · 1 internal anchor

  1. [1]

    Classical Dissipative St¨ ormer Problem - CDSP (friction, and no Brownian Motion) 7 C

    Numerical results 7 B. Classical Dissipative St¨ ormer Problem - CDSP (friction, and no Brownian Motion) 7 C. Brownian motion in the Stochastic-Dissipative St¨ ormer Problem - SDSP-dissipation and stochastic effects 8

  2. [2]

    Escape rates 11

    Numerical results 9 D. Escape rates 11

  3. [3]

    The numerical algorithm for the escape rate. 12 IV. Conclusions 12 Acknowledgements 14 Conflict of interest 14 ∗ tiberiu.harko@aira.astro.ro † gabriela.mocanu@aira.astro.ro Keywords 14 References 14 I. INTRODUCTION Magnetic fields are a common occurrence on all as- trophysical, and even cosmological scales, their presence being detected from the galactic or...

  4. [4]

    The solution for two CSP cases was obtained and it was verified by the K1 − 0 method of [49] that one case is not chaotic, and one is chaotic

    Numerical results In our numerical simulations the particle is injected into the dipole field from an initial position ⃗R0 = (X0,Y 0,Z 0), with initial velocities ⃗V0 = (Vx0,V y0,V z0). The solution for two CSP cases was obtained and it was verified by the K1 − 0 method of [49] that one case is not chaotic, and one is chaotic. The solutions differ only by th...

  5. [5]

    Numerical results In the case of the periodic motion, the combined ef- fects of the presence of dissipative and stochastic forces is represented in Figs. 5-7. The effects of the increase of the dissipation coefficient, generally leading to a decrease in the radiation intensity, and periodic motion patterns are compensated by the presence of the stochastic fo...

  6. [6]

    The algorithm H used is as follows: We ran each tra- jectory for L = 150000 timesteps

    The numerical algorithm for the escape rate. The algorithm H used is as follows: We ran each tra- jectory for L = 150000 timesteps. We considered the last L1 = 10000 steps of each trajectory and asked if this ar- ray is well fitted by a straight line distance vs. time. If the R2 of the fit is ≥ 0. 9, then the trajectory is being counted an escape trajectory...

  7. [7]

    Klein, A

    U. Klein, A. Fletcher, Galactic and Intergalactic Mag- netic Fields , Springer International Publishing, Cham, Switzerland, 2015

  8. [8]

    Wielebinski, R

    R. Wielebinski, R. Beck R. (Eds.), Cosmic Magnetic Fields, Lect. Notes Phys. 664, Springer, Berlin, Heidel- berg, 2005

  9. [9]

    L¨ uhr, J

    H. L¨ uhr, J. Wicht, S. A. Gilder, and M. Holschneider (Eds)., Magnetic Fields in the Solar System, Planets, Moons, and Solar Wind Interactions , Astrophysics and Space Science Library, Springer International Publishing AG, Cham, Switzerland, 2018

  10. [10]

    C. C. Finlay et al., Geophysical Journal International 2010 183, 1216

  11. [11]

    St¨ ormer,Arch

    C. St¨ ormer,Arch. Sci. Phys. Nat. 1907, 24, 5

  12. [12]

    St¨ ormer,Arch

    C. St¨ ormer,Arch. Sci. Phys. Nat. 1907, 24, 113

  13. [13]

    St¨ ormer,Arch

    C. St¨ ormer,Arch. Sci. Phys. Nat. 1907, 24, 221

  14. [14]

    St¨ ormer,Astrophysical Journal 1913, 38, 311

    C. St¨ ormer,Astrophysical Journal 1913, 38, 311

  15. [15]

    St”ormer, Geofys ´ ısk Publikationer1921, 1, 269

    C. St”ormer, Geofys ´ ısk Publikationer1921, 1, 269

  16. [16]

    St¨ ormer, Terrestrial Magnetism and Atmospheric Electricity 1917, 22, 97

    C. St¨ ormer, Terrestrial Magnetism and Atmospheric Electricity 1917, 22, 97

  17. [17]

    St¨ ormer,Astrophysica Norvegica 1934, 1, 1

    C. St¨ ormer,Astrophysica Norvegica 1934, 1, 1

  18. [18]

    St¨ ormer,The Polar Aurora, Clarendon Press, Oxford, UK, 1955

    C. St¨ ormer,The Polar Aurora, Clarendon Press, Oxford, UK, 1955

  19. [19]

    Dragt, Rev

    A. Dragt, Rev. Geophys. 1965, 3, 255

  20. [20]

    Dragt, J

    A. Dragt, J. Geophys. Res. 1976, 81, 2327

  21. [21]

    M. A. Almeida, I. C. Moreira, and H. Yoshida, J. Phys. A Math. Gen. 1992, 25, L227

  22. [22]

    F. J. Pav´ on-Carrasco and A. De Santis, Frontiers in Earth Science 2016, 4, 40

  23. [23]

    Schuster and K

    R. Schuster and K. O. Thielheim, J. Phys. A: Math. Gen. 1987, 20, 5511

  24. [24]

    J. E. Howard, M. Hor´ anyi, and G. R. Stewart, Phys. Rev. Lett. 1999, 83, 3993

  25. [25]

    H. R. Dullin, M. Hor´ anyi, J. E. Howard, Physica D 2002, 171, 178

  26. [26]

    et al., Physica D 2004, 197, 242

    I˜ narrea M. et al., Physica D 2004, 197, 242. 15

  27. [27]

    et al., Chaos, Solitons and Fractals 2009, 42, 155

    I˜ narrea M. et al., Chaos, Solitons and Fractals 2009, 42, 155

  28. [28]

    V. Epp, M. A. Masterova, Astrophysics and Space Sci- ence 2014, 353, 473

  29. [29]

    V. Epp, O. N. Pervukhina, Monthly Notices of the Royal Astronomical Society 2018, 474, 5330

  30. [30]

    Pi˜ na, E

    E. Pi˜ na, E. Cort´ es,European Journal of Physics 2016, 37, 065009

  31. [31]

    E. K. Kolesnikov, Geomagnetism and Aeronomy 2017, 57, 137

  32. [32]

    Leghmouche, N

    A. Leghmouche, N. Mebarki, A. Benslama, New Astron- omy 2023, 98, 101931

  33. [33]

    Ershkov, E

    S. Ershkov, E. Prosviryakov, D. Leshchenko, N. Bur- masheva, Mathematical Methods in the Applied Sciences 2023, 46, 19364

  34. [34]

    S. V. Ershkov, J. Appl. Comput. Mech. 2025, in press, https://doi.org/10.22055/jacm.2025.48054.4946

  35. [35]

    Balescu, Statistical dynamics: Matter out of equilib- rium, Imperial College Press, London, 1997

    R. Balescu, Statistical dynamics: Matter out of equilib- rium, Imperial College Press, London, 1997

  36. [36]

    S Lemons, D

    D. S Lemons, D. L Kaufman, IEEE Transactions on Plasma Science 1999, 27, 1288

  37. [37]

    Balescu, Aspects of Anomalous Transport in Plasmas , Series in Plasma Physics, Institute of Physics Publishing, Bristol and Philadelphia, 2005

    R. Balescu, Aspects of Anomalous Transport in Plasmas , Series in Plasma Physics, Institute of Physics Publishing, Bristol and Philadelphia, 2005

  38. [38]

    W. T. Coffey, Yu. P. Kalmykov, J. T. Waldron, The Langevin Equation, with Applications to Stochastic Prob- lems in Physics, Chemistry, and Electrical Engineering , World Scientific, Singapore, 2005

  39. [39]

    Harko, G

    T. Harko, G. R. Mocanu, Monthly Notices of the Royal Astronomical Society 2012, 421, 3102

  40. [40]

    Harko, C

    T. Harko, C. S. Leung, and G. R. Mocanu, Eur. Phys. J. C 2014, 74, 2900

  41. [41]

    Mocanu, N

    G. Mocanu, N. Magyar, A. Pardi, A. Marcu, Monthly No- tices of the Royal Astronomical Society 2014, 439, 3790

  42. [42]

    Harko, G

    T. Harko, G. Mocanu, Eur. Phys. J. C 2016, 76, 160

  43. [43]

    G. R. Mocanu, Romanian Astronomical Journal 2019, 29, 41

  44. [44]

    G. R. Mocanu, Rom. Rep. Phys. 2020, 72, 105

  45. [45]

    G. R. Mocanu, IEEE Transactions on Plasma Science 2021, 49, 1527

  46. [46]

    I. N. Toptygin, G. D. Fleishman, Astrophys. Space Sci- ence 1987, 132, 213

  47. [47]

    Reville, J

    B. Reville, J. G. Kirk, Astrophys. J. 2010, 724, 1283

  48. [48]

    S. R. Kelner, F. A. Aharonian, D. Khangulyan, Astro- phys. J. 2013, 774, 61

  49. [49]

    Dil˜ ao, R

    R. Dil˜ ao, R. Alves-Pires R., Chaos in the St¨ ormer prob- lem, in Progress in Nonlinear Differential Equations and Their Applications , Vol. 75, 175-194, Birkh¨ auser Verlag, Basel, Switzerland, 2007

  50. [50]

    L. D. Landau, E. M. Lifshitz, The classical theory of fields, Oxford, Pergamon Press, 1971

  51. [51]

    Huang, Lectures on Statistical Physics and Protein Folding, World Scientific Publishing Company Pte

    K. Huang, Lectures on Statistical Physics and Protein Folding, World Scientific Publishing Company Pte. Ltd., Singapore, 2005

  52. [52]

    Giles, A

    L. Giles, A. Keller, S. Heinrich, H. Niederreiter (Eds) ., Monte Carlo and Quasi-Monte Carlo Methods , 343-358, Springer Berlin Heidelberg, Berlin, Heidelberg, 2006

  53. [53]

    J. Larsen, Correlation Functions and Power Spectra, Sec- tion for Cognitive Systems, Informatics and Mathemati- cal Modelling , Technical University of Denmark, Copen- hagen, Eighth Edition, 2009

  54. [54]

    S. V. Vaseghi, Advanced Digital Signal Processing and Noise Reduction , Second Edition, John Wiley & Sons Ltd, New York, Singapore, 2000

  55. [55]

    G. A. Gottwald, I. Melbourne, Proceedings of the Royal Society A 2004, 460, 603

  56. [56]

    Chandrasekhar, Reviews of Modern Physics 1943, 15, 1

    S. Chandrasekhar, Reviews of Modern Physics 1943, 15, 1

  57. [57]

    A. A. Tikhonov, K. G. Petrov, Cosmic Research 2002, 40, 203

  58. [58]

    W. D. Gonzalez et al., J. Geophys. Res. 1994, 99, 5771

  59. [59]

    Zhao et al., Solar Physics 2021, 296, 66

    M.-X. Zhao et al., Solar Physics 2021, 296, 66

  60. [60]

    Yokoyama, H

    T. Yokoyama, H. Shinagawa, H. Jin, J. Geophys. Res. Space Physics 2014, 119, 10474

  61. [61]

    Walach, A

    M.-T. Walach, A. Grocott, S. E. Milan, Geophysical Re- search: Space Physics 2021, 126, e2020JA028512

  62. [62]

    K. S. Thorne, Astrophysical Journal Supplement 1963, 8, 1

  63. [63]

    D. B. Papadopoulos, I. Contopoulos, K. D. Kokkotas, N. Stergioulas, General Relativity and Gravitation 2015, 47, 49

  64. [64]

    T. M. Burinskaya, M. M. Shevelev, Plasma Physics Re- ports 2016, 42, 929

  65. [65]

    T. M. Burinskaya, M. M. Shevelev, Plasma Physics Re- ports 2017, 43, 910