![\"plasma-oscillations\"](\"http:\/\/underthehood.blogwyrm.com\/wp-content\/uploads\/2016\/09\/Plasma-Oscillations.png\")
<\/a><\/p>\nHow does the plasma respond to this out of equilibrium condition? The first approximation one can make is that the electrons are free to move but the ions are fixed in space.\u00a0 This corresponds to the physical approximation that the ions have an infinite mass \u2013 an approximation that is well supported by the ~2000 \u2013 1 ratio between the mass of a proton and that of an electron.<\/p>\n
The electrons will feel a strong restoring force (red arrow) and they will move accordingly.\u00a0 Since the ions are fixed, they will feel an equal but opposite restoring force (blue arrow) but they will not be able to respond.<\/p>\n
The electric field that the electrons and ions will feel, which is most easily obtained from Gauss\u2019s law, is non-zero between the strips of positive and negative charge in an exactly analogous fashion as to the textbook problem on parallel plate capacitors.\u00a0 The density of ions in the red region is $$n_0$$ and the volume is $$A \\Delta x$$.\u00a0\u00a0 Since $$\\Delta x$$ is small compared to the other dimensions, the electric field is closely approximated as being entirely in the $$x$$-direction and the total flux through a Gaussian surface whose boundaries coincide with the ion region is due strictly to the left and right faces.\u00a0 The right face flux is exactly zero since the flux from the electrons exactly cancels the flux from the ions.\u00a0 The left face flux, which is equal to the total flux, is then<\/p>\n
\\[ \\Phi_E = E A \\; , \\]<\/p>\n
which is proportional to the total charge enclosed<\/p>\n
\\[ \\Phi_E = \\frac{e n_0 A \\Delta x}{\\epsilon_0} \\; .\\]<\/p>\n
Equating gives<\/p>\n
\\[ E= \\frac{e n_0 }{\\epsilon_0 } \\Delta x \\; ,\\]<\/p>\n
which is the usual result for the field between a parallel plate capacitor ($$n_0 \\Delta x$$ is the surface charge density).<\/p>\n
As the electrons were all displaced uniformly and are experiencing a uniform force, they will move together as a unit, that is to say cooperatively.\u00a0 The mass of the unit is simply the mass density $$m_e n_0$$ times the volume consumed $$V = A \\ell_0$$.\u00a0 The center of mass of the slab is at $$\\ell\/2- \\Delta x$$ with the corresponding acceleration $$\\Delta {\\ddot x}$$.\u00a0 The force on the slab is<\/p>\n
\\[ F = q_{tot} E = -e n_0 V \\frac{e n_0}{ \\epsilon_0} \\Delta x\\]<\/p>\n
and the equation of motion of the entire slab is given by:<\/p>\n
\\[ m_e n_0 V \\frac{d^2}{dt^2} \\Delta x = -e n_0 V \\frac{e n_0} {\\epsilon_0} \\Delta x \\; .\\]<\/p>\n
Dividing out the common factors gives<\/p>\n
\\[ \\frac{d^2}{dt^2} \\Delta x = -\\frac{e^2 n_0}{m_e \\epsilon_0} \\Delta x \\; .\\]<\/p>\n
From this equation we would conclude that, for small displacements, the electron slab would oscillate at the angular frequency<\/p>\n
\\[ \\omega_{pe} = \\sqrt{ \\frac{e^2 n_0}{m_e \\epsilon_0} } \\; , \\]<\/p>\n
which is called the electron plasma frequency (note it is really an angular frequency).<\/p>\n
Relaxing the assumption that the ions don\u2019t move makes the analysis somewhat more complicated but only changes the result in a small but crucial way.\u00a0 The location of both electron and ion centers-of-mass now matter and must be tracked separately.\u00a0 The picture of their motion is now<\/p>\n
Since both species are free to move, the thickness of the exposed electron or ion regions (red and blue) is not determined solely by the movement of their respective centers-of-mass (the ion region pictured above is much larger than $$\\Delta x_i$$) .\u00a0 The total thickness of the exposed regions depends on the relative motion between both centers-of-mass<\/p>\n \\[ t = \\Delta x_i – \\Delta x_e \\; .\\]<\/p>\n An initial analysis may suggest that there are 8 separate cases: positive or negative variations for each species (2×2) as well as whether the ions are more displaced versus the electrons.\u00a0 In fact, there are only two cases based on whether $$t$$ is positive (ion strip is on the right) or whether $$t$$ is negative (ion strip is to the left).\u00a0 The following table shows the 8 variations for $$\\ell = 1$$.<\/p>\n The parameters with the \u2018min\u2019 designate the leftmost extent of the electron slab ($$e_{min}$$) and the ion slab ($$i_{min}$$).\u00a0 The parameter $$t_{low} = i_{min} – e_{min}$$.\u00a0 A negative value, such as in case 1, means that the ion slab has shifted more to the left than the electrons and that the exposed ion region (red) is now on the left while the exposed electron region (blue) is now on the right.\u00a0 A positive value of $$t_{low}$$ means the converse.\u00a0 The corresponding parameters at the rightmost side are defined in an analogous fashion.\u00a0 The color coding assists in visualizing which species is in which location.\u00a0 Note that the entire table can be summarized by $$\\Delta x_i – \\Delta x_e$$, whose magnitude and sign match $$t_{low}$$ and $$t_{high}$$.\u00a0 Also note that the displacements in the table are purely geometric and that the actual dynamic displacements would not shift the total center-of-mass since the only forces in the problem are internal.<\/p>\n The electric field in between the two uncovered charge regions can also be obtained by Gauss\u2019s law, but it is instructive to use the expression for the electric field due to a sheet of charge that and superposition. The electric field due to the ions will have a magnitude of<\/p>\n \\[ E_i = \\frac{e n_0}{2 \\epsilon_0} \\left(\\Delta x_i – \\Delta x_e \\right) \\]<\/p>\n using the standard result from elementary E&M. <\/p>\n Likewise the electric field due to the electrons will also have the same magnitude $$E_e = E_i$$. Their directions will be the same in the region between them (the neutral plasma) and opposite outside, so that a non-zero field with magnitude $$2 E_i$$ will result only in the middle region. The direction of the field depends on which one is on the right versus the left but it doesn\u2019t matter for the sake of the dynamical analysis. \\[ m_i \\frac{d^2}{dt^2} \\Delta x_i = -\\frac{e^2 n_0}{\\epsilon_0} \\left(\\Delta x_i – \\Delta x_e\\right) \\]<\/p>\n and<\/p>\n \\[ m_e \\frac{d^2}{dt^2} \\Delta x_e = \\frac{e^2 n_0}{\\epsilon_0} \\left(\\Delta x_i – \\Delta x_e\\right) \\; . \\] <\/p>\n Combining these two equations gives<\/p>\n \\[ \\frac{d^2}{dt^2} \\left( \\Delta x_i – \\Delta x_e \\right) = -\\frac{e^2 n_0}{\\epsilon_0}\\left(\\frac{1}{m_e} + \\frac{1}{m_i} \\right) \\left(\\Delta x_i – \\Delta x_e\\right) \\]<\/p>\n from which we conclude that the frequency is now<\/p>\n \\[ \\omega_p^2 = \\frac{e^2 n_0}{\\epsilon_0}\\left(\\frac{1}{m_e} + \\frac{1}{m_i} \\right) = \\frac{e^2 n_0}{\\epsilon_0 m_e} + \\frac{e^2 n_0}{\\epsilon_0 m_i} \\equiv \\omega_{pe}^2 + \\omega_{pi}^2 \\; .\\] <\/p>\n The above analysis has a direct analog in a mechanical system of two unequal masses $$m_1$$ and $$M_2$$ connected by a spring of spring constant $$k$$ and equilibrium length $$\\ell_0$$.<\/p>\n Defining the deviations of each mass from the equilibrium positions as<\/p>\n \\[ \\Delta x_1 = x_1 – x_1^{(0)} \\]<\/p>\n and<\/p>\n \\[ \\Delta x_2 = x_2 – x_2^{(0)} \\] <\/p>\n the distance between them takes the form<\/p>\n \\[ x_2 – x_1 = x_2^{(0)} + \\Delta x_2 – x_1^{(0)} – \\Delta x_1 \\equiv \\ell \\; .\\]<\/p>\n Combining these two definitions gives<\/p>\n \\[ x_2 – x_1 = \\Delta x_2 – \\Delta x_1 + \\ell_0 \\equiv \\delta \\ell + \\ell_0 \\]<\/p>\n from which the potential follows as<\/p>\n \\[ V = \\frac{1}{2} k (\\ell – \\ell_0)^2 = \\frac{1}{2} k (\\Delta x_2 – \\Delta x_1)^2 \\; .\\]<\/p>\n The Lagrangian takes the form <\/p>\n \\[ L = \\frac{1}{2} m_1 \\Delta \\dot x_1^2 + \\frac{1}{2} M_2 \\Delta \\dot x_2^2 – \\frac{1}{2} k (\\Delta x_2 – \\Delta x_1)^2 \\; , \\]<\/p>\n with the equations of motion becoming<\/p>\n \\[ \\frac{d}{dt} \\frac{\\partial L}{\\partial \\dot x_1} – \\frac{\\partial L}{\\partial x_1} = m_1 \\Delta \\ddot x_1 + k(\\Delta x_2 – \\Delta x_1 ) = 0 \\]<\/p>\n and<\/p>\n \\[ \\frac{d}{dt} \\frac{\\partial L}{\\partial \\dot x_2} – \\frac{\\partial L}{\\partial x_2} = M_2 \\Delta \\ddot x_2 – k(\\Delta x_2 – \\Delta x_1 ) = 0 \\]<\/p>\n These equations combine nicely into one expression for the dynamics of the change in the separation between the two masses<\/p>\n \\[ \\delta \\ddot \\ell = -\\left( \\frac{k}{m_1} + \\frac{k}{M_2} \\right) \\delta \\ell \\; , \\]<\/p>\n which immediately tells us that the frequency of the oscillations is given by <\/p>\n \\[ \\omega^2 = k \\left(\\frac{1}{m_1} + \\frac{1}{M_2} \\right) = \\omega_1^2 + \\omega_2^2 \\; .\\]<\/p>\n In the limit as $$ M_2 \\rightarrow \\infty $$, the equation of motion simplifies to <\/p>\n \\[ \\delta \\ddot \\ell + \\frac{k}{m_1} \\delta \\ell = 0 \\; ,\\]<\/p>\n with frequency of <\/p>\n \\[ \\omega^2 = \\frac{k}{m_1} = \\omega_1^2 \\; .\\]<\/p>\n Thus there is a perfect analogy between the two pictures and the cooperative motion of the electrons and ions within the plasma can be interpreted in strictly mechanical terms.<\/p>\n","protected":false},"excerpt":{"rendered":" This month\u2019s installment focuses on that aspect of plasma dynamics that is strongly analogous with mechanical vibration: plasma oscillations.\u00a0 Plasma oscillations were first predicted by Paul Langmuir in the 1920s… Read more ><\/a><\/p>\n","protected":false},"author":2,"featured_media":0,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[1],"tags":[],"class_list":["post-886","post","type-post","status-publish","format-standard","hentry","category-uncategorized"],"_links":{"self":[{"href":"https:\/\/underthehood.blogwyrm.com\/index.php?rest_route=\/wp\/v2\/posts\/886","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/underthehood.blogwyrm.com\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/underthehood.blogwyrm.com\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/underthehood.blogwyrm.com\/index.php?rest_route=\/wp\/v2\/users\/2"}],"replies":[{"embeddable":true,"href":"https:\/\/underthehood.blogwyrm.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=886"}],"version-history":[{"count":9,"href":"https:\/\/underthehood.blogwyrm.com\/index.php?rest_route=\/wp\/v2\/posts\/886\/revisions"}],"predecessor-version":[{"id":895,"href":"https:\/\/underthehood.blogwyrm.com\/index.php?rest_route=\/wp\/v2\/posts\/886\/revisions\/895"}],"wp:attachment":[{"href":"https:\/\/underthehood.blogwyrm.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=886"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/underthehood.blogwyrm.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=886"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/underthehood.blogwyrm.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=886"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}<\/a><\/p>\n<\/a><\/p>\n
\nThe equation of motion for each slab is then given by<\/p>\n<\/a><\/p>\n