pith. sign in

arxiv: 2604.12523 · v1 · submitted 2026-04-14 · ⚛️ physics.plasm-ph

Nonlinear Energy Transfer Analysis in Developing Plasma Turbulence

Pith reviewed 2026-05-10 14:26 UTC · model grok-4.3

classification ⚛️ physics.plasm-ph
keywords plasma turbulencenonlinear energy transferRayleigh-Taylor modesdrift-wave modesRitz methodKim methoddensity fluctuationsspectral analysis
0
0 comments X

The pith

Nonlinear interactions transfer energy from Rayleigh-Taylor modes to low-frequency drift-wave modes in plasma turbulence.

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

The paper examines nonlinear energy transfer among spectral components of plasma density fluctuations measured in the IMPED device, where signatures of Rayleigh-Taylor and drift-wave instabilities appear. It applies the existing Ritz and Kim computational methods, developed for single-field turbulence, to both simulations and experimental data to quantify transfers arising from quadratic coupling. The analysis shows energy moving from RT modes to a lower-frequency DW mode at different radial positions. These methods succeed in capturing the transport but only when the fluctuations satisfy statistical conditions such as adequate kurtosis and spatial stationarity. Understanding this cascade helps explain how turbulence develops and organizes in confined plasmas.

Core claim

Energy transfer analysis at different radial locations using the Ritz and Kim methods reveals the transfer of energy from RT modes to a comparatively low-frequency DW mode through quadratic coupling processes, demonstrating the capability of the methods to quantify spectral energy transport in the plasma turbulence.

What carries the argument

The Ritz and Kim methods, which compute nonlinear energy transfer rates among frequency components of a fluctuating field based on quadratic coupling in single-field turbulence models.

If this is right

  • The Ritz and Kim methods quantify spectral energy transport when applied to developing plasma turbulence.
  • Energy moves from higher-frequency RT modes to lower-frequency DW mode at multiple radial locations in the experiment.
  • Method validity requires checking higher-order moments and spatial stationarity of the measured fluctuations.

Where Pith is reading between the lines

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

  • The approach could be used to track how instabilities saturate in other plasma confinement experiments.
  • It offers a direct test for whether reduced turbulence models correctly reproduce quadratic energy cascades.
  • Similar spectral transfer diagnostics might apply to fluid or atmospheric turbulence where quadratic interactions dominate.

Load-bearing premise

The data must possess spatial stationarity and statistical properties such as sufficiently high kurtosis for the quadratic coupling calculations to apply reliably.

What would settle it

Applying the methods to data lacking nonlinear interactions, or to stationary data with low kurtosis, would produce zero or inconsistent transfer rates between the RT and DW modes.

Figures

Figures reproduced from arXiv: 2604.12523 by Daniel Raju, Kunal Singha, Lavkesh Lachhvani, Prabal Chattopadhyay, Rosh Roy, Sandip Das, Tanmay Karmakar.

