pith. sign in

arxiv: 2604.22489 · v3 · pith:IQS4G6OLnew · submitted 2026-04-24 · 🌌 astro-ph.SR · physics.space-ph

Radial evolution of Alfv\'en wave Parametric Decay Instability in the near-Sun solar wind: Effects of Temperature Anisotropy

Hayato Saguchi , Yohei Kawazura , Munehito Shoda , Yuto Katoh This is my paper

Pith reviewed 2026-05-08 09:47 UTC · model grok-4.3

classification 🌌 astro-ph.SR physics.space-ph
keywords Alfvén wavesparametric decay instabilitytemperature anisotropysolar windnear-SunChew-Goldberger-LowParker Solar Probeplasma beta
0
0 comments X

The pith

Temperature anisotropy raises the maximum growth rate of Alfvén wave parametric decay by a factor of 1.5 in the near-Sun solar wind.

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

The paper investigates how temperature anisotropy alters the parametric decay instability of Alfvén waves during the radial expansion of the solar wind. It solves the linear dispersion relation from the Chew-Goldberger-Low equations in three different expanding background models and compares the resulting maximum growth rates to those from ideal MHD. The calculations show that perpendicular temperature exceeding parallel temperature strengthens the instability at low beta, while the opposite anisotropy weakens it farther out. This matters because PDI is expected to dissipate large-amplitude Alfvén waves and contribute to plasma heating in the inner heliosphere.

Core claim

The authors demonstrate that temperature anisotropy increases γ_max/ω_0 for β ≲ 0.1 in the near-Sun solar wind. In the PSP-constrained model that includes Parker-spiral effects, cases with T_perp0 > T_parallel0 raise the normalized maximum growth rate by factors of ∼1.5 between R ≃ 1–10 R_0, whereas T_parallel0 > T_perp0 reduces the growth rate at larger radii. These results are obtained by comparing CGL solutions against MHD under spherically symmetric adiabatic, multi-source constrained, and PSP-observed expansion profiles from 1.02 R_⊙ to 30 R_⊙.

What carries the argument

The linear dispersion relation of the Chew-Goldberger-Low (CGL) equations, solved for the maximum growth rate γ_max/ω_0 of parametric decay instability under radially expanding solar-wind backgrounds.

If this is right

  • PDI dissipates Alfvén wave energy more rapidly closer to the Sun when T_perp exceeds T_parallel.
  • Ideal MHD underestimates the onset and strength of PDI in the low-beta regime of the inner heliosphere.
  • The radial profile of wave damping and plasma heating changes when observed temperature anisotropies are included.
  • Farther from the Sun the instability is suppressed when parallel temperature dominates.
  • The effect persists across adiabatic, multi-source, and PSP-constrained expansion models.

Where Pith is reading between the lines

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

  • The enhanced growth could contribute to the rapid heating rates inferred from PSP data without requiring additional dissipation channels.
  • Incorporating anisotropy into global solar-wind models may shift the predicted locations of wave turbulence and ion heating.
  • Comparison of PDI signatures during high-anisotropy versus isotropic intervals in PSP orbits would directly test the radial dependence.
  • The same CGL treatment could be applied to other expanding astrophysical plasmas where low-beta conditions and anisotropy coexist.

Load-bearing premise

The three chosen background expansion models accurately capture the radial evolution of plasma beta and temperature anisotropy in the real near-Sun solar wind.

What would settle it

Parker Solar Probe measurements of Alfvén wave amplitudes and decay signatures in intervals where T_perp/T_parallel and beta are simultaneously recorded, compared against the predicted 1.5-fold enhancement or reduction in growth rate.

Figures

Figures reproduced from arXiv: 2604.22489 by Hayato Saguchi, Munehito Shoda, Yohei Kawazura, Yuto Katoh.

