pith. machine review for the scientific record. sign in

arxiv: 2605.11305 · v1 · submitted 2026-05-11 · ⚛️ physics.plasm-ph · astro-ph.SR· physics.space-ph

Recognition: 2 theorem links

· Lean Theorem

Dynamic Alignment: A Fragile Survival Effect

Amir Jafari

Pith reviewed 2026-05-13 00:51 UTC · model grok-4.3

classification ⚛️ physics.plasm-ph astro-ph.SRphysics.space-ph
keywords dynamic alignmentMHD turbulenceElsasser incrementssurvival effectsolar windturbulence cascadeangular correlation
0
0 comments X

The pith

Dynamic alignment in magnetohydrodynamic turbulence is primarily a selective survival effect favoring intense small-angle events.

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

The paper argues that what is measured as dynamic alignment in MHD turbulence mainly reflects the longer survival of intense events with small angles between Elsasser increments. In simulations, the typical angle stays close to random while only the strongest events show much smaller angles. Tests that shuffle the data or track how events move between amplitude-angle sectors confirm a genuine negative covariance between amplitude and misalignment. This pattern holds in solar wind observations as well. If correct, it means standard diagnostics are not capturing a volume-filling cascade ordering but a bias toward persistent intense aligned fluctuations.

Core claim

These results indicate that dynamic alignment, as measured by conventional weighted diagnostics, is best understood as selective sampling of longer-lived intense small-angle events, not as a cascade-wide alignment of typical MHD fluctuations.

What carries the argument

A finite-time state-retention analysis showing that high-amplitude large-angle events depart their sector faster than high-amplitude small-angle events, supported by shuffled-null tests that isolate the negative covariance between amplitude and angular misalignment.

Load-bearing premise

The shuffled-null tests and finite-time state-retention analysis isolate a genuine survival mechanism without confounding effects from other unmeasured properties of the turbulence or sampling.

What would settle it

A simulation in which high-amplitude events with large angles persist in their amplitude-angle sector at the same rate as those with small angles would falsify the proposed survival mechanism.

Figures

Figures reproduced from arXiv: 2605.11305 by Amir Jafari.

