Nonlinear dispersion relation information

This article has multiple issues. Please help improve it or discuss these issues on the talk page. (Learn how and when to remove these template messages) The topic of this article may not meet Wikipedia's general notability guideline. Please help to demonstrate the notability of the topic by citing reliable secondary sources that are independent of the topic and provide significant coverage of it beyond a mere trivial mention. If notability cannot be shown, the article is likely to be merged, redirected, or deleted. Find sources: "Nonlinear dispersion relation" – news · newspapers · books · scholar · JSTOR (May 2024) (Learn how and when to remove this message) This article needs additional citations for verification. Please help improve this article by adding citations to reliable sources. Unsourced material may be challenged and removed. Find sources: "Nonlinear dispersion relation" – news · newspapers · books · scholar · JSTOR (May 2024) (Learn how and when to remove this message) This article may be too technical for most readers to understand. Please help improve it to make it understandable to non-experts, without removing the technical details. (May 2024) (Learn how and when to remove this message) This article has no lead section. Please improve this article by adding one in your own words. (June 2024) (Learn how and when to remove this message) (Learn how and when to remove this message)

A nonlinear dispersion relation (NDR) is a relation that assigns the correct phase velocity ${\displaystyle v_{0}}$ to a nonlinear wave structure. As an example of how diverse and intricate the underlying description can be, we deal with plane electrostatic wave structures ${\displaystyle \phi (x-v_{0}t)}$ which propagate with ${\displaystyle v_{0}}$ in a collisionless plasma. Such structures are ubiquitous, for example in the magnetosphere of the Earth, in fusion reactors or in the laboratory. Correct means that this must be done according to the governing equations, in this case the Vlasov-Poisson system, and the conditions prevailing in the plasma during the wave formation process. This means that special attention must be paid to the particle trapping processes acting on the resonant electrons and ions, which requires phase space analyses. Since the latter is stochastic, transient and rather filamentary in nature, the entire dynamic trapping process eludes mathematical treatment, so that it can be adequately taken into account “only” in the asymptotic, quiet regime of wave generation, when the structure is close to equilibrium.

This is where the pseudo-potential method in the version of Schamel,^[1][2] also known as S method, comes into play, which is an alternative to the method described by Bernstein, Greene and Kruskal.^[3] In the S method the Vlasov equations for the species involved are first solved and only in the second step Poisson's equation to ensure self-consistency. The S method is generally considered the preferred method because it is best suited to describing the immense diversity of electrostatic structures, including their phase velocities. These structures are also known under Bernstein–Greene–Kruskal modes or phase space electron and ion holes, or double layers, respectively. With the S method, the distribution functions for electrons ${\displaystyle f_{e}(x,v)=F_{e}(\epsilon _{e},\sigma _{e})}$ and ions ${\displaystyle f_{i}(x,u)=F_{i}(\epsilon _{i},\sigma _{i})}$, which solve the corresponding time-independent Vlasov equation,

${\displaystyle (v\partial _{x}+\phi '(x)\partial _{v})f_{e}(x,v)=0}$ and ${\displaystyle (u\partial _{x}-\theta \phi '(x)\partial _{u})f_{i}(x,u)=0}$, respectively,

are described as functions of the constants of motion. Thereby the unperturbed plasma conditions are adequately taken into account. Here ${\displaystyle \epsilon _{e}={\frac {v^{2}}{2}}-\phi }$ and ${\displaystyle \epsilon _{i}={\frac {u^{2}}{2}}-\theta (\psi -\phi )}$ are the single particle energies of electrons and ions, respectively ; ${\displaystyle \sigma _{e}={\frac {v}{|v|}}}$ and ${\displaystyle \sigma _{i}={\frac {u}{|u|}}}$ are the signs of the velocity of untrapped electrons and ions, respectively. Normalized quantities have been used, ${\displaystyle \theta ={\frac {T_{e}}{T_{i}}}}$ and it is assumed that ${\displaystyle 0\leq \phi \leq \psi }$, where ${\displaystyle \psi }$ is the amplitude of the structure.

Of particular interest is the range in phase space where particles are trapped in the potential wave trough. This area is (partially) filled via the stochastic particle dynamics during the previous creation process and is expressed and parameterized by so-called trapping scenarios.

