arxiv: v1 [physics.plasm-ph] 5 Aug 2012

Transcription

1 Classial mirosopi derivation of the relativisti hydrodynamis equations arxiv: v1 [physis.plasm-ph] 5 Aug 2012 P. A. Andreev Department of General Physis, Physis Faulty, Mosow State University, Mosow, Russian Federation. (Dated: June 29, 2018) We present mirosopi derivation of the relativisti hydrodynamis (RHD) equations diretly from mehanis omitting derivation of kineti equation. We derive ontinuity equation energy-momentum onservation law. We also derive equation of evolution of partiles urrent. In non-relativisti hydrodynamis equation of partiles urrent evolution oinide with the equation of momentum evolution, the Maxwell s equations ontain onentration eletri urrent (whih proportional to the partiles urrent), so, to get a lose set of equations we should have equations of evolution of the onentration the partiles urrent. Evolution of the partiles urrent depends on the eletrial magneti fields. Thus, we obtain the set of the RHD equations as the set of the ontinuity equation, the equation of partiles urrent the Maxwell equations. This approximation does not require to inlude the evolution of momentum allows to onsider physial problems. Certanly, partiles urrent evolution equation ontains some new funtions whih we an express via onentration partiles urrent or we an derive equation for this funtions,, thus, get to more general approximation. This approximation also developed disussed in this paper. I. INTRODUCTION Relativisti kinetis hydrodynamis have been in a enter of attention. A lot of papers books dediated to this topi have been published in last several years. Wide reent review of relativisti kinetis has been presented in R. Hakim s book [1], but we also should mention a paper about relativisti kinetis [2], where was suggested a method of equations derivation whih we suppose to use for derivation of hydrodynamis equations. When one onsiders the relativisti hydrodynamis equations he usually uses onservation laws: onservation of partiles number momentum energy, in Lorentz ovariant form, along with the Maxwell equations. To get a lose set of the hydrodynami equations we should present a onnetions between the momentum density P the veloity field v. In the relativisti hydrodynamis (RHD) we have two equations instead of the Euler equation in the non-relativisti hydrodynamis, where the Euler equation is both the momentum balane equation the equation of evolution of the partiles urrent j (whih emerges in the ontinuity equation proportional to the veloity field j = nv, where n is the onentration of partiles). In the RHD s these equations are different. We an use partiles urrent evolution instead of evolution of momentum density. So, we do not need to find onnetion between P v = j/n. In the Newton s mehanis it s relativisti generalization the law of momentum evolution is the basi dynamial law. It gives us mirosopi dynamial piture. When we are going to derive marosopi dynamial equations we should find one ontaining information keeping in the Newton s seond law, the partiles urrent evolution equation is one of them, as well as the momentum balane equation, these two equations math in the nonrelativisti ase. When one has relativisti momentum balane equation it is a hard job to find a onnetion between momentum density P(r,t) partiles urrent j(r,t) = n(r,t)v(r,t) to obtain the lose set of the hydrodynamis equations for the system of partiles with the temperature. Where is the onnetion between the momentum p i (t) the veloity v i (t) of the one partile p i (t) = m iv i (t), (1) 1 v2 i (t) 2 v i (t) = p i(t) p 0i (t), (2) where is the speed of light, is proportional to the energy of the partile (t) = E i(t) = p 2 i(t) 2 +m 2 i 4. (3) If we onsider system of partiles, it is no simple matter to suggest a onnetion between P = i mv i(t)/ 1 v 2 i(t)/ 2 v = i v i(t). When we onsider the relativisti dynamis the Newton s equation still an be used, but we must onsider the relativisti onnetion (1) between p i (t) v i (t), so d dt p i(t) = ee(r i (t),t)+ e v i(t) B(r i (t),t), (4) where E(r i (t),t) B(r i (t),t) are the eletri magneti fields ating on the i-th partile. Otherwise, for partile s aeleration we have d dt v i(t) = e m 1 v2 i(t) 2 ( E(r i (t),t)+ 1 [v i(t), B(r i (t),t)] 1 ) v i(t)(v 2 i (t)e(r i (t),t)). (5)

