In a previous article we established an approximate estimate for the frequency of oscillations in a “cold plasma” because of small disturbances in the arrangement of the electrons and ions. In this article we will establish a much better model for oscillations in a plasma, which takes into consideration the fact that Plasmas tend to be hot (the thermal factor). To do this we establish a differential equation using facts we can take as a given. Firstly note that we model our plasma as a fixed lattice of ions, with moving electrons. These moving electrons cause the plasma oscillations we are trying to model. The frequency at which a hot plasma oscillates is called the * Langmuir frequency*, and the method detailed in this article attempts to find an estimate of this frequency.

We start by looking at the conservation of charge equation:

Where is the number density of particles in the plasma, and is the charge on the electron. We notice that the charge of the electron is clearly a constant, so it can be factored out of the derivative with respect to , and we can divide through by this electron charge to get the following:

It is clear that we have the derivative of a product with respect to . We can therefore expand this using the product rule:

We next assume that the variation in velocities and number density is small. Doing so we can state the following:

and

Where is the constant number density when the plasma is still, and is the velocity of the electrons when the plasma is still (0). The terms and represent small changes in their respective quantities. We now have a general expression for both the number density and the velocity of the electrons in terms of a constant *“still state”* and small a small variation in each. We can therefore substitute into the conservation of charge equation for both and (remembering that ):

We can now neglect the derivatives in the form and , as these are second order terms and will be negligible due to the assumptions underlying this approximation of the Langmuir frequency:

We now move on to the equation of motion of the fluid. To derive this for our plasma we take the general equation of motion for a small volume of fluid:

This is simply derived from Newton’s second law, , where is the mass density, thus the term corresponds to , per unit volume. This means that the right hand side corresponds to the force on the fluid, and it does. The term corresponds to the force on a small volume per unit volume. The term (which means the vector gradient of the pressure on the fluid) is an extra force the fluid feels due to changes in pressure at the edges of the small volume . We can simplify the equation of motion like so:

We now only work in the direction so we can drop vector notation. To get to this equation (above), I substituted for the mass density, in terms of the number density and the mass, simply called the term for simplicity, and replaced pressure with as this is the pressure in a gas. We can therefore take the vector grad (remembering we are only working in the direction) to get .

We can now make the same substitutions as before (with number density and velocity in terms of small changes in such quantities) to get the following (remembering that derivative of constants is 0):

We can again neglect second order terms, like before leaving us with:

We now have the following two equations, with some similar terms, we can therefore look to make a substitution:

and

If we take the temporal derivative of the modified conservation of charge equation , an interesting substitution becomes possible:

The term now appears in both the equations we previously derived, allowing us to make a substitution. We first get make the subject of the motion equation:

We can now substitute for in the other equation:

We now have a differential equation in that we can solve. To do so we have to consider what the term represents once again. We defined to be equal to where is the external force on the charged particles in our plasma. When considering this external force we can use simple electric fields to determine its magnitude as with being the charge on the electron. Since thus . Since , and we can simplify the differential equation. Also because change in force per particle is equal to the external force multiplied by the small change in electron number density. This is because the change in electron number density is what creates the electric field, thus force. As the plasma is modelled as a static lattice of ions with moving electrons. So:

The derivative of the electric field strength with respect to is given by the following- (where is the charge density of the particles in the plasma), thus we can replace this in our equation:

We know that charge density () is the charge per unit volume thus (where is the electron number density) allowing us to find an interesting result:

The term is the square of the cold plasma oscillation frequency we worked out here () we can therefore substitute this in, and once finished we will have an approximation of the Langmuir frequency, in terms of the cold plasma frequency, allowing us to see how the thermal energy of the plasma affect the oscillating frequency:

We now need to solve the differential equation, which ends up being very simple. We just have to assume that the changes in electron number density can be modelled as a plane wave, (where is the wave vector). We also recognise that the other changes in number densities in the equation also must correspond to changes in electron number densities, as explained above- it’s the changes in electron number density that matter. We can therefore find the required partial derivatives of the plane wave to substitute in for the and terms, as these are simple to work out:

Then we can substitute this all back into our differential equation:

The plane wave term() is now common to all the terms so we can cancel it:

We now have a function for the Langmuir frequency we are trying to derive, but we can make a few further simplifications. The kinetic energy of the electrons in a hot plasma is equal to where is the thermal velocity of the electrons (velocity due to thermal energy). This means that if we substitute the kinetic energy in terms of the thermal energy, the mass of the electron cancels. Furthermore , which was the number density of the electrons at an equilibrium point, is equal to where is the atomic number of the ions in the plasma, and is their number density. This can be said because at equilibrium the number of ions is equal to the number of electrons in the plasma. Thus:

Therefore an approximation to the Langmuir frequency of oscillation of a plasma can be given by the following:

***Disclaimer- This represents my own work as part of a physics challenge. I cannot be sure of the validity of each step of my derivation, so I encourage anybody reading this to give me any tips if something isn’t quite correct. Thank you.*

## Comments 1