View source for Draft:Schamel method

{{AFC submission|d|v|u=SaschaXX|ns=118|decliner=Ldm1954|declinets=20240609112747|reason2=nn|ts=20240418102048}} <!-- Do not remove this line! -->

{{AFC comment|1=You are getting closer, but 
*At the moment almost all the refs include Schabel, which are primary sources. You need fully independent secondary sources by others.
*Many sections cite no sources 
*The lead is far too technical. You need a general intro for a high school student.
*WP is not a textbook or guide, see [[WP:!]]
*Why is this notable? You have not explained relevance. [[User:Ldm1954|Ldm1954]] ([[User talk:Ldm1954|talk]]) 11:27, 9 June 2024 (UTC)}}

----

{{Short description|An appropriate method for finding electrostatic equilibrium solutions}}
{{Draft topics|physics}}
{{AfC topic|stem}}

{{Multiple issues|{{Notability|date=March 2024}}
{{Additional citations|date=March 2024}}
{{Technical|date=March 2024}}}}
The '''Schamel method''' (S method) <ref name="myref1"/> is a method for finding [[Electrostatics|electrostatic]] [[equilibrium solution]]s  <math>\phi(x-v_0t)</math> for the Vlasov-Poisson system. It is an alternative to the method described by Bernstein, Greene and Kruskal <ref name="myref2"/>, as it first solves the Vlasov equations for the species involved and only in a 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 variety of electrostatic structures, including their phase velocities. These structures, also known under [[Bernstein–Greene–Kruskal modes]] or phase space [[electron]] and [[ion]] holes, or double layers, respectively, are ubiquitously found in collisionless [[Plasma (physics)|plasma]]s such as in the [[Earth]]'s [[magnetosphere]], in fusion machines, or in the laboratory.

With the S method, the distribution functions for electrons <math>f_e(x,v)=F_e(\epsilon_e,\sigma_e)</math> and ions <math>f_i(x,u)=F_i(\epsilon_i,\sigma_i)</math>, which solve the corresponding time-independent Vlasov equation,

<math>(v\partial_x +\phi'(x)\partial_v)f_e(x,v)=0</math> and  <math>(u\partial_x -\theta\phi'(x)\partial_u)f_i(x,u)=0</math>, respectively,

are described as functions of the [[constants of motion]]. Thereby the unperturbed plasma conditions are adequately taken into account. Here <math>\epsilon_e=\frac{v^2}{2}-\phi</math> and
<math>\epsilon_i=\frac{u^2}{2}-\theta(\psi-\phi)</math> are the single particle energies of electrons and ions, respectively ; <math>\sigma_e=\frac{v}{|v|}</math>  and  <math>\sigma_i=\frac{u}{|u|}</math> are the signs of the velocity of untrapped electrons and ions, respectively. Normalized quantities have been used, <math>\theta=\frac{T_e}{T_i}</math> and it is assumed that <math>0 \le\phi \le\psi</math>, where <math>\psi</math> 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, <math>\phi_{xx} = n_e - n_i=: - \mathcal {V} '(\phi;...) </math>, where <math>n_e,n_i</math> are the corresponding densities obtained by velocity integration of <math>F_e(\epsilon_e,\sigma_e)</math> and  <math>F_i(\epsilon_i,\sigma_i)</math>, is integrated whereby  the pseudo-potential <math>\mathcal{V}(\phi;...)</math> is introduced. The result is <math>\frac{\phi_x^2}{2} + \mathcal{V(\phi;...)} = 0</math>, which represents the pseudo-energy. In
<math>\mathcal{V}(\phi;...)</math> the points stand for the different trapping parameters <math>B_s, D_{1s}, D_{2s}, C_s, s=e,i</math>. Integration of the pseudo-energy results in <math>x(\phi)= \int_\phi^\psi \frac{d\xi}{\sqrt{-2\mathcal{V}(\xi;...)}} </math>, which yields by inversion the desired <math>\phi(x)</math>. In these expressions the canonical form of <math>\mathcal{V}(\phi;...)</math> is used already.

There are, however,  two further trapping parameters  
<math>\Gamma_s,s=e,i</math>, which are missing in the canonical pseudo-potential. In its extended previous version <math>\mathcal{V}_0(\phi;...,\Gamma_e,\Gamma_i,v_0)</math> they fell victim to the necessary constraint that the gradient of the potential <math>\phi(x)</math> vanishes at its maximum. This requirement is

<math>\mathcal{V}_0(\psi;...,\Gamma_e,\Gamma_i,v_0)=0</math> 

