About this Journal Submit a Manuscript Table of Contents
International Journal of Optics
Volume 2012 (2012), Article ID 648741, 8 pages
Research Article

Quantum Character of Electromagnetic Langmuir Oscillations in Conventional Electron-Ion Plasma

Institute for High Energy of Joint Institute for High Temperature, Russian Academy of Science, Ijorskaja 13/19, Moscow 127412, Russia

Received 22 January 2012; Revised 16 April 2012; Accepted 17 April 2012

Academic Editor: Adam Miranowicz

Copyright © 2012 Boris Alexandrovich Veklenko. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


It is shown that the low-temperature plasma near-thermodynamic equilibrium cannot be classical because of a quantum nature of the longitudinal electromagnetic field and electron interaction with Rayleigh-Jeans distribution of Langmuir waves. The theory requires introduction of a dimensionless quantum charge whose value is greater than unity leading to a liquid-like behavior of the plasma.

1. Introduction

Plasma is called classical and the low-temperature one, if electrons and ions in it are governed by equations of classical electrodynamics and obey Maxwell distribution at thermodynamic equilibrium. The latter condition implies that the de-Brogle temperature wave length is much less than the mean distance between the particles in plasma. Here, is the concentration of particles. The inequality written above does not mean that the electromagnetic field in plasma could be described also by classical physics methods. The latter is widely used as a traditional approach [15]. Describing the electromagnetic field with the help of the strength , we suppose that it is in a quantum coherent state [6] and the occupation numbers are much greater than unity. Although the latter is meaningful, the former has no grounds. Similarly, one could suppose that the field in plasma is in Fock quantum state with the mean value . This is more preferable since there are the Fock states, which are generated by a thermal source. In any case, a classical electromagnetic field description of the plasma should be justified. The procedure of field quantization in plasma leads to the appearance of the quantum of Langmuir oscillations energy with frequency . We use the rationalized Gauss unit system. Here, is electron charge, is its mass, and is Planck constant. The quanta do not play any role only if the classical energy of two-particle interaction at the distance of Debye radius , where is the temperature, satisfying the inequality:

After substitution the explicit expressions for and , we get where is the mean electron velocity in plasma. For , it occurs that and . At first glance, one could think that the inequality (2) is satisfied and the quantization procedure could be omitted. But, surprisingly, a dimensionless charge (2) appears in the theory, with the value greater than one. The analogous dimensionless charge is known in quantum electrodynamics which uses transversal electromagnetic waves as fine structure constant or Sommerfeld constant. Its small value permits there to use a perturbation technique. Basically, investigating the plasma, we deal with the same electrodynamics but dealing with the longitudinal waves. In this theory, the parameter denoting the vacuum light velocity is absent. It is replaced now by . That is why the dimensionless charge (2) occurs being great and there appears a question about the applicability of the perturbation technique, which is the base of modern theory of plasma. We state that, in our case, the interaction constant is big and plasma has liquid-like but not a gaze structure. Moreover, the value of dimensionless charge (2) arises at and the limiting transition to the conventional classical theory of plasma is not evident any more.

To resolve the problem, we develop a quantum field theory of plasma and investigate the conditions of the applicability of the classical theory.

2. Principle Equations

Let us consider the plasma in thermodynamic equilibrium in which the ions are uniformly distributed over the space forming a certain background, which, along with electrons, is neutral. We use the Heisenberg representation to describe the quantum electron field , where is the spatial coordinate and is an arbitrary moment of time. Spin effects are omitted. Since the plasma electrons obey Maxwell distribution by definition, the results do not depend on the kind of statistics. Thus, we can use Bose-Einstein commutation relations:

The field operator is described by the Schrodinger equation:

The electromagnetic field operator obeys the Maxwell equation: and the conventional commutation relation:

We assume Einstein's summation rule over repeating indices . The current density is defined by expression:

If the classical current exists in the system, it is necessary to add it to . The system of (4), (5), (7) is written in gauge condition with a zero scalar potential. It permits us to write the quantum Maxwell equations in a simple form (5)-(6). We are interested in the longitudinal electromagnetic field which is the potential field. That is why instead of zero scalar potential we introduce the pseudo-potential by definition: or after Fourier transformation:

Further, it is convenient to exclude the vector potential from (4) in order to deal with the linear term . Such a problem can be solved by using the auxiliary operator:

The integration is performed over any curve connecting points and , where is any arbitrary point. By changing the variable: one can write (4) in the form: and write instead of (7).

We omit tildes below. If is a classical function, we come to the conventional description of the classical plasma (see [7]). Below, we investigate how quantum longitudinal plasma waves affect the electrons. In a classical plasma, there are undamped Langmuir oscillations [1, 4] at the wave numbers . It means in connection with commutations relations (6) that

Here, and are the annihilation and the creation Langmuir’s wave operators with vector and polarization index . For the transversal waves , for the longitudinal one, , is the Heaviside step function. The longitudinal part of vector potential satisfies the equation:

Now, from (9), one gets