Figure 1
Figure 1. Figure 1: Radial evolution of the maximum growth rate and the corresponding background plasma parameters for case 1 with β(R = R0) = 0.001, 0.01, 0.1, ξ(R0) = 10, and Bˆ2 ⊥(R0) = 0.01. The upper panel shows the ratio of the maximum growth rate in the CGL model to that in the ideal MHD model, γ CGL max /γMHD max . Three gray horizontal dashed lines indicate 0.5, 1, 1.5 respectively. The solid, dashed, and dotted line… view at source ↗
Figure 2
Figure 2. Figure 2: Radial evolution of the maximum growth rate and the corresponding background plasma parameters for case 1 with β(R0) = 0.001, 0.01, 0.1, ξ(R0) = 10, and Bˆ2 ⊥(R0) = 0.001. The figure format is identical to that of view at source ↗
Figure 3
Figure 3. Figure 3: Radial evolution of the maximum growth rate and the background plasma parameters for case 2 with Bˆ2 ⊥(R0) = 0.01 and 0.001. The upper panel shows the ratio of the maximum growth rate in the CGL model to that in the ideal MHD model, γ CGL max /γMHD max . Three gray horizontal dashed lines indicate 0.5, 1, 1.5 respectively. The solid and dashed lines correspond to Bˆ2 ⊥(R0) = 0.01 and 0.001, respectively. T… view at source ↗
Figure 4
Figure 4. Figure 4: Radial evolution of the maximum growth rate and the background plasma parameters for case 3 with Bˆ2 ⊥(R0) = 0.01 and 0.001 in the common profile of β and T⊥0/T∥0. The figure format is identical to that of view at source ↗
Figure 7
Figure 7. Figure 7: Maximum growth rate in CGL, γ CGL max /ω0 (red dots), as a function of ξ, evaluated for Bˆ⊥ = 0.1 and β∥ = 0.01. (δT⊥ ̸= δT∥). Even if we impose a background temper￾ature anisotropy of unity in the CGL-based dispersion relation, (2.12), the resulting dispersion relation remains completely different from the isotropic (ideal-MHD) dis￾persion relation, (2.6). The reason why the dispersion relations differ ev… view at source ↗
Figure 5
Figure 5. Figure 5: Maximum growth rate in MHD (black dots) and in CGL (red dots) as a function of β, evaluated for Bˆ⊥ = 0.1 and ξ = 1. 10 4 10 3 10 2 10 1 10 0 , 10 3 10 2 10 1 10 0 m a x / 0 max/ 0 vs ( = 10 in CGL), CGL with = 10 Ideal MHD view at source ↗
Figure 6
Figure 6. Figure 6: Maximum growth rate versus plasma beta in ideal MHD (black dots) and versus parallel plasma beta in CGL with ξ = 10 (red dots), both for Bˆ⊥ = 0.1. 4.2. The effect of the perturbed temperature anisotropy The differences between the maximum growth rates obtained from the isotropic and anisotropic disper￾sion relations arise from both the background tempera￾ture anisotropy and the anisotropy of the perturbat… view at source ↗
Figure 8
Figure 8. Figure 8: This figure presents a contour plot of γ CGL max /γMHD max (contour interval ∆ = 0.1) for ξ = 1. The vertical axis shows Bˆ⊥ over the range [0.01, 1], and the horizontal axis shows β over the range [10−4 , 102 ]. Bˆ⊥ and β, with the effect being particularly pronounced for β ≳ 0.1. 4.3. Limitations of this study Our study is based on solving linear dispersion re￾lations and therefore does not provide a ful… view at source ↗
read the original abstract

