pith. sign in

arxiv: 2604.17113 · v1 · submitted 2026-04-18 · ⚛️ physics.plasm-ph

Thermal Effects on Buneman Instability: A Vlasov-Poisson Study

Chingangbam Amudon , Sanjeev Kumar Pandey , Rajaraman Ganesh This is my paper

Pith reviewed 2026-05-10 06:28 UTC · model grok-4.3

classification ⚛️ physics.plasm-ph
keywords Buneman instabilityVlasov-Poisson simulationthermal effectsgrowth rateplasma heatingdensity inhomogeneitybeam-plasma interaction
0
0 comments X

The pith

Maximum growth rate of Buneman instability remains independent of temperature ratio

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

The paper investigates thermal effects on the Buneman instability through Vlasov-Poisson simulations. It establishes that the maximum growth rate follows the familiar dependence on the mass ratio but shows little to no variation with the temperature ratio between ions and electrons. This result differs from both fluid models and linearized kinetic theory. The simulations further reveal that ion density inhomogeneities control the conversion of electron beam energy into plasma heating, with efficiency dropping as the plasma warms and density variations diminish.

Core claim

Numerical integration of the Vlasov-Poisson system shows that the peak growth rate of the Buneman instability recovers the (m/M)^{1/3} scaling yet stays essentially constant across different temperature ratios of the constituent species. The amplitude of the resulting ion density inhomogeneity regulates the generation of sidebands and thereby the efficiency with which beam energy is transferred into bulk plasma temperature.

What carries the argument

Self-consistent Vlasov-Poisson evolution of distribution functions for electrons and ions, initialized with a drifting electron beam and thermal velocities for both species.

If this is right

  • Predictions of instability growth in plasmas can proceed without precise temperature ratio inputs in the kinetic regime.
  • Reduced ion density inhomogeneity in warm plasmas leads to lower sideband activity and less efficient beam energy deposition.
  • The energy transfer mechanism is self-regulated by the nonlinear density structures that develop.
  • Kinetic models are necessary to capture the temperature-independent growth unlike fluid approximations.

Where Pith is reading between the lines

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

  • This temperature independence could allow simpler parametrizations in large-scale plasma simulations for space weather or fusion research.
  • Experimental measurements of growth rates in varying temperature conditions might confirm or refute the numerical findings.
  • The control by ion density inhomogeneity points to potential ways to manipulate energy transfer through initial density perturbations.

Load-bearing premise

The Vlasov-Poisson numerical method and the specific initial conditions faithfully represent the thermal physics without numerical artifacts suppressing the dependence on temperature ratio.

What would settle it

Observing or simulating the growth rate over a range of electron-to-ion temperature ratios from 0.01 to 100 and finding whether the maximum value stays constant or varies systematically.

Figures

Figures reproduced from arXiv: 2604.17113 by Chingangbam Amudon, Rajaraman Ganesh, Sanjeev Kumar Pandey.

