Copy_____ THE UNIVERSITY OF MICHIGAN COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRICAL ENGINEERING Radiation Laboratory EVOLUTION OF AN INHOMOGENEOUS RAREFIED PLASMA AND ASSOCIATED HIGH FREQUENCY ELECTROMAGNETIC RADIATION by R. K. Ritt 2764-10-T - RL-2078 April 1962 The work described in this report was partially supported by the ADVANCED RESEARCH PROJECTS AGENCY, ARPA Order Nr. 120-61, Project Code Nr. 7400. ARPA Order Nr. 120-61, Project Code Nr. 7400 Contract DA 36-039 SC-75041 Department of the Army Project Nr. 3A99-23-001-01 Contract With: U. S. Army Signal Research and Development Laboratories Fort Monmouth, New Jersey Administered through: OFFICE OF RESEARCH ADMINISTRATION * ANN ARBOR

THE UNIVERSITY OF MICHIGAN 2764-10-T ASTIA AVAILABILITY NOTICE: Qualified Requestors may obtain copies of this report from ASTIA.

THE UNIVERSITY OF MICHIGAN 2764-10-T EVOLUTION OF AN INHOMOGENEOUS RAREFIED PLASMA AND ASSOCIATED HIGH FREQUENCY ELECTROMAGNETIC RADIATION by R. K. Ritt 2764-10-T April 1962 ARPA Order Nr. 120-61, Project Code Nr. 7400 Contract Nr. DA 36-039 SC-75041 Department of the Army Project Nr. 3A99-23-001-01 The work described in this report was partially supported by the ADVANCED RESEARCH PROJECTS AGENCY ARPA Order Nr. 120-61, Project Code Nr. 7400 Prepared for THE ADVANCED RESEARCH PROJECTS AGENCY and the U. S. Army Signal Research and Development Laboratories Ft. Monmouth, N. J. I I i

THE UNIVERSITY OF MICHIGAN 2 764-10-T TABLE OF CONTENTS ABSTRACT iii I INTRODUCTION 1 II THE TERMINAL DISTRIBUTION 3 III THE PERTURBATION 9 IV THE OSCILLATION FREQUENCIES 16 V CONCLUSIONS 25 APPENDIX A APPENDIX B APPENDIX C REFERENCES ---- I —__________ J

THE UNIVERSITY OF MICHIGAN 2764-10-T ABSTRACT If N(x) 0= a-< x a } No x r x > a is the distribution function of the positive ions in an electron-ion gas, if the ions are held fixed, and if the initial distribution of the electrons is v2Ii N(x) e 22 we first calculate the limiting electron distribution as t -4 oo. Assuming that the true electron distribution for finite t is a perturbation from this limiting distribution, we find that the perturbation equations can be satisfied by transvers electromagnetic waves with certain discrete frequencies, is. These frequencies fall into two classes, their real parts having, respectively, the form R(i sn) 2, and R(is ) - -- n 2a n h e3/2 where c = velocity of light in vacuum, and h = -, where w is the plasma frequency. plasma frequency. I iii

THE UNIVERSITY OF MICHIGAN 2 764-10-T INTRODUCTION This report presents a highly simplified model of the electron distribution in the rear of, and slightly behind, an object traveling through a rarefied, totally ionized gas; its velocity is much greater than the rms velocities of the ions and much less than the rms velocities of the electrons. It has been shown, by Dolph and Weil (1959), and by Sawchuk (1962), that behind such an object there is a 'hole' in which the ion density is very low in comparison to the ion density outside the hole and that the electron-ion configuration is electrically neutral. Our first simplification is to replace this hole by a well defined vacuum, and to make the problem essentially one-dimensional by letting the hole be the region between two parallel planes, 2a units apart. We then investigate the manner in which the hole 'fills up'. The time interval is one in which only the electron motion is considered, the ions being assumed not to move. We first assume that the only forces on the electrons are the Coulomb forces of the ions and other electrons, and that the phenomenon is governed by Vlasov's equations. In Appendices A and B. we extend a result of Iordaneskii ( 1951 to show that these equations have a unique solution, and that the solutions have a limit as t-o co. In Section II, we find the limiting electron distribution and electric field. The limiting distribution we call the terminal distribution. In Section III, we assume that the forces on the electrons are electromagnetic forces, and regard the electron distribution and the electromagnetic field as perturbations of the terminal distribution and the electrostatic field found in Section II. Writing the appropriate linearized equations, and assuming that the electromagnetic field is transverse, we secure a partial differential equation for the electric field. When we take the Laplace transform of this equation (as a function of the complex value, s ), we find that for large s, the poles of this 1

THE UNIVERSITY OF MICHIGAN 2764-10-T Laplace transform can be evaluated by fairly elementary methods. Appendix C provides a mathematical justification for the use of these methods. In Section III, we represent the electron distribution by approximating, linearly, the true electron distribution, which enables us to find, routinely, asymptotic representations of the poles of the Laplace transform. These are of two types. If h = e 3/2X, where X = a/! V w, a being the rms velocity of the electrons, w being the plasma frequency (we call X the Debye length, and it differs by a numerical factor from the usual Debye length), and if a/h ~>> 1, then one set of poles have the form s [ -C log n7rc + nir ] n a a 2 j and another set has the form c c 3 h +n.1 (c is the velocity of light in a vacuum). We hesitate to attempt a heuristic explanation for these poles, which correspond to very high frequency, damped, electromagnetic radiation, except to say that they are a consequence of the collision-free Boltzmann and Maxwell equations. The electron-ion configuration appears to behave as some sort of resonant cavity, which selects from the frequency spectrum of the electron motion certain privileged ones which are allowed to propagate. m 2

THE UNIVERSITY OF MICHIGAN 2764-10-T II THE TERMINAL DISTRIBUTION We consider a gas, consisting of ions and electrons, which at t = 0 has the following properties: a) The ion number distribution is [No, /xl >I a N(x) = 0, lxl < a b) The electron distribution function is N(x) V(v), where v2 V(v)= --- e 2-2O < v< O 1 -ooav<oo We assume that the evolution of this system for t > 0, is governed by the Vlasov equations f + 8f q _. E Sf t0 at ax m ' av (2-1) ax K \f(x, v, t)dv-N(x). -00 In (2-1) the functions f and E represent the electron distribution function and the electric field, respectively. The constants, q, m, e0 represent respectively the charge and mass of the electron, and the dielectric constant of free space. The use of (2-1) carries with it the tacit assumption that the ions do not move, that the forces on the electrons are electric Coulomb forces, and that the electron and ion densities are large enough to permit the use of a distribution function. Equation (2-1) must be solved subject to the initial condition f(x, v, 0) = N(x)V(v). Since the initial data is not differentiable, there\are, in general, no i - 3

