arxiv: v1 [physics.plasm-ph] 31 Jul 2018

Transcription

1 arxv: v1 [physcs.plasm-ph] 31 Jul 218 Consstency n Drft-ordered Flud Equatons Jakob Gath and Matthas Wesenberger Department of Physcs, Techncal Unversty of Denmark, DK-28 Kgs. Lyngby, Denmark jagath@fysk.dtu.dk, mattw@fysk.dtu.dk Abstract We address several concerns related to the dervaton of drft-ordered flud equatons. Startng from a fully Gallean nvarant flud system, we show how consstent sets of perturbatve drft-flud equatons n the case of a sothermal collsonless flud can be obtaned. Treatng all the dynamcal felds on equal footng n the sngular-drft expanson, we show under what condtons a set of perturbatve equatons can have a non-trval quas-neutral lmt. We gve a sutable perturbatve setup where we provde the full set of perturbatve equatons for obtanng the frst-order corrected felds and show that all the constants of moton are preserved at each order. Wth the dynamcal feld varables under perturbatve control, we subsequently provde a quanttatve analyss by means of numercal smulatons. Wth drect access to frst-order correctons the convergence propertes are addressed for dfferent regmes of parameter space and the valdty of the frst-order approxmaton s dscussed n the three settngs: cold ons, hot ons and fnte charge densty.

2 Contents 1 Introducton 1 2 Non-relatvstc flud equatons Perturbatve setup Smulaton Intal condtons Results Dscusson 16 A Non-relatvstc lmt and Gallean nvarance 18 B Conserved charges 19 C Alternatve perturbatve setups 2 References 21 1 Introducton The scrape-off layer n magnetcally confned plasmas s typcally treated as a cold and collson-rch regon where the turbulence s characterstcally low frequency and long wavelength n nature. In ths regme, flud models are arguably wdely applcable [1 4], as they are sutable for studyng the stablty and transport propertes of the plasma on scales much longer than the on gyroradus and much lower than the cyclotron frequency. In fact, many key characterstcs of the scrape-off layer are captured by the drft-flud theory (see e.g. [5]). For drft-flud models, t s common to treat the flud-momentum equaton n a sngular perturbaton expanson that makes t possble to solve explctly for the perpendcular components of the flud velocty. Together wth the addtonal flud conservaton equatons ths consttutes the so-called drft-ordered flud equatons. Perturbatve flud equatons derved n ths way are used n many plasma physcs smulaton codes [6 8]. However, there are several fundamental concerns wth respect to the consstency of such perturbatve schemes: () In the standard dervaton of drft-ordered flud equatons, the dynamcal feld varables are (questonably) rarely expanded on equal footng. It s therefore mpossble to systematcally check how well the approxmaton represents the true underlyng system as there s no noton of order. To control the expanson n a coupled non-lnear system, all dynamcal felds should n general be treated on equal footng such that the correctons to the dynamcs appear n an ordered scheme. In ths way, one can be montor the correctons and estmate whether a gven order of approxmaton has converged to a suffcent degree. () A onzed plasma s prmarly governed by the electromagnetc force. In many cases, one s therefore nterested n systems where (at least a part of) the electromagnetc felds are selfconsstently obtaned alongsde the non-relatvstc flud varables, e.g. the electrc feld sustaned by a plasma n a background magnetc feld. However, electromagnetsm s a Lorentz nvarant theory and therefore a sutable non-relatvstc lmt of the Maxwell equatons s requred to obtan a coupled system (wth the non-relatvstc flud equatons) that s fully nvarant under Gallean transformatons. () One of the defnng crtera for an onzed gas to be characterzed as a plasma s to have a suffcently hgh densty such that local concentratons of charge s shelded out at very short dstances. In such cases, the plasma can be sad to exst n a quas-neutral state, where the denstes of the charged speces are equlbrated. However, all nterestng electromagnetc phenomena ultmately 1

3 reles on separaton of charge. There exsts therefore a subtle nterplay between the quas-neutral property and the dynamcs of the plasma. For ths reason, t s mportant that the quas-neutral assumpton s mplemented on the same footng along wth any other assumpton such that t becomes a natural consequence of the perturbatve setup and comes under perturbatve control. In general, f an assumpton s mplemented before treatng the full underlyng system of equatons, the perturbatve nature of the assumpton s lost. (v) Treatng the flud-momentum equaton as a sngular perturbaton problem allows for a separaton between the fast and slow dynamcs and the constructon of the drft expanson. However, snce the momentum equaton s now degenerate t can not possble satsfy all ntal condtons that can be specfed for the full underlyng flud system (.e. the soluton space s naturally restrcted). The ntal condtons wll nstead be completely determned by the ntal condtons for the other dynamcal felds. One should therefore make sure that the perturbaton problem under consderaton s not ll-posed. Before proceedng, we start wth a general dscusson about the quas-neutral assumpton n drft-ordered flud models. In the general drft-ordered setup, the perpendcular components of the velocty are decomposed nto three prmary terms; the E B, the damagnetc, and the polarzatondrft velocty. The latter two manfest the separaton of charge n the system. However, under a strct quas-neutral assumpton no charge separaton n the system s n prncple allowed whch gves rse to some nterestng constrants on the possble orderng of the drft veloctes to ensure the consstency of the system of equatons. Typcally, the E B-drft velocty s taken to be the man contrbuton to the velocty feld (and advecton of densty) whle the polarzaton-drft velocty by constructon of the sngular perturbaton expanson must appear at hgher order. Now, n many drft-flud models the damagnetc-drft velocty s typcally assumed to be of the same magntude as the E B-drft velocty. However, n a proper perturbatve setup ths would requre the charge separaton to be treated at the same order as the E B-drft velocty whch s n drect contradcton wth the assumpton that the separaton of charge s neglgble. Such cases can only make sense for suffcently dlute plasmas where the quas-neutral assumpton can be lfted (or f the damagnetcdrft velocty does no convecton as n a constant magnetc feld confguraton [2]). Otherwse, t wll have to appear at the succeedng order where t n conjuncton wth the polarzaton-drft velocty provdes a tme evoluton for the electrc potental. In fact, the quas-neutral assumpton reduces the charge-current conservaton equaton to the dvergenceless current constrant and only due to the polarzaton-drft velocty ths appears as a dynamcal constrant equaton. Alternatvely, t s natural to attempt to push the convecton due to the damagnetc-drft velocty to the next order as t s arguably small (n a nearly constant magnetc feld). However, t turns out that ths wll push the polarzaton-drft velocty to an even hgher order reducng the current conservaton equaton to a purely spatal constrant that can only be satsfed by a constant magnetc feld or zero pressure gradents. We go through the steps of obtanng a consstent perturbatve expanson of a smple collsonless sothermal flud governed by the Lorentz force n Sec. 2 that address all of the above lsted ponts. Subsequently, we fnd solutons to the perturbatve system usng the numercal tools avalable n the Feltor lbrary [9] n the case of a two speces flud n Sec. 3. Through the numercal solutons we perform a qunttatve analyss and estmate the physcal regme of valdty of the presented perturbatve system. We nclude three appendces. In App. A, we show how a consstent nonrelatvstc lmt of the Maxwell equatons can be obtaned. In App. B, we provde some detals on the conserved charges of the forced flud system presented n Sec. 2. In App. C, two alternatve perturbatve setups vald for dlute plasmas (n contrast wth the equatons gven n Sec. 2.1) are gven. 2 Non-relatvstc flud equatons A charged non-relatvstc flud n (3+1)-dmensons s consttuted by the contnuty, momentum, energy and and charge-current conservaton equatons. Here, we wll consder a collsonless sother- 2