Figure 1
Figure 1. Figure 1: FIG. 1: Space-averaged Electron and Ion Velocity Distribution i.e., [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Temporal Evolution of Fourier Electric mode (fundamental and higher harmonics) [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Phase Space Information of electron (first row) and ion (third row) with their [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: Temporal Evolution of the Ratio of Potential Energy with initial Kinetic Energy, [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Growth Rate (in Base 10) of Buneman instability for the different plasma limit – [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: Temporal Evolution of dominant Fourier Electric mode, [PITH_FULL_IMAGE:figures/full_fig_p011_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7: Variation of Maximum Growth Rate, [PITH_FULL_IMAGE:figures/full_fig_p012_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8: Temporal Evolution of dominant Fourier Electric mode for White Noise [PITH_FULL_IMAGE:figures/full_fig_p012_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9: Temporal Evolution of Fourier Electric mode (fundamental and higher harmonics) [PITH_FULL_IMAGE:figures/full_fig_p013_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10: The Phase Space of Electron (first row) and Ion (third row) with their density [PITH_FULL_IMAGE:figures/full_fig_p014_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11: The Temporal Evolution of Ratio of Electrostatic Energy with Initial Kinetic [PITH_FULL_IMAGE:figures/full_fig_p015_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12: Numerical Solution of the Fluid Dispersion Relations - Cold Plasma and Warm [PITH_FULL_IMAGE:figures/full_fig_p019_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: FIG. 13: Information on the Domain Decomposition and Sub-Domain. [PITH_FULL_IMAGE:figures/full_fig_p020_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: FIG. 14: Electron Distribution Function: (a) the space-averaged [PITH_FULL_IMAGE:figures/full_fig_p021_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: FIG. 15: Temporal Evolution of (a) dominant Fourier Electric mode, [PITH_FULL_IMAGE:figures/full_fig_p021_15.png] view at source ↗
Figure 15
Figure 15. Figure 15: Initially, the wave gets Landau damped; as the system evolves, the wave-particle [PITH_FULL_IMAGE:figures/full_fig_p022_15.png] view at source ↗
read the original abstract

Buneman instability has been extensively studied, and related aspects, namely anomalous resistivity, have been explored in detail using analytical theory as well as numerical simulations based on Particle-in-Cell and Vlasov solvers. Most numerical studies have focused on understanding the nonlinear evolution of the instability. In the present study, the growth rate of the Buneman instability in the presence of thermal effects of the constituent species (i.e., ions and electrons) is investigated. It is observed that the growth rate differs significantly from that obtained using fluid models (both cold and warm) as well as from linearized kinetic models. While the well-known result of $(m/M)^{1/3}$ dependence of the maximum growth rate is recovered, it is shown that the maximum growth rate is essentially independent of the temperature ratio of the constituent species. It is further demonstrated numerically that the amplitude of ion density inhomogeneity self-consistently controls the transfer of electron beam energy into the bulk plasma temperature. In particular, as one moves from the cold to the warm plasma limit, the decrease in ion density inhomogeneity reduces the generation of sidebands and thus lowers the transfer efficiency.

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 presents a Vlasov-Poisson numerical study of the Buneman instability incorporating thermal effects for both electrons and ions. It reports that the maximum growth rate recovers the expected (m/M)^{1/3} scaling with mass ratio but is essentially independent of the electron-to-ion temperature ratio, in contrast to predictions from both cold/warm fluid models and linearized kinetic theory. The work further claims that the amplitude of self-consistent ion density inhomogeneity controls the efficiency of energy transfer from the electron beam into bulk plasma heating, with reduced inhomogeneity in the warm limit suppressing sideband generation.

Significance. If the numerical results hold under proper validation, the reported temperature-ratio independence would be significant, as it contradicts the expected damping role of thermal velocities through the plasma dispersion function and could revise models of anomalous resistivity and beam-plasma energy coupling. The nonlinear demonstration that density inhomogeneity self-consistently modulates sideband generation and transfer efficiency provides a concrete mechanistic insight. The direct use of the Vlasov-Poisson system without ad-hoc closures is a methodological strength.

major comments (2)
  1. [Abstract] Abstract: The central claim that the maximum growth rate 'is essentially independent of the temperature ratio' is presented without any resolution studies, error bars on extracted growth rates, or explicit comparison of simulated ω_i(k) against numerical roots of the warm-plasma dispersion relation. This is load-bearing, as the skeptic concern of velocity-space dissipation preferentially damping high-k modes at higher T_e/T_i cannot be ruled out from the given information.
  2. [Abstract] Abstract: No details are provided on the Vlasov-Poisson solver parameters (velocity-grid spacing, time step, or any Landau-damping filter) or on the precise diagnostic used to extract growth rates (e.g., early-time exponential fit window or specific k-mode selection). Without these, the reported deviation from linearized kinetic predictions cannot be assessed for numerical fidelity.
minor comments (1)
  1. The abstract states that results 'differ significantly' from fluid and kinetic models but does not reference the specific figures or sections where the growth-rate curves or energy-transfer diagnostics are shown.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for highlighting important points regarding validation and reproducibility. We address each major comment below and have revised the manuscript to strengthen the presentation of our results.

read point-by-point responses

  1. Referee: [Abstract] Abstract: The central claim that the maximum growth rate 'is essentially independent of the temperature ratio' is presented without any resolution studies, error bars on extracted growth rates, or explicit comparison of simulated ω_i(k) against numerical roots of the warm-plasma dispersion relation. This is load-bearing, as the skeptic concern of velocity-space dissipation preferentially damping high-k modes at higher T_e/T_i cannot be ruled out from the given information.

    Authors: We agree that the abstract, being a concise summary, does not include supporting validation details. The main text presents growth-rate results across temperature ratios via direct Vlasov-Poisson simulations, but we acknowledge the need for explicit checks. In the revised manuscript we add a dedicated subsection with resolution studies (doubling velocity-space resolution and confirming growth rates converge within a few percent), error bars obtained from the linear-fit procedure, and a direct overlay of simulated ω_i(k) against numerically solved roots of the warm-plasma dispersion relation. These additions allow readers to assess whether the reported temperature-ratio independence is physical or an artifact of numerical dissipation. revision: yes

  2. Referee: [Abstract] Abstract: No details are provided on the Vlasov-Poisson solver parameters (velocity-grid spacing, time step, or any Landau-damping filter) or on the precise diagnostic used to extract growth rates (e.g., early-time exponential fit window or specific k-mode selection). Without these, the reported deviation from linearized kinetic predictions cannot be assessed for numerical fidelity.

    Authors: We accept that the original submission omitted these implementation specifics. The revised manuscript now contains an expanded Numerical Methods section that reports the velocity-grid spacing, time-step size, absence of any explicit Landau-damping filter, and the exact diagnostic: growth rates are obtained from a least-squares exponential fit to the early-time electric-field energy (linear phase) for the single k-mode that exhibits the largest growth. These details enable independent assessment of the deviation from linearized kinetic theory. revision: yes

Circularity Check

0 steps flagged

No circularity: results are direct numerical outputs from Vlasov-Poisson integration

full rationale

The paper reports growth rates obtained by solving the Vlasov-Poisson system numerically for varying temperature ratios. It recovers the known (m/M)^{1/3} scaling for the cold case and observes that the maximum growth rate shows little dependence on T_e/T_i within the simulated parameter space. These are empirical findings from the solver, not an analytical derivation that reduces to a fitted parameter, self-citation, or ansatz by construction. No load-bearing step equates a prediction to its own input; the work is a standard kinetic simulation study whose central claim rests on the fidelity of the numerical method rather than on any self-referential logic.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claims rest on the standard Vlasov-Poisson system plus numerical discretization; no new entities or ad-hoc parameters are introduced in the abstract.

axioms (2)
  • standard math The Vlasov-Poisson equations govern the evolution of the electron and ion distribution functions.
    Invoked implicitly as the simulation framework throughout the abstract.
  • domain assumption Initial conditions consist of a drifting electron beam and stationary ions with finite thermal spreads.
    Required to set up the Buneman configuration but not detailed in the abstract.

pith-pipeline@v0.9.0 · 5508 in / 1340 out tokens · 26228 ms · 2026-05-10T06:28:26.232241+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

41 extracted references · 41 canonical work pages

  1. [1]

    This reduced steepening is due to the finite thermal spread of the considered species. In the present study, the electron hole-ion soliton coupling is captured with high precision and resolution (see the electron phase space, first row and ion density, fourth row in column iv and v of Figure 3). Furthermore, the electron’s streaming/beam energy is entirel...

    work page

  2. [2]

    dispersion relations are used to validate the results of the simulation. From Figure 5, it is observed that the introduction of thermal effects into the mathe- matical model introduces correction on the growth rate of the instability (refer Appendix A for details), thus deviates from the cold plasma fluid results; the maximum growth rate has decreased and...

    work page

  3. [3]

    3 and Eqn

    generally overestimates the growth rate of the instability while the kinetic dispersion (Eqn. 3 and Eqn. 4) provides accurate growth rates for small values of Doppler shift but tends to underestimate them at higher values in warm plasma limit. The overlapping of the two Kinetic Dispersion Solvers (shown in Figure 5) with regards to the value of growth rat...

    work page

  4. [4]

    A. Wong, R. Motley, and N. d’Angelo, Physical Review133, A436 (1964)

    work page 1964

  5. [5]

    Buneman, Phys

    O. Buneman, Phys. Rev. Lett.1, 8 (1958)

    work page 1958

  6. [6]

    Hirose, Plasma Physics20, 481 (1978)

    A. Hirose, Plasma Physics20, 481 (1978)

    work page 1978

  7. [7]

    Yoon and T

    P. Yoon and T. Umeda, Physics of Plasmas17(2010)

    work page 2010

  8. [8]

    C. Watt, R. Horne, and M. Freeman, Geophysical Research Letters29, 4 (2002)

    work page 2002

  9. [9]

    Omura, W

    Y. Omura, W. Heikkila, T. Umeda, K. Ninomiya, and H. Matsumoto, J. Geophys. Res.: Space Physics108(2003)

    work page 2003

  10. [10]

    Ishihara, A

    O. Ishihara, A. Hirose, and A. Langdon, Physical Review Letters44, 1404 (1980)

    work page 1980

  11. [11]

    Ishihara, A

    O. Ishihara, A. Hirose, and A. Langdon, The Physics of Fluids24, 452 (1981)

    work page 1981

  12. [12]

    Hirose, O

    A. Hirose, O. Ishihara, and A. Langdon, The Physics of Fluids25, 610 (1982)

    work page 1982

  13. [13]

    Niknam, D

    A. Niknam, D. Komaizi, and M. Hashemzadeh, Physics of Plasmas18(2011)

    work page 2011

  14. [14]

    Niknam, H

    A. Niknam, H. Roozbahani, M. Hashemzadeh, and D. Komaizi, Physics of Plasmas21(2014)

    work page 2014

  15. [15]

    Shokri and A

    B. Shokri and A. Niknam, Physics of plasmas12(2005)

    work page 2005

  16. [16]

    Hashemzadeh, Physica A: Statistical Mechanics and its Applications459, 68 (2016)

    M. Hashemzadeh, Physica A: Statistical Mechanics and its Applications459, 68 (2016)

    work page 2016

  17. [17]

    Buneman, Phys

    O. Buneman, Phys. Rev.115, 503 (1959)

    work page 1959

  18. [18]

    R. S. Rajawat and S. Sengupta, Physics of Plasmas24(2017)

    work page 2017

  19. [19]

    Raghunathan and R

    M. Raghunathan and R. Ganesh, Physics of Plasmas20(2013)

    work page 2013

  20. [20]

    S. K. Pandey, Ph.D. thesis, Homi Bhabha National Institute (2023)

    work page 2023

  21. [21]

    Trivedi, Ph.D

    P. Trivedi, Ph.D. thesis, Homi Bhabha National Institute (2019)

    work page 2019

  22. [22]

    C. Z. Cheng and G. Knorr, J. Comput. Phys.22, 330 (1976)

    work page 1976

  23. [23]

    Collela and P

    P. Collela and P. R. Woodward, J. Comput. Phys.54, 174 (1984)

    work page 1984

  24. [24]

    PlasmaPy,https://www.plasmapy.org(2025)

    work page 2025

  25. [25]

    Goldston,Introduction to Plasma Physics(CRC Press, 2020)

    R. Goldston,Introduction to Plasma Physics(CRC Press, 2020). 22

    work page 2020

  26. [26]

    Krall and A

    N. Krall and A. Trivelpiece,Principles of Plasma Physics(McGraw-Hill, 1973)

    work page 1973

  27. [27]

    Petkaki, C

    P. Petkaki, C. Watt, R. Horne, and M. Freeman, J. Geophys. Res.: Space Physics108(2003)

    work page 2003

  28. [28]

    Petkaki, M

    P. Petkaki, M. Freeman, T. Kirk, C. Watt, and R. Horne, J. Geophys. Res.: Space Physics 111(2006)

    work page 2006

  29. [29]

    I. H. Hutchinson, (2024), arXiv:2411.12821 [physics.plasm-ph]

    work page arXiv 2024

  30. [30]

    Pavan, P

    J. Pavan, P. Yoon, and T. Umeda, Physics of Plasmas18(2011)

    work page 2011

  31. [31]

    Tavassoli, O

    A. Tavassoli, O. Chapurin, M. Jimenez, M. Papahn Zadeh, T. Zintel, M. Sengupta, L. Cou¨ edel, R. J. Spiteri, M. Shoucri, and A. Smolyakov, Physics of Plasmas28(2021)

    work page 2021

  32. [32]

    P. K. Kaw, A. Lin, and J. Dawson, The Physics of Fluids16, 1967 (1973)

    work page 1967

  33. [33]

    S. K. Pandey, J. Mahapatra, and R. Ganesh, Physica Scripta97, 105602 (2022)

    work page 2022

  34. [34]

    Tappert, The Physics of Fluids15, 2446 (1972)

    F. Tappert, The Physics of Fluids15, 2446 (1972)

    work page 1972

  35. [35]

    Sakanaka, The Physics of Fluids15, 304 (1972)

    P. Sakanaka, The Physics of Fluids15, 304 (1972)

    work page 1972

  36. [36]

    Tagare, Plasma Physics15, 1247 (1973)

    S. Tagare, Plasma Physics15, 1247 (1973)

    work page 1973

  37. [37]

    S. K. Pandey and R. Ganesh, AIP Advances11(2021)

    work page 2021

  38. [38]

    S. K. Pandey and R. Ganesh, Physica Scripta96, 125616 (2021)

    work page 2021

  39. [39]

    S. K. Pandey and R. Ganesh, Physica Scripta96, 125615 (2021)

    work page 2021

  40. [40]

    S. K. Pandey and R. Ganesh, Physica Scripta99, 125608 (2024)

    work page 2024

  41. [41]

    Manfredi, Physical Review Letters79, 2815 (1997)

    G. Manfredi, Physical Review Letters79, 2815 (1997). 23

    work page 1997