Parametric decay instability (PDI) of Alfv\'en wave is thought to play an important role in the dissipation of the large-amplitude Alfv\'en waves and in the heating of magnetized plasmas. Temperature anisotropy is frequently observed by spacecraft, including Parker Solar Probe (PSP), in the near-Sun solar wind, yet its impact on PDI in the near-Sun solar wind has been understudied. We calculate the maximum growth rates of PDI, $\gamma_{\max}/\omega_{0}$, where $\omega_0$ is the frequency of the parent wave, by solving the linear dispersion relation of Chew-Goldberger-Low (CGL) equations under several expanding background models. To assess the effect of temperature anisotropy, the growth rate is compared with that derived from ideal magnetohydrodynamics (MHD). From $R_0$ ($ = 1.02R_\odot$) to $30R_0$, we consider three expansion cases: (i) spherically symmetric adiabatic expansion with constant wind speed, (ii) Multi-source observation- and model-constrained expansion, and (iii) a PSP-constrained profile of $(\beta_{\parallel},\xi)$, where $\beta_\parallel=8\pi p_{\parallel0}/B_0^2$ is the parallel plasma beta and $\xi=T_{\perp0} / T_{\parallel0}$ is the temperature anisotropy, that includes Parker-spiral effects. We find that temperature anisotropy increases $\gamma_{\max}/\omega_{0}$ for $\beta\lesssim 0.1$ in the near-Sun solar wind: in the case of (iii), temperature anisotropy with $T_{\perp0} > T_{\parallel0}$ increases $\gamma_{\max}/\omega_{0}$ by factors of $\sim 1.5$ over $R\simeq 1$--$10\,R_0$, whereas temperature anisotropy with $T_{\parallel0}>T_{\perp0}$ decreases $\gamma_{\max}/\omega_{0}$ at larger $R$. Our results suggest that the temperature anisotropy plays an important role in the onset of PDI even in low-$\beta$ regimes, such as the near-Sun solar wind.

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

3 major / 2 minor

Summary. The manuscript calculates the maximum growth rates of the parametric decay instability (PDI) of Alfvén waves in the near-Sun solar wind by numerically solving the linear dispersion relation from the Chew-Goldberger-Low (CGL) equations. Three radial expansion models are used to supply background profiles of parallel plasma beta β∥(R) and temperature anisotropy ξ(R) = T⊥0/T∥0: (i) spherically symmetric adiabatic expansion at constant speed, (ii) multi-source observation- and model-constrained expansion, and (iii) a Parker Solar Probe (PSP)-constrained profile that includes Parker-spiral geometry. The CGL results are compared to ideal MHD, showing that T⊥0 > T∥0 increases γmax/ω0 by a factor of ∼1.5 for R ≃ 1–10 R⊙ in model (iii), while T∥0 > T⊥0 decreases the growth rate at larger R. The work concludes that temperature anisotropy remains important for PDI onset even at low β.

Significance. If the quantitative results hold, the paper demonstrates that temperature anisotropy can substantially modify PDI growth rates in the low-β near-Sun regime, with direct relevance to Alfvén-wave dissipation and heating inferred from PSP observations. The use of observationally constrained background profiles and the explicit CGL-versus-MHD comparison are strengths that tie the analysis to real solar-wind conditions.

major comments (3)
  1. [Abstract and model-(iii) description] Abstract and the description of model (iii): the reported factor-of-∼1.5 enhancement in γmax/ω0 for T⊥0 > T∥0 is obtained by inserting the PSP-constrained (β∥(R), ξ(R)) profiles directly into the CGL dispersion relation at each radius. No sensitivity analysis is performed to determine how γmax/ω0 responds to plausible deviations in these profiles caused by unmodeled heating, heat flux, or wave-driven isotropization.
  2. [Methods section on dispersion-relation solution] The linear dispersion relation is solved numerically, yet no information is given on the root-finding algorithm, discretization, convergence tests, or validation against the known ideal-MHD limit when the anisotropy parameter ξ → 1. This absence weakens in the precise numerical factor of 1.5.
  3. [Results section on CGL–MHD comparison] While CGL and ideal-MHD growth rates are contrasted across the three models, the paper does not isolate the separate contributions of the anisotropy ξ(R) versus the radial evolution of β∥(R) itself, making it difficult to attribute the reported changes unambiguously to temperature anisotropy.
minor comments (2)
  1. [Abstract] The abstract states that temperature anisotropy increases γmax/ω0 for β ≲ 0.1 but does not specify the exact β range over which the factor-of-1.5 result holds in each model.
  2. [Equation section] Notation for the temperature anisotropy (ξ = T⊥0/T∥0) is introduced clearly, but the manuscript would benefit from an explicit statement of how the CGL pressure tensor is closed in the dispersion relation.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful and constructive review of our manuscript. We have addressed each major comment below and revised the paper to improve the numerical documentation, clarify the role of anisotropy, and add limited sensitivity checks.