2 2 In this paper we pay attention to the lassi mehanis hydrodynamis, but where are also papers dediated to the quantum hydrodynamis [3], [4], [5], the relativisti quantum hydrodynamis [6], [7]. This paper is organized as follows. In Se. II we present a brief review of the relativisti hydrodynamis for partiles system with the small thermal veloity spread. In Se. III we onsider existing in the literature methods of obtaining of the RHD equations for the system of relativisti partiles with the temperature. In Se. IV we derive the RHD equations from the piture of mirosopi motion of partiles desribed by the Newton s law. In Se. V we onsider temporal evolution of new funtion appeared at derivation desribed in the Se. IV. In Se. VI we derive the equation of evolution for the energy-momentum density to show that our treatment give the well-known equation. In Se. VII we disuss ways to obtain a lose set of equations. In Se. VIII we onsider dispersion of waves in the relativisti plasma. In Se. IX we present the brief summary of our results. II. RELATIVISTIC HYDRODYNAMICS OF PARTICLES SYSTEM AT SMALL TEMPERATURE As a simple example of relativisti many partile system we onsider a relativisti eletron beam. If it s beam is monoenergeti or it has small thermal spread of partile veloities when all partiles have (near) equal veloities u, so momentum of the system P = N mu/ 1 u 2 / 2 = Nmu/ 1 u 2 / 2 is proportional to the momentum of single partile, it onnets with the partile veloity by well-known relativisti formula. Here below, we use u as veloity field for approximately monoenergeti system of partiles. Here below, for simpliity, we onsider one speies. Following by Ref.s [8], [9], [10], [11], we an write an example of using in literature set of the RHE t n+ (nv) = 0, (6) t p+(v )p = ee+ e v B, (7) where momentum of medium is presented in the form p = γmv, where γ = (1 v 2 / 2 ) 1/2, analogously to single partile, to lose the set of equations. In equations (6) (7)nis the onentration of partiles, v is the veloity field, p = P/n, where P is the density of momentum, E B are the eletri magneti fields. III. SHORT REVIEW OF BASIC POINTS OF RELATIVISTIC HYDRODYNAMICS In well-known ourse of theoretial physis written by L. D. Lau E. M. Lifshitz [12], [13] we an find a marosopi derivation of the relativisti hydrodynamis for the system of neutral partiles. This derivation based on treatment of the marosopi bit of medium. They start their onsideration from definition of the stress-energy tensor ˆT in the rest frame. In this ase they get ˆT = e p p p, (8) whereeis the inner energy density,pis the pressure, these quantities are defined in the rest frame, i.e. the frame (may be loal) where is no marosopi motion of the medium. Having (8) they use the Lorentz transformation to get in an arbitrary inertial frame. In the result they find T ab = hu a u b / 2 +pg ab, (9) where g ab is the metri sign onvention suh that g ab = diag(1, 1, 1, 1), h = e+pis the enthalpy density, Latin indexesa,b = 0, 1, 2, 3. The stress-energy tensor ontains information about density of dynamial quantities desribing the system, onsequently, its an be used to write equation evolution of the medium whih is an analog of the non-relativisti Euler equation, we also should write the ontinuity equation whih we present here in the terms of the four-dimensional variables j a x a = 0 (10) is the number of partiles onservation law, wherej a is the four-urrent. T ab x b = 0 (11) is the stress-energy onservation law. In formulas (10) (11) x a = (t,x,y,z), where t,x,y,z are the independent variables, whereas in the mirosopi mehanis we have (x(t), y(t), z(t), t). To write (8) we should onsider a loal frame where this piee of medium does not move. Generally speaking this piee an be involved in omplex motion, moving with aeleration relatively to some global inertial frame. Moreover, if the piee is not moving marosopially in hosen inertial frame, its partiles in their thermal motion an move with relativisti veloities. Thus, we an not use nonrelativisti formulas in this frame too. Presented in the next setion desription is free of these restritions. More general approximation in ompare with the desribed above was developed by S. M. Mahajan in 2003 [14], effetively used in many papers, for example we present several of them [15], [16], [17]. This approximation allows to get minimal oupling for the system of relativisti harged partiles. It inludes thermal distribution of partiles founds on the Lorentz invariane of equation of motion. When we onsider system of harged