4 mal flud system wthout dsspatve terms of charged partcles governed by the Lorentz force, t n a + (n a v a ) =, m a d t v a = T a log(n a )+q a ( φ+v a B). (1) The latn subscrpt a labels the speces of the partcle and bold symbols are spatal vectors. The total dervatvesdefnedbyd t t +v a. Thephyscalquanttes assocated wththendvdualspeces are the partcle mass m a and charge q a whle the felds are densty n a, pressure p a, temperature T a (n unts of k B ) and velocty v a. Each speces s assumed to satsfy the equaton of state p a = n a Ta. The partcle dynamcs s governed by the exstence of a gven external background magnetc feld B B(x) created by some external steady current and satsfyng B =. The dynamcal electromagnetc feld satsfes the non-relatvstc (electrc lmt) of the Maxwell equatons (see App. A for detals). In ths lmt there s no magnetc back-reacton on the partcle dynamcs. The magnetc feld appearng n the Lorentz force of Eq. (1) s therefore gven entrely by the external background feld B. Although the magnetc part of the dynamcal electromagnetc feld exsts n the non-relatvstc lmt (see Eq. (2)), the feld equaton s decoupled from the rest of the dynamcal system. The only couplng between the partcle dynamcs Eq.(1) and the dynamcal electromagnetc feld s therefore through the Posson equaton for the electrc potental, 2 φ = ρ ε, (2) where ε s the electrc (vacuum) permttvty. The electrc potental appearng n Eq. (2) s the leadng order n the non-relatvstc expanson. Wth the addton of ntal and boundary condtons Eq. (1) and Eq. (2) provdes a coupled set of equatons for the dynamcal felds n a,v a and φ whch are drven by the external magnetc feld B. For clarty of the presentaton, we wll n ths work present the above system n the case of two speces of charged partcles, electrons and ons wth atomc number Z = 1, labeled by a = {e,}, respectvely. The system of equatons we consder s however straght-forwardly generalzable to the case of multple speces. The charge densty and current densty are gven by ρ = q n +q e n e, J = n q v +n e q e v e. (3) Charge-current conservaton s therefore drectly satsfed by the sum of the contnuty equatons. We wll consder the system (1)-(3) n a perturbatve setup where the felds are assumed to approxmate the soluton n the form of perturbatve expansons. It s therefore convenent to frst put the system n dmensonless form. Introducng the characterstc length l, frequency ω, and velocty v = ωl together wth the scale of the magnetc feld B, densty n, and temperatures T a wth T e = m v 2, the system n dmensonless form s ελ 2 φ = n e n, t n a + (n a v a ) =, εµ a d t v a = ετ a log(n a )+z a ( φ+v a B), (4) where all quanttes are now dmensonless. The parameter z a specfes the sgn and multple of the elementary charge for the speces and the dmensonless parameters are µ a = ma m, τ a = T a T e, ε = ω ω c, and λ = ω2 c, = ε ωp, 2 B 2 m wth the on cyclotron frequency ω n c, = e B m and on plasma frequency ωp, 2 = ne2 ε m. The physcal parameters are the mass rato µ a and the temperature rato τ a together wth the ntal and boundary condtons, whle the remanng parameters ε and λ smply set the unts. In partcular, f ω s chosen to be ω c,, the characterstc length s l = ρ s m Te e B that the characterstc velocty v s the cold on acoustc speed and λ = λ2 D ρ 2 wth the Debye length s λ 2 D = ε T e. The dynamcal system (4) descrbes a non-relatvstc flud of charged partcles. It ne 2 therefore has a thermodynamcal descrpton that naturally gves rse to a number of conserved quanttes. In App. B, we provde the assocated conserved quanttes and show how the total energy can be obtaned from the flud energy conservaton equaton. such 3