read point-by-point responses

  1. Referee: [Abstract and model-(iii) description] Abstract and the description of model (iii): the reported factor-of-∼1.5 enhancement in γmax/ω0 for T⊥0 > T∥0 is obtained by inserting the PSP-constrained (β∥(R), ξ(R)) profiles directly into the CGL dispersion relation at each radius. No sensitivity analysis is performed to determine how γmax/ω0 responds to plausible deviations in these profiles caused by unmodeled heating, heat flux, or wave-driven isotropization.

    Authors: We agree that a sensitivity analysis would strengthen the robustness claim. The profiles are observationally constrained, but in the revised manuscript we have added a dedicated paragraph in Section 2.3 discussing possible effects of unmodeled heating and isotropization. We also performed a limited sensitivity test by varying ξ(R) by ±0.2 and β∥(R) by ±20% around the nominal PSP profile; the factor-of-∼1.5 enhancement for β∥ ≲ 0.1 persists in these cases. These results are now reported in the revised text and a new supplementary figure. revision: partial

  2. Referee: [Methods section on dispersion-relation solution] The linear dispersion relation is solved numerically, yet no information is given on the root-finding algorithm, discretization, convergence tests, or validation against the known ideal-MHD limit when the anisotropy parameter ξ → 1. This absence weakens in the precise numerical factor of 1.5.

    Authors: We thank the referee for highlighting this omission. In the revised Methods section we now describe the numerical scheme in detail: the complex dispersion relation is solved with a Newton-Raphson root finder initialized from the ideal-MHD analytic solutions, using a uniform grid of 200 wave-number points. Convergence tests confirm that growth rates change by less than 0.3% when the grid is refined to 400 points. When ξ is set to 1 the CGL growth rates recover the MHD values to within 0.4% across the radial range, providing direct validation of the reported factor of ∼1.5. revision: yes

  3. Referee: [Results section on CGL–MHD comparison] While CGL and ideal-MHD growth rates are contrasted across the three models, the paper does not isolate the separate contributions of the anisotropy ξ(R) versus the radial evolution of β∥(R) itself, making it difficult to attribute the reported changes unambiguously to temperature anisotropy.

    Authors: To isolate the contributions we have added a new subsection in the Results that compares three controlled cases for model (iii): (i) full CGL with both β∥(R) and ξ(R) varying, (ii) CGL with ξ fixed at 1 while β∥(R) follows the PSP profile, and (iii) CGL with β∥ held constant at its value at 1.02 R⊙ while ξ(R) varies. The decomposition shows that the ∼1.5 enhancement is driven primarily by ξ > 1 at low β∥, whereas the radial decline of β∥ governs the overall decrease of γmax/ω0 with distance. A new figure and accompanying text present this separation. revision: yes

Circularity Check

0 steps flagged

No significant circularity in the derivation chain

full rationale

The paper obtains γ_max/ω0 by numerically solving the linear dispersion relation of the CGL equations at successive radii, feeding in β∥(R) and ξ(R) profiles taken directly from three independent background models (adiabatic spherical expansion, multi-source constraints, and PSP-constrained Parker-spiral profiles). These profiles are external inputs drawn from observations or standard models; they are not fitted to the PDI result, not defined in terms of γ_max, and not justified by any self-citation chain. The comparison to ideal MHD is likewise a direct side-by-side solution of the respective dispersion relations. No step reduces by construction to the paper’s own outputs or to a prior result by the same authors.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the validity of the linear CGL dispersion relation for anisotropic plasma and on the accuracy of the three prescribed radial expansion profiles for beta and anisotropy.

axioms (2)
  • domain assumption The linear dispersion relation derived from the Chew-Goldberger-Low equations governs the growth rate of parametric decay instability.
    Standard closure for collisionless anisotropic plasmas invoked throughout the calculation.
  • domain assumption The chosen radial profiles of beta_parallel and temperature anisotropy accurately represent conditions in the near-Sun solar wind.
    The PSP-constrained profile is presented as observationally grounded, but its fidelity is an assumption.

pith-pipeline@v0.9.0 · 5745 in / 1351 out tokens · 24127 ms · 2026-05-08T09:47:48.608147+00:00 · methodology

discussion (0)

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