3 3 partiles, in right-h side of the equation (11) the eletromagneti field emerges. Therefore, instead of (11) we have T ab x b = F ab j b, (12) where F ab = a A b b A a is the tensor of the eletromagneti field A a is the four-potential of the eletromagneti field. Introduing p (e/)a instead of p we an rewrite (12) as T ab x b = 0, (13) where T ab is the sum of the stress-energy tensor T ab the eletromagneti stress-energy tensortem ab, where T ab EM = 1 (F a g d F bd 14 ) 4π gab F d F d In some ases, see for example set. II of this paper, we an presentt ab via four-veloityu a = (γ,γu/) as T ab = a U b b U a. (14) To aount the statistial information of the partiles system in this formula in Ref. [14] a new tensor S ab = a (fu b ) b (fu a ) was introdued. It ontains parameter f, whih is the funtion of the temperaturet. Thus, in Ref. [14] a following equation was suggested U b M ab = 0, (15) where M ab is the tensor that ouples the eletromagneti the fluid fields,m ab = F ab +(m 2 /e)s ab. Following by the Ref.s [16] we an represent the spaelike omponents of equation (15) ( t + u )(fγu) = e m (E+ u B) 1 p. (16) mnγ We also represent orresponding ontinuity equation t (γn)+ (γnu) = 0. (17) Usually equation (17) onsidered as the ontinuity equation in the arbitrary frame. From this point of view equation (6) is the ontinuity equation in rest frame. However, below we show that at mirosopi derivation of the ontinuity equation in arbitrary inertial frame we get (6). Let us repeat the problem we disuss in the paper: we want to find the form of the RHD equations for the system of partiles with the temperature, so, where there is distribution of partiles by veloities,, in general, partiles an move with relativisti thermal veloities. One of the ways to solve with problem we just desribe following to the Ref. [14], but where was suggested the minimal oupling model. It inludes the ontribution of the thermal motion in the RHD s, it does not hange number of equations ( variables) in the model. However, it is interesting to see a way of widening of the model, a way of it s further development beyond of the minimal oupling. Besides, one more ruial moment, we do not want to make onnetion between the density of the momentum P(r, t) the veloity field v(r,t) for the system of relativisti partiles with the temperature, but we want to be able to solve desribed problem. Another way to find equation evolution of the fourmomentum is a use of a kineti equation. In the physial kinetis the distribution funtion f(r, p,t) is defined in the six dimensional spae of possible value of the oordinate r momentum p of individual partiles. An individual partile s momentum veloity onnet by formula (1). Introduing the distribution funtion as number of partiles in the viinity of a point r of the physial spae having momentum in the viinity of a point p of the momentum spae, following to the Ref. [2], we an write f(r, p,t) = dη δ(r+ξ r i (t)) δ(p+η p i (t)). r p From this formula we see that p onneted with the momentum of individual partiles p i (t). Thus, in the kineti equation we an write p = γmv. It is well-known that the hydrodynami funtions an be found using the distribution funtion, for example: n(r,t) = dpf(r, p,t), j(r,t) = dp p f(r, p,t). p 0 In this ase, one should use the kineti equation to derive the hydrodynami equations (see for example Ref. [18]). IV. MICROSCOPIC DERIVATION OF RELATIVISTIC HYDRODYNAMICS EQUATIONS FOR CHARGED PARTICLES To start the derivation of the RHD equations we hoose the inertial frame. Next, we onsider a sphere around eah point of spae. Eah moment of time we an alulate a number of partiles (or total mass of partiles) in eah sphere. To find the onentration of partiles in the viinity of a point of the three dimensional physial spae we should ount the number of the partiles divide it on the volume of the viinity ρ(r,t) = 1 N(r,t) m i. (18) Mirosopi number of partiles (total mass) in the viinity of the point r hanges during the time. This sum also