Figure 1
Figure 1. Figure 1: Non-Linear Black Box Model The analytical linear and quadratic transfer functions defined in the following way [10]: Lf = 1.0 − 0.4 f 2 f 2 Nyq + i 0.8 f fNyq , (28) 8 [PITH_FULL_IMAGE:figures/full_fig_p008_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Power Spectrum of Input and Output signals after 5 [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: (a) Auto Bispectrum and (b) Cross Bispectrum after 5 [PITH_FULL_IMAGE:figures/full_fig_p010_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Convergence plot of (a) Auto and Cross Power Spectra (b) Auto and Cross Bispectra. [PITH_FULL_IMAGE:figures/full_fig_p010_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Benchmarking of Lf with Analytical Function [PITH_FULL_IMAGE:figures/full_fig_p011_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Benchmarking of Q f1,f2 f with Analytical Function. As discussed earlier, we use the input and output data at iteration = 5, for which the estimated transfer functions show good agreement with the analytical results. A small deviation is observed in the linear transfer function near zero frequency, and the contour plot of the estimated quadratic transfer function appears slightly noisier than its analytica… view at source ↗
Figure 7
Figure 7. Figure 7: (a)Probability distribution of input time series having [PITH_FULL_IMAGE:figures/full_fig_p012_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Skewness and Kurtosis plot with number of Iteration [PITH_FULL_IMAGE:figures/full_fig_p013_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Plot of Lf with Analytical Function at iteration = 7. 13 [PITH_FULL_IMAGE:figures/full_fig_p013_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Plot of Q f1,f2 f with Analytical Function at iteration = 7. In this study, using artificial data, we observed that the Ritz method fails to correctly estimate the transfer functions once the excess kurtosis exceeds approximately 0.40, which occurs after iteration = 5. As shown earlier, both the linear and quadratic transfer functions exhibit good agreement with the analytical functions at iteration = 5 i… view at source ↗
Figure 11
Figure 11. Figure 11: Simulation validation of the Kim method. Panels (a) and (b) show the real and imaginary [PITH_FULL_IMAGE:figures/full_fig_p017_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Panels (a) and (b) show the real and imaginary parts of the linear transfer function [PITH_FULL_IMAGE:figures/full_fig_p018_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Growth rate γf and nonlinear energy transfer Tf /2Pf as a function of frequency. The growth rate and energy transfer function are plotted in [PITH_FULL_IMAGE:figures/full_fig_p019_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Cross-correlation function (CCF) at a radial position of 5.76 cm showing a positive time delay [PITH_FULL_IMAGE:figures/full_fig_p021_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Radial variation of statistical parameters (mean, standard deviation, skewness, and kurtosis) [PITH_FULL_IMAGE:figures/full_fig_p021_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: Probability distribution functions of the measured density fluctuations (˜n [PITH_FULL_IMAGE:figures/full_fig_p022_16.png] view at source ↗
Figure 17
Figure 17. Figure 17: Comparison of power spectra of input(X) and output(Y) signals at two radial locations. [PITH_FULL_IMAGE:figures/full_fig_p023_17.png] view at source ↗
Figure 18
Figure 18. Figure 18: (a) Auto bicoherence analysis. (b) Density fluctuation spectra for ˜n [PITH_FULL_IMAGE:figures/full_fig_p024_18.png] view at source ↗
Figure 19
Figure 19. Figure 19: (a) Power transfer function (Tf (f1, f2)). (b) Energy transfer function (W f NL) for ˜n of DW and RT modes. From the two-dimensional contour plot of Tf (f1, f2) as shown in [PITH_FULL_IMAGE:figures/full_fig_p025_19.png] view at source ↗
Figure 20
Figure 20. Figure 20: (a) Power transfer function (Tf (f1, f2)). (b) Energy transfer function (W f NL) for ˜n of DW and RT modes. The accuracy of the nonlinear energy transfer estimation is assessed using the energy mismatch parameter Wmis (shown in the legends of [PITH_FULL_IMAGE:figures/full_fig_p026_20.png] view at source ↗
Figure 21
Figure 21. Figure 21: (a) Auto bicoherence analysis. (b) Density fluctuation spectra for ˜n [PITH_FULL_IMAGE:figures/full_fig_p027_21.png] view at source ↗
Figure 22
Figure 22. Figure 22: (a) Power transfer function (Tf (f1, f2)). (b) Energy transfer function (W f NL) for ˜n of DW and RT modes. From the two-dimensional contour plot of Tf (f1, f2) shown in [PITH_FULL_IMAGE:figures/full_fig_p028_22.png] view at source ↗
Figure 23
Figure 23. Figure 23: (a) Power transfer function (Tf (f1, f2)). (b) Energy transfer function (W f NL) for ˜n of DW and RT modes. and validity of the Ritz method at this radial position. In summary, at this radial location the dominant interaction shows that the RT mode at 11.3 kHz transfers energy to the DW mode at 2.5 kHz as well as to the RT mode at 8.8 kHz. In the next section, we present the growth-rate analysis at the ra… view at source ↗
Figure 24
Figure 24. Figure 24: Growth rate spectrum γf as a function of frequency. The dashed line indicates γf = 0. Here, we provide a qualitative discussion of the growth rate plot rather than a quantitative analysis. This approach allows us to focus on the overall trends, such as which modes are growing or damping and how their behavior correlates with the energy transfer between modes, without delving into exact numerical values or… view at source ↗
read the original abstract

Energy transfer among various spectral components of fluctuating physical parameters in plasma occurs due to the nonlinear interactions, but these effects are typically not captured by the traditional linear spectral methods. Plasma density fluctuations measured in the Inverse Mirror Plasma Experimental Device (IMPED) have signatures of nonlinear mode interactions among various instability modes, i.e. Rayleigh-Taylor (RT) and Drift-Wave (DW) modes. In this paper, the energy transfer among these modes as a result of nonlinear wave interactions (through the quadratic coupling processes) have been investigated in detail. The existing computational methods for single field turbulence model such as Ritz method and Kim method have been explored to understand the turbulence dynamics. Both methods are applied and validated in simulation as well as experimental data from IMPED for developing plasma turbulence. We find that the validity and applicability of the methods depend on the statistical nature of the data, particularly higher-order moments such as kurtosis, and on spatial stationarity. Energy transfer analysis at different radial locations using these methods reveals the transfer of energy from RT modes to a comparatively low-frequency DW mode, demonstrating the capability of the method to quantify spectral energy transport in the plasma turbulence.

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 paper applies the Ritz and Kim quadratic-coupling methods (from single-field turbulence theory) to density fluctuation measurements in the IMPED device to quantify nonlinear energy transfer among Rayleigh-Taylor (RT) and drift-wave (DW) modes in developing plasma turbulence. Both methods are validated on simulations and on the experimental data; the central result is directional energy transfer from RT modes to a lower-frequency DW mode at multiple radial locations. The abstract explicitly conditions applicability on data kurtosis and spatial stationarity.

Significance. If the stationarity and statistical assumptions are satisfied, the work supplies a concrete, quantitative example of spectral energy transport in a laboratory plasma, extending existing methods to a developing-turbulence regime. The dual validation on simulation and experiment is a positive feature.

major comments (2)
  1. [Abstract / Results] Abstract and Results section: the manuscript states that 'the validity and applicability of the methods depend on the statistical nature of the data, particularly higher-order moments such as kurtosis, and on spatial stationarity,' yet reports neither measured kurtosis values nor any stationarity tests (e.g., time-series segmentation or cross-spectral consistency checks) for the IMPED density time series at the different radial locations. Because the directional energy-transfer claim rests on the quantitative applicability of the Ritz/Kim estimators, this omission is load-bearing.
  2. [Methods / Validation] Methods / Validation section: while the abstract asserts that both methods 'were validated on simulations and experiment,' the text provides no explicit equations for the transfer-rate estimators, no error-propagation formulas, and no data-selection criteria (e.g., stationarity windows or kurtosis thresholds) used to accept or reject individual time series. Without these, the experimental transfer directions cannot be independently verified.
minor comments (1)
  1. [Methods] Notation for the quadratic-coupling coefficients and the definition of the energy-transfer direction (sign convention) should be stated once in a dedicated subsection for clarity.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful and constructive review. The comments correctly identify gaps in the presentation of statistical validation and methodological details that are important for assessing the applicability of the Ritz and Kim methods. We address each point below and have revised the manuscript to incorporate the requested information.

read point-by-point responses
  1. Referee: [Abstract / Results] Abstract and Results section: the manuscript states that 'the validity and applicability of the methods depend on the statistical nature of the data, particularly higher-order moments such as kurtosis, and on spatial stationarity,' yet reports neither measured kurtosis values nor any stationarity tests (e.g., time-series segmentation or cross-spectral consistency checks) for the IMPED density time series at the different radial locations. Because the directional energy-transfer claim rests on the quantitative applicability of the Ritz/Kim estimators, this omission is load-bearing.

    Authors: We agree that the specific kurtosis values and stationarity test results should have been reported to allow readers to evaluate the applicability of the methods. Although the abstract notes the dependence on these factors, the numerical values and test outcomes were omitted from the original text. In the revised manuscript we have added a dedicated paragraph in the Results section that reports the measured kurtosis for the density fluctuations at each radial location (all values lie within the range consistent with the method assumptions) together with the outcomes of stationarity tests performed by segmenting the time series and verifying cross-spectral consistency. These additions directly support the directional energy-transfer conclusions. revision: yes

  2. Referee: [Methods / Validation] Methods / Validation section: while the abstract asserts that both methods 'were validated on simulations and experiment,' the text provides no explicit equations for the transfer-rate estimators, no error-propagation formulas, and no data-selection criteria (e.g., stationarity windows or kurtosis thresholds) used to accept or reject individual time series. Without these, the experimental transfer directions cannot be independently verified.

    Authors: We acknowledge that the original submission did not include the explicit mathematical expressions for the Ritz and Kim transfer-rate estimators, the associated error-propagation formulas, or the precise data-selection criteria. This omission limits independent verification. The revised Methods section now contains the full equations for the quadratic-coupling estimators, the derivation of the error estimates, and the quantitative thresholds (kurtosis bounds and stationarity-window lengths) applied to accept or reject time series. The validation subsections for both the simulations and the IMPED data have been expanded to reference these criteria explicitly. revision: yes

Circularity Check

0 steps flagged

No circularity: established methods applied to independent experimental data

full rationale

The paper applies the pre-existing Ritz and Kim methods (described as established computational techniques for single-field turbulence) to compute quadratic-coupling energy transfers from measured density fluctuations in the IMPED device. The reported directional transfer from RT modes to low-frequency DW modes is obtained directly as the output of these methods on the time-series data, after separate validation on simulations. No step in the provided derivation reduces by construction to a fitted parameter, self-definition, or renaming of the input data; the central claim does not depend on a load-bearing self-citation chain or imported uniqueness theorem from the authors' prior work. The analysis remains self-contained as an application to external measurements.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

Analysis rests on standard spectral assumptions for quadratic nonlinear coupling and on the data satisfying the statistical conditions (low kurtosis, spatial stationarity) explicitly stated as necessary for method validity. No new entities or free parameters are introduced in the abstract.

axioms (2)
  • domain assumption Nonlinear interactions occur through quadratic coupling processes among spectral components
    Invoked when applying Ritz and Kim methods to compute energy transfer.
  • domain assumption Data statistics (kurtosis) and spatial stationarity determine method applicability
    Explicitly stated as conditions that affect validity in the abstract.

pith-pipeline@v0.9.0 · 5521 in / 1291 out tokens · 31176 ms · 2026-05-10T14:26:11.408245+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

37 extracted references · 37 canonical work pages

  1. [1]

    W. Horton. Drift waves and transport in magnetized plasmas.Reviews of Modern Physics, 71(3):735–778, 1999. 30 S. Das, L. Lachhvani and Others Nonlinear energy transfer in plasma turbulence

  2. [2]

    C. M. Ritz, E. J. Powers, and R. D. Bengtson. Experimental measurement of nonlinear energy transfer in a turbulent plasma.Phys. Fluids B: Plasma Phys., 1(1):153–159, 1989

  3. [3]

    D. E. Smith, E. J. Powers, and G. S. Caldwell. Fast-fourier-transform spectral-analysis techniques as a plasma fluctuation diagnostic tool.IEEE Transactions on Plasma Science, PS-2:261–272, 1974

  4. [4]

    D. E. Smith, E. J. Powers, and G. S. Caldwell. Nonlinear three-wave interactions in plasma turbu- lence.IEEE Transactions on Plasma Science, PS-2:261–267, 1974

  5. [5]

    D. E. Smith and E. J. Powers. Experimental investigation of nonlinear wave interactions in a turbulent plasma.Physics of Fluids, 16:1373–1380, 1973

  6. [6]

    Elgar and R

    S. Elgar and R. T. Guza. Observations of bispectra of ocean surface gravity waves.Journal of Fluid Mechanics, 161:425–448, 1985

  7. [7]

    Kim and Edward J

    Young C. Kim and Edward J. Powers. Digital bispectral analysis and its applications to nonlinear wave interactions.IEEE Trans. Plasma Sci., PS-7(2):120–131, 1979

  8. [8]

    R. W. Miksad, F. L. Jones, and E. J. Powers. Energy transfer in drift-wave instabilities.Physics of Fluids, 26:1402–1410, 1983

  9. [9]

    Y. C. Kim, J. M. Beall, E. J. Powers, and R. W. Miksad. Bispectrum and nonlinear wave coupling. Physics of Fluids, 23:258–263, 1980

  10. [10]

    Ch. P. Ritz and E. J. Powers. Estimation of nonlinear transfer functions for fully developed turbu- lence.Physica D: Nonlinear Phenomena, 20(2-3):320–334, 1986

  11. [11]

    Ch. P. Ritz, E. J. Powers, R. W. Miksad, and R. S. Solis. Nonlinear spectral dynamics of a transitioning flow.Physics of Fluids, 31(12):3577–3588, 1988

  12. [12]

    An improved method for bispectral analysis of nonlinear wave interaction system.Physica Scripta, 95(5):055202, 2020

    Yu-Hang Shen, Jia-Yin Li, Tian Li, and Jia Li. An improved method for bispectral analysis of nonlinear wave interaction system.Physica Scripta, 95(5):055202, 2020

  13. [13]

    J. S. Kim, R. D. Durst, R. J. Fonck, E. Fernandez, A. Ware, and P. W. Terry. Technique for the experimental estimation of nonlinear energy transfer in fully developed turbulence.Physics of Plasmas, 3(11):3998–4009, 1996

  14. [14]

    K. I. Kim, E. J. Powers, C. P. Ritz, R. W. Miksad, and F. J. Fischer. Bispectral analysis of nonlinear drift oscillations of moored vessels.IEEE Journal of Oceanic Engineering, OE-12(4):568–575, 1987

  15. [15]

    Anderson, E-j Kim, B

    J. Anderson, E-j Kim, B. Hnat, and T. Rafiq. Elucidating plasma dynamics in hasegawa-wakatani turbulence by information geometry.Physics of Plasmas, 27(2):022307, 2020

  16. [16]

    Plasma edge turbulence.Physical Review Letters, 50(9):682, 1983

    Akira Hasegawa and Masahiro Wakatani. Plasma edge turbulence.Physical Review Letters, 50(9):682, 1983

  17. [17]

    Krall and Alvin W

    Nicholas A. Krall and Alvin W. Trivelpiece.Principles of Plasma Physics. McGraw-Hill, New York, 1973

  18. [18]

    P. W. Terry and P. H. Diamond. Resistivity-gradient-driven turbulence.Physics of Fluids, 28:1419– 1430, 1985

  19. [19]

    D. E. Smith, E. J. Powers, and G. S. Caldwell. Fast-fourier-transform spectral-analysis techniques as a plasma fluctuation diagnostic tool.IEEE Transactions on Plasma Science, PS-2:261–270, December 1974

  20. [20]

    C. L. Nikias and A. P. Petropulu.Higher-Order Spectra Analysis: A Nonlinear Signal Processing Framework. Prentice Hall, Englewood Cliffs, New Jersey, 1993

  21. [21]

    Pseudo-three-dimensional turbulence in magnetized nonuniform plasma.Physical Fluids, 21(1):87–92, 1978

    Akira Hasegawa and Kunioki Mima. Pseudo-three-dimensional turbulence in magnetized nonuniform plasma.Physical Fluids, 21(1):87–92, 1978

  22. [22]

    M. H. Freilich and R. T. Guza. Nonlinear wave interactions in ocean surface gravity waves.Philo- sophical Transactions of the Royal Society of London. Series A, 311:1–26, 1984. 31 S. Das, L. Lachhvani and Others Nonlinear energy transfer in plasma turbulence

  23. [23]

    Hasegawa, C

    A. Hasegawa, C. G. Maclennan, and Y. Kodama. Nonlinear behavior and turbulence spectra of drift waves.Physics of Fluids, 22:2122–2129, 1979

  24. [24]

    M. D. Millionshchikov. On the theory of homogeneous isotropic turbulence.Doklady Akademii Nauk SSSR, 32:611–614, 1941

  25. [25]

    P. W. Terry, P. H. Diamond, K. C. Shaing, L. Garcia, and B. A. Carreras. The spectrum of resistivity gradient driven turbulence. Technical Report IFSR-174, Institute for Fusion Studies, Austin, Texas, 1986

  26. [26]

    Xia and M

    H. Xia and M. G. Shats. Spectral energy transfer and generation of turbulent structures in toroidal plasma.Physics of Plasmas, 11(2):561–571, 2004

  27. [27]

    Bendat and Allan G

    Julius S. Bendat and Allan G. Piersol.Random Data: Analysis and Measurement Procedures. Wiley, 4th edition, 2010

  28. [28]

    Spiegel and Larry J

    Murray R. Spiegel and Larry J. Stephens.Schaum’s Outline of Statistics. McGraw-Hill Education, 6th edition, 2018

  29. [29]

    Inverse mirror plasma experimental device (imped)–a magnetized linear plasma device for wave studies.Journal of Plasma Physics, 81(2):345810203, 2015

    Sayak Bose, PK Chattopadhyay, J Ghosh, S Sengupta, YC Saxena, and R Pal. Inverse mirror plasma experimental device (imped)–a magnetized linear plasma device for wave studies.Journal of Plasma Physics, 81(2):345810203, 2015

  30. [30]

    Inverse mirror plasma experimental device—a new magnetized linear plasma device with wide operating range.Review of Scientific Instruments, 86(6), 2015

    Sayak Bose, Manjit Kaur, PK Chattopadhyay, J Ghosh, and YC Saxena. Inverse mirror plasma experimental device—a new magnetized linear plasma device with wide operating range.Review of Scientific Instruments, 86(6), 2015

  31. [31]

    Rosh Roy, Tanmay Karmakar, Lavkesh Lachhvani, Prabal Chattopadhyay, Abhijit Sen, and Sayak Bose. An experimental study of the existence regions and non-linear interactions of drift wave and kelvin–helmholtz instabilities in a linear magnetized plasma.Physics of Plasmas, 32(2), 2025

  32. [32]

    Rosh Roy, Tanmay Karmakar, Lavkesh Lachhvani, Bhoomi Khodiyar, Prabal Chattopadhyay, and Abhijit Sen. An experimental study of the rayleigh–taylor and kelvin–helmholtz instabilities, and their nonlinear interactions in a linear magnetized plasma column.Physics of Plasmas, 32(12), 2025

  33. [33]

    On the transition to drift turbulence in a magnetized plasma column.Physics of Plasmas, 12(5), 2005

    MJ Burin, GR Tynan, GY Antar, NA Crocker, and C Holland. On the transition to drift turbulence in a magnetized plasma column.Physics of Plasmas, 12(5), 2005

  34. [34]

    M. Xu, G. R. Tynan, C. Holland, Z. Yan, S. H. Muller, and J. H. Yu. Study of nonlinear spectral energy transfer in frequency domain.Physics of Plasmas, 16(4):042312, 2009

  35. [35]

    M. Xu, G. R. Tynan, C. Holland, Z. Yan, S. H. Muller, and J. H. Yu. Fourier-domain study of drift turbulence driven sheared flow in a laboratory plasma.Physics of Plasmas, 17(3):032311, 2010

  36. [36]

    P. Manz, M. Ramisch, U. Stroth, V. Naulin, and B. D. Scott. Bispectral experimental estimation of the nonlinear energy transfer in two-dimensional plasma turbulence.Plasma Physics and Controlled Fusion, 50(3):035008, 2008

  37. [37]

    S.-I. Itoh, K. Itoh, Y. Nagashima, and Y. Kosuga. On the application of cross bispectrum and cross bicoherence.Plasma and Fusion Research, 12:1101003, 2017. 32