THE UNIVERSITY OF MICHIGAN 2764-10-T solutions to (2-1), unless we enlarge our definition of a solution to mean a pair f(x, v, t), E(x, t), where aE/ax exists almost everywhere and satisfies the second equation in (2-1), and f(x, v, t), instead of satisfying the first equation of (2-1), remains constant on the characteristics of the first equation of (2-1). In Appendix A to this report, we show that with the additional restriction, | E(-oo, t) = 0, a unique solution of the type just described exists, and in Appendix B, we show that as t-oo, the pair f, E have limits, g, G, where g and G satisfy vag - Gg = vx m av (2-2) ax_= _ e 0 g(x, v)dv -N(x) ax C o It is further shown, in Appendix B, that N No 2 2 cr g= e e 2r (2-3) =- G ax m We shall call g the terminal electron distribution and W(x) = g(x,v)dv= N e -o0 the terminal electron density. In Appendix B it is observed that /(x) = V(-x) and G(x) = -G(-x), and that /(co) = G(co) = 0. It is then sufficient to study the equations only for 0 l x a oo, with G(0) = 0. Using (2-3) and the second equation in (2-2), and letting, w 4

THE UNIVERSITY OF MICHIGAN 2764-10-T q P (x) = - m G, we obtain the pair of equations dx = dx } (2-4) /q N0 In (2-4), =7 --, which is the "plasma frequency" and; (oo) = (0) = O(oo)=0. o To solve this system we consider separately the regions 0 < x < a, a x, and use the fact that both i and I are continuous (but not a / a x ). The region 0 < x < a. In this region, N(x) = 0, so dx TXL d 2 2 = c e dx A first integral is easy to obtain, and using the fact that i (0) = 0, we have fr2 0? 2 =2 2 c2 e -e - o in which o M = ~(0), and is yet to be determined. If z2 = e, it is elemer 0 0 to derive WZ = =v cu z tan rX x =u2lo z 2 se z o ] = a2 log z sec 2 a (2-4') (2-5) itary (2-6) - -- -- -- 5

I THE UNIVERSITY OF MICHIGAN 2764-10-T The region a < x. Here, N(x) = No, and dx - (2-4") e(y2 Using the fact that i (oo) = /(oo) = 0, we obtain a first integral I- o 2 = 2 2C2 e - 1 - 2-. (2-7) The solutions to (2-4") can then be obtained by quadrature, but the integrals are non-elementary. However, using (2-5), (2-7), and the continuity of i and /, we obtain 2 O(a) z = 1+ O a, or, from (2-6), - wz a) z2 l+log z2 sec2 0~ (2-8) Equation (2-8) can be used to determine z, if it is required that b and i be bounded for 0 1 x < a. In this case, the solution of (2-8) corresponds to the first intersection of the graphs of the functions 2 2 z0- 1 -logz0 wz a and logsec2 { Jf~- U wa which occurs for 0 < z < 1. It is easy to verify that as -- oo, z —0. 0 o0 Let us define the Debye length, X, by the equation X -, (2-9) 6

THE UNIVERSITY OF MICHIGAN 2764-10-T and let us proceed under the assumption that X /a is very small. Now for 0 < x < a, the electron density function can be given explicitly, 2 2 Z0X W(x) = No o sec2 2 Using equation (2-8), we have zo -1 W(a) = N e, and o Wt(a) = W(a) e -zo Therefore, approximating W(x) linearly, for 0 < x < a, we obtain a graph (Rigure 2. 1).,N ez2-1 No o e X [ez2-1 21- 1/2 e el 0. 1 FIGURE 2. 1 The dotted line in Fig. 2. 1 is a more accurate sketch of the curve. 7

- THE UNIVERSITY OF MICHIGAN 2764-10-T For x > a, W(x) ->N. This approach, for large x can easily seem to be exponential. But approximating W(x) linearly to the right of x = a, we obtain the graph (Figure 2. 2) 1/2.(e-l)e' X,A /i I a FIGURE 2.2 On the basis of the above reasoning we shall approximate the terminal electron distribution by the graph (Figure 2. 3) h =e X FIGURE 2.3 -" 8

THE UNIVERSITY OF MICHIGAN 2764-10-T III THE PERTURBATION The computation in the preceding section assumed that the forces on the electrons are Coulomb forces. A more accurate set of equations would have been the Boltzmann-Maxwell system: at af - (E -+v x H)af iy- (E'+/ E ovxH)-y0 at +V_ ax - -- aEH 7xE+ p ^ - O - o at V7 —E =- f d In (3-1), f is the electron distribution function; E, H are the electric and magnetic components of the electromagnetic field. Instead of attempting to solve these equations we shall employ the following device: At t = 0, we shall assume that the fields and the distribution function have the following form: E= G+ eE - O Hf = gH (3-2)+ f -o f=g+ef In (3-2), e is the perturbation parameter, G is the electrostatic field obtained by taking G of the preceding section and allowing it to be the x-component of a vector field, and g is the terminal electron distribution, modified by replacing v2 by 9 ---

THE UNIVERSITY OF MICHIGAN 2764-10-T v- and changing the normalization from 1 1 to V27 a (27)3/2 a3 If the terminal distribution is truly determined by the methods of the preceding section, and if the distribution is truly determined by the system (3-1), then for some to< oo, the form (3-2) is justified for t > to. As a matter of labeling we set the value of to = 0. We then insert (3-2) into (3-1), and disregard all terms which are not linear in e. We then obtain the linearized system: af af q &f -+v - G -E vxH - -0 at - ax m- av mE +O- av aH x Eo+o at — 0 aE 7xH -~ -o Ff d3 7x Ho- E - - f d (3-3) q 3 d - E = - f d v V- H =0 -0O Since ag /av is proportional to v, the term involving H, in the first equation of (3-3), vanishes. Because of the linearity of the system (3-3),we can now seek a solution of (3-3) in which the electromagnetic field is transverse; that is, E has - O only the y-component, E, and H has only the z-component, H. The equation then becomes becomes 10

THE UNIVERSITY OF MICHIGAN 2764-10-T Ofo fof q ag + G(x) - E Oz 3z _E aH 0 3 ax o at - aH aE f 3 Ox + eO - = qv v f d - Ox o at y o - E q C 3 ay f d vfd aH = qJV f d3 3 v V f3v 0 ay i x o o If we multiply the first equation in (3-4) by v and v, respectively, integrate with respect to v, and make use of the remaining equations in (3-4), we obtain the equations O (OE 1 q E a2E 1 a2 E wo 2 C a2 1 ai i2 E/ 2 E t) n = f0 O0 q Ofi0 q O 3x2- c2 at2 E e E=/~o ~ ~ G(x) E at ax' x av m av X y aE In addition, E = E(x,y,t), since ~ = 0. We shall now use the third equation of (3-5) to determine the right hand members of the first two equations in (3-5). 11

THE UNIVERSITY OF MICHIGAN 2764-10-T At time t = 0, let f f (x, y, z, v,v,v ). Let x(x, v, t), v(x,v, t) be the o o x y z o x o x solutions of the equations: dx dt x dvx q d t- G(x); x = x oVx=v when t = 0; i. e. the equations of the electron trajectories in the unperturbed system. Let x (x, v, t), v (x, v, t) be the inverse of these functions (see Appendix A). Then, taking into account the fact that vx ~gX- q G(x) Vx= 0 and applying the method of characteristics to the third equation in (3-5), we obtain the equation f (, )=f (x (x, v%, t), y-vyt, z-v t;v(x, vx, t), vy, vz ) + ag(x,vxv,vz) (t (3-6) y 0 Let us designate the second member of (3-6) as f. It is easy to see that r af(l) 2 r t 0 G(x) _ I E(x,(x,v,t-r),y-v (t-/')dZ~ - m- a a J at xo(x t-, y-vy(t- d Hence, the second equation in (-5) becomes Hence, the second equation in (3-5) becomes 12

THE UNIVERSITY OF MICHIGAN 2764-10-T 82E 1 a2E w2 0m a aX2 Q C2 ae C2 In ey av ' t * at E(xo(x,vxt-),y-vy(t —?),)d d v= (3-7) 2 0 ~ Aoq2 ao,3 = Vv d v, m y ' - x in which fO has the arguments given in (3-6). An analogous equation can be obtained O in place of the first equation in (3-5), and in principal, these two equations can be — used to find E. However, we shall concern ourselves onlywith (3-7), whose right hand member we shall designate as p. If we take the Laplace transform of (3-7) we obtain: r. 2 2 22 2O,2f q dX2 /(E)- 2I o(E) + m (3-8) o(p) +s E(x,y, O)+ E(x, y, ). t Now, we shall concern outselves with the ratio E(xo(xvxt-T),y-Vy(t-r't),')d --. (3-9) 4(E) 13 ---

THE UNIVERSITY OF MICHIGAN - 2764-10-T If we examine the numerator we see that it can be written as 02 [a E(xo(x,vx,t-?'),y-Vy(t-r),?P)d T-E(x, y, t) at 0 =s h -st E(xi(xVX,t-X), y-v (t-h)e, ri )ddtec fe 0 0:ing the substitutions tu=st, X = sZ?, the ratio becomes,(E) By mak J e- E(x,y, ) 0o Now E(xo(x, vx, (M-X)z), y-vy(p-X)z, Xz) can be expanded in powers of z, giving the series E(x, y, 0)+ X - at x V aty a... t x + ay -o Similarly, E(x,y, pz) = E(x,y, 0) +p at at=0 z+... If, in these expressions, z is replaced by, and it is recalled that (3-9) must be inserted into the integral in (3-8), so that the contribution of the term involving /aE aE \ Vx iX+Vy 4gyt=o will be zero, it is easy to see that the ratio (3-9), inserted into (3-8), is O( s ). s 14

THE UNIVERSITY OF MICHIGAN 2764-10-T Therefore, as s-+oo, o ay aat-E(xo(x,Vxt-7),y-vy(t-?.r)drd3 m vy t x y 0 = 2 - e Z(E)Oo(( ) The analytic continuation of the Laplace transform^ into the left half plane will satisfy the same relationship. Therefore, we can write (3-8) in the form: d2 — 2 1 - >(E):(E)= (P)+ s E(x, y,0) + c (3-10) a + a E(x,y,0). We wish to find the poles of/(E) associated with large values of I s. These poles will correspond to high frequency components of E. In order to find these poles, we must find two linearly independent solutions of the homogeneous equation associated with (3-10), YI, YII, which respectively are asymptotic to + Ir2+2J x Cz e as x-++ oo (see Appendix C ); the poles of '(E) will then be the zeros of the Wronskian of I and yII. Because our evaluation is asymptotic, we can neglect the 0( ) term in (3-10), and the problem of finding high frequency components of E is reduced to the problem of finding, asymptotically, the zeros of the Wronskian of the functions yI, yII, associated with the differential operator d - s2 + e *'y02,,dx2 -2 Y 15

THE UNIVERSITY OF MICHIGAN 2764-10-T IV THE OSCILLATION FREQUENCIES We have now to compute the zeros of the Wronskian of the two linearly independent solutions YI, YII of the equation yTT (s2+w2e )yO y, s 0 i x for which yi e as x- oo, and y e as x- -oo, where 1/2 si = (s 2+2) the function formed by putting a cut from -i U to i ao, and choosing the branch which is positive for positive real s. Because i(-x) = / (x), we can let yII(x) = yI(-x), and W(yI, y) = 2y (O)y(0). Therefore we need only to find y (x), and only to consider non-negative values of x. Instead of the exact function e, we shall use the approximate function found in Section II: 1 x >x2 /2 X - (x-xl) x < x < X2 (4-1) 0 0<x <xi, 3/2 1/2 where h = e / X (X being the Debye length), xi = a -e /X, x2 = xl +h. In the course of finding yi, we shall retain only the dominant terms in the asymptotic expressions, and we shall discard factors which will not affect the computation of those large values of s which annihilate the product y,(0)yi (0). For x > x2, from (4-1), we see that /s e, 16

THE UNIVERSITY OF 2764-10-T MICHIGAN - so that dYI s (X2) = - dx c (4-2) If the new variable - c2h a = ( -j (jr 2/3 Li t (4-3) is introduced, the equation to be:solved, \for xl < x < x2, is dz2 -y= 0 Z1 Z< Z2, where 2 as1 Z = 2 2 2 as2 Z1 =-2 c Further, dy_ = /2 d dz dx For this equation we choose the two solutions 2/3 z P(z) 1/4 { -2/3 z3/2 Q(z)- e 1+ 5 +3 +.. 48z3/2 48 z3/2 I these representations being asymptotic for large I z. Retaining just the dominant terms, we obtain 17

THE UNIVERSITY OF 2 764-1O-T dP~ 1/2p~ 1/4 e2/3 z / MICHIGAN -n dP -z/ dz N+ z dz ~Ql z/2 dz p 1 -3/2 1/4 e2/3 z 11 4 1 -3/2 1/4 e -2/3 3/ Q 1i..z/4 e-2/3 3/2 (4-4) Now for x1 <x -x2 =, AP+ BQ, where A and B are constants. Since the values of y I(z2), ~~Z da (X2) are respectively 1 and -a from (4-2), we have IV AP(Z2)+ BQ(Z2) = 1 1 (/2 AP?(Z2)+ BQY(Z2) = - These equations have the solution A-Q'(z9)+ Qz z A= W(P, Q) B = - P I(Z2)+ 42 P(Z2)-.1 ~ WCP, Q) so that., using (4-4),, we obtain 1 -3/2 -2/3 2 zI e P~) 2/3 z3/2e P~z)+ e 2 Q(Z),$ Thus Y(l 1 2"2-3/2 e-2C/3j Z2 I e P(zi)+e2/ ` Q(z1) (4-5) dy1 1 -/ -2/ 3/2 q O 1l/2 Iy(X )"'-z e/z Q"2 P)1+e Q'( zl) 18

THE UNIVERSITY OF MICHIGAN 2764-10-T S~ S Now, for 0 < x x, = C e c + D e c thus, the quantitites S S S S C e X + D e- cX z1 (C ec l-D e c X) are respectively asymptotic to the right hand members of the equations in (4-5). Before solving for C and D, we observe that 2s 3/2 3/2 2 ok f(2+ 23/2 23 h (2+) /2 -s - +w) 31 23hC23/231ssh 3 c3 - 3 c -s ~ Then solving for C, D, and discarding factors, as before, we obtain: 1 ( -sh/c -3/2 _esh/c -3/2 )eS(X-X )/c + -(xe z )e 1 e-sh -3/2 -3/2 +eahc s(x-x)/c +(64e Z1 Z1 + )e Then, expressing zl and z2 in terms of s, we secure the fact that the values of s, for which the product YI(O) y, (0) vanishes,, are obtained, asymptotically, by solving the equations 1 w2c -sh/c sh/c -sxl/c +(4c2 esh/+esh/c sxl / -)e _h256 e +e )e =0. 8 hs3 64 h (4-6) We shall now give a discussion of approximating solutions of (4-6) for certain - ranges of values of s. The equations (4-6) were derived by using certain asymptotic forms of the solutions of the Airy equation; therefore we are tacitly assuming that we are dealing with values of s which cause the variable z, in the above discussion, to be large. Now Izj- l 2-''( C2hc )2/3 12 19 --

THE UNIVERSITY OF MICHIGAN 2764-10-T If S slI I / z 4/3 Ihi ' hw This parameter is large if the plasma frequency, w, is small in comparison with the frequency -, or equivalently, since hcvcr, if the r.m. s. velocity of the electrons is small in comparison with c. This is an assumption we shall make. If I I c I I1 / c I4z3 (h )2 Isi -1 -1 x zi) h Since we are assuming\that ~<< 1, whether [z) will be large, for such values of s, is dependent upon the exact values of the parameters, and cannot be decided on the basis of the general assumption which we made in the case of I sI| v h Hence, for I sI,,,, the above analysis is invalid unless the corresponding value of I|z is large. Therefore, in the contrary case, a separate analysis must be given. We can now proceed. CasI Is c ( c )4/3 h 2 Case I: Is], I h ( (-~x >>1. Here, s | Ax ' << 1, so we approximate the first parenthesis in (4-6) by -2sh/c, and the second parenthesis by (1 + 4w c /64 h s ). Let us examine the term c /h s. If 4 2 4 6 4 4 C | c S6c - he = X w h x 6 ( v 6 h2s6 h2 c4 But from ( )4/3 h)2>> 1 hwo x we have 4 xl 6 c h ~<< h4w4 20

THE UNIVERSITY OF MICHIGAN 2764-10-T so 4 2 w c 26 << 1 h s Thus we replace (4-6), by the equation esx1/c + 2 (4-7) sx!4 Letting z =, this last equation is equivalent to the four equations c z + t ez +i xl (4-7.) e=+ - e =+ (4-7') 2cz ' 2cz Now, letting z = -x + iy, the solutions of (4-7'), with Re z < 0, correspond to the intersection, for x > 0, of the curve 2 2 W2 X 2x 2 y = e -x 4 c2 with either of the two curves x tan y = y, x = -y tan y. A simple sketch of the graphs of these curves shows that the intersections take place for + nwT yn~ 2 cn7r x log -, n=l, 2,... n oWx1 Therefore the associated complex oscillation frequencies are c f cnr- + n'r i n x1 x1 x 2 J Because xl a, we prefer to write this as 21 I

THE UNIVERSITY OF MICHIGAN 2764-10-T a n-log + n- (4-8) n a ucowa 2 It is interesting to make the following observation. Suppose that instead of the approximation (4-1), we had used instead the approximation e e/2= 1, (xl > a (4-9) =0, Ixl < a This, in effect, assumes that h/xl is so small as to have no influence, and is consistent with the negation of ( ) 4/3(J ) >> 1 hw Xi -S1x/C In this case, it is easy to verify, by letting I = e for x > a and YI = A e+ Bes for 0 < x < a, that the condition for y (O)yI(O) = 0 is 2sa/c + W2 e -4s2 which is the same as (4-7). Thus, the same estimates (4-8), are obtained. It is reasonable to conjecture that regardless of the value of c 4/3 h )2 hw Xas long as h/x << 1, the values (4-8) will be obtained. We have not atte mpted this analysis, because in this report we are primarily concerned with those values of s for which Isl h, which case we shall now consider. Case II: IsI v -. We shall rewrite (4-6) as X1 1 e —T ^C z e +i(h) z e 22

THE UNIVERSITY OF MICHIGAN 2 764-10-T CZ by setting s 2h Since we are assuming s = 0 ( ), z = 0(1), and since Re(z) 0, the second term in the numerator may be neglected. Thus, we arrive at 14-10) e =3 (4-10) c )2 3 z l+(-h) z e To estimate the zeros of (4-10) we first set z = 2rni and observe that the left numbe is 1 and the right number is -1 [ c 2 3 27rnix,/h 1 +( )(27r ni) e r / However, if z is given a small negative real part, the left number remains close to 1. However, if z is given a small negative real part, the left members remains close to 1, and the term c 2 3 t-z (-) z e xi -- is greatly diminished, because - >> 1. Thus we expect to find solutions in the neighborhood of 27rni. To see this more clearly, setting z = -3y +iy, and equating the phase and magnitude of both sides of (4-10), we obtain _1 1 sin( y-3tan -) e'y = sin( hy-3 tan-1 ) (4-11) e2jy-2e/ycosy-1= (h)4y6(1+32)3e2 (y Let / > 0, and consider the second equation of (4-11). When y = 0 both sides are zero. If we call the left side f(y) and the right side g(y), we see that f(y)- oo as y-+oo, and that g(y) has a maximum at y = 3h/xlx, and for y > 3h/xl, g(y) decreases 23 --

THE UNIVERSITY OF MICHIGAN 2764-10-T monotonically to zero. Further, 3h c 4 h 6 1+32) (1+132)3 (= ( ()6 >> (e)6 since 2/3 On the other hand, 3 6bh 3h 3h f(13) eX1 -2e cos(3) -1 so that for /3 sufficiently small, 3h 3h g(1) > f(-) ~ But since f(y)-*oo as y- oo and g(y)- 0 as y-coo, the second equation of (4-11), 3h has a branch y(3) > - for sufficiently small.3 Now, examining the first equation of (4-11), we see that for /3 = 0 the values of y which satisfy the equation are the solutions of sin( y — ) = sin (y- 3) or, since x2 - xl = h, n = 2n7r. Thus, this equation has branches y n(3) for which Yn(0) = 2n7r. Further it is an elementary calculation to show that dyn(3) d/3 = Yn(0), so these branches actually extend into the region / > 0. These branches intersect 3h 3h the branch y(3) > x-l for n ^xy(0); the corresponding complex frequencies, are then given by the formula s - 3 h i] +( c -2 - ni (4-12) nh L a J 24

I -- THE UNIVERSITY OF MICHIGAN 2764-10-T V CONCLUSIONS We find that as the terminal distribution is approached, transverse electromagnetic fields can be generated, having the frequencies given by (4-8) and (4-12). Whether these attenuated oscillations are detectable depends upon, of course, the validity of our model, and upon the energy which is radiated. The energy available for the production of these radiations is certainly no greater than the energy necessary to create the initial ion-electron distribution of Section II. In order to find the amount of this energy which goes into the electromagnetic fields described above, it is necessary to solve the equations (3-1), which is just what we have managed to avoid doing. However, the nature of the attenuation in (4-12) gives hope that these oscillations are detectable. In regard to the model used, we have reproduced our calculations for a cylindrical geometry; the material in Section II is more complicated, and we have not been able to obtain a value for h, short of numerical computation for fixed parameter values. Aside from this, however, the work proceeds as above, and identical formulas are secured. 25

THE UNIVERSITY OF MICHIGAN 2764-10-T APPENDIX A We are concerned with the existence and uniqueness of solutions to the system v +A(x, t)-=O, - =+ f(x,v,t),dv-N(x), at ax av ax -- -00 subject to the initial conditions f(x, v, o), =N(x)V(v), and the boundary condition A(-co, t)=0. Iordaneskii (1959) sketched a proof of such a theorem; our proof follows in the main the outline given by Iordaneskii. The arguments are sufficiently delicate to warrant their detailed exposition, and since our hypotheses are slightly different, we give a complete proof. The hypotheses on N(x) is that it satisfy 0< N(x) < N0, and that[N(x)-N] be summable on -oo < x < co. The hypotheses on V(v) is that it be continuously differentiable, that it be positive, monotone decreasing for increasing Ivl, and that it be summable and have a second moment on -0o < v < 0o. The problem is replaced by the one of finding an A(x, t) having the properties that if x (x, v, t), v (x, v, t) are the characteristics of the system dt dx dv 1 v A ' then 1 A(x,t)= = _ N(xo(r,v,t)V(vo(Q,v,t) )dv-N(r) d. We shall first assume that N(x) is continuous, which is part of Iordaneskii's hypothesis, and then pass to the more general situation. Let A(x, t) have the following properties: - -- - - -- m 26

THE UNIVERSITY OF MICHIGAN 2764-10-T (a) II Al < P(t), P(t) continuous [IIAl = SUPJA (b) a~ is continuous in x, and I | | <p(t), p(t) continuous ax a x _____ (A) (c) lim A(x,t) = 0 x- -oo (d) A(x,o)=0 For brevity, when A(x, t) satisfies all the properties (A), we shall write A(x,t) e(A). Lemma 1 Let A(x, t) e (A). Let 0(xo,v,,Vo (x, t), t) be the unique solutions of the system dt (C) _ = A(0,t) dt which have, respettively, the values xo, vo at t = 0. Then 0(xo, vo, t), (xO, vo, t) have continuous partial derivatives with respect to xo, vo; 3A - =^, forallt; axo avo a Vo axo and 5t I0(xo, vo, t)-(Xo+vot)l 1 (t-P)p(r)d?' 0 (xo Vot)-Vo, P( dr0 0 27

THE UNIVERSITY OF MICHIGAN 2764-10-T Proof: The inequalities are trivial, and the rest of the statements are well-known. Lemma 2 Let 0(xo, vo, t), O(xO, vo, t) be the same as in Lemma 1. Then the mapping x = 0(xo, vo, t) v = V(xo, vo, t) is one-one and onto the (x, v) plane from the (xo, vo) plane. The inverse functions xo = xO(x, v, t) v = v(x,v, t) satisfy the inequalities x (x, v,t)- (x-v t) r )d 0t [v(x vt)-vl P()d?' O 0 and also satisfy the equations av av av o 0 0 - +v + A(x,t)- =0 at ax av 8x ax ax o 0 0 +v +A(x,t) =0 at ax av (x(x, v, t), v (x, v, t) are the characteristics of the equations dt dx _ dv 1 v A(x, t) and shall be referred to as the characteristics.) 28

THE UNIVERSITY OF MICHIGAN - 2764-10-T Proof: x (x, v, t), v(x, v, t) are the initial values of 0, 4 in (C), such that 0, ~ have, respectively, the values x, v at t. If 0 (x, v, ), i(x, v, a) are the unique solutions of do- * (C*) =-A(0, O t-cx) dc having, respectively, the values (x, v) at a =, then x (x, v, t), v (x, v, t) are, respectively, 0 (x, v, t), * (x, v, t). The proof follows, trivially. Lemma 3 Let v (x, v, t) be as in Lemma 1. Then for every pair (v, t), lim v (x,v,t) = v. x-)-o o Proof: From (C+), t 0* (c) x-va + (T-t+ ) P(q)d T?, -CT l (ca)-v I (A(0 (a't), t-c') d'. 0 From the second of these inequalities, t t Iv (x,v,t)-vI I5 LA(0c (),t-cr)l da= IA(0(t-?), T)l d'r, and from the first t 0(t-?)_! x-v(t-') +, (t'-2')P(T")d'. r 29

THE UNIVERSITY OF MICHIGAN 2746-10-T For fixed t, v, T, lim x- -oo 0*(t-T)=-oo, so that, using (A) - (c) lim IA (0 (t-?),T )I = 0. x -- -0o AI (0*(t —'),7)[ $ P(r); But therefore, by the bounded convergence theorem lim Ivo(x,v,t)-v ]- 0. x — -oo Lemma 4 Let y1(o), y2(0r) be the solution of y" -p(t-or)y=0 for which yi(0) = y (0) = 1, Y1 (0) = Y2(0) = 0 Then and ax ax av ' Y2(t), ax yi1(t) ax at — IvI yl(t) + y2(t) P(t) Proof: Let us examine (C*). From the first equation av a av av \ 4 / Since the left member of this equation and a v we have are continuous functions of (x, v), ao\v v av av 30 --

I THE UNIVERSITY OF MICHIGAN 2764-10-T Similarly, 8 80 au eav) A _0 t-) a* ~ —W $ av I Also, when ac = 0, =0, av= 1 -I av Then a av ) au 3v Pao* = - - - 1 %.O av a,i a__0-.A a I + a aA (-(X)auav -1) — a $0 ' -[av 0 Thus - 1 < av 5 p(t-c') {I-av + 0 + 7x do'u I D+cIff S p(t-cx')I- +c' 0 0 A(cx) =50 p(t-cx') a' dcx' dc"t. 00 do' do" +.A(c), where Letting H(a)= 5p(t-0') av 0 0 +a' I dc' dx ", we obtain H(0) = H'(0) = 0, and H"(cr) < p(t-cx)H(ac) - p(t-a)JA(a). 31

THE UNIVERSITY OF MICHIGAN — 2764-10-T Now, let X (a, c') = y2(cr)yl(r')-yl(o)y2('), X (c, r)=O, - (a,o )=l, and da d2X 2 = p(t-o)X>/0 if X>O. Therefore, for cr;', X is an increasing function, and in particular X /- 0. Therefore H(a) S \0 X(a, )p(t-a')A(CT-)dcy' = A(c ') y2(a)dy(lr')-yli()dy2(c')} = -jAo ) +J ' y2(cr)dy(a') -y (a)dy2( o'] = -A ()- f+ y2(). Therefore, + ay2( ), so + t y2(t) t Thus axC < y2(t), and in a similar fashion, | < yl(t) * vI 3x From Lemma 2, al < Ivl yl(t) + y2(t) P(t) Lemma 5 Let Al(x, t)e(A), A2(x, t)e(A), the p(t) being the same for both, and i\ Al(x, t)-A2(x, t)ll $ P, 2(t), a continuous function. Let x (x,v, t),v (x,v,t), ()(x,v, t), v )(x,v, t) be the characteristics corresponding to A and A2, as in Lemma 2. Let and be as in Lemma 4. responding to A1 and A2, as in Lemma 2. Let yi and y2 be as in Lemma 4. 32 -- --

THE UNIVERSITY OF MICHIGAN 2764-10-T If k (t,2 )=y2(t)yl(t-T)-yl(t)y2(t-t), X'(t, )=y2(t)yl(t-'r')-y{l(t)y2(t-T'), then (1) (2)I (5 k(t,') P1 2(r)d x -x 2d 2 0 0 0 |vo -vo I X '(t,?')P 2 (r) d?. O0 Proof: We have from (C ), d (01-02) = - -( 2) dcr( )= [Al(1, t-a)- (0, t-u)] + (01.,t-()-A2 (0 2,t-a)] } The proof then proceeds in a fashion similar to that of Lemma 4. 00 Lemma 6 Let V(v) be non-negative, S V(v)dv=l, monotonic decreasing 0 for v > 0, and monotonic increasing for v O. co If A(x, t)e(A) and v (x, v, t) is the corresponding characteristic, then OD t V(v) (x, v, t) )dv < 1+2 V(O)S P(?)d. -00 0 Proof: Let = P(2')d?. From Lemma 2, v-f3 v v+, so if v /3, 05 v-3 v, andtthus V(v ) t V(v-3). Also if v -,, V(v ) V(v+). 33

THE UNIVERSITY OF MICHIGAN 2764-10-T Therefore, V(5 ) V(v+13)dv+ V(v-3)dv+ V(v )dv -0oo -o = 1+ V(v )dv, 1+2 V(O)3. 0 Lemma 7 Let A(x, t) e(A). Let x (xvt, v, t) v (x, v, t) be the characteristics. Let V(v) be as in Lemma 6, and in addition be continuously differentiable and have a second moment. Let N(x) be continuous, 0- N(x) < N with 0 OD (N -N(x) )dx co. O -00 Then (A) = S N(x (,v,t) )V(vo(r,v,t))dv-N({) d{ -00 -OD __/ exists and ^(A)(A). Proof: From Lemma 6,, N(x (r,v,t) )V(v (,v, t) )dv -c exists for all,. If 13= P(')dr2, o - I 34

THE UNIVERSITY OF MICHIGAN 2764-10-T the integrand is bounded by N V(v+3) for vZ -3, by N V() for -3 vg13, and by N V(v-3) for v >-; hence O 00 lim N( (o,,t) )V(vo( v,t))dv = limr N(x (M,v,t))V(v(,v,tdv -00 0 =S ( N(x (V(Vot))(v))d)(r-S)v S X-oo X -co The last integral has a limit as X —, because of the inteeveryability of N(x)-N 0 exists. This last integral can be written as y; (N(xO)-N )V(VO)dvd + N S (V(v )-V(v) )dv dU +D; [No-N(q)] dU. The last integral has a limit as X-P-Go, because of the integrability of N(x)-N | The second integral can be examined in the following fashion: if P(X, x, t)= ) (V(vo)-V(v) ) d_, I 35

THE UNIVERSITY OF MICHIGAN 2764-10-T then X 00 00 [V(vo)-V(v)]dvd= 5 P(X, x, t)dv, X -00o oo0 P (,x^ 0) -0, a (X, x, t)= -,( v -b ) +A(, t) a( ] d = at i ab v J - X and.C=x - -v V(v _ 0 ~ O IX — -v(v) d s~At t V' ~av This follows from the fact that V has a continuous derivative, and Lemma 2. Thus X 00 X -00 00 L(vo)-V(v)] vtv(vo(x,v, 2)) - -00 -- -V(v(X, v, )] dv d 2 Now since V(v) has a second moment, v V(v) is integrable. For v >3, v V(vO) v V(v-f3) = (v-3)V(v-3)+o3 V(v-3), where 3 = P(T')dZ. A similar 0 result is obtained for v - -3, so that the integrand in this last integral is bounded by an integrable functions, uniformly in X. From Lemma 3, and the bounded Iconvergence theorem, 36

THE UNIVERSITY OF MICHIGAN 2764-10-T x N5 5 {V(vo)-V(v)} dv co 00 Ot oo = N t vf V(vo(x, vSZ)-V(v) dv dr 0 -O0 Now, the integral 5X 5 (N(xo)-No)V(vo)dv d, X -oo by the transformation of coordinates x=Xo:Xb(lt),v,= V t), oVo(,v, t) becomes S(N(xo)-N) V(vo) d o, "5 where, is the region in the (xo, vo) plane bounded by the two curves 0(xo, vo, t)= x, 0(xo, vo, t) = X. By Lemma 1, this region lies between the two straight lines t xo+Vot = x +5 (t-')P(2)dr= x 0 xo+vot = X t-)P(')d = x2 0 Thus X -x 55(N(xo)-N)Vvo)dV(vo)doo-N( V(vo)dvo dx -00 x -X -00 37

THE UNIVERSITY OF MICHIGAN 2764-10-T therefore, [No-N(xo V(vO)dvodxo X -oo remains bounded as X-+ -oo, so x OD Q5 3 (N(xo)-No)V(vo)dvd l -00 -00 exists. This shows thatf(A) exists. It is obvious that (A) satisfies (A)-(c) and (A)-(d). a(A) fN(xo)V(vo)dv-N(x), -0o which is bounded in absolute value by No 1+2 V(O) 5 P(2)di, as we have seen. Hence, (A)-(b) is satisfied. Now, let us examine the integral SV(o)do, in which the region of integration, in the (xo, vo) plane lies between the curve o(xo, vo, t)=x, and the straight line xo+vot=x, with the convention that the integrand is to be taken positive in those simply connected parts of this region which lie above the straight line, and negative in those which lie below. Because of the inequality\ in Lemma 1, the vertical distance between any two points in this region is bounded by t 2 (t-T) P(T)dT, JQ ] m 38 - --

I I THE UNIVERSITY OF MICHIGAN. 2764-10-T so the integral obviously exists. By Green's theorem, this integral can be represented as,co rVo(x,r7, t) I(x, t)=- t) V(~)d~ a x>(x,, t) dr, j x-xo(x, n, t) t I(x, 0) = 0, and oo a rI (vo)a V(X toI aX dra at at a d7 - V V(0)d a ( at -00 X-Xo t Because of Lemmas 3 and 4, the fact that V(v) has a second moment, this last integral can be integrated by parts; we obtain a I - [(v{ ao a + V(-( _ ) dr at at a a t r 7 V t2 a_ j 7 - - r V(vo)dn+f Xv(vo)( r))^ drj, -00 -0O From above, we have already seen that this is the same as a 0x 00 + at [V(Vo,)-V(r7) dcr -00 -00 I 39 --

THE UNIVERSITY OF MICHIGAN 2764-10-T Therefore we have x O0 JfV(vo)dJI __ '[V(v~(x,vt) )-V(v)] dv. A -00 0-sO Now, by rewriting the representation for (A), we obtain X-Xn J(A) = ~ (N(xo)-No)V(vo)dvodxo+5 o-N()] d w -oo c-o 00 + (N(xo)-No)V(vo)dj+ No5 A-O -OD [V(vo)-V(v) dv. By using the preceding equality, (J) where 3(A) = Ao(x, t)+i N(xo)V(vo)d o, Ao(x,t)= N-\vt)V(v)dv-N() d, -0o O which is easily seen to be the same as the sum of the first two integrals in the immediately preceding formula. It is clear that, uniformly in x, J (A)]; [No-N( )] d+ N 1+2 V(O) P( )dj This shows that A) satisfies (A)-, and completes the proof of thelemma This shows that J(A) satisfies (A)-(a), and completes the proof of the\lemma. w 40

THE UNIVERSITY OF MICHIGAN 2764-10-T Corallary: Let Ale(A), and An+l= 6(An), n=1,2,.... Then Ane(A) for all n, and there exists a pji(t), P(t) such that pn(t) p(t), Pn(t) 4 P (t), all n; p(t), P(t) depend only on Al. Proof: The first statement is immediate. From the proof of the Lemma, it is clear that Pn+l, Pn+1 can be chosen to be p S [No-N(1 db+Nfl+2 V(O)j P(})di} It Pn+ = No 1+2 V(O) Pn (t)dt Thus, P(t), ap(t) can be obtained by a Picard type iteration. Lemma 8 Let Al(x, t), A2(x, t), be as in Lemma 5. Then there exists continuous functions Kl(t,?), K2(t,?) dependent upon P, p above, such that t?' IY(A1)- (A2)|| 0 [Kl(t.T) A1-A2 " +.K2(t, ) IAi-A21 d 3 d?. Proof: From the preceding lemma, o(A1)-j(A2)=j (N(xo)-No)V(vo)d 4Jo x + NO ) (v)(, v, t) )-V(vo ),v, t)] dvdwhere A12 is the region between the two curves 0(1) (xo, v t), 0(2) (xo,t)= x, with the same convention in sign as before. By Green's theorem, the first integral 41

THE UNIVERSITY THEUNIERSTYOF MICHIGAN 2764-10-T can be written N (1) - [ 2NX )-Nj (1: (2) -00 a(1) a(2)1 (1) (1) (Vo 1)a7 () V(~)d~ + (2) (2) +50 C (.4 (N(x C (1) X0 -00 -OSw (1) 92) V0 (1) -00 (1) x2o (N()NO (2) -00 V(~)d~} d-q. ~Thus 5(N(x O)V(vO)I dJ0- rvo -v~2Po V1O S(N()N0d 'A O I(1) (2) +N01 x. -x 0 1. ~As for the integral x OD OD F (1) (2)1 -- LV(VO )-V(v0 )j dv~d~, 42

THE UNIVERSITY OF MICHIGAN 2764-10-T from the preceding lemma, this is precisely t Go - v S [v(ol)-V(v() dv t. O -00 But 00 00 ) vV(vo )-V(v 0j dv = v V'(0)(v () ) dv G-oo -o0 where v - P()d' 0-v + P(t?)d'.t v ) #0 But then CV (1))-V(v 2) )dv I v(1) (2)l v IVT(2) IdO )v(-(v ))-V(vo )dv v - v ) and, in a fashion similar to that in which an upper bound on co _) V(v0)dv -GO was found, this last integral is bounbed by a function of? dependent upon P () alone. These results, together with Lemma 5, complete the proof. Theorem: Let N(x), V(v) be as in Lemma 7. Then there exists a unique A(x, t)e(A) for which Yf(A)=A. Proof: The Corollary of Lemma 8, and Lemma 9, show that starting with any function, A1 e (A), the sequence of iterates, A = (A ), converges uniformly on every bounded inerval. Lemma 9 also gives niqueness. every bounded t interval. Lemma 9 also gives uniqueness. l 43

THE UNIVERSITY OF MICHIGAN 2764-10-T We now pass to the less restrictive hypothesis - namely, that N(x), instead of being continuous, be merely measurable; the above techniques of proof fail because the partial derivatives of the characteristics, with respect to x and v, can no longer be assumed to exist. Let N(x) be the limit, therefore, of a sequence of continuous functions Nm(x), in which no generality is lost in assuming that the N (x) satisfy the condition 0, Nm(x) < N(x). Let A(m)(x,t) be the corresponding solution to the equation m(A ) = A(m) Each A()(x, t) e(A), and can be obtained by the iteration process. If the same first interate, say A '(x, t) 0 is used for all m, then the corollary to Lemma 7 sthe _Im) (m) shows that the P (t) and p (t) are all the same, having the common value P(t) and p(t), respectively. The only modification is that in the proof of the corollary, the iteration formula for P must be changed by replacing N(P) by inf N (). n m m Now, using (c), in Lemma 7, A(m x,t)=A(m(x,t)+ Jff / m V(v)dd Am where oo x-xo (m) ( T m (xt)= JL N m (xo)-N V(vo)dvodxo+ 9 N0-N( )J d~, - J-00 -00 o o and A corresponds to the region between the curves x +v t=x, 0(mxo, v t)=x, (m) o o(i) 0 0 having its obvious meaning. How, by the Lebegue theorem, A(m)(x, t) has the 0 limit Aj x, t)= \ [; 0N(S-vt)V(v)dv-N( ) d, and the integrals and the integrals 44 -m - -

THE UNIVERSITY OF MICHIGAN 2764-10-T JJ N(m xo)V(vo)dJo m all have the same bound, as in Lemma 7, determined by P(t). Therefore, a subsequence of the A (x, t), say A )(x, t), converges to some function A(x, t). Obviously A(x, t) is bounded by P(t), and although not necessarily differentiable with respect to x, it satisfies the Lipshitz condition A(xl, t)-A(x2, t) | p(t) I xi-x2 since this condition is satisfied by all the A(m)(x, t). Hence, the functions x j), v (j) arising from the solution of (C*), will converge to functions x, v which are characteristics for A(x,t). But the curve 0(mj)(xo, vo, t)=x has the parameterization x =x( j )(xt) v=v (x, r, t), so that the regions Am O O ' 0 M converge to the region A, corresponding to the curve x = x (x, r, t), v v (x,t). Thus by the Lebesque theorem A(x,t)=A (x,t)+/J N(x )V(v )dJ 0 O 6 0 0 0 But this is equivalent to A(x, t)= ) N(xv (, v, t))V(v (,v,t)dv-N(r), -00 -00 and x, v are the characteristics for A(x, t), we have established the existence of a solution. If Al(x, t), A2(x, t) are two solutions, the estimates of Lemma 5, remain valid. Since Al(x, t) and A2(x, t) are absolutely continuous, if we impose the additional hypotheses that N(x) is almost everywhere continuous, the characteristics x( v, x() v( will be almost everywhere differentiable, and proof 0 0 0 0 of Lemma 8 is easy to modify to give the conclusion of that lemma. Hence, uniqueness is established. We have then the theorem: -_45

THE UNIVERSITY OF MICHIG 2764-10-T Theorem: Let N(x) be measurable, 0 ~ N(x). N, and 00 co f | N -N( I)| d U < oD. Let V(v) > 0, having a continuous derivative, second moment, and co \ V(v)dv = 1. Further, let V(v) be decreasing for Iv| increasing. Then there e lim A(x, t) = 0 x — oo rAN -,xists an A(x, t) such that rl x 00 A(x, t) = N(x0(.* v., t) -OD 00 )V(v (, v, 0 where x (, v,, v v (, v, t) are the characteristics for A. If N(x) is almost everyO O where continuous, A(x, t) is unique. 46

THE UNIVERSITY OF MICHIGAN 2764-10-T APPENDIX B We shall, in this appendix, sketch the proof that if N(x) = N(-x), V(v) = V(-v), the solution A(x,t), of Appendix A, and the distribution function f(x,v, t) = N [ (x, v, t)IV [(x, v, t) have limiting values as t - co, and in part, we can find the functional equations which must be satisfied by these limiting values. The first observation is that, due to the symmetry of the functions N(x), V(v), that x (x,v,t) = -x (-x,-v,t); v (x,v,t) = -v (-x-v, t) so that 0 0 0 0 f(x, v, t) = f)-x, -v, t) and A(-x, t) = -A(x, t). From what has gone before, this shows that x 0o0 A(x,t) = f J ff(M,v,t)dv - N()} d, (B-1) o - 00 and oo (x, t) = - v f (x,v,t) dv. (B-2) at J -00 From (B-l) it is possible to conclude that A(x,t), is bounded, and A(x, t) > o for x > o. The details of the argument are based on the following heuristic idea. If A(x, t) were unbounded, thinking of f (x, v, t) as the distribution function of the electrons, an increasing number of electrons would have to enter the region between o and x, and remain there. On the other hand, electrons entering - -- 47

THE UNIVERSITY OF MICHIGAN 2764-10-T this region are accelerated positively, this acceleration increasing as A(x, t) increases. This decreases the number of electrons which can enter. The only alternative is that A(x, t), being unbounded, is oscillating with large amplitude. But (B-l) shows that A(x, t) is bounded from below, which is a final contradiction. By similar arguments, it can be shown that A(x, t) > o for x > o. But then this implies that x (x,v,t) -e oo for x> o, v < o and x (x,v,t) -> -oo for O O x > o as t -- oo. By considering the equations for the characteristics, it is easy to show that t x [v (xv,t)j2 2= o 0 (X, v, t), V (x,v, t),'t where 0 [x (x,v,t), v (x,v,t),] is the x coordinate, at time, of the particle which arrives at x with velocity v, at time t. But then the remarks above show that as t -- co, x v(xvt) -2 A(, t)d -c0 Then, from (B-2), A 4 -_ vN V - 2 A(, t) d dv at 0 -o -00 and since the integrand is an odd function of v, and since the integrand is an odd function of v, - o as t -oo, which, at 1 __48

THE UNIVERSITY OF MICHIGAN 2764-10-T together with the fact that A is bounded, shows that A (x, t) has a limit, A (x), x and F(x,v,t) has alimit F*(x,v) = NoVfv2-2 J A*(")dJ. Thus we have the functional equation x ~0 A*(x) = 0 N V Iv2-2 J A*( )d] dv - N() d 0 00 -00 which is the content of equation (2-2) and the remarks at the bottom of page 4. 49

THE UNIVERSITY OF MICHIGAN 2764-10-T APPENDIX C Suppose that the function y(x, s) is the Laplace transform of some function of (x, t), and that d - [s2 + W(x)] y =h (x, s), (C-l) where h(x, s) = 0( I si ), 0, W(x), and W(x)-;Wo as x-+t oo. The solutions of the homogeneous equation dx2 - [s+ W(x)] = (C-2) will be a linear combination of the two functions yl(x, s), y2(x, s), which, for both large I sl and large x have the asymptotic forms [s ] V2 e yl(x, s) e. s2 +W(x)] y2(X, S) e 2 +W(] -4 [s2+W(x)J Since any two solutions of (C-1) differ by a linear combination of Y and Y2, and since any non-zero linear combination of y1 and Y2 must increase exponentially in s for some value of x, there is only one solution of (C-1) which is 0 (1 ) for all x, and this must be the solution y(x, s) which we seek, because y(x, s) is, a priori, a Laplace transform. I 50

THE UNIVERSITY OF MICHIGAN 2764-10-T Now if we form the Green's function: 1 Yl(x)y2(x') x< ' G(x, x')= (C-4) W(Y1, Y2) Yl((x')y2(x) x >X and consider S G(x,x')h(x',s)dx', (C-5) this is the solution of (C-1) which vanishes as x — + o. But from the general theory \2 of such operators, it is known that the L norm of such operators is 0 (i) and since h(x', s) = O( Is I), the solution given by (C-4) is O (l-), at least when h(x, s) is in L2(-oo, oo). But by approximating h(x, s) by functions in L2(-oo, oo ), we see that the appropriate form for the solution of y(x, s) which is O (, or equivalently, which is a Laplace transform, is given by (C-4), so that the poles of y(x, s) are given as the zeros of W(yl, y2). I I m 51

THE UNIVERSITY OF MICHIGAN --- m 2764-10-T REFERENCES Dolph, C. L., and Weil, H., 1959, The University of Michigan Radiation Laboratory Report 2778-2-F, Studies in Radar Cross Sections XXXVII. Iordaneskii,S. V, 1959, A Solution of the Cauchy Problem for the Kinetic Equation of Electron Plasma, Dok. Akad. Nauk. USSR, Vol. 127, 509-512. Sawchuk, W., 1962, The University of Michigan Radiation Laboratory Report 2764-9-T. - -- 52