4 4 hanges from one point of spae to another. It is not suitable to work with the sum whih up limit of summation depend on r t. Using the Dira s delta funtion we an rewriteρ(r,t) in the following way ρ(r,t) = 1 m i δ(r+ξ r i (t)). (19) If we onsider system of partiles with equal masses we an write ρ(r,t) = mn(r,t), wheren(r,t) is the onentration of partiles. Differentiating the partiles onentration with respet to time we obtain the ontinuity equation t n+ j = 0. (20) We do not hose any partiular inertial frame. One might onsider the rest frame, but we start from the mirosopi desription, the notion rest frame does not lear from this point of view. Separate partiles have no information about: do they part of the marosopially motionless system or not? This question might be answered on marosopi sale only! Current of partiles appears in ontinuity equation, it s evident form is j(r,t) = 1 v i δ(r+ξ r i (t)) = 1 p i p 0i δ(r+ξ r i (t)). (21) For the following it is suitable to introdue the veloity field v(r,t) = j(r,t) n(r,t). (22) Differentiation quantity (21) with respet to time we get simple, but very important equation: equation of partiles urrent evolution. One has very interesting form t j α + β Π αβ = eηe α + e εαβγ η β B γ e ηαβ E β, (23) where Greek indexes α, β are used for 1, 2, 3; Π αβ is the urrent of j or the urrent of partiles urrent Π αβ (r,t) = 1 v α i vβ i δ(r+ξ r i (t)). (24) In the non-relativisti ase Π αβ oinides with the urrent of momentum. In the equation of partiles urrent evolution (23) three new funtion appear. Evident form of them are η(r,t) = 1 δ(r+ξ r i (t)), (25) p 0 i η α (r,t) = 1 η αβ (r,t) = 1 v α i v α i v β i δ(r+ξ r i (t)), (26) δ(r+ξ r i (t)). (27) In the non-relativisti limit we find that η n/m, η α j α /m, η αβ Π αβ /(m). Comparing the first the third terms in the right-h side of the equation (23), inluding the evident form of the onentration n Π αβ, we an see that the third term is proportional tov 2 / 2, exists in the semi-relativisti approximation only. So, we should neglet it in the non-relativisti limit, in ompare with the first term. In the result we have the usual nonrelativisti Euler equation m( t j α + β Π αβ ) = ene α + e εαβγ j β B γ. We make our alulations in some inertial frame, but we do not onsider transition from one frame to another. For suh transition we an use the Lorentz transformation for the mirosopi quantities, using it we an find lows of transformation of marosopi quantities, but this topi lay of the paper. Equation (23) other dynamial equations are written in the self-onsistent field approximation, we also neglet by ontribution of the eletri dipole other moments of the medium. To do it we suppose that E α i (r i,t) = E α i (r+ξ,t) E α i (r,t), B α i (r i,t) = B α i (r+ξ,t) B α i (r,t). A way of appearing of them is desribed in Ref.s [19], [20] in the non-relativisti ase. To make the paper not to overload we left this topi for next one. The eletri E the magneti B fields arising in equation (23) ( in other dynamial equations presented in this paper) satisfy to the Maxwell s equations B = 0, (28) B = 1 te+ 4πenv, (29) E = 4πen, (30) E = 1 tb. (31)

5 5 V. EVOLUTION OF THE NEW FUNCTIONS Below we derive the equations of the momentum the energy evolution, but stritly speaking we do not need them to work with the RHD. In fat we already have equation of motion -it is equation (23). Equation of the urrent evolution (23) introdue to us the three new funtions (η, η α, η αβ ) we should find a way to lose the set of the hydrodynamis equations to solve partiular problems. One of possible ways to do it is to find equations of evolution of η, η α, η αβ. Certanly we an expet that in this ase some new funtion will appear. We have evident form ofη η α, so, we an differentiate them with respet to time obtain equations of these quantities evolution. In the result we have t η + α η α = eζ α E α, (32) t η α + β η αβ = eζe α + e εαβγ ζ β B γ 2 e ζαβ E β. (33) As we expeted several new funtions appear. They evident form is presented in this setion below. But, from these equations (32) (33) we find onnetion betweenη,η α, η αβ. From (32) we see thatη α is the urrent ofη,, equation (33) shows us thatη αβ is the urrent ofη α. In general, the urrent of a quantity f, whih we designate asf α, might be presented in the form of f α = f v α, (34) where v α is the veloity field introdued above. However, we neglet by the thermal motion in this formula, but it is very useful approximation. If we inlude the thermal motion in the formula (34) it assume the formf α = f v α +f α, where f α present ontribution of the thermal motion. For example, we onsider the non-relativisti urrent of the momentum Π αβ, whih oinides with the urrent of partile urrent. In this ase Π αβ j α v β, knowing that j α = nv α, we haveπ αβ nv α v β. Inluding ontribution of the thermal motion we an write Π αβ nv α v β +p αβ, where p αβ is the tensor of pressure. Usually one onsider the salar pressure p, whih onnets with the tensor in following form p αβ = pδ αβ, where δ αβ is the Kroneker delta. Thus, we an approximately writeη α = ηv α η αβ = η α v β / = ηv α v β /. In the result, to underst meaning of η, η α, η αβ we need to underst meaning of η only. To do it we should onsider its semi-relativisti approximation. Here, we present evident form of three funtions whih arise in equations (32) (33) ζ(r,t) = 1 ( ) 2 δ(r+ξ r i (t)), (35) ζ α (r,t) = 1 v α i ( ) 2 δ(r+ξ r i (t)), (36) ζ αβ (r,t) = 1 ( ) 2 v α i v β i δ(r+ξ r (t)). i (37) It an be shown that ζ, ζ α, ζ αβ onnet with eah other in the same way as η, η α, η αβ. Thus we have ζ α = ζv α ζ αβ = ζv α v β /. VI. STRESS-ENERGY TENSOR Usually, at the RHD s desription one onsider tensor of energy-momentum. It appears at onsideration of evolution of the energy-momentum four-vetor density. To show our desription oinides with well-known we derive the equation of the energy-momentum four-vetor density evolution. The stress-energy tensor appears with. Let s introdue the density of the momentum P(r,t) = 1 p i δ(r+ξ r i (t)) = 1 m i v i 1 v2 i 2 δ(r+ξ r i (t)). (38) This quantity is a part of the stress-energy tensor T ab. Their onnetion is T 0α = P. Time part of the four-momentum vetor is the energy. So, we present the density of the energy T 00 (r,t) = 1 δ(r+ξ r i(t)). (39) As in the ases desribed above we should differentiate the definition of a quantity to find an equation of it s evolution. So, we present the equation of the momentum density evolution 1 tt 0α + β T αβ = ene α + e εαβγ j β B γ. (40) In this equation only one new quantity appears, it is the tensor of momentum urrent T αβ (r,t) = 1 p α i p β i δ(r+ξ r i (t)). (41) In the same way we find the energy density evolution equation 1 tt 00 + α T α0 = e je, (42)