The function deals with transversal waves fall out of calculations. We shall use below some auxiliary constructions. strait calculation yields

The averaging in formulae (17) performs both in quantum and statistical senses over Gibbs distribution, and if , then (Rayleigh-Jeans approximation)

3. Green Functions and

For investigation of the equation system (5)–(16) it is convenient to use the method of Green functions. Let by definition

Using (12), we get

Using (12) again, we have for the right term in (20)

In order to get the closed equation system, we use the following conventional procedure of breaking up the correlators in (21):

Such breaking is certainly correct if the fields do not interact or the interaction is sufficiently weak. In practice, the region of its validity is much greater. Its accuracy is greater than the accuracy of a conventional Hartree-Fock approximation when one supposes the correlator braking in the first (20) equation of system. But now, one breaks the correlators in the second equation (21) of the system. In some cases by such a way, one can develop a theory of phase transitions [8]. If , one can suppose that . Now, (21)-(22) could be easily solved using Fourier transformation:

Finally, from (20), one gets

In accordance with (16), we took into account that

The function is analytical in upper half plane of complex by definition due to Heaviside step function . That is why the term appears in (24).

In accordance with (24), the interaction of plasma electrons with thermal quantum electromagnetic fields is described by a pole leading to qualitative deformation of their properties. Formula (24), in accordance with (18), can be rewritten as where and is the Bohr radius of Hydrogen atom. So, the procedure of quantization of the longitudinal electromagnetic field leads to the emergence of the dimensionless charge introduced in the introduction qualitatively. Formula (26) can be rewritten in a form:

The limit corresponds to free electrons without field . If , the free electrons are replaced by two types of quasi-particles described by the two last terms of (28). It is the first result of present paper. Namely, these quasi particles describe electrodynamics, thermodynamic, and especially kinetic properties of plasma at . Since such inequality is always satisfied for low-temperature plasma, the quantum parameter should be present in all formulae. Nevertheless, one should stress that the used formalism and, in particular, the procedure of correlator breaking does not guarantee the accuracy of quantitative results at . For this reason, the extrapolation of the results to the region is expected to lead as in theory of liquids to semiquantitative predictions. But at the same time, this way proves that theories which use limit (Vlasov equation) are not correct.

It follows from formula (28) that at there exists any “camouflage effect” that sometimes hides the Plank's constant . For instance, in such region and Plank's constant disappeared since using explicit form (27) for , we get

That is why, in this region, the erroneous picture of electron clouds as an ensemble of free electrons sometimes leads to the pseudo reasonable results following, in particular, from Vlasov equation [2].

4. Chemical Potential of Electron Gas

On can find the correlation function of electron gas with the help of Green functions and fluctuation-dissipation theorem [9]:

For nondegenerate gas,

One gets from (28) that

Since one has the equation for chemical potential which is easily solved where is the chemical potential of ideal electron gas. Formula (36) possesses a remarkable property. The difference in both regions as and according to equality , that is, does not depend on Planck constant. So only for this difference depends on . We remark that in plasma always . For typical values of parameters, and , it turns out that , , . If the concentration diminishes then increases but the ratio diminishes. So at higher concentration the procedure of electromagnetic field quantization according to (36) leads to the violation of the Sacha formula: where is the ionization potential. The corrector term unclosed in brackets depends on only if is of the order of unity.

5. Plasma Dielectric Permittivity for Longitudinal Waves

Let the classical current   generating the classical potential exists in plasma. Now, (12) takes the view:

This equation can be rewritten in integral form: where is the solution of (40) at . We use approximation:

The electron density at first approximation over the potential is where

According to (33) and (34), one gets

Now, the increment of the electron concentration in plasma as a function of due to their interaction with field can be written as where

The plasma dielectric permittivity can be found from (5) and equation of discontinuity:

On the other hand, from the mean value equality , it follows that . After Fourier transformation of (5), multiplying the result by and taking into account the presence of classical current and longitudinal waves, one gets

After comparing this expression with the definition of dielectric permittivity : one gets


H.c. means the complex conjugate values and simultaneously changing and .

The explicit form of the dielectric permittivity one can obtain from (52) and (28) if simultaneously one takes into account that for not too large with respect to the inequalities and take place:

From this formula, it follows that there are three resonance regions in plasma: —the Vlasov-Landau region at small frequencies and the new resonance regions in vicinities of quantum frequencies and . The existence of these regions manifests themselves in a characteristic way in the imaginary part of the dielectric permittivity. From (53), one gets :

From this equality at , as expected, one finds the Landau attenuation. As the parameter goes up, the plasma starts to absorb at frequencies . Then, this effect disappears. And if ,

In other words, the plasma absorbs at not only at low frequencies, but at high quantum frequencies as well. At and low frequencies, is a half of the imagine part of permittivity predicted by Landau theory [3] at . The intensity of absorption at the frequencies is, respectively, larger than intensity of low frequencies absorption.

Since where is the conductivity of substance one states that the parameter plays an important role in kinetic properties of plasma.