Figure 1
Figure 1. Figure 1: shows the covariance mechanism in the main DNS ensemble. The amplitude-weighted proxy is smaller than the unweighted sine average at the plotted separa￾tions, and shuffling the amplitudes relative to the angles removes the difference. The covariance is negative across the same range and collapses to zero under the shuffled null. This is the first indication that strong apparent alignment is not a volume-fi… view at source ↗
Figure 2
Figure 2. Figure 2: previews one of the central numerical re￾sults of the paper. The unweighted mean angle stays only modestly below the random three-dimensional ex￾pectation and shows no simple monotone decrease over the plotted separations, indicating at most weak aver￾age alignment. By contrast, much smaller angles appear when one conditions on large Els¨asser-increment am￾plitudes, whereas conditioning on large current de… view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4 [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5 [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6 [PITH_FULL_IMAGE:figures/full_fig_p012_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7 [PITH_FULL_IMAGE:figures/full_fig_p012_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8 [PITH_FULL_IMAGE:figures/full_fig_p014_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9 [PITH_FULL_IMAGE:figures/full_fig_p014_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10 [PITH_FULL_IMAGE:figures/full_fig_p023_10.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12 [PITH_FULL_IMAGE:figures/full_fig_p024_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: summarizes the result over the fifteen 3203 subvolumes. The weighted total-variation residuals are 0.065±0.004, 0.046±0.002, 0.039±0.001, 0.054±0.003 for the triplets (128, 96, 64), (96, 64, 48), (64, 48, 32), and (128, 64, 32), respectively, with uncertainties given as SEM across cubes. Thus the composed transition through an intermediate scale reproduces the directly measured transition to within a few … view at source ↗
Figure 14
Figure 14. Figure 14: FIG. 14 [PITH_FULL_IMAGE:figures/full_fig_p035_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: FIG. 15 [PITH_FULL_IMAGE:figures/full_fig_p035_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: FIG. 16 [PITH_FULL_IMAGE:figures/full_fig_p036_16.png] view at source ↗
Figure 17
Figure 17. Figure 17: FIG. 17 [PITH_FULL_IMAGE:figures/full_fig_p037_17.png] view at source ↗
Figure 18
Figure 18. Figure 18: FIG. 18 [PITH_FULL_IMAGE:figures/full_fig_p038_18.png] view at source ↗
Figure 19
Figure 19. Figure 19: FIG. 19 [PITH_FULL_IMAGE:figures/full_fig_p039_19.png] view at source ↗
Figure 20
Figure 20. Figure 20: FIG. 20 [PITH_FULL_IMAGE:figures/full_fig_p039_20.png] view at source ↗
read the original abstract

Dynamic alignment in magnetohydrodynamic (MHD) turbulence is usually interpreted as a cascade-wide tendency of Elsasser increments to become increasingly collinear at smaller scales. We argue instead that the standard measurements mainly detect a conditional survival effect of intense events. In high-resolution Johns Hopkins MHD simulations, the typical folded Elsasser-increment angle remains only modestly below the random-orientation baseline and shows no evidence for a rigid, monotone, volume-filling ordering of the cascade. Much smaller angles appear primarily in the strongest Elsasser-amplitude events, while conditioning on current density leaves the angle close to its unweighted behavior. Shuffled-null tests show that this reduction is caused by a genuine negative covariance between event amplitude and angular misalignment, not by weighting alone. Cross-scale angular correlations are measurable but decaying, indicating partial and non-rigid persistence of the local alignment field. A finite-time state-retention test directly supports the proposed mechanism: high-amplitude large-angle events leave their amplitude--angle sector faster than high-amplitude small-angle events, while incoming transitions continually replenish the large-angle sector. NASA Wind solar-wind observations show the same angle--amplitude hierarchy and negative covariance in Taylor-sampled Elsasser increments. These results indicate that dynamic alignment, as measured by conventional weighted diagnostics, is best understood as selective sampling of longer-lived intense small-angle events, not as a cascade-wide alignment of typical MHD fluctuations.

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 / 2 minor

Summary. The paper claims that dynamic alignment in MHD turbulence, as conventionally measured by weighted diagnostics on Elsasser increments, is not evidence of a cascade-wide, volume-filling tendency toward collinearity at small scales. Instead, it primarily reflects a selective survival effect in which intense, small-angle events persist longer than large-angle ones of comparable amplitude. This is supported by Johns Hopkins high-resolution simulations showing only modest deviations of typical angles from random baselines, a negative amplitude-angle covariance confirmed by shuffled-null tests, decaying cross-scale correlations, a finite-time state-retention analysis demonstrating differential exit rates from amplitude-angle sectors, and consistent patterns in NASA Wind solar-wind data. Conditioning on current density is reported to leave the angle distribution nearly unchanged.

Significance. If the central interpretation holds, the result would substantially revise the understanding of dynamic alignment away from a fundamental cascade property toward a fragile, conditional sampling bias affecting only the strongest events. This carries implications for theories of MHD energy transfer and dissipation. The manuscript earns credit for its use of high-resolution simulations, independent observational confirmation, multiple null tests (shuffled, cross-scale, retention), and direct comparison to external datasets without reliance on fitted parameters or circular definitions.

major comments (2)
  1. [finite-time state-retention test] Finite-time state-retention analysis: the test establishes faster departure of high-amplitude large-angle events and replenishment of the large-angle sector, but does not condition on local current density |j|, vorticity, or dissipation rate. Although the abstract states that conditioning on |j| leaves the angle distribution nearly unchanged, this does not demonstrate that |j| (or correlated quantities) is independent of retention time; the survival effect could still be driven by these confounders rather than angle itself.
  2. [results and discussion of retention and covariance tests] Interpretation of the mechanism: the negative covariance and differential retention are presented as evidence that alignment drives selective survival, yet the tests do not isolate angle from other properties known to correlate with both amplitude and geometry in MHD turbulence. A direct test conditioning the retention analysis on |j| (or dissipation) would be required to support the causal claim that alignment itself produces the observed lifetime difference.
minor comments (2)
  1. Clarify the precise definition and folding procedure for the 'Elsasser-increment angle' used in the weighted and unweighted diagnostics, especially how it differs from the standard angle between increments.
  2. The abstract refers to 'folded' angles and 'state-retention'; ensure these terms are defined at first use in the main text with explicit formulas or pseudocode for reproducibility.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive report and positive assessment of the manuscript's methodology and implications. We address the two major comments on the finite-time state-retention analysis below, focusing on the potential role of current density |j| as a confounder.

read point-by-point responses

  1. Referee: [finite-time state-retention test] Finite-time state-retention analysis: the test establishes faster departure of high-amplitude large-angle events and replenishment of the large-angle sector, but does not condition on local current density |j|, vorticity, or dissipation rate. Although the abstract states that conditioning on |j| leaves the angle distribution nearly unchanged, this does not demonstrate that |j| (or correlated quantities) is independent of retention time; the survival effect could still be driven by these confounders rather than angle itself.

    Authors: We agree that the retention test does not explicitly condition on |j|. However, the analysis is already restricted to high-amplitude events, which are strongly correlated with elevated |j|. Within this population, the angle distribution remains nearly unchanged upon conditioning on |j|, indicating that |j| does not preferentially select for small angles. The differential retention rates are therefore observed as a function of angle at comparable amplitudes (and thus comparable |j| levels). We will add a clarifying paragraph in the revised manuscript explaining this reasoning and why |j| is unlikely to be the dominant driver. revision: partial

  2. Referee: [results and discussion of retention and covariance tests] Interpretation of the mechanism: the negative covariance and differential retention are presented as evidence that alignment drives selective survival, yet the tests do not isolate angle from other properties known to correlate with both amplitude and geometry in MHD turbulence. A direct test conditioning the retention analysis on |j| (or dissipation) would be required to support the causal claim that alignment itself produces the observed lifetime difference.

    Authors: The shuffled-null tests already isolate the amplitude-angle covariance from weighting artifacts. The retention test then demonstrates longer persistence of small-angle events at fixed high amplitude. Given the demonstrated insensitivity of the angle distribution to |j| conditioning, we maintain that the survival effect is attributable to angle rather than |j| or dissipation. A full re-analysis conditioning retention times on |j| would be computationally intensive and, in our view, unlikely to change the conclusions. We will expand the discussion to explicitly address potential confounders and the limits of the current isolation. revision: partial

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper's claims rest on direct empirical measurements from external Johns Hopkins MHD simulations and NASA Wind observations, including shuffled-null tests for covariance and finite-time state-retention analysis for differential survival rates. These steps compare observed angle-amplitude hierarchies and transition rates against baselines without reducing to self-defined quantities, fitted parameters renamed as predictions, or load-bearing self-citations. The reinterpretation of dynamic alignment as selective sampling follows from the data patterns rather than by construction, and conditioning statements (e.g., on current density) are presented as independent checks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The paper introduces no new free parameters or invented entities. It reanalyzes existing data under standard MHD and statistical frameworks.

axioms (2)
  • domain assumption Turbulence obeys the MHD equations and Elsasser increments are well-defined
    Underpins all angle and amplitude measurements in both simulations and observations.
  • standard math Shuffling data preserves amplitude distributions while destroying correlations
    Central to the null test establishing genuine negative covariance between amplitude and misalignment.

pith-pipeline@v0.9.0 · 5545 in / 1290 out tokens · 78630 ms · 2026-05-13T00:51:59.077132+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.

Reference graph

Works this paper leans on

30 extracted references · 30 canonical work pages

  1. [1]

    Larger-cube and component-weighting checks As an additional consistency check, we repeated the one-scale angle diagnostics on a substantially larger cutout from the same MHD dataset, using a single 448 3 cube att= 57 with spatial rangex, y, z= 289–736. We do not include this run in the main 15-cube ensemble averages, since it is a separate larger-volume r...

    work page

  2. [2]

    These diagnostics use only the velocity and magnetic fields; pressure is not required

    Numerical implementation of angle and cross-scale correlation diagnostics This appendix describes the numerical implementation of the local alignment-angle and cross-scale correlation diagnostics. These diagnostics use only the velocity and magnetic fields; pressure is not required. 25 For each subvolume we form the Els¨ asser fields z± =u±B. For each sep...

    work page

  3. [3]

    transition

    Chapman–Kolmogorov consistency check As a consistency check on the Markov-in-scale reduc- tion, we tested the Chapman–Kolmogorov relation for the signed cosine field cr(x) =q r(x). Here “transition” is meant in an empirical scale-space sense, not as a physical-time trajectory of a material ob- ject. At each sampled base point and separation direc- tion, w...

    work page

  4. [4]

    For this part of the calculation we use a smooth compact-support separable kernel, not a Gaussian filter

    Compact-support filtering for the geometric diagnostic The geometric diagnostic requires coarse-grained Els¨ asser fields and their effective accelerations. For this part of the calculation we use a smooth compact-support separable kernel, not a Gaussian filter. In one dimension the unnormalised kernel is ϕ(s) =    exp − 1 1−s 2 ,|s|<1, 0,|s| ≥1. For a...

    work page

  5. [5]

    The subscale stress tensor is computed componentwise as τ ± ℓ,ij = z± i z∓ j ℓ −z ± i,ℓz∓ j,ℓ, where the overbar denotes compact-support filtering

    Discrete Els¨ asser acceleration and control coefficient For each scaleℓ, we compute filtered fieldsz ± ℓ and Πℓ. The subscale stress tensor is computed componentwise as τ ± ℓ,ij = z± i z∓ j ℓ −z ± i,ℓz∓ j,ℓ, where the overbar denotes compact-support filtering. The resolved nonideal contribution is N ± ℓ =−∇ ·τ ± ℓ +ν∇ 2z± ℓ +F ℓ, with ν=η= 1.1×10 −4. The...

    work page

  6. [6]

    This part uses the resolved Els¨ asser increments and does not require pressure

    Local-perpendicular angle and amplitude diagnostics The stochastic angle–amplitude diagnostic tests whether amplitude-weighted alignment is caused by a genuine coupling between local angle and event strength. This part uses the resolved Els¨ asser increments and does not require pressure. Unlike the compact-support geo- metric diagnostic above, the local ...

    work page

  7. [7]

    Full definitions of thec r,s r, Pearson, Spear- man, and surrogate diagnostics are given in Appendix B; here we summarize only the implementation

    Cross-scale correlation and surrogate calculations The cross-scale calculations use the same local-B L- perpendicular centered increments as the one-scale angle diagnostics. Full definitions of thec r,s r, Pearson, Spear- man, and surrogate diagnostics are given in Appendix B; here we summarize only the implementation. For each pair of separations (r i, r...

    work page

  8. [8]

    At a given separationr, each valid sampled point is assigned to an amplitude–angle state

    Finite-time state-retention diagnostic The finite-time state-retention diagnostic tests the residence-time interpretation of the amplitude–angle statistics. At a given separationr, each valid sampled point is assigned to an amplitude–angle state. We then ask whether a sample that is initially in one state remains in that same state after a finite simulati...

    work page

  9. [9]

    The fifteen 320 3 subvolumes provide the main statistical ensemble

    Robustness ensemble and reproducibility The 448 3 cube provides the larger-volume reference calculation used in several geometric and stochastic di- agnostic figures. The fifteen 320 3 subvolumes provide the main statistical ensemble. The geometric calcula- tion uses the compact-support filter and theℓ-ranges de- scribed above. The angle, stochastic, and ...

    work page 1998

  10. [10]

    To check that the result is not tied to this particular threshold, Fig

    Amplitude-threshold hierarchy The main comparison uses the top 10% of events ranked by Aτ =|δ τz+| |δτz−|. To check that the result is not tied to this particular threshold, Fig. 16 compares the full sample with the top 20%, top 10%, and top 5% amplitude populations. The folded angle decreases monotonically as the amplitude threshold is made more selectiv...

    work page

  11. [11]

    Figure 17 shows the leave-one-out stability of both the folded-angle hierarchy and the normalized angle–amplitude covari- ance

    Leave-one-out stability We test whether any single 24-hour interval controls the WIND50 ensemble by recomputing the ensemble means after removing each interval in turn. Figure 17 shows the leave-one-out stability of both the folded-angle hierarchy and the normalized angle–amplitude covari- ance. The high-amplitude population remains substan- tially more a...

    work page

  12. [12]

    For each 24-hour interval we compute the mean solar-wind speed Vsw =|⟨v(t)⟩ t| and an interval-level normalized cross helicity, σc = E+ −E − E+ +E − , E ± = D z±(t)− ⟨z ±⟩t 2E t

    Solar-wind speed and Alfv´ enicity splits We next test whether the angle–amplitude hierarchy is confined to a particular broad solar-wind regime. For each 24-hour interval we compute the mean solar-wind speed Vsw =|⟨v(t)⟩ t| and an interval-level normalized cross helicity, σc = E+ −E − E+ +E − , E ± = D z±(t)− ⟨z ±⟩t 2E t . Herez ± =v±b A, and⟨·⟩ t denote...

    work page

  13. [13]

    To test whether the result is spe- cific to that stream sequence, we also analyze a sepa- rate diverse validation ensemble

    Diverse-interval validation ensemble The primary WIND50 ensemble is a clean contigu- ous 1998 sequence. To test whether the result is spe- cific to that stream sequence, we also analyze a sepa- rate diverse validation ensemble. Fifty additional 24- hour Wind intervals were selected independently of the 1998 sequence and distributed across multiple years a...

    work page 1998

  14. [14]

    Boldyrev, Spectrum of magnetohydrodynamic turbu- lence, Phys

    S. Boldyrev, Spectrum of magnetohydrodynamic turbu- lence, Phys. Rev. Lett.96, 115002 (2006)

    work page 2006

  15. [15]

    Mason, F

    J. Mason, F. Cattaneo, and S. Boldyrev, Dynamic align- ment in driven magnetohydrodynamic turbulence, Phys. Rev. Lett.97, 255002 (2006)

    work page 2006

  16. [16]

    Biskamp,Magnetohydrodynamic Turbulence(Cam- bridge University Press, 2003)

    D. Biskamp,Magnetohydrodynamic Turbulence(Cam- bridge University Press, 2003)

    work page 2003

  17. [17]

    A. A. Schekochihin, S. C. Cowley, W. Dorland, G. W. Hammett, G. G. Howes, E. Quataert, and T. Tatsuno, Astrophysical gyrokinetics: Kinetic and fluid turbu- lent cascades in magnetized weakly collisional plasmas, The Astrophysical Journal Supplement Series182, 310 (2009)

    work page 2009

  18. [18]

    A. A. Schekochihin, Mhd turbulence: A biased review, Journal of Plasma Physics88, 905880501 (2022)

    work page 2022

  19. [19]

    Beresnyak, Spectral slope and kolmogorov constant of mhd turbulence, Phys

    A. Beresnyak, Spectral slope and kolmogorov constant of mhd turbulence, Phys. Rev. Lett.106, 075001 (2011)

    work page 2011

  20. [20]

    Lazarian, G

    A. Lazarian, G. L. Eyink, A. Jafari, G. Kowal, H. Li, S. Xu, and E. T. Vishniac, 3d turbulent reconnection: Theory, tests, and astrophysical im- plications, Physics of Plasmas27, 012305 (2020), https://doi.org/10.1063/1.5110603

    work page doi:10.1063/1.5110603 2020

  21. [21]

    Y. Li, E. Perlman, M. Wan, Y. Yang, C. Meneveau, R. Burns, S. Chen, A. Szalay, and G. Eyink, A public turbulence database cluster and applications to study La- grangian evolution of velocity increments in turbulence, Journal of Turbulence9, N31 (2008), arXiv:0804.1703 [physics.flu-dyn]

    work page arXiv 2008

  22. [22]

    nasa.gov/, accessed 2026-05-10

    NASA Space Physics Data Facility Coordinated Data Analysis Web (CDAWeb),https://cdaweb.gsfc. nasa.gov/, accessed 2026-05-10

    work page 2026

  23. [23]

    R. P. Lepping, M. H. Acu˜ na, L. F. Burlaga, W. M. Far- rell, J. A. Slavin, K. H. Schatten, F. Mariani, N. F. Ness, F. M. Neubauer, Y. C. Whang, J. B. Byrnes, R. S. Ken- non, P. V. Panetta, J. Scheifele, and E. M. Worley, The Wind Magnetic Field Investigation, Space Science Re- views71, 207 (1995)

    work page 1995

  24. [24]

    R. P. Lin, K. A. Anderson, S. Ashford, C. Carlson, D. Curtis, R. Ergun, D. Larson, J. McFadden, M. Mc- Carthy, G. K. Parks, H. R` eme, J. M. Bosqued, J. Coute- lier, F. Cotin, C. D’Uston, K.-P. Wenzel, T. R. Sander- son, J. Henrion, J. C. Ronnet, and G. Paschmann, A Three-Dimensional Plasma and Energetic Particle Inves- tigation for the Wind Spacecraft, S...

    work page 1995

  25. [25]

    Jafari, Lagrangian approach to reconnection and topology change, Phys

    A. Jafari, Lagrangian approach to reconnection and topology change, Phys. Rev. E111, 065212 (2025)

    work page 2025

  26. [26]

    G. I. Taylor, The spectrum of turbulence, Proceed- ings of the Royal Society of London. A. Math- ematical and Physical Sciences164, 476 (1938), https://royalsocietypublishing.org/rspa/article- pdf/164/919/476/40796/rspa.1938.0032.pdf

    work page arXiv 1938

  27. [27]

    Friedrich and J

    R. Friedrich and J. Peinke, Description of a turbulent cascade by a fokker-planck equation, Phys. Rev. Lett. 78, 863 (1997)

    work page 1997

  28. [28]

    Renner, J

    C. Renner, J. Peinke, and R. Friedrich, Experimental in- dications for markov properties of small-scale turbulence, Journal of Fluid Mechanics433, 383–409 (2001)

    work page 2001

  29. [29]

    Zwanzig, Memory effects in irreversible thermodynam- ics, Phys

    R. Zwanzig, Memory effects in irreversible thermodynam- ics, Phys. Rev.124, 983 (1961)

    work page 1961

  30. [30]

    Mori, Transport, collective motion, and brown- ian motion*), Progress of Theoretical Physics33, 423 (1965), https://academic.oup.com/ptp/article- pdf/33/3/423/5428510/33-3-423.pdf

    H. Mori, Transport, collective motion, and brown- ian motion*), Progress of Theoretical Physics33, 423 (1965), https://academic.oup.com/ptp/article- pdf/33/3/423/5428510/33-3-423.pdf

    work page 1965