6 6 wheret α0 is the urrent of energy, whih simply onnets with the momentum density T α0 = T 0α. Analogously to the non-relativisti ase we an write T αβ = P α v β + p αβ, where p αβ pδ αβ is the pressure tensor pis the isotropi pressure. Thus we an rewrite equation (40) in the following form t P α + β (P α v β +pδ αβ ) = ene α + e εαβγ j β B γ. (43) Evolution of the energy densityt 00 the momentum density P depends on the onentration n, the partiles urrent j, the eletri E the magneti B fields. However, evolution of n, j, E, B do not depend on P T 00. Therefore, to have the losed set of the RHD equations we do not need to onsider P T 00 evolution, put our attention forn, j η. If we put p = 0 P = np in equation (43) inlude equation (6) or (20) (whih is the same) we find that equation (43) mah with the (7). At studying of the relativisti plasma an approximation is used for the pressure p. In order to lose the set of the RHD equation one use the equation of state for an ideal gas: p = nk B T, where k B is the Boltzmann onstant T is the temperature. It we onsider evident form of Π αβ T αβ we will see that they have the same tensor struture, i.e. they both depend on v α i vβ i. Thus, we an use for π αβ, where Π αβ = nv α v β + π αβ, the same approximation as for p αβ. Consequently, we an put π αβ = φ δ αβ φ = φ(n,t). VII. CLOSING OF SET OF THE HYDRODYNAMICS EQUATIONS During the paper we attain that we an desribe relativisti plasma by means ofn,j α,η α,e α, B α. In this ase we inlude some of effets aused by temperature, but some of them we lost. To study relativisti plasma with the large temperatures we should onsider ontribution of thermal motion at least in η α η αβ. However, we suppose not to onsider it in the paper. We already introdued new quantity η, we should underst it s meaning. For this purpose we onsider semi-relativisti limit for η, find out it s ontribution in dispersion of waves try to approximately express it via n η. Even after getting approximate onnetion between n, v, η, we suppose to onsider η as an independent variable, along with onentration n veloity field v. Thus, the semi-relativisti approximation of η has form η = n/m Π αα /2m 2, where Π αα is the trae of the tensor of the urrent of the partiles urrent (24), Π αα = nv 2 +3φ, whereφis the urrent of the partiles urrent on the thermal veloities. We have got it using the formula 1 = 1 p 0i m i (1 v2 2 2 ), whih is the semi-relativisti approximation for the inverse time omponent of the four-momentum. VIII. DISPERSION OF LONGITUDINAL WAVES IN RELATIVISTIC PLASMA We suggest that in equilibrium state the relativisti eletron plasma (we suppose that ions motionless) is desribed dy following parametersn 0, η 0, v 0 = 0. To study dynamis of small perturbation we use equations (20), (23), (30). For the first step we suggest that φ depends on onentration n only, but below we will aount that φ depends on n η. We notie that in non-relativisti limit φ beomes pressure p, η beomes onentration n, dependene on n η redues to dependene on n. Considering evolution of small perturbations around the equilibrium state we an find it s dispersion dependenies, whih has form of ( ) φ ω 2 = 4πe 2 η 0 + k 2, (44) n 0 we an see that η 0 appears instead of the equilibrium onentration n 0, φ emerges instead of the pressure p. In the semi-relativisti limit from (44) we get ω 2 = 4πe2 n 0 m 6πe2 φ 0 m 2 + ( ) φ k 2, n 0 where we use that Π αα 0 0+3φ 0, sine v 0 = 0. Consequently, in the non-relativisti limit we have ω 2 = 4πe2 n 0 m +3k BT m k2, (45) in the seond term the equation of state of the ideal gas was used at the adiabati ondition with the rate of adiabat equals 3. Inluding that φ depends on both onentration n η we still get the formula (44). However if we onsider the wave of partile onentration in an eletron beam we obtain dispersion equation (ω ku) 2 1 m k2 n φ+ 4πe2 η0 2Uk mn 0 (ω ku) ηφ 4πe 2 η 0 = 0, (46) whih desribes the longitudinal waves in the eletron beam. It is the equation of the third degree in ontrast to (44) whih is equation of the seond degree. Thus we an see that an aount of φ(η) dependene ould lead to some new effets. However, to get solution we should present equation of state φ = φ(n,η,t), whih has not found yet. Presented equations it s onsequene give us a lot of new open questions, but they also give another view on the development of the RHD.