As for the real part of dielectric permittivity at , (52) shows the conventional results, and in particular the Vlasov dispersion formula [2]:

We remind that really Z>>1. If grows, but still it is small, the conventional results of Vlasov theory change considerably. Let for the beginning . The explicit form of equation Re describing oscillation spectrum of the system in linear approximation over and by has the form: where and slightly depend on the frequency. In the present conditions, it occurs that . That is why can be taken into account by the perturbation technique. The only one solution of (58) satisfies inequalities (59). As a result, we have

So, taking into account inequality at small , we state that the Vlasov spectrum (57) changes significantly.

The more important fact is the depending of the result through on Planck constant . The larger is , the larger is the region of validity of the formula (60).

Let now . The equation describing the oscillation spectrum has the form:

In this equation, the last term can be taken into account using perturbation technique. It yields

With respect to (60), it is another brunch of spectrum. The symbol means the principal value of the integral. The expression (62) is written in the form, which allows comparison with Vlasov formula (57). It is evident that in the region under consideration the Vlasov formula is correct but the presence of parameter means that all the oscillation modes through imagine part are damped. Moreover, the dispersion condition (57) can hardly be implemented in practice. In real electron-ion low-temperature plasma condition is always .

Now, we consider the more important case when . The equation describing the oscillation spectrum reads where

Also, (63) describes two spectrum branches: as a solution of biquadrate equation. In conditions under consideration, the following inequality is correct . The only one result coinciding with inequalities (64) follows:

This formula differs significantly from Vlasov’s formula (57). We state that for the real conditions the presence (63) of the second term dramatically changes the Vlasov result. This branch of spectrum depends principally on since .

6. Debye Screen Radius

Equation (50) can be cast in the form: where is the density of external classical charge present in plasma and

If the point charge is placed at point and does not depend on time, then . Using (68) at , one gets the conventional result:

According to (53) at arbitrary and ,

The term multiplied by describes the Debye radius renormalization. The last term takes part in radius and charge renormalization. If , then using (36): and as a consequence, we return for to the formula (69) in which now instead of the value is placed where

This formula does not depend on . Such a fact once again manifests the “camouflage effect” explaining in some cases the pseudoreasonable results following from conventional classical theory of plasma.

7. Conclusion

The electromagnetic field is always quantized. The question is whether the quantum properties of the field are significant or not in a particular system.

In the present work, we show that as quantization procedure of longitudinal (Langmuir) electromagnetic field and taking into account the electron interaction with Rayleigh-Jeans distribution of Langmuir waves at finite temperature change dramatically the structure of the electron cloud in a low-temperature plasma. The dimensionless quantum parameter is proportional to the charge squared whose value is much larger than the one appears in theory. It makes impossible to use a perturbation technique. We tried to construct the closed system of equations, the solution of which using perturbation technique (small ) returns to the consequences of Vlasov equation. In real conditions of the low-temperature electron-ion plasma, that is why, in any case, the plasma optical properties could not be investigated using Vlasov’s equation. The analysis of the obtained system of equations at large leads to the conclusion that the quantum properties of plasmas manifest themselves dramatically in kinetic and hydrodynamic properties.


  1. L. Tonks and I. Langmuir, “Oscillations in ionized gases,” Physical Review, vol. 33, no. 2, pp. 195–210, 1929. View at Publisher · View at Google Scholar · View at Scopus
  2. A. A. Vlasov, “On vibration properties of electron gas,” Journal of Experimental and Theoretical Physics, vol. 8, no. 3, p. 291, 1938 (Russian).
  3. L. D. Landau, “O kolebaniiakh elektronnoi plazmy,” Journal of Experimental and Theoretical Physics, vol. 16, p. 574, 1946 (Russian).
  4. A. F. Alexandrov, L. S. Bogdankevich, and A. A. Rukhadze, Principles of Plasma Electrodynamics, Springer, Berlin, Germany, 1984.
  5. P. K. Shukla and B. Eliasson, “Nonlinear aspects of quantum plasma physics,” Uspekhi Fizicheskikh Nauk, vol. 180, pp. 55–82, 2010.
  6. R. Glauber, in Quantum Optics and Electronics, C. DeWitt, A. Blandin, and C. Cohen-Tannoudji, Eds., Universite deGrenoble. Ecole d’ete Physique Theorique 1964.Gordon and Breach, Science Publishers, New York, NY, USA, 1965.
  7. M. V. Kuselev and A. A. Rukhadze, “On the quantum description of the linear kinetics of a collisionless plasma,” Physics-Uspekhi, vol. 42, no. 6, p. 603, 1999.
  8. V. L. Bonch-Bruevich and S. V. Tjablikov, The Method of Green Functions in Statistical Mechanics, URSS, Moscow, Russian, 1961.
  9. H. B. Callen and T. A. Welton, “Irreversibility and generalized noise,” Physical Review, vol. 83, no. 1, pp. 34–40, 1951. View at Publisher · View at Google Scholar · View at Scopus