5 2.1 Perturbatve setup We wll n the followng assume non-zero postve flud denstes and effectvely consder two spatal dmensons by assumng that the drecton of the external magnetc feld b = B B s a preferred constant drecton. Here, we have defned B B 2. We wll for smplcty neglect the dynamcs n the drecton of the external magnetc feld. The parallel component v a, = b v a of the velocty s therefore assumed to vansh. In the case where the dynamcs of the denstes n a s much slower than the dynamcs of the velocty felds v a, the system (4) exhbts two tme scales. It s therefore natural to splt the system nto two subsystems, one for the fast dynamcs and another for the slow dynamcs. Ths s possble n the framework of sngular perturbaton methods. For many practcal purposes, we do not requre to resolve the fast partcle dynamcs,.e. the physcs on the scale of the on gyro-frequency. We thus consder a long-wavelength expanson where all dervatves come wth a factor of ε. In addton, we are manly nterested n the systems where the denstes are suffcently large for the system to be n a near neutral state. Note that n the lmt λ =, the system of equatons (4) are nvarant under scalng transformatons of the densty. We thus consder the two-flud system (4) n a perturbaton n ε whle keepng λ, µ a and τ a fnte. 1 The perturbatve equatons are equvalently obtaned by choosng ω = ω c, and scalng the ndependent varables by t ε 1 t and x ε 1 x, together wth the scalng of the electrc potental φ ε 1 φ such that the small scale structures of the electrcal potental are suppressed. Other choces of perturbaton setups are dscussed n App. C. The partcular placement of the expanson parameter ε n Eq. (4) puts both Posson s equaton and the momentum equaton nto the form of a sngular perturbaton problem, where the zerothorder equaton has a lower dfferental order than the unperturbed equaton. Ths means that the degenerate system can not satsfy all ntal and boundary condtons n general. In partcular, the velocty felds v a are no longer dynamcal and ther ntal condtons are specfed by the ntal condtons for the denstes n a and electrc potental φ. We nevertheless expect that a soluton to the full system (4) wll approach the soluton of the perturbatve system much n the same way ordnary dynamcal systems behave [1] (at least for some ntal condtons and for some tme nterval). Under ths assumpton, we can therefore restrct to consderng the sngular perturbaton problem and consder the dynamcal felds n a formal power seres, n a = k= ε k n (k) a, φ = ε k φ (k), v a = k= k= ε k v (k) a. (5) Snce we can not hope to prove convergence, we wll smply assume that the expanson behaves as a dervatve expanson and consder the frst two terms. The best we could do would be to estmate the radus of convergence by extrapolaton usng a fnte number of terms. In the end of the computaton, we of course set ε to unty and keep n mnd that the expressons only hold as a long-wavelength, low frequency expanson.e. for suffcently slowly varyng confguratons. We thus consder the system (4) order-by-order n ε. Due to the sngular perturbatve nature of the momentum equaton, we can solve for the components of the veloctes perpendcular to the drecton of the external magnetc feld. The frst couple of orders of the perpendcular components of the flud veloctes are 2 v () a, = v() E, v(1) a, = v(1) E +v() Da +v() Pa, v(2) a, = v(2) E +v(1) Da +v(1) Pa. (6) wth the perpendcular drft veloctes defned by v (k) E b φ(k) B v (1) Da, v () Da τ b logn () a z a B, v () Pa µ a b d () t v () E z a B, (7) ). τ ( a z a n () B b n (1) a n (1) a logn()), v (1) Pa µ ( a z a B b d () t v (1) a, +v(1) a, v() E 1 For example, for B = 1T and n = 1 19 m 3, one has λ µ e. 2 In App. C, we dscuss perturbatve setups where the damagnetc-drft velocty appear at zeroth order. 4

6 One mght expect that the electron-polarzaton drft s neglgble compared to the on-polarzaton drft as t s suppressed wth a factor of µ e. However, for completeness we wll keep t throughout. We wll drop the perpendcular subscrpt n the followng. At zeroth order, the Posson equaton s smply ρ () = and the system s therefore neutral. In fact, ths mples that any charge separaton n the system that mght occur, e.g. polarzaton or damagnetc velocty drfts must enter at a hgher order. Consequental, the veloctes must be equal v () = v e () and therefore J () =. We can now defne a sngle flud densty n () = n () = n () e (at ths order) and snce the parallel component of the velocty felds s neglected wrte the contnuty equaton, t n () + Γ () =. (8) where we have ntroduced the densty flux Γ () = n () v () E. The system of zeroth order equatons s closed by the frst-order Posson equaton λ 2 φ () = ρ (1) and the correspondng charge-current densty conservaton equaton, t ρ (1) + J (1) =. (9) The frst-order correcton to the current s J (1) = n () (v () D + v() P ) + ρ(1) v () E damagnetc and polarzaton-drft veloctes v () D = v() D v () De and v() P = v() P v () We note n passng that n the case of multple speces, the on densty n () by the total densty of (postve) ons such that factors of v (1) wth the combned Pe, respectvely. can smply be replaced E does not appear n the current. Due to the sngular drft expanson of the perpendcular veloctes (6), tme dervatves are ntroduced n the current J (1) tself by the polarzaton drft. In fact, ths s the reason that the seemngly strct neutral approxmaton λ =, n whch Posson s equaton n (k) = n (k) e, can make sense even though t s equvalent to mposng J (k) =. Ths s exactly what s meant by the quas-neutral approxmaton. To evolve the charge densty ρ (1) both the electrc potental φ () and ts tme dervatve t φ () have to be evaluated at equal tmes. Usng the Posson equaton and the charge current conservaton equatons (9) one obtans the ellptc equaton D [ t φ ()] + J (1) =. (1) where we have ntroduced the modfed current J (1) = n () (v () D + v() S ) + ρ(1) v () E wth v() ( ) ( ) S = (µ +µ e ) 1 B v () E φ () B and the dfferental operator D (µ +µ e ) n() +λ B 2 usng that φ = φ wth = b ( b). The boundary value problem gven by Eq. (1) requres an ndependent set of boundary condtons for the tme-dervatve of the φ () to be specfed. The two equatons (8) and (1) provde a consstent system for the felds n () φ (), t φ () through the correspondng source J (1). The dynamcs at the lowest order s thus drven by the charge separatng terms v () Da and v() Pa even though they are n prncple frst-order quanttes. At ths pont t s worth commentng on a couple of detals of how the perturbatve equatons comes out. Frst, we note that t s suffcent to solve the ellptc equaton (1) and smply treat φ () as a dynamcal feld (n place of solvng the Posson equaton). Ths comes at the cost of exchangng the boundary condtons for φ () wth ntal condtons. Proceedng to the next order, where the second-order tme dervatve of the electrc potental s requred, we note that even t φ () can be treated as a dynamcal feld. In ths case the boundary condtons for t φ () are traded wth an addtonal set of (dependent) ntal condtons. In ths way, each order only ntroduces one addtonal ellptc equaton. Secondly, the set of equatons at a gven order s only closed by usng nformaton about the velocty feld from the succeedng order. Ths s generc for all orders of ths partcular perturbatve setup and makes the damagnetc and polarzaton-drft veloctes as well as the charge and current denstes of the next order avalable n terms of the lower order quanttes. In ths sense, Eq. (9) s redundant for obtanng the zeroth-order felds. Fnally, there exsts an alternatve (but equvalent) set of equatons that takes the place of Eq. (1) at the cost of 5

7 ntroducng an auxlary feld (see e.g. [11] for the zeroth order). The presented set of equatons (8) (1) have the advantages that they have a straght-forward generalzaton to systems wth multple partcpatng speces. The frst-order correctons to the denstes can be computed from the contnuty equatons t n (1) a + Γ (1) a =, (11) where the correcton to the flux densty s Γ (1) a = n () v a (1) +n (1) a v () E. Note that J(1) = Γ (1) Γ (1) e. The system of frst-order equatons s closed by the Posson equaton λ 2 φ (1) = ρ (2) and the second-order charge-current densty conservaton t ρ (2) + J (2) =. (12) Thesecond-ordercorrectontothecurrentsJ (2) = n () (v (1) D +v(1) P )+n(1) v (1) n (1) e v e (1) +ρ (2) v () E wth the combned damagnetc and polarzaton-drft veloctes v (1) D = v(1) D v(1) De and v(1) P = v(1) P v (1) Pe, respectvely. As a consequence of the sngular scheme, the second-order correcton to the current J (2) now nvolves both the tme dervatve of the correcton to the electrc potental t φ (1) and the second-order dervatve t 2φ(). An equaton for the latter s obtaned drectly from Eq. (1) whle an equaton for the frst s obtaned n a smlar way to Eq. (1),.e. we take the combnaton of the frst-order Posson and charge-current equaton (12), [ D t φ (1)] + J (2) =, D [ 2 tφ ()] + ( t J(1) ) = (13) wth the modfed current J (2) = n () (v (1) D +v(1) S ) +n(1) v (1) n (1) e v e (1) + ρ (2) v () E wth v(1) S = b B ( t (µ (v () D + v () P ) + µ e(v () De + v() Pe )) + v() E (µ v (1) + µ e v e (1) ) + (µ v (1) + µ e v e (1) ) v () E ). The second equaton s obtaned from Eq. (1), t s therefore an equaton n terms of zeroth-order felds t J(1) = n () t (v () D +v() S )+ρ(1) t v () E J(1) v () E Γ() (v () D +v() P ). Wth Eq.(13), t s possble to treat φ (1) and t φ () as dynamcal felds at the cost of exchangng ther boundary condtons wth consstent ntal condtons. The complete perturbatve system of equatons to frst order s thus gven by Eq. (8)-(9), Eq. (11) and Eq. (12) together wth the two ellptc problems gven by Eq. (13) and the dynamc evoluton of the felds t φ (), φ (), and φ (1). Notce that smlarly to the precedng order, the second-order correctons to the charge and current denstes are avalable through Eq. (12) n terms of only frst-order quanttes. Conserved quanttes A partcular placement of the perturbaton parameter ε can potentally break the conservaton of a conserved charge of the full system, 3 such that t s conserved only f all orders are ncluded. We therefore check f the perturbaton equatons retans the conservaton of the constants of moton at each order. The conserved charges of the full system are gven n App. B. If contrbutons from the boundary vansh (.e. vanshng fluxes at the boundary), the conservatve form of Eq. (8)-(9) and Eq. (11)-(12) ensures that mass and charge conservaton s preserved at each order. The ntegrated felds n (), n (1) a, ρ (1), ρ (2) are therefore constants of moton. In ths case, the total energy E s therefore also preserved at each order. Assumng that the felds vansh at the boundary, the frst couple of terms n the expanson of the total energy densty are gven by ( E () = λ ( φ ()) 2 + n () τ a logn () + µ ( a φ () ) 2) 2 2 B 2, (14) a ( ) τ a n (1) a logn () +1 + E (1) = λ φ () φ (1) + a ( ( µ a 1 B 2 φ() 2 n(1) a φ() +n () φ (1) + τ a logn () + µ )) a d () t v () E z a z a 3 For example by splttng the convecton term n a conservaton equaton.. (15) 6

8 These can be obtaned ether drectly by nsertng the expanson nto Eq. (28) or by usng the equatons of moton. The lowest contrbuton E () gven by Eq. (14) s obtaned by multplyng Eq. (1) by φ (), ntegratng over space and usng ntegraton by parts. Ths wll leave two resdual terms, one can be exchanged by the consderng the volume ntegral over τ a t (n () logn () ) and usng the contnuty equaton (8), the other by usng the Posson equaton. In a smlar fashon, the frst-order correcton E (1) gven by Eq. (15) can be found by consderng Eq. (13) and usng the remanng equatons of moton. Dscusson Some general comments on the presented perturbatve equatons can now be made. In many cases the gradent of the densty s substantal and t would seemngly be preferable to treat the damagnetc-drft velocty at the same order as the E B-drft velocty. In ths regard, the velocty expanson gven by Eq. (6) s seemngly an unconventonal drft-orderng (see e.g. [4]). However, ths requres a perturbatve setup where charge separaton can be accommodated at zeroth order. Ths means that the separaton of charge s assumed to be suffcently large as to be treated as a zeroth-order quantty. But n such cases, a quas-neutral lmt λ = can not make sense. It s therefore partcularly nterestng to consder the dfferent lmts of the parameter λ whch sets the scale of the densty,.e. λ B 2 n. The strct quas-neutralty lmt λ s reached for suffcently large denstes. In ths lmt, the system of equatons (4) becomes nvarant under scale transformatons of the densty. The quas-neutral system therefore predcts equvalent dynamcs for all magntudes of denstes. On the other hand, for small denstes λ O(1), the small scale structure of the electrc potental becomes ncreasngly mportant and charge separaton can consstently be ncluded at the lowest order. Ths always comes at the cost of removng the possblty of a quas-neutral lmt. In App. C, we dscuss two such dlute perturbatve setups. Even though the charge separatng terms enters at frst order n the perturbatve setup presented n Sec. 2.1, the damagnetc and polarzaton-drft correctons for the next order are stll avalable at any gven order. It s therefore possble to obtan a better approxmaton for the velocty feld by v a k ṽ a (k) +O(v (k+1) E ) wth ṽ (k) a = v (k) E +v(k) Da +v(k) Pa. (16) Naturally, the partcle speces appear n the system (4) on equal footng; ther dynamcs dstngushed by mass, charge and n the sothermal system by ther constant temperature. As the sgn of the carred charge and the mass rato µ e are fxed, the asymmetry between the dynamcs of the speces s solely controlled externally by the on-temperature rato τ. However, at zeroth order gven by Eq. (8)-(1) the dynamcs s fully determned by the combned parameters µ + µ e and τ +τ e multplyng the combned-drft veloctes v () P and v() D, respectvely. At ths order, the system therefore does not dstngush between the speces. The asymmetry between the speces appears at frst order through the mass rato µ e as the frst-order combned-drft polarzaton velocty v (1) P now ntroduces cross-terms wth the polarzaton-drft velocty, µ 2 µ 2 e µ 2 and the damagnetc-drft velocty µ τ µ e τ e µ τ. One therefore expects that non-zero on temperature plays a sgnfcant role for the frst-order correctons. 3 Smulaton We perform a quanttatve analyss of the perturbatve equatons gven by Eqs. (8)-(13) usng numercal methods. In partcular, we are nterested n what parameter regmes the perturbaton seres converges to a degree where we can trust that t represents the dynamcs of the full system gven by Eq. (4). Ths s possble, snce we have drect access to the frst-order correctons of the dynamcal felds n a numercal smulaton. We can thus montor when the correctons become too mportant to say that the lowest order represent a converged seres on ts own. In order to control the dfferent dynamcal features of the two-flud we consder the smulated dynamcs of a seeded blob 7

9 usng a numercal mplementaton of the perturbaton equatons. Here, we are manly nterested n the valdty of the perturbatve equatons and are not partcular concerned wth the blob dynamcs tself. We wll therefore restrct to consderng the felds themselves and not the center of mass quanttes. A great deal of work can be found n the lterature on the dynamcs of seeded blobs (see e.g. [12] and references wthn). The numercal mplementaton and convergence of the numercal methods have been verfed usng the method of manufactured solutons. For actual smulatons, we have used the conserved charges,.e. mass, charge and energy of the system to montor f the numercal soluton stays wthn physcal constrants. Let us now defne the operator C p (n) (f) f(n) p wth n,p N f (n 1) + that gves the L p -norm rato p of the n th order feld coeffcent to the (n 1) th order feld coeffcent. We wll use ths operator to estmate the mportance of the hgher-order correctons. Furthermore, we ntroduce the vortcty ω a v a wth the parallel component ω a b ω a. 3.1 Intal condtons We choose a Cartesan coordnate system {t,x,y,z} wth the drecton of the magnetc feld set to b = e z. The magntude of the external magnetc feld s set to B 1 (x) = 1 + κ(x x ), where the curvature s parametrzed by κ satsfyng the postvty condton κx < 1. The dynamcal felds are put on a dscrete two-dmensonal grd spanned by the x and y coordnates. We run smulatons wth a seeded blob on a doman D wth perodc boundary condtons n the y-drecton and Drchlet condtons n the x-drecton. The blob s ntalzed n the zeroth-order densty n () by a two-dmensonal Gaussan functon n () = 1+Aexp ( (x x ) 2 ) 2σ 2, n (k) a =, (17) wth ampltude A, equal wdth σ n both drectons and ts ntal center of mass at x = x D. We note that the system s by constructon neutral at zeroth order and the ntal condtons are compatble wth the perturbatve assumpton p (n a ) 1. The boundary condtons for the denstes are n () D 1 and n (k) a D =. If the total area of the smulaton doman s denoted by V, the total densty n the system s dvn n(v +2πAσ 2 ). The electrc potental φ (to all orders) must be zero ntally because of the ntal condtons for the denstes. The perpendcular components of the veloctes are gven completely n terms of the densty and electrc potental due to the (sngular) drft expanson (6). Ther ntal and boundary condtons are therefore also determned by the ntal and boundary condtons of the dynamcal felds. The choce of ntal condtons must therefore be chosen wth care. Whle the above specfed condtons nnocuously sets the blob at rest(to lowest order) ntally, the zeroth-order velocty can not possbly resemble the true velocty feld at ntal tmes t (as the correcton completely domnates the zeroth order,.e. p (v a ) 1). However, snce the system s nvarant under tme-reversal, a smulaton can be contnued from any pont n tme. Thus, f the system farly quckly orders tself, a zerothorder smulaton ntalzed from rest can be thought of startng from a non-zero ntal velocty at a later tme t 1 where the frst-order correctons are suffcently suppressed and the ntal perod can be dscarded. The queston then becomes whether the evolved ntal condton s of nterest. Thus f we consder a smulaton on some fnte nterval t t T, we expect rapdly changng velocty felds n a neghborhood of t,.e. t t t 1 and a transton to slowly changng felds (compared to the densty and electrc potental) n the remanng nterval t 1 t T. 3.2 Results The dynamcs of a gven soluton s determned by the external parameter κ, the physcal parameter τ representng the rato between on and electron temperature, the shape of the blob gven by the parameters A and σ and fnally the scale λ B 2 n. The parameter λ sets the scale of the densty. Recall, that when the system s not strctly quas-neutral (λ ), the system equatons are not 8

10 nvarant under scalng transformatons of the denstes. We can therefore study the dynamcs at dfferent scales of the densty. We consder the followng regmes of parameter space () cold ons τ = n the quas-neutral lmt λ =, () hot ons τ > n the quas-neutral lmt λ =, and () the system at fnte denstes λ > (non-zero charge densty). The frst case vares the ntal condtons and shows the sgnfcance of the spatal scales, whle the latter two shows the sgnfcance of temporal scales by varyng a physcal parameter. In all cases we set the curvature parameter to κ =.15. Cold ons n the quas-neutral lmt We start by consderng the dynamcal case consstng of cold ons τ = n the quas-neutral lmt λ = (nfnte densty lmt). In ths regme, the system of equatons s densty scale nvarant and n (k) e = n (k). The dynamcal felds for dfferent parameter sets of {A,σ} can therefore to some extent be related by scalng transformatons even though the system s non-lnear. The transformatons can be estmated by usng the ntal condtons (17) and the equatons of moton. They should therefore at least hold at ntal tmes. We omt the detals on how these can be obtaned and smply lst the estmated scalng transformatons for the operator C appled on the felds n Tab. 1. Feld Tme Value 2 (φ) σ 1 2à 1 σ 1 2à 1 2 (n) σ 1 2à 1 σ 1 2A (v e), 2 (ṽ e) σ 1 2 à 1 2 σ 3 2 à (v ), 2 (ṽ ) σ 1 2 σ 1 2 Table 1: Scalng transformatons of the operator C appled on the felds n the quas-neutral (λ = ) cold on (τ = ) lmt wth respect to {A,σ}. Here à A 1+A. The scalng relatons hold to good approxmaton wth respect to σ whle the relatons for the ampltude A only holds for early tmes as shown n Fg. 2. Note that for the densty, the value scales wth A and not Ã. Wenowconsdernumercalsmulatonsforarangeofntalcondtons{ 1 4 A 4,1 σ 7} whch sets the ntal length scale of the densty and thus determnes the ntal strength of the electron damagnetc-drft velocty. In Fg. 1, we show a typcal evoluton of the densty feld n n () +n (1) and the separate expanson coeffcents. Fgure 1: The densty feld n the quas-neutral lmt λ = and cold on lmt τ = for the parameters A = 1 and σ = 2. From left to rght, the feld n (), n (1), and the frst-order approxmaton of the densty n n () +n (1) are shown (overlapped) at three dfferent tmes {,245,377}ω c,. The densty s advected towards the rght whle ncreasng (decreasng) n the upper (lower) 9

11 half plane,.e. the blob moves slghtly upwards. Ths s generally the case n the cold on lmt as the cross terms n J (2) between the polarzaton and damagnetc-drft veloctes are domnated by the electrons. For suffcent non-zero on temperature the roles are reversed and the blob moves downwards. 2 1 () σ (v e,y) σ (v e,x) ω c, tσ 1 2 () Ã (v e,y) 1 2 Ã (v e,x) ω c, tã () σ (v,x) σ (v,y) ω c, tσ 1 2 (v) 2 (v,x) 2 (v,y) ω c, t 1 3 σ = 1 σ = 2 σ = 3 σ = 4 σ = 5 σ = 6 σ = 7 A = 1 4 A = 1 2 A = 3 4 A = 1 A = 2 A = 3 A = 4 Fgure 2: The scalng transformatons lsted n Tab. 1 for the operator C appled on the component of the electron (left) and on velocty (rght) for (, ) fxed {λ =,τ =,A = 1} wth varyng wdth σ, and (, v) fxed {λ =,τ =,σ = 4} wth varyng ampltude A. Here, Ã A 1+A. In Fg. 2, the velocty felds for dfferent smulatons are shown after applyng the scalng transformatons lsted n Tab. 1. For a fxed ampltude A, the scaled felds agree wthn good approxmaton for all values of σ on the entre tme nterval (see Fg. 2(,)). For varyng ampltude A, the two velocty felds evolves on dfferent tme scales (see Fg. 2(,v)). So whle the scalng relatons provde almost perfect agreement for the electron velocty (as the damagnetc-drft velocty domnates on the entre tme nterval), the on velocty ntally evolves on the same tmescale for all values of the ampltude, but wll after the ntal acceleraton grow at a faster rate for larger ampltudes. The densty and electrc potental are therefore also affected at later tmes albet to a smaller degree (as the relatve magntude between the correctons for the two speces are several orders). 1

12 .5.25 () 2 (v e,x) 2 (v e,y) 2 (v,x) 2 (v,y) C (2) 2 (v e,x) C (2) 2 (v e,y) C (2) 2 (v,x) C (2) 2 (v,y) () ω c, t ω c, t 1 3 Fgure 3: The operator C appled on the electron and on velocty components. Both fgures show the absolute values for the partcular smulaton parameters {λ =,τ =,A = 1,σ = 5}. The second-order correcton to the velocty feld (6) s approxmated by neglectng v (2) E. The operator C appled on the components of the velocty felds for both the electron and on speces s shown n Fg. 3 for the partcular set {A = 1,σ = 5}. The frst-order correctons to the velocty felds decreases farly rapdly over an ntal perod after whch t becomes slowly varyng. Snce there are densty gradents at all tmes, the frst-order correcton to the electron velocty (whch s domnated by the damagnetc-drft velocty) s orders of magntude larger than the correcton to the on velocty. For suffcently small scales (.e. small σ or large ampltude A), the electron-velocty feld s therefore not very well represented by the lowest order as s especally evdent from the zeroth-order vortcty feld ω (see Fg. 5). However, ths s not explctly seen from the frst-order correcton to the densty (Fg. 4()) as the convecton s to a large extent drven by the E B-drft velocty. Snce the E B-convecton n turn s drven by the sub-leadng charge separatng terms through the electrc potental, t becomes necessary to consder the second-order correctons to the velocty n order to check the convergence. Fg. 3() show the C operator appled on the second-order correctons for the veloctes. It shows that the second-order correctons are almost of the same order as the frst-order correcton to the densty. Thus, ncludng the frst-order correcton n the electron velocty we see the same order of approxmaton (n) 2 (φ) () (ṽ e,x) 2 (ṽ e,y) 2 (ṽ,x) 2 (ṽ,y) () ω c, t ω c, t 1 3 Fgure 4: The operator C appled on () the densty and electrc potental, and () the corrected electron and on velocty gven by Eq. (16). Both fgures show the absolute values for the partcular smulaton parameters {λ =,τ =,A = 1,σ = 5}. Note that n the quas-neutral lmt we have n = n e = n. 11

13 Strctly speakng, t s necessary to check one order above the truncaton-order to estmate the convergence of the gven order of approxmaton. However, n the gven perturbaton setup we have access to hgher order correctons to the velocty feld at a gven order. It therefore makes sense to nclude the terms n the velocty feld at the gven order (post smulaton). Ths s exactly the corrected velocty feld ṽ a gven by Eq. (16). In Fg. 4, the corrected velocty felds, densty and electrc potental are shown. The frst-order correctons to ṽ a s not surprsngly magntudes smaller than the true frst-order velocty correcton v a (1). In some sense, we can therefore use the corrected velocty to gan an estmate of how well the overall approxmaton s at a gven order. However, t should be stressed that whle the corrected velocty feld gves a better approxmaton at a gven order, t does not consttute the flux Γ at that same order. InFg. 4, all frst-ordercorrectons areshowntogrow wthtme. Oneway tonterpretthssthat when the densty blob begns to spread out, small scales starts to develop as s partcular promnent when the blob tal begns to spral. Comparng the densty feld (Fg. 1) wth the magntude of the frst-order correcton to the vortcty (Fg. 5) by scalng tme approprately, we see that the onset of ths local spnnng s pcked up by the vortcty feld. In ths case, the vortcty therefore provdes a measure for the emergence of small scale structures. We also note that n the cold on lmt the relatve vortcty can wthn good approxmaton be estmated n terms of lowest-order quanttes by θ a ω() Da 2 ω () E 2. Ths quantty s sometmes also denoted the fnte Larmor radus (FLR) strength parameter n the case of ons as t s used as a measure for the strength of the charge separaton (ω e) 2 (ω ) θ e ω c, t 1 3 Fgure 5: The operator C appled on the electron and on vortcty along the magnetc feld drecton n comparson wth the FLR parameter θ e defned n the text. The fgure show the absolute values for the partcular smulaton parameters {λ =,τ =,A = 1,σ = 5}. The numercal smulatons and estmated scalng transformatons (lsted n Tab. 1) support that all hgher-order correctons are suppressed at large spatal scales,.e. wth ncreasng σ and decreasng ampltude A. As soon as the dynamcs s consdered on suffcently small scales the correctons becomes too mportant to assert that the soluton to the degenerate system has converged suffcently to represent the true dynamcs. The perturbaton equatons are therefore ndeed vald as a long-wavelength expanson. In partcular, we note that the zeroth-order system gven by Eqs. (8)- (1) can provde a converged soluton n the cold on lmt (at least on some tme nterval) f the system s consdered for suffcently large scales and f the corrected velocty s used(post smulaton) to approxmate the velocty feld. Hot ons n the quas-neutral lmt We consder numercal smulatons n the case of non-zero on temperature τ n the quas-neutral lmt (λ = ). In partcular, we consder two parameter seres, () one wth varyng on temperature {A = 1,σ = 4, τ 4}, and () one wth varyng wdth {A = 1,τ = 1,1 σ 7}. Smlarly to the cold on case, we lst n Tab. 2 scalng transformatons based on the ntal condtons (17) 12

14 and the equatons of moton. The transformatons relates the dynamcal felds for dfferent σ and τ and holds for suffcently large on-temperature rato τ. Feld Tme Value 2 (φ) σ 1 2 (1+τ ) 1 2 σ τ 4 2 (n) σ 1 2 (1+τ ) 1 2 σ τ 4 2 (v e), 2 (ṽ e) σ 1 2(1+τ ) 1 σ (v ), 2 (ṽ ) σ 1 2 (1+τ ) 1 2 σ τ 4 Table 2: Scalng relatons for the operator C appled on the felds n the regme where the correcton to the on velocty s domnated by the damagnetc-drft velocty,.e. at suffcently large ontemperature rato τ. For non-zero on temperature the asymmetry between the speces becomes much more apparent n the dynamcs as the cross terms between the polarzaton and damagnetc-drft veloctes are much larger due to the on mass rato µ beng unty. We therefore observe a transton pont n the dynamcs wth ncreasng on temperature when the cross terms are n balance. At suffcently large on-temperature, ths balance s broken and we leave the cold on lmt. In ths case, the blob has a downwards movement. An example of ths s shown n Fg. 6. In ths regme, the dynamcs happens on a faster tme scale (compared to the cold on lmt, see e.g. Fg. 1). The faster movement can therefore result n the blob leavng behnd a depleton n the densty whch can only be flled out by ncludng hgher-order correctons. For the frst-order system ths explctly shows up when the correctons to the densty becomes as large as the zeroth order and a negatve densty s obtaned. In Fg. 6, ths s seen at the last tme stamp. Fgure 6: The densty feld n the quas-neutral lmt λ = and non-zero on temperature τ = 1 for the parameters A = 1 and σ = 2. From left to rght, the feld n (), n (1), and the frstorder approxmaton of the densty n n () +n (1) are shown (overlapped) at three dfferent tmes {,133,2261}ω c,. The color scale for the frst-order correcton has been cropped to show more of the structure. The faster dynamcs s a drect consequence of the larger charge separaton n the system at non-zero on temperature whch generates a larger electrc potental. At some pont n tme, the E B-drft velocty therefore grows at a much faster rate whch s reflected drectly n the densty movement. Whether the frst-order soluton gves such large correctons depends on the parameters as can be seen from Fg. 7. In partcular, there exst a crtcal parameter doman for whch the frst-order system nevtably wll loose ts predctve power as the frst-order correcton becomes comparable wth the zeroth order. In ths parameter doman hgher-order correctons are necessary for the seres expanson to gve a meanngful approxmaton. We note that ths statement refers to 13

15 the unform convergence of the soluton. Lookng at the overall densty feld, the correctons are very localzed n space and the correctons can therefore be suppressed by ncreasng the blob wdth. C (1) (n) () τ = τ = 1 4 τ = 1 2 τ = 3 4 τ = 1 τ = 2 τ = 3 τ = 4 () σ = 1 σ = 2 σ = 3 σ = 4 σ = 5 σ = 6 σ = ω c, t ω c,t 1 3 Fgure 7: The coeffcent (n) n the quas-neutral lmt λ = wth A = 1 for () fxed wdth σ = 4 at varous on-temperature ratos τ, and () for fxed τ = 1 at dfferent values of σ. For some doman of the on-temperature rato and wdth the frst-order correctons wll be larger than the zeroth-order densty after some pont n tme. In Fg. 8, the frst-order coeffcent 2 s shown for the densty and corrected velocty felds for the partcular parameters {σ = 4,τ = 1 2 }. It s clear that the correcton to the velocty feld grows wth the on temperature(comparng wth Fg. 4). Nevertheless, for low on temperatures and blob wdth, the zeroth-order approxmaton can provde a farly converged soluton. For moderate on temperatures and blob wdth (see Fg. 8()), the frst-order correctons becomes necessary as the zeroth-order approxmaton can not capture the dynamcs sngle-handed. Fnally, for hgh on temperatures, the correctons to the velocty feld becomes comparable wth the zeroth-order approxmaton and hgher-order correctons are requred. Also, t s nterestng to note that due to the large E B-drft velocty, the corrected velocty feld ṽ a almost algn for both speces when the on temperature rato s near unty even though the damagnetc and polarzaton-drft veloctes come wth opposte contrbutons for the two speces (n) 2 (φ) () (ṽ e,x) 2 (ṽ e,y) 2 (ṽ,x) 2 (ṽ,y) () ω c, t ω c, t 1 3 Fgure 8: The operator C appled on () the densty and electrc potental, and () the corrected electron and on velocty gven by Eq. (16). Both fgures show the absolute values for the partcular smulaton parameters {λ =,τ = 1 2,A = 1,σ = 4}. Note that n the quas-neutral lmt we have n = n e = n. 14

16 Naturally, the correcton to both vortcty felds also ncreases wth the on temperature wth the electrons stll havng much faster local spnnng. In partcular, θ e decreases wth ncreasng on temperature and provdes a poor estmate of the relatve electron vortcty. On the other hand, the FLR parameter for the ons θ agrees to a much larger degree wth the relatve on vortcty at moderate temperatures The numercal smulatons for varous on temperature together wth the scalng transformatons gven n Tab. 2 show that the perturbaton equatons to frst order are vald for low frequences and second or hgher-order correctons are necessary to resolve the faster dynamcs. In addton, the dependence on σ reveals (as for the cold on lmt) that the perturbaton equatons are vald as a long-wavelength expanson. Fnte densty lmt We consder numercal smulatons for fnte denstes by consderng non-zero postve λ for a fxed set of ntal condtons {A = 1,σ = 2}. We consder two parameter seres, () one for vanshng on temperature τ =, and () one for unty on-electron temperature rato τ = 1. For a non-zero λ, the Posson equaton enters the system of equatons and the scale nvarance of the system s broken. The parameter λ s thus a measure for the sgnfcance of the Posson equaton. For fnte λ, the perturbaton equatons stll requres the system to be neutral at the lowest order (by constructon),.e. ρ () =. The charge densty therefore evolves due to the splt n frst-order densty correctons. Wth ncreasng λ, the electrc potental s determned progressvely more by the charge densty nstead of the charge separatng damagnetc and polarzaton-drft veloctes. As a consequence, the dynamcs happens on a much slower tme scale wth ncreasng λ. C (2) 2 (ρ) λ = 1 λ = 5 λ = 1 λ = 5 λ = 1 () () ω c, t ω c, t 1 3 Fgure 9: The C 2 operator appled on the charge densty ρ for () fxed {τ =,A = 1,σ = 2} wth varyng lambda λ, and () fxed {τ = 1,A = 1,σ = 2} wth varyng lambda λ. The legend apples for both cases. In Fg. 9(), the relatve sze between the second-order correcton and the frst-order correcton to the charge densty s shown n the cold on lmt τ = for dfferent values of λ. The sgnfcance of the second-order correcton ncreases as the densty becomes more dlute. Already for λ = 1 correspondng to n 1 15 m 3 at B = 1T, we observe qute sgnfcant correctons. For suffcently dlute systems, the frst-order corrected charge densty (.e. ρ = ρ () + ρ (1) ) therefore fals to converge. For non-zero on temperature, there s a balance between the contrbuton from the charge separatng veloctes and the charge densty. Whle the dynamcs s faster compared to the cold on lmt, as shown for the case of τ = 1 n Fg. 9(), the charge densty helps wth smoothng out the flud denstes such that small scales do not develop to the same extent. At fnte denstes, the correctons to the felds are therefore suppressed to a much larger degree than n the correspondng quas-neutral lmt λ = where negatve denstes are encountered. However, the second-order correcton to the charge densty s stll substantal. In all cases, t therefore seems 15

17 more approprate to use the perturbaton equatons gven n App. C n the dlute regme. Here, the charge separaton s approprately assumed to be large and can appear at the lowest-order. 4 Dscusson We have consdered a charged non-relatvstc flud system governed by the Lorentz force where n partcular, the electrc feld were treated self-consstently through the back-reacton on the partcle dynamcs. The system was consdered n a perturbatve setup sutable for studyng the long wavelength and low frequency dynamcs of the system. Expandng the dynamcal felds n a formal seres expanson a set of perturbatve equatons for the frst-order corrected felds were obtaned. We showed that the perturbaton equatons preserved the conserved charges of the full underlyng system at each order of the expanson. Furthermore, we showed that a non-trval quas-neutral lmt was possble only because t does not allow charge separatng terms at the lowest order. In partcular, ths excluded the damagnetc-drft velocty to appear at the lowest order n the expanson. Perturbatve schemes that allows the damagnetc-drft velocty at the lowest order were gven n App. C. These schemes do not allow a non-trval quas-neutral lmt as the polarzaton-drft velocty does not appear at the same order and thus does not ntroduce tme dependence n the electrcal current. Fnally, the perturbatve setup consdered n the man text also has a natural generalzaton to multple speces. Subsequently, usng numercal smulatons gvng explct access to the frst-order correctons to the dynamcal feld varables, we have provded a quanttatve analyss showng that the perturbaton equatons are ndeed vald at long wavelengths and low frequences. In partcular, we have consdered three dfferent parameter regmes from whch we summares: () For zero on temperature, the lowest-order approxmaton can provde a farly converged soluton to the dynamcs at suffcent large scales (f the corrected velocty feld s used). () For non-zero on temperature, the dynamcs generally happens on faster tme scales. When the on temperature s non-zero, but suffcently low the frst-order corrected felds provde a relatvely converged approxmaton, whle above a certan temperature hgher-order correctons are requred n order to resolve the faster dynamcs. () For fnte denstes, the Posson equaton gves rse to a fnte charge densty (.e. a splt between the densty for the speces). In ths case the dynamcs happens on a slower tme scale. However, the small scale structures n the electrc potental becomes ncreasngly more mportant. In ths regme, the dlute perturbaton equatons gven n App. C becomes applcable. Furthermore, we have provded some approxmatve scalng transformatons that relates the soluton between dfferent set of parameters n the cases where the system s densty scale nvarant (.e. n the quas-neutral lmt). We have verfed that the scalng transformatons are n farly good agreement wth the numercal smulatons. In concluson, we have addressed all the ponts ()-(v) lsted n the ntroducton. Ths work establshes the very frst step n developng perturbatve drft-flud equatons where all dynamcal felds are under full perturbatve control. There are therefore several natural extensons: Frstly, t would be nterestng to lft the assumpton of a vanshng parallel velocty component v a, generalzng the equatons to a full (3+1)-dmensonal system. In partcular, ths poses two challenges. Snce the parallel component can not be treated n the same sngular perturbaton as the perpendcular components of the velocty, t must be treated dynamcally. Its evoluton can therefore happen on a dfferent tme scale. Also, the relatve magntudes between the components of the electrcal current could potentally be large, thus posng a challenge f the system s requred to be neutral at the zeroth order. Secondly, t would be desrable to lft the sothermal assumpton on the flud thereby makng the nternal energy dynamcal. Besdes addng addtonal degrees of freedom to the system, ths should n prncple a straght-forward extenson. Fnally, t would be nterestng to nclude dsspatve correctons, lke vscosty or other resstve terms and furthermore check how these modfy the conserved currents. Another drecton that hold nterest s to nclude a self-consstent magnetc back-reacton on the partcle dynamcs. Ths can be done by ether ncludng hgher-order correctons from the non-relatvstc expanson of the Maxwell and flud equatons, or by usng the magnetc lmt of 16