7 7 IX. CONCLUSION We have presented mirosopi derivation of the relativisti hydrodynamis equation. We have derive as well as well-known new equations. Among well-known equation we an mention ontinuity equation, momentum balane equation energy balane equation. We have presented the partiles urrent evolution equation, sine the partiles urrent simply related with the veloity field j = nv. At this equation derivation new funtions have appeared. To underst their meaning to onsider their influene on partiles dynamis we have derive equations of theirs evolution. We have suggested the losed set of the RHD equations used its to onsider dispersion of the olletive exitations in the relativisti plasma. reevpa@physis.msu.ru [1] R. Hakim. Introdution to Relativisti Statistial Mehanis Classial Quantum (WS, 2011) [2] L. S. Kuz menkov, Theoretial Mathematial Physis 86, 159 (1991). [3] P. A. Andreev, L. S. Kuz menkov, M. I. Trukhanova, Phys. Rev. B 84, (2008). [4] P. A. Andreev, L. S. Kuz menkov, Phys. Rev. A 78, (2008). [5] P. K. Shukla, B. Eliasson, Rev. Mod. Phys. 83, 885 (2011). [6] F. A. Asenjo, V. Munoz, J. A. Valdivia, S. M. Mahajan, Phys. Plasmas 18, (2011). [7] F. Haas, B. Eliasson, P. K. Shukla, Phys. Rev. E 85, (2012). [8] R. D. Hazeltine S. M. Mahajan, ApJ, 567, 1262 (2002). [9] A. Bret, ApJ, 699, 990 (2009). [10] A. Bret, C. Deutsh, Physis of Plasmas, 13, (2006). [11] A. Bret, M. E. Diekmann, Physis of Plasmas, 15, (2008). [12] L. D. Lau, E.M. Lifshitz Hydrodynamis, Vol. 6 of Course of Theoretial Physis (Pergamon, London, 1959). [13] L. D. Lau, E.M. Lifshitz Field theory, Vol. 2 of Course of Theoretial Physis (Pergamon, London, 19????). [14] S. M. Mahajan, Phys. Rev. Lett. 90, (2003). [15] S. M. Mahajan, Phys. Rev. Lett. 100, (2008). [16] F. A. Asenjo, V. Munoz, J. A.Valdivia, T. Hada, Physis of Plasmas, 16, (2009). [17] F. A. Asenjo, F. A. Borotto, A. C.-L. Chian, V. Munoz, J. A.Valdivia, E. L. Rempel, Phys. Rev. E, 85, (2012). [18] L. D. Lau, E.M. Lifshitz Physial kinetis, Vol. 10 of Course of Theoretial Physis (Pergamon, London, 1981). [19] M. A. Drofa, L. S. Kuz menkov, Theoretial Mathematial Physis 108, 849 (1996). [20] L. S. Kuz menkov P. A. Andreev, will be presented in PIERS Proeedings, Augoust 19-23, Mosow, Russia 2012.