In the second step, Poisson's equation, ${\displaystyle \phi _{xx}=n_{e}-n_{i}=:-{\mathcal {V}}'(\phi ;...)}$, where ${\displaystyle n_{e},n_{i}}$ are the corresponding densities obtained by velocity integration of ${\displaystyle F_{e}(\epsilon _{e},\sigma _{e})}$ and ${\displaystyle F_{i}(\epsilon _{i},\sigma _{i})}$, is integrated whereby the pseudo-potential ${\displaystyle {\mathcal {V}}(\phi ;...)}$ is introduced. The result is ${\displaystyle {\frac {\phi _{x}^{2}}{2}}+{\mathcal {V(\phi ;...)}}=0}$, which represents the pseudo-energy. In ${\displaystyle {\mathcal {V}}(\phi ;...)}$ the points stand for the different trapping parameters ${\displaystyle B_{s},D_{1s},D_{2s},C_{s},s=e,i}$. Integration of the pseudo-energy results in ${\displaystyle x(\phi )=\int _{\phi }^{\psi }{\frac {d\xi }{\sqrt {-2{\mathcal {V}}(\xi ;...)}}}}$, which yields by inversion the desired ${\displaystyle \phi (x)}$. In these expressions the canonical form of ${\displaystyle {\mathcal {V}}(\phi ;...)}$ is used already.

There are, however, two further trapping parameters ${\displaystyle \Gamma _{s},s=e,i}$, which are missing in the canonical pseudo-potential. In its extended previous version ${\displaystyle {\mathcal {V}}_{0}(\phi ;...,\Gamma _{e},\Gamma _{i},v_{0})}$ they fell victim to the necessary constraint that the gradient of the potential ${\displaystyle \phi (x)}$ vanishes at its maximum. This requirement is

${\displaystyle {\mathcal {V}}_{0}(\psi ;...,\Gamma _{e},\Gamma _{i},v_{0})=0}$

and leads to the term in question, the nonlinear dispersion relation (NDR). It allows the phase speed ${\displaystyle v_{0}}$ of the structure to be determined in terms of the other parameters and to eliminate ${\displaystyle v_{0}}$ from ${\displaystyle {\mathcal {V}}_{0}(\phi ;...,\Gamma _{e},\Gamma _{i},v_{0})}$ to obtain the canonical ${\displaystyle {\mathcal {V}}(\phi ;...)}$. Due to the central role, the two expressions ${\displaystyle {\mathcal {V}}(\phi ;...)}$ and ${\displaystyle {\mathcal {V}}_{0}(\phi ;...,\Gamma _{e},\Gamma _{i},v_{0})}$ are playing, the S method is sometimes also called Schamel's pseudo-potential method (SPP method). For more information, see references 1^[1] and 2.^[2]

A typical example for the canonical pseudo-potential ${\displaystyle {\mathcal {V}}(\phi ;...)}$ and the NDR is given by^[1] ${\displaystyle -{\mathcal {V}}(\phi ;B_{e},B_{i})/\psi ^{2}={\frac {k_{0}^{2}}{2}}\varphi (1-\varphi )+B_{e}{\frac {\varphi ^{2}}{2}}(1-{\sqrt {\varphi }})+B_{i}{\frac {\theta ^{3/2}}{2}}(1-(1-\varphi )^{5/2}-{\frac {1}{2}}\varphi (5-3\varphi ))}$

and by

${\displaystyle k_{0}^{2}-{\frac {1}{2}}Z_{r}'({\frac {|v_{D}-v_{0}|}{\sqrt {2}}})-{\frac {\theta }{2}}Z_{r}'({\sqrt {\frac {\theta }{2\delta }}}v_{0})=B_{e}+{\frac {3}{2}}\theta ^{3/2}B_{i}-\Gamma _{e}-\Gamma _{i}}$

where ${\displaystyle \varphi :=\phi /\psi }$ and ${\displaystyle \delta :=m_{e}/m_{i}}$, respectively.

These expressions are valid for a current-carrying, thermal background plasma described by Maxwellians (with a drift ${\displaystyle v_{D}}$ between electrons and ions) and for the presence of the perturbative trapping scenarios ${\displaystyle (\Gamma _{s},B_{s})}$, s=e,i.

A well known example is the Thumb-Teardrop dispersion relation, which is valid for single harmonic waves and is given as a simplified version of the NDR above with zero trapping parameters and a vanishing drift. It reads

${\displaystyle k_{0}^{2}-{\frac {1}{2}}Z_{r}'({\frac {v_{0}}{\sqrt {2}}})-{\frac {\theta }{2}}Z_{r}'({\sqrt {\frac {\theta }{2\delta }}}v_{0})=0}$.

It has been thoroughly discussed in ^[4] but mistakenly^[5] as a linear dispersion relation.

A plot of ${\displaystyle \omega _{0}:=k_{0}v_{0}}$ for a more general type of structures, the periodic cnoidal electron holes,^[6] is presented in Fig 1. for the case ${\displaystyle B_{e}=B,B_{i}=0}$ and immobile ions (${\displaystyle \theta =0}$) showing the effect of the electron trapping parameter ${\displaystyle B}$ for positive and negative values.

Finally, it should be mentioned that an NDR has the nice property that it remains valid even if the potential ${\displaystyle \phi (x)}$ has an undisclosed form, i.e. can no longer be described by mathematically known functions.^[1]