and leads to what is called the [[nonlinear dispersion relation (NDR)]]. It allows the phase speed <math>v_0 </math> of the structure to be determined in terms of the other parameters and to eliminate <math>v_0 </math> from <math>\mathcal{V}_0(\phi;...,\Gamma_e,\Gamma_i,v_0)</math> to obtain the canonical <math>\mathcal{V}(\phi;...)</math>. Due to the central role, the two expressions <math>\mathcal {V}(\phi;...)</math> and <math>\mathcal {V}_0(\phi;...,\Gamma_e,\Gamma_i ,v_0)</math>  are playing, the S method is sometimes also called Schamel's pseudo-potential method (SPP method). For more information, see references 1.<ref name="myref1"/> and 3<ref name="myref3"/>

A typical example for the canonical pseudo-potential and the NDR is given by<ref name="myref3"/>
<math>-\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))</math> 

and by 

<math>k_0^2 - \frac{1}{2}Z_r'(\frac{|v_D-v_0|}{\sqrt2}) - \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</math>

where <math>\varphi:=\phi/\psi</math> and <math>\delta:=m_e/m_i </math>, respectively.

These expressions are valid for a current-carrying, thermal background plasma described by Maxwellians (with a drift <math>v_D</math> between electrons and ions) and for the presence of the perturbative trapping scenarios <math>(\Gamma_s, B_s)</math>, s=e,i.

==Innovations achieved through the S method<ref name="myref3"/>==

===Failure of conventional wave theories===
Probably the most important achievement of the S method is the demonstration of the failure of linear wave theories in connection with stationary solutions. The reason for this is that coherence and trapping are interdependent and remain relevant even in the limit range of vanishingly small amplitudes.<ref name="myref3"/>

===Unlimited variety of potential structures <math>\phi(x)</math>===
The theory based on the S method <ref name="myref3"/> predicts an unlimited variety of potential structures of the solitary and cnoidal types, not just <math>sech^4(x) </math> or <math>e^{-x^2} </math> solitary structures or single harmonic waves extending the few discrete modes known from linear wave theories appreciably. Most of them appear in an undisclosed form, that is, <math>x(\phi)</math> cannot be represented by mathematically known functions and therefore neither can its inverse <math>\phi(x)</math>.

===Function space of solutions===
The function space of Vlasov-Poisson (VP) solutions is completely detached from the corresponding linear Vlasov-Poisson (lVP) space. There is no overlap and neither the linear Landau theory nor the Van Kampen theory can bridge this gap.

===Trapping scenarios - resonance broadening===
The chaotic behavior of the particles near resonance when captured by the wave is addressed through so-called trapping scenarios. Various expressions are used with corresponding parameters that stand for the different trapping channels. Two of them (<math>\Gamma_e</math>, <math>\Gamma_i</math>) have no effect on the structure because they are missing in the canonical form of <math>\mathcal{V}(\phi)</math>. However, they have an impact on the NDR as they  allow <math>v_0</math> to be shifted. This results in “resonance broadening” as predicted but not approved in the 1960s on the basis of [[Airy wave theory|linear wave theory]].

===Nonlinear dispersion relation - constraints===
As said, <math>\mathcal{V}_0(\psi;...,\Gamma_e,\Gamma_i,v_0)=0</math>, represents a closed connection between the various parameters and is therefore suitable for a unique determination the phase speed <math>v_0</math> in dependence on the other parameters. 

The conditions for the existence of these structures can moreover only be accessed and derived from the NDR, such as
<math>\theta > 3.5</math> for long-wavelength ion-acoustic structures or <math>\theta < 3.5</math> for solitary ion holes <ref name="myref4"/>,<ref name="myref5"/>
It should be noted also that in the special case of a single wave the NDR takes the form of a "thumb-teardrop" dispersion relation, which is therefore strictly nonlinear. 

===Slow electron  and ion acoustic modes===
In a two-component, Maxwellian background plasma there are two acoustic modes, the slow electron acoustic <math>v_0=1.307</math> and the slow ion acoustic
<math>u_0=1.307 \sqrt{\theta/\delta}</math> mode, where <math>\delta=m_e/m_i</math> is the mass ratio. They are strictly nonlinear, complement the linear wave spectrum and are the basis for the propagation of solitary electron and ion holes.

===Privileged solitary electron hole and Schamel equation===
The showpiece of the S method is the privileged solitary electron hole because it has the smoothest of all trapped particle distributions. It therefore has the best chance of surviving collisions in a weakly dissipative plasma environment <ref name="myref6"/>,<ref name="myref7"/>,<ref name="myref8"/>,<ref name="myref9"/>. It reads

