pith. machine review for the scientific record. sign in

arxiv: 2604.20963 · v1 · submitted 2026-04-22 · ⚛️ physics.plasm-ph · astro-ph.HE

Recognition: unknown

Turbulence Mode Decomposition and Anisotropy in Magnetically Dominated Collisionless Plasmas

Samuel T. Sebastian , Siyao Xu , Yue Hu , Luca Comisso , Saikat Das , Joonas N\"attil\"a
Authors on Pith no claims yet

Pith reviewed 2026-05-09 22:42 UTC · model grok-4.3

classification ⚛️ physics.plasm-ph astro-ph.HE
keywords relativistic plasma turbulencecollisionless plasmasmode decompositionturbulence anisotropyAlfven modeskinetic simulationsmagnetically dominated plasmas
0
0 comments X

The pith

Alfvén and slow modes in relativistic collisionless plasma turbulence follow Goldreich-Sridhar anisotropy scaling while fast modes are isotropic with higher kinetic energy.

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

This paper uses three-dimensional fully kinetic simulations to examine how turbulence breaks down into different modes in relativistic, magnetically dominated collisionless plasmas. It adapts an existing decomposition technique to separate the turbulence into Alfvén, fast, and slow modes and measures their properties. Alfvén and slow modes remain anisotropic and follow the expected scaling law, fast modes stay isotropic, and fast modes carry a larger share of kinetic energy than in non-relativistic cases. These results indicate stronger interactions between Alfvén and fast modes under relativistic collisionless conditions. The study also reports weaker scale dependence for dynamic alignment and shows that thermal fluctuations near kinetic scales flatten the velocity structure function while reducing anisotropy.

Core claim

By extending the Cho and Lazarian decomposition to relativistic collisionless plasmas in 3D fully kinetic simulations, the work shows Alfvén and slow modes are anisotropic with Goldreich-Sridhar scaling, fast modes are isotropic, the kinetic energy fraction in fast modes exceeds non-relativistic MHD values suggesting stronger Alfvén-fast coupling, dynamic alignment exhibits weaker scale dependence, and thermal fluctuations flatten the turbulent velocity structure function while weakening anisotropy and alignment near kinetic scales.

What carries the argument

The Cho and Lazarian decomposition method extended to relativistic collisionless plasmas, which separates turbulence into Alfvén, fast, and slow modes to quantify their anisotropy, energy fractions, and scale-dependent alignment.

Load-bearing premise

The Cho and Lazarian decomposition technique developed for non-relativistic MHD turbulence applies to relativistic collisionless plasmas after only minor adjustments and still correctly identifies the physical modes.

What would settle it

A simulation result in which the kinetic energy fraction of fast modes equals the non-relativistic MHD value, or in which fast modes display clear anisotropy instead of isotropy, would contradict the central findings.

Figures

Figures reproduced from arXiv: 2604.20963 by Joonas N\"attil\"a, Luca Comisso, Saikat Das, Samuel T. Sebastian, Siyao Xu, Yue Hu.

