Information Propagation in Predator-Prey Dynamics of Turbulent Plasma

TL;DR

The paper addresses the mechanism and causality of predator–prey–like oscillations in magnetically confined plasmas by constructing a minimal stochastic model with intrinsic noise for zonal flow $Z$ and drift‑wave turbulence $W$, and analyzing it via system‑size expansion. It derives deterministic rate equations and linear Langevin fluctuations around the coexistence fixed point, demonstrating that quasi‑cycles arise from stochastic amplification and that the steady state is non‑equilibrium Gaussian with nonzero probability current. It then quantifies predictive causality using information flow, proving that information about $Z$ propagates to $W$ in steady state across parameter regimes, even without quasi‑cycles, which has implications for turbulence control guided by zonal‑flow characteristics. The work provides a principled information‑theoretic framework for turbulence regulation in fusion plasmas and offers insights applicable to general predator–prey dynamics in complex systems.

Abstract

Magnetically confined fusion plasmas exhibit predator-prey-like cyclic oscillations through the self-regulating interaction between drift-wave turbulence and zonal flow. To elucidate the detailed mechanism and causality underlying this phenomenon, we construct a simple stochastic predator-prey model that incorporates intrinsic fluctuations and analyze its statistical properties from an information-theoretic perspective. We first show that the model exhibits persistent fluctuating cyclic oscillations called quasi-cycles due to amplification of intrinsic noise. This result suggests the possibility that the previously observed periodic oscillations in a toroidal plasma are not limit cycles but quasi-cycles, and that such quasi-cycles may be widely observed under various conditions. For this model, we further prove that information of the zonal flow is propagated to turbulence. This result suggests that turbulence behavior may be predictable to a certain extent based on zonal flow characteristics.

Figures (4)

• Figure 1: (a) Schematic of the model. (b) Vector fields for parameter values $p=0.5, d=b=\gamma=1$ (left) and $p=2, d=b=\gamma=1$ (right). The color bar denotes the magnitude of the vector fields. The red markers indicate the stable fixed points, and the blue markers indicate the unstable fixed points. (c) The trajectory of the Langevin equation. The parameter values are $\epsilon=1/2, \beta=1, D^W=1/2$. The characteristic period is given by $2\pi/\omega_c\simeq10$.
• Figure 2: (a) Schematic of information flow $\dot{I}^{Z\rightarrow W}_\lambda$. (b) $(\epsilon,\beta)$-dependence of $\dot{I}^{Z\rightarrow W}_{\lambda=1/2}$ with $D^W=1/2$. The dashed line denotes the resonance condition $2\epsilon \beta-(1-\epsilon)^2=0$. The stochastic amplification occurs in the parameter region above this line.
• Figure 3: $\tau$ dependence of $C^{WZ}(\tau)$, $C^{ZW}(\tau)$, and $|C^{WZ}(\tau)|-|C^{ZW}(\tau)|$. The parameter values are $\epsilon=D^W=1/2$, $\beta=1$. The characteristic period is given by $T=2\pi/\sqrt{\epsilon \beta-(1-\epsilon)^2/2}\simeq10$.
• Figure 4: $\tau$ dependence of $C^{WZ}(\tau)$, $C^{ZW}(\tau)$, and $|C^{WZ}(\tau)|-|C^{ZW}(\tau)|$. The parameter values are $\epsilon=0.2$, $D^W=1/2$, $\beta=1$. In this case, the resonance condition $\epsilon \beta-(1-\epsilon)^2/2>0$ is not satisfied.