<math>-\mathcal{V(\phi)}= B_e \frac{\phi^2}{2}(1-\sqrt{\frac{\phi}{\psi}}) ; \phi(x)= \psi sech^4(\frac{\sqrt B_e}{4} x) ; v_0= 1.307( 1-B_e+ \Gamma_e)</math>

Note that these expressions only contain the two parameters <math>B_e,\Gamma_e</math>, which are assumed to be small. It is a solution to the Schamel equation derived in <ref name="myref10"/> by the reduced scaling method using appropriate stretched coordinates.

===Solitary holes of different polarity===
By the appropriate selection of trapping scenarios the existence of solitary electron holes of positive <math>\phi(x)>0</math> and negative <math>\phi(x)<0</math> polarity was confirmed <ref name="myref3"/>,<ref name="myref11"/> . This meant that corresponding observations in the magnetosphere <ref name="myref12"/>could be explained theoretically for the first time.

===Marginal stability of single waves===
Contrary to some statements in the literature, there is currently no reliable linear stability theory for these structures, especially solitary structures.
An exception, however, are harmonic single wave structures, for which marginal stability is predicted regardless of the drift in a current-carrying plasma <ref name="myref3"/>. This statement  contradicts Landau's theory for corresponding linear single waves.

===Negative wave energy holes - anomalous transport===
The energy of a structured plasma associated with these nonlinear waves can typically be negative. These negative energy holes therefore act as attractors in this dynamical system, providing an abundance of secondary wave excitations that contribute to anomalous resistivity and heating in intermittent turbulence and anomalous transport.<ref name="myref7"/>,<ref name="myref8"/>,<ref name="myref9"/>

===Extensions===
Based on the S method, various extensions have been presented in the past.
In <ref name="myref13"/>  a generalization to relativistic electrons and the inclusion of magnetic field effects were presented. The structures, observed ubiquitously in circular particle accelerators and storage rings, could be explained in  <ref name="myref9"/> by use of the S method considering the focusing property of wall impedances. Quantum corrections of holes, investigated by <ref name="myref9"/> <ref name="myref14"/>, revealed the tendency to bring this extreme non-equilibrium state closer to thermodynamic equilibrium, an effect that is due to the tunneling of particles through the potential barrier. Finally reference <ref name="myref15"/> dealt with the nonlinear shielding of planar test charges and pointed out the important role of particle trapping.

== References ==
{
{reflist |
refs=<ref name="myref1">{{cite journal|first=H. |last=Schamel|journal=Plasma Physics|volume=14|year=1972|page=905|title=Stationary solitary, snoidal and sinusoidal ion acoustic waves|issue=10|doi=10.1088/0032-1028/14/10/002|bibcode=1972PlPh...14..905S}}</ref>
<ref name ="myref2">{{Cite journal|last1=Bernstein|first1=Ira B.|last2=Greene|first2=John M.|last3=Kruskal|first3=Martin D.|date=1957|title=Exact Nonlinear Plasma Oscillations|journal=Physical Review|volume=108|issue=3|pages=546–550|doi=10.1103/PhysRev.108.546|bibcode=1957PhRv..108..546B|hdl=2027/mdp.39015095115203|hdl-access=free}}</ref>
<ref name ="myref3">{{Cite journal |last=Schamel |first=H. |date=2023 |title=Pattern formation in Vlasov–Poisson plasmas beyond Landau caused by the continuous spectra of electron and ion hole equilibria |url=https://doi.org/10.1007/s41614-022-00109-w |journal=Reviews of Modern Plasma Physics |language=en |volume=7 |issue=1 |pages=11 |doi=10.1007/s41614-022-00109-w |issn=2367-3192|arxiv=2110.01433 |bibcode=2023RvMPP...7...11S }}</ref>
<ref name ="myref4">{{Cite journal|last1=Chen|first1=L.-J.|last2=Thouless|first2=D.|last3=Tang|first3=J.-M.|date=2004|title=Bernstein Greene Kruskal solitary waves in three-dimensional magnetized plasma|journal=Physical Review E|volume=69|issue=5|pages=55401|doi=10.1103/PhysRevE.69.055401|pmid=15244872 |arxiv=physics/0303021 |bibcode=2004PhRvE..69e5401C }}</ref>
<ref name ="myref5">{{Cite journal|last1=Schamel|first1=H.|last2=Chakrabarti|first2=N.|date=2023|title=Response to Comment on "On the Evolution Equations of Nonlinearly Permissible, Coherent Hole Structures Propagating Persistently in Collisionless Plasmas [Ann. Phys.(Berlin) 2023,2300102]|journal=Annalen der Physik |volume=2023|issue=3 |pages=2300441|doi=10.1002/andp.202300441 }}</ref>
<ref name ="myref6">{{Cite journal|last1=Korn|first1=J.|last2=Schamel|first2=H.|date=1996|title=Electron holes and their roles in the dynamics of current-carrying weakly collisional plasmas. Part 1. Immobile ions|journal=Journal of Plasma Physics |volume=56|issue=2 |pages=307–337|doi=10.1017/S0022377800019280 |bibcode=1996JPlPh..56..307K }}</ref>
<ref name ="myref7">{{Cite journal|last1=Korn|first1=J.|last2=Schamel|first2=H.|date=1996|title=Electron holes and their roles in the dynamics of current-carrying weakly collisional plasmas. Part 2. Mobile ions|journal=Journal of Plasma Physics |volume=56|issue=2 |pages=339–359|doi=10.1017/S0022377800019292 |bibcode=1996JPlPh..56..339K }}</ref>
<ref name ="myref8">{{Cite journal|last1=Schamel|first1=H.|last2=Korn|first2=J.|date=1996|title=Plasma transport in the presence of self-generated phase space structures|journal=Physica Scripta |volume=T63|pages=63–68|doi=10.1088/0031-8949/1996/T63/009 |bibcode=1996PhST...63...63S }}</ref>
<ref name ="myref9">{{Cite journal|last1=Luque|first1=A.|last2=Schamel|first2=H.|date=2005|title=Electrostatic trapping as a key to the dynamics of plasmas, fluids and other collective systems|journal=Physics Reports |volume=415|issue=5–6 |pages=261–359|doi=10.1016/j.physrep.2005.05.002 }}</ref>
<ref name="myref10">{{cite journal | last=Schamel | first=Hans | title=A modified Korteweg-de Vries equation for ion acoustic wavess due to resonant electrons | journal=Journal of Plasma Physics | publisher=Cambridge University Press (CUP) | volume=9 | issue=3 | year=1973 | issn=0022-3778 | doi=10.1017/s002237780000756x | pages=377–387| bibcode=1973JPlPh...9..377S | s2cid=124961361 }}</ref>
<ref name ="myref11">{{Cite journal|last1=Schamel|first1=H.|last2=Chakrabarti|first2=N.|date=2023|title=On the evolution equations of Nonlinearly Permissible, Coherent Hole Structures Propagating Persistently in Collisionless Plasmas|journal=Annals of Physics |volume=535|issue=8 |pages=2300102|doi=10.1002/andp.202300102 |bibcode=2023AnP...53500102S }}</ref>
<ref name ="myref12">{{Cite journal|last1=Holmes|first1=J.C.|last2=Ergun|first2=R.E.|last3=Newman|first3=D.L.|last4=Wilder|first4=F.D.|last5=Sturner|first5=A.P.|last6=Goodrich|first6=K.A.|last7=Torbert|first7=B.L.|last8=Giles|first8=B.L.|last9=Strangeway|first9=R.J.|last10=Burch|first10=J.L.|date=2018|title=Negative potential solitary structures in the magnetosheath with large parallel width|journal=Journal of Geophysical Research: Space Physics |volume=123|pages=132–145|doi=10.1002/2017/JA024890
|doi-broken-date=2024-04-19 }}</ref> 
<ref name ="myref13">{{Cite journal|last1=Eliasson|first1=B.|last2=Shukla|first2=P.K.|date=2006|title=Formation and dynamics of coherent structures involving phase-space vortices in plasmas|journal=Physics Reports |volume=422|issue=6 |pages=225–290|doi=10.1016/j.physrep.2005.10.003 |bibcode=2006PhR...422..225E }}</ref>
<ref name ="myref14">{{Cite journal|last1=Luque|first1=A.|last2=Schamel|first2=H.|last3=Fedele|first3=R.|date=2004|title=Quantum corrected electron holes|journal=Physics Letters A |volume=324|issue=2–3 |pages=185–192|doi=10.1016/j.physleta.2004.02.049 |arxiv=physics/0311126|bibcode=2004PhLA..324..185L }}</ref>
<ref name ="myref15">{{Cite journal|last1=Das|first1=N.|last2=Schamel|first2=H.|date=2005|title=Nonlinear shielding of planar test charge in one-dimensional Vlasov-Poisson plasma|journal=Journal of Plasma Physics |volume=71|issue=6 |pages=769–784 |doi=10.1017/S0022377805003776 |bibcode=2005JPlPh..71..769D }}</ref>
}}