Figure 1
Figure 1. Figure 1: Velocity field (normalized by the rms velocity vrms) in the mid-z plane from the PIC (upper panel) and MHD (lower panel) simulations. to model the plasma. When the turbulence is fully devel￾oped, heating causes increase in plasma β and the plasma skin depth de, with β ≈ 0.3 and de ≈ 5 cells corresponding to the data analyzed in this work. 2.2. MHD Simulation The 3D isothermal MHD turbulence simulation anal… view at source ↗
Figure 2
Figure 2. Figure 2: SFv(r∥) and SFv(r⊥) (normalized by v 2 rms) measured in the PIC (upper panel) and MHD (lower panel) simulations. The GS95 anisotropic scalings are indicated by the dashed and dash￾dotted lines. The vertical dotted lines indicate de in the PIC simula￾tion and the numerical dissipation scale ld in the MHD simulation. where v is the fluid velocity, cs is the sound speed, ρ is the fluid density, B is the magne… view at source ↗
Figure 3
Figure 3. Figure 3: SFv(r∥) and SFv(r⊥) of different turbulence modes for the PIC (left) and MHD (right) simulations. In [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Kinetic energy fractions of different turbulence modes measured from the PIC and MHD simulations. fast, and slow modes measured in the PIC simulation, and fA ≈ 0.60, ff ≈ 0.11, and fs ≈ 0.29 for the MHD sim￾ulation. With the solenoidal driving adopted in the MHD simulation, the small ff indicates inefficient generation of fast modes due to the weak coupling between Alfven and fast ´ modes (Cho & Lazarian 2… view at source ↗
Figure 5
Figure 5. Figure 5: SFv(r∥, rξ, rλ) for the PIC (left panel) and MHD (right panel) simulations. on MHD simulations of Alfvenic turbulence with resolutions ´ up to 40963 (Beresnyak 2015). In the regime of relativistic turbulence in collisionless plasmas, for a systematic compar￾ison with the theory of dynamic alignment, higher-resolution PIC simulations of Alfvenic turbulence will be needed for the ´ convergence study. In Appe… view at source ↗
Figure 6
Figure 6. Figure 6: is seen. 4. CONCLUSIONS We perform turbulence mode decomposition and investi￾gate the turbulence anisotropy of relativistic turbulence in strongly magnetized collisionless pair plasma with the 3D PIC simulation. We focus on the dynamics of turbulence and thus the turbulent velocity fluctuations. We find the SF of total velocity fluctuations are consis￾tent with the GS95 anisotropic scaling in the inertial … view at source ↗
Figure 7
Figure 7. Figure 7: Alignment angle θv,b(r) (left panel) and alignment slope α (right panel) measured in PIC simulations with different spatial resolu￾tions. APPENDIX We perform a lower-resolution 5123 PIC simulation with the same initial conditions as detailed in Section 2. We analyze the simulation when turbulence is fully developed, corresponding to β ≈ 0.6, and de ≈ 5 cells [PITH_FULL_IMAGE:figures/full_fig_p009_7.png] view at source ↗
read the original abstract

We use the 3D fully kinetic simulation to study different turbulence modes and turbulence anisotropy of relativistic turbulence in magnetically dominated collisionless plasmas. We extend the method developed by Cho & Lazarian (2002) for decomposing non-relativistic magnetohydrodynamic (MHD) turbulence into Alfv\'en, fast, and slow modes to the regime of collisionless plasmas. We find that Alfv\'en and slow modes are anisotropic, following the Goldreich & Sridhar (1995) scaling, while fast modes are isotropic. We observe a larger kinetic energy fraction of fast modes compared to that in the non-relativistic MHD turbulence, suggesting a stronger coupling of Alfv\'en and fast modes in relativistic magnetized turbulence in collisionless plasmas. We further examine the dynamic alignment and find a weaker scale dependence of the alignment angle than previously proposed. The dominant thermal fluctuations in the kinetic range can cause flattening of the turbulent velocity structure function and weakening of the turbulence anisotropy and dynamic alignment near the kinetic scales.

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

No significant circularity; results extracted from direct simulation diagnostics

full rationale

The paper extends the Cho & Lazarian (2002) decomposition to relativistic collisionless plasmas and applies it to 3D kinetic simulation outputs to measure mode energies, GS95 anisotropy scalings, and dynamic alignment. No equations reduce the reported kinetic energy fractions or anisotropy measurements to quantities defined by fitted parameters or self-referential inputs. The central claims are obtained from Fourier-space projections and structure functions computed on the simulation data itself, with the extension treated as a methodological choice rather than a derivation that collapses to its own assumptions by construction. This keeps the analysis self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the assumption that the non-relativistic mode decomposition remains valid in the relativistic collisionless limit and on the interpretation of simulation outputs as cleanly separated Alfvén, fast, and slow modes.

axioms (1)
  • domain assumption The Cho & Lazarian (2002) decomposition method applies to relativistic collisionless plasmas
    The paper states that the method is extended but does not detail relativistic or kinetic corrections required for the decomposition.

pith-pipeline@v0.9.0 · 5497 in / 1479 out tokens · 31996 ms · 2026-05-09T22:42:29.321781+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

39 extracted references · 39 canonical work pages

  1. [1]

    , keywords =

    Armstrong, J. W., Rickett, B. J., & Spangler, S. R. 1995, ApJ, 443, 209, doi: 10.1086/175515

    work page doi:10.1086/175515 1995

  2. [2]

    2025, The spectrum of magnetized turbulence in the interstellar medium

    Bhattacharjee, A. 2025, The spectrum of magnetized turbulence in the interstellar medium. https://doi.org/10.1038/s41550-025-02551-5

    work page doi:10.1038/s41550-025-02551-5 2025

  3. [3]

    2012 , month = apr, journal =

    Beresnyak, A. 2012, Monthly Notices of the Royal Astronomical Society, 422, 3495, doi: 10.1111/j.1365-2966.2012.20859.x —. 2014, The Astrophysical Journal Letters, 784, L20, doi: 10.1088/2041-8205/784/2/L20 —. 2015, The Astrophysical Journal Letters, 801, L9, doi: 10.1088/2041-8205/801/1/L9

    work page doi:10.1111/j.1365-2966.2012.20859.x 2012

  4. [4]

    2006, PhRvL, 96, 115002, doi: 10.1103/PhysRevLett.96.115002

    Boldyrev, S. 2006, Phys. Rev. Lett., 96, 115002, doi: 10.1103/PhysRevLett.96.115002

    work page doi:10.1103/physrevlett.96.115002 2006

  5. [5]

    Chen, C. H. K., Bale, S. D., Bonnell, J. W., et al. 2020, The Astrophysical Journal Supplement Series, 246, 53, doi: 10.3847/1538-4365/ab60a3

    work page doi:10.3847/1538-4365/ab60a3 2020

  6. [6]

    2021, ApJL, 923, L13, doi: 10.3847/2041-8213/ac3afa

    Chernoglazov, A., Ripperda, B., & Philippov, A. 2021, ApJL, 923, L13, doi: 10.3847/2041-8213/ac3afa

    work page doi:10.3847/2041-8213/ac3afa 2021

  7. [7]

    2002, Phys

    Cho, J., & Lazarian, A. 2002, Phys. Rev. Lett., 88, 245001, doi: 10.1103/PhysRevLett.88.245001 —. 2003, Monthly Notices of the Royal Astronomical Society, 345, 325, doi: 10.1046/j.1365-8711.2003.06941.x

    work page doi:10.1103/physrevlett.88.245001 2002

  8. [8]

    Cho, J., Lazarian, A., & Vishniac, E. T. 2002, The Astrophysical Journal, 564, 291, doi: 10.1086/324186

    work page doi:10.1086/324186 2002

  9. [9]

    Cho, J., & Vishniac, E. T. 2000, The Astrophysical Journal, 539, 273, doi: 10.1086/309213

    work page doi:10.1086/309213 2000

  10. [10]

    R., & Muzio, M

    Comisso, L., Farrar, G. R., & Muzio, M. S. 2024, ApJL, 977, L18, doi: 10.3847/2041-8213/ad955f

    work page doi:10.3847/2041-8213/ad955f 2024

  11. [11]

    2018, Phys

    Comisso, L., & Sironi, L. 2018, Phys. Rev. Lett., 121, 255101, doi: 10.1103/PhysRevLett.121.255101 —. 2019, The Astrophysical Journal, 886, 122, doi: 10.3847/1538-4357/ab4c33

    work page doi:10.1103/physrevlett.121.255101 2018

  12. [12]

    2021, PhRvL, 127, 255102, doi: 10.1103/PhysRevLett.127.255102

    Comisso, L., & Sironi, L. 2021, PhRvL, 127, 255102, doi: 10.1103/PhysRevLett.127.255102

    work page doi:10.1103/physrevlett.127.255102 2021

  13. [13]

    Crutcher, R. M. 2012, Annual Review of Astronomy and Astrophysics, 50, 29, doi: 10.1146/annurev-astro-081811-125514 Del Zanna, L., Bucciantini, N., & Landi, S. 2025, A&A, 702, A171, doi: 10.1051/0004-6361/202556255

    work page doi:10.1146/annurev-astro-081811-125514 2012

  14. [14]

    2022, Science Advances, 8, eabn7627, doi: 10.1126/sciadv.abn7627

    Dong, C., Wang, L., Huang, Y .-M., et al. 2022, Science Advances, 8, eabn7627, doi: 10.1126/sciadv.abn7627

    work page doi:10.1126/sciadv.abn7627 2022

  15. [15]

    Eyink, G., Bandak, D., Goldenfeld, N., & Mailybaev, A. A. 2021, Dissipation-Range Fluid Turbulence and Thermal Noise. https://arxiv.org/abs/2107.13954

    work page arXiv 2021

  16. [16]

    1949, Physical Review, 75, 1169, doi: 10.1103/PhysRev.75.1169 21

    Fermi, E. 1949, Phys. Rev., 75, 1169, doi: 10.1103/PhysRev.75.1169

    work page doi:10.1103/physrev.75.1169 1949

  17. [17]

    2: Strong alfvenic turbulence

    Goldreich, P., & Sridhar, H. 1995, The Astrophysical Journal, 438, 763, doi: 10.1086/175121

    work page doi:10.1086/175121 1995

  18. [18]

    Superdiffusion of cosmic rays in compressible magnetized turbulence.MNRAS2022,512, 2111–2124, [arXiv:astro-ph.GA/2111.15066]

    Hu, Y ., Lazarian, A., & Xu, S. 2022, MNRAS, 512, 2111, doi: 10.1093/mnras/stac319 8

    work page doi:10.1093/mnras/stac319 2022

  19. [19]

    Damping of MHD turbulence in a partially ionized medium.MNRAS 2024,527, 3945–3961, [arXiv:astro-ph.GA/2306.10010]

    Hu, Y ., Xu, S., Arzamasskiy, L., Stone, J. M., & Lazarian, A. 2024, MNRAS, 527, 3945, doi: 10.1093/mnras/stad3493

    work page doi:10.1093/mnras/stad3493 2024

  20. [20]

    2021, ApJ, 911, 37, doi: 10.3847/1538-4357/abea18

    Hu, Y ., Xu, S., & Lazarian, A. 2021, ApJ, 911, 37, doi: 10.3847/1538-4357/abea18

    work page doi:10.3847/1538-4357/abea18 2021

  21. [21]

    M., & Hopkins, P

    Hu, Y ., Xu, S., Lazarian, A., Stone, J. M., & Hopkins, P. F. 2025, The Astrophysical Journal, 994, 142, doi: 10.3847/1538-4357/ae1127

    work page doi:10.3847/1538-4357/ae1127 2025

  22. [22]

    , keywords =

    Inoue, T., Asano, K., & Ioka, K. 2011, The Astrophysical Journal, 734, 77, doi: 10.1088/0004-637X/734/2/77

    work page doi:10.1088/0004-637x/734/2/77 2011

  23. [23]

    2010, The Astrophysical Journal, 720, 742, doi: 10.1088/0004-637X/720/1/742

    Kowal, G., & Lazarian, A. 2010, The Astrophysical Journal, 720, 742, doi: 10.1088/0004-637X/720/1/742

    work page doi:10.1088/0004-637x/720/1/742 2010

  24. [24]

    Larson, R. B. 1981, Monthly Notices of the Royal Astronomical Society, 194, 809, doi: 10.1093/mnras/194.4.809

    work page doi:10.1093/mnras/194.4.809 1981

  25. [25]

    Lazarian, A., & Vishniac, E. T. 1999, The Astrophysical Journal, 517, 700, doi: 10.1086/307174

    work page doi:10.1086/307174 1999

  26. [26]

    2001, The Astrophysical Journal, 562, 279, doi: 10.1086/323470

    Lithwick, Y ., & Goldreich, P. 2001, The Astrophysical Journal, 562, 279, doi: 10.1086/323470

    work page doi:10.1086/323470 2001

  27. [27]

    2026, Acta Mechanica Sinica, 42, 325230, doi: 10.1007/s10409-025-25230-x

    Liu, X., Zhao, C., Zhou, Z., Cheng, W., & Si, T. 2026, Acta Mechanica Sinica, 42, 325230, doi: 10.1007/s10409-025-25230-x

    work page doi:10.1007/s10409-025-25230-x 2026

  28. [28]

    2023, Advances in Aerodynamics, 5, 3, doi: 10.1186/s42774-022-00134-w

    Ma, Q., Yang, C., Chen, S., Feng, K., & Zhang, J. 2023, Advances in Aerodynamics, 5, 3, doi: 10.1186/s42774-022-00134-w

    work page doi:10.1186/s42774-022-00134-w 2023

  29. [29]

    2001, Astrophys

    Maron, J., & Goldreich, P. 2001, Astrophys. J., 554, 1175, doi: 10.1086/321413

    work page doi:10.1086/321413 2001

  30. [30]

    2006, Phys

    Mason, J., Cattaneo, F., & Boldyrev, S. 2006, Phys. Rev. Lett., 97, 255002, doi: 10.1103/PhysRevLett.97.255002

    work page doi:10.1103/physrevlett.97.255002 2006

  31. [31]

    F., & Ostriker, E

    McKee, C. F., & Ostriker, E. C. 2007, Annual Review of Astronomy and Astrophysics, 45, 565, doi: 10.1146/annurev.astro.45.051806.110602 N¨attil¨a, J. 2022, A&A, 664, A68, doi: 10.1051/0004-6361/201937402 —. 2024, Nature Communications, 15, 7026, doi: 10.1038/s41467-024-51257-1 N¨attil¨a, J., & Beloborodov, A. M. 2021, The Astrophysical Journal, 921, 87, d...

    work page doi:10.1146/annurev.astro.45.051806.110602 2007

  32. [32]

    T., & Comisso, L

    Sebastian, S. T., & Comisso, L. 2025, ApJL, 994, L1, doi: 10.3847/2041-8213/ae1696

    work page doi:10.3847/2041-8213/ae1696 2025

  33. [33]

    M., Mullen, P

    Stone, J. M., Mullen, P. D., Fielding, D., et al. 2024, arXiv e-prints, arXiv:2409.16053, doi: 10.48550/arXiv.2409.16053

    work page doi:10.48550/arxiv.2409.16053 2024

  34. [34]

    2016, The Astrophysical Journal Letters, 831, L11, doi: 10.3847/2041-8205/831/2/L11 —

    Takamoto, M., & Lazarian, A. 2016, The Astrophysical Journal Letters, 831, L11, doi: 10.3847/2041-8205/831/2/L11 —. 2017, Monthly Notices of the Royal Astronomical Society, 472, 4542, doi: 10.1093/mnras/stx2292

    work page doi:10.3847/2041-8205/831/2/l11 2016

  35. [35]

    M., Howes, G

    TenBarge, J. M., Howes, G. G., Dorland, W., & Hammett, G. W. 2014, Computer Physics Communications, 185, 578, doi: 10.1016/j.cpc.2013.10.022

    work page doi:10.1016/j.cpc.2013.10.022 2014

  36. [36]

    2023 Review of Solar Energetic Particle Prediction Models

    Whitman, K., Egeland, R., Richardson, I. G., et al. 2023, Advances in Space Research, 72, 5161, doi: https://doi.org/10.1016/j.asr.2022.08.006

    work page doi:10.1016/j.asr.2022.08.006 2023

  37. [37]

    Cosmic-Ray Propagation: Nonlinear Diffusion Parallel and Perpendicular to Mean Magnetic Field.ApJ 2008,673, 942–953, [arXiv:astro-ph/0710.2617]

    Yan, H., & Lazarian, A. 2008, The Astrophysical Journal, 673, 942, doi: 10.1086/524771

    work page doi:10.1086/524771 2008

  38. [38]

    A., Werner, G

    Zhdankin, V ., Uzdensky, D. A., Werner, G. R., & Begelman, M. C. 2018, MNRAS, 474, 2514, doi: 10.1093/mnras/stx2883

    work page doi:10.1093/mnras/stx2883 2018

  39. [39]

    Zrake, J., & MacFadyen, A. I. 2012, The Astrophysical Journal, 744, 32, doi: 10.1088/0004-637X/744/1/32 9 Figure 7.Alignment angleθ v,b(r)(left panel) and alignment slopeα(right panel) measured in PIC simulations with different spatial resolu- tions. APPENDIX We perform a lower-resolution512 3 PIC simulation with the same initial conditions as detailed in...

    work page doi:10.1088/0004-637x/744/1/32 2012