arxiv: v1 [physics.plasm-ph] 9 Dec 2018

Transcription

1 Quantized Plasmon Excitations of Electron Gas in Potential Well M. Akbari-Moghanjoughi 1 1 Faculty of Sciences, Department of Physics, arxiv: v1 [physics.plasm-ph] 9 Dec 2018 Azarbaijan Shahid Madani University, Tabriz, Iran (Dated: December 18, 2018) Abstract Using the Schrödinger-Poisson system in this paper the basic quantum features of plasmon excitations in a free noninteracting electron gas with arbitrary degeneracy is investigated. The standing wave solution of the free electron gas is derived from the corresponding linearized pseudoforce system with appropriate boundary conditions. It is shown that the plasmon excitation energies for electron gas confined in an infinite potential well are quantized eigenvalues of which are obtained. It is found that any arbitrary degenerate quantum electron gas possesses two different characteristic length scales, unlike the classical dilute electron gas, with the smaller length scale corresponding to the single particle oscillation and the larger one due to the collective Langmuir excitations. The probability density of the free electron gas in a box contains fine structures which are modulated over a larger pattern. The envelope probability density profile for the electron Fermi gas confined in an impenetrable well in different energy states are found to be quite similar to that of the free electron confined to an infinite potential well. However, the illustrative features of the plasmon theory presented in this research can be further elaborated in order to illuminate a wide range of interesting physical phenomenon involving both single particle as well as the collective features. PACS numbers: q,71.10.Ca, d 1

2 I. HISTORICAL BACKGROUND Quantum behavior of atomic scale phenomena has fundamentally altered all our intuition of physical world, since the pioneering discoveries of Max Planck [1], Erwin Schrödinger [2], Louis de Broglie [3] and many others. Almost all of the modern scientific developments in the fields like solid state [4, 5], optoelectronics [6] and many technological achievements concerning miniaturized semiconductor integrated circuitry [7, 8], nano- and microelectronic device fabrication [9], plasmonics [10], etc. are indebted to our extensive knowledge of the quantum theory. Quantum effects appear when the characteristic length of the system compares to the thermal de Broglie wavelength of the particle leading to effective quantum interactions. Such interactions may be either due to the confinement of the particle in a potential well of length l comparable to this wavelength or due to increase in the particle number density to the extreme limit where the interparticle distances compares to this quantum length so that the single particle wave functions in the system start to overlap [11]. However, there is fundamental difference between the later two cases as follows. The single particle Schrödinger wave function possesses only one characteristic length determined by the width of the interacting potential well, whereas, the many particle wave function has two characteristic lengths with the smaller one corresponding to the single particle diffraction effect and the larger one due to the collective interactions. In plasmas the larger length scale is closely related to the well known Thomas-Fermi length. While the celebrated Schrödinger equation satisfactorily deals with many natural phenomena regarding the noninteracting Bosonic particle systems, it may not be well suited for dense fermionic systems like quantum plasmas in which collective quantum electrostatic interactions are natural. The field of quantum plasmas itself has a long development tradition. A big list of pioneering contributions to this field includes those of Fermi [12], Madelung [13], Hoyle and Fowler [14], Chandrasekhar [15], Bohm [16], Pines [17], Levine [18], Klimontovich and Silin [19] and many others which have led to discovery of many interesting collective quantum properties of dense and degenerate ionized environments. Quantum plasmas are relevant not only to metals, semiconductors and nanostructures but also are quite relevant to the astrophysical environments like big planetary cores, stellar media, white dwarfs and neutron star crusts. There are many recent literature which explore many interesting aspects of quantum plasmas [20 43]. In this paper we develop a method in order to solve for the 2

3 quantized eigenvalues of the Schrödinger-Poisson system for the electron gas of arbitrary degeneracy via the pseudoforce concept [44]. We show that such a plasmonic system has two distinct characteristic lengths and when confined to the infinite potential well, similar to the quantum particle problem, leads to quantization of the plasmon excitation energies. Current results show akin similarities to the problem of particle confined to the infinite potential well and further reveals the fundamental features of collective effects in quantization of plasmon energy levels and properties of electron wave function in a free electron gas. The paper is organized as follows. The equation of state (EoS) of the finite temperature degenerate free electron gas is described in Sec. II. In Sec. III the hydrodynamic formulation is presented. The equivalent nonlinear Schrödinger-Poisson (SP) formulation is given in Sec. IV. The coupled pseudoforce model [45] and its standing wave solution is given in Sec. V. The numerical analysis and discussion is presented in Sec. VI and conclusions are drawn in Sec. VII. II. EQUATION OF STATE OF THE FREE ELECTRON GAS Let us consider a free electron gas with an arbitrary degree of degeneracy characterized by the number density and statistical quantum pressure for electron fluid as follows n(µ,t) = 21/2 m 3/2 π 2 h 3 P(µ,T) = 23/2 m 3/2 3π 2 h EdE e β(e µ) +1, E 3/2 de e β(e µ) +1. (1a) (1b) where β = 1/(k B T) and E is the energy, µ is the chemical potential of electron gas and other parameters in (4) have their conventional physical meanings. The number density and pressure in (4) can also be rewritten in terms of familiar analytic Li functions n(µ,t) = NLi 3/2 [ exp( βµ)], P(µ,T) = N β Li 5/2[ exp( βµ)], (2) where Li ν (z) is the polylogarithm function of the order ν with the argument z. Also, the prefactor N appearing in (2) is defined as below N = 2 Λ 3 = 2 ( mkb T 2π h 2 ) 3/2, (3) 3

4 with Λ being the electron de Broglie thermal wavelength. The EoS of isothermal free electron gas with arbitrary degeneracy is then written as [46] P(µ,T) = n(µ,t) Li 5/2 [ exp( βµ)] β Li 3/2 [ exp( βµ)]. (4) Note that in the classical limit βµ 1 we have Li ν [ exp( βµ)] exp( βµ) and the classical isothermal EoS, P = nk B T in (4) is retained. In the opposite fully degenerate extreme limit (βµ 1), we have lim T 0 µ = ǫ F in which ǫ F is the corresponding Fermi temperature and the polylog function is approximated as lim x Li ν ( e x ) = x ν /Γ(ν+1). Therefore, in the fully degenerate isothermal free electron gas we have P = (2/5)nk B T F in which T F = ǫ F /k B is the Fermi temperature of the (zero temperature T T F ) electron gas. The Fermi energy is related to the electron number density via ǫ F = h 2 kf 2 /2m with k F = (3π 2 n) 1/3 being the Fermi wavenumber. The typical value of the Fermi energy for metals range from 1 to 12eV. Using (2) it is confirmed that the identity P/ µ = n(µ,t) generally holds for the isothermal electron gas. In what follows we consider different models which describe collective quantum aspects of electron gas via the hydrodynamic coupling of quantum Bohm and electrostatic potentials. Some clarification and distinction between current models and the Hartree many body approximation, which assumes a description based on a single Slater determinant, usually considered in solid state literature, is in order. The complete many body description of a degenerate electron gas must be an idealistic formulation based on the N-electron wave function, ψ(x 1,x 2,...,x N,t) with all interactions between electron taken into account. As it is quite obvious this is a formidable task requiring extensive numerical work. However, in the weak correlation regime the electron gas is considered as essentially noninteracting (except the strong short range Pauli exclusion interactions dominating the system) where the many electron wave function can be decoupled into many single electron wave functions with one Schrödinger equation governing each and with the common electrostatic potential of the gas coupling them to each other and to the Poisson s equation. Therefore, the following hydrodynamic and Schrödinger-Poisson models are based on the Dawson s multistrem model [20]. Some recent developments based on the Wigner-Poisson and quantum hydrodynamic models appear in Refs. [47, 48] (see also the references 4

5 therein). III. THE HYDRODYNAMIC MODEL The dynamic properties of an electron fluid can be explored via the hydrodynamic model consisting of continuity, momentum balance and Poisson s equations, as follows n + (nu) = 0, t (5a) [ ] u m t +(u )u = (eφ µ+eφ ext )+ h2 2m φ = 4πen, ( ) n, (5b) n in which we have used the identity P(µ,T) = n(µ,t) µ as an isothermal electron fluid EoS, described above. Unlike the electrostatic potential φ which is due to the electron gas itself, the component φ ext can be any externally applied potential in the environment. The second term in the right hand side of the momentum balance equation denotes the quantum Bohm potential characterizing the electron quantum diffraction effect. The model (5) consists of a closed equation set to be used for full description of dynamics and evolution of electron plasma oscillations in an unmagnetized arbitrary degenerate fermionic gas. Note thatinthezero temperature, T T F, limit thefermi energyreplaces thechemical potential in (5). At finite temperatures, the temperature dependence of the chemical potential may be obtained using the following Sommerfeld expansion in the limit T < T F [ ( ) 2 ( ) 4 µ(t) = ǫ F 1 π2 T π4 T + ]. (6) 12 T F 80 T F (5c) IV. THE SCHRÖDINGER-POISSON MODEL Let us now consider the 1D analogue of (5) model in the absence of external electrostatic potential φ ext = 0. Thus, using the well-known Madelung transformation ψ(x,t) = R(x,t)exp[iS(x,t)/ h] in which R(x,t) = n(x,t) and u(x,t) = (1/m) S(x,t)/ x the 1D hydrodynamic model (5) can be cast into the following effective Schrödinger-Poisson model 5

6 [22] with a full nonlinear contribution to the field-density perturbations i h ψ t = h2 ψ 2 eφψ +µ(r,t)ψ, 2m x (7a) 2 φ 2 x = 4πen 0R 2. 2 (7b) Notehowever that theappearance ofn 0 which isthe equilibrium number density ofarbitrary degenerate electron gas in Poisson s equation indicates the normalization of R with R 0 in the model. Note also that (7) has a potential beyond the original Schrödinger equation by making contribution to the electrostatic interactions and well describes the plasmon wave packet evolution as will be consequently revealed. V. PSEUDOFORCE MODEL AND STANDING WAVE SOLUTION We are now interested in obtaining the small amplitude standing wave solution of the free electron arbitrary degenerate gas confined in an infinite potential well of length l. To this end, let us consider the solutions of type ψ(x,t) R(x)exp[iS(t)/ h] and φ φ(x). This ensures stationary wave types with standing u = 0 character inside the potential well, as desired. The model (7) is then decomposed into the following coupled pseudoforce model h 2 d 2 R(x) +eφ(x)r(x) µ(r,t)r(x) = ǫr(x), 2m dx 2 (8a) i h ds(t) = ǫs(t), dt (8b) d 2 φ(x) = 4πen dx 2 0 R 2 (x), (8c) in which ǫ is the energy eigenvalue of the coupled Schrödinger-Poisson system. The equation (8b) immediately yields S(t) = exp( iǫt/ h). In order for obtaining the small amplitude solution to the system (8) let us study the linearized system in which µ(r,t) µ 0 (R 0,T). Then we have the simplified following coupled differential equation system to solve with appropriate initial values d 2 R dx + 2m(ǫ µ 0) 2 h 2 R 2me h 2 φ = 0, d 2 φ dx 8πen 0R = 0. 2 (9a) (9b) 6

7 The coupled pseudoforce system (9) has the following simple solutions with the tentative initial conditions R(0) = R (0) = φ (0) = 0 and φ(0) = φ 0 R(x) = A R [cos(k + x) cos(k x)], φ(x) = A φ [ cos(k x) (k /k + ) 2 cos(k + x) ], (10a) (10b) where k + and k are the characteristic high and low plasmon wavenumbers, respectively. These parameters together with the corresponding solution coefficients are given below k + = αk p A R = ǫ2 p δǫ 2, k = αk p eφ 0 /δǫ 1 4ǫ 2p /δǫ2, A φ = 1 1 4ǫ2 p δǫ 2, α = δǫ 2ǫ p, (11a) (k + /k 0 ) 2 φ 0, (11b) (δǫ/ǫ p ) 1 4ǫ 2p /δǫ2 where δǫ = ǫ µ 0, k p = 2mǫ p / h and ǫ p = hω p is the plasmon energy quanta with ω p = 4πe 2 n 0 /m being the electron plasma frequency and n 0 the unperturbed electron gas number density. Note that we use µ 0 for the unperturbed chemical potential of the electron gas. We may now proceed with quantizing the energy eigenvalues by confining the free electron gas in the infinite potential well of width l. It is then required that the real function R(x) representing the number density of the electron fluid to vanish at the boundaries of the well. However, we use a more general periodic boundary condition as (k + ± k )l = 2πn with an arbitrary quantum integer number n, not to be confused with the electron gas number density. We then obtain the following quantized energy eigenvalue ǫ, ( ǫ n = µ 0 +2ǫ p 1+ λ2 pn 2 ), λ 2l 2 p = 2π. (12) k p At zero temperature, T T F, we have the limiting value of µ 0 = ǫ F. The quantized values of the characteristic wavenumbers then read k + = k p 1+ λ2 pn l 2 ( 1+λ 2 p n 2 /2l 2) 2, (13a) k = k p 1+ λ2 pn 2 2l ( 1+λ 2 p n 2 /2l 2) 2, (13b) 7

8 It should be noted that for any given value of n, the high an low wavenumbers k + and k are reciprocal of each other, i.e., k + k = 1. In the normalized scheme l = l/λ p, x = x/λ p, k ± = k ± /k p, ǫ = ǫ/ǫ p and µ 0 = µ 0 /ǫ p and using the normalization 1 0 R2 (x)dx = 1 to define the parameter φ 0 in (10), one arrives at following simplified results for energy values of the free electron gas δǫ n = ǫ n µ 0 = 2+ n2 l, k 2 ± = 1+ n2 1 1± 1 2l 2 (1+n 2 /2l 2 ) 2. (14) The normalized solutions corresponding to (10) are then given as follows ( ) ( ) k+ k k+ +k R n (x) = [cos(k + x) cos(k x)] = 2sin x sin x, (15a) 2 2 [ (1+n φ n (x) = 2 /2l 2) + (1+n 2 /2l 2 ) 1] [ ( ) 2 2 k cos(k x) cos(k + x)]. (15b) It is clearly rmarked that solutions ψ n (x,t) = R n (x)exp(iǫ n t/ h) have two characteristic length scales 2π/k ± corresponding to high frequency single-particle and low frequency collective excitations. It is also confirmed that the free electron energy dispersion in the Fermi gas follows the familiar relation δǫ n = h 2 k 2 n /2m with k2 n = k2 + k2 + k + being the total quantized momentum. We then have formulated the double length-scale theory of the arbitrary degenerate free electron gas confined in a potential well exactly analogous to the standard quantum theory of particle in a box. Current theory may be easily extended to 3D case with the isotropic electrostatic potential of the form φ(x) related to the common electron number density by Poisson s equation for all three independent dimensions. By using the standard separation of variables procedure (before the lumbarization) with the assumption ψ(x,y,z) = ψ 1 (x)ψ 2 (y)ψ 3 (z) in (8), three similar independent coupled linear equations for each dimension can be obtained as follows d 2 R i (x i ) dx i 2 d 2 φ i (x i ) dx i 2 + 2m(ǫ i µ) h 2 R i (x i ) 2me h 2 φ i(x i ) = 0, 8πen 0 R(x i ) = 0, (i = 1,2,3). (16a) (16b) where x 1 = x, x 2 = y and x 3 = z. Therefore, for an electron gas confined in a cubic box of sides l one obtains the eigenvalues given in (14) with new definitions n 2 = n 2 x +n2 y +n2 z and δǫ = ǫ i = h 2 ki/2m 2 with ki 2 = ki+ 2 +ki. 2 The normalized wavefunction ψ(x,y,z) i i 8

9 assuming φ 01 = φ 02 = φ 02 = φ 0 and taking other initial values to be zero is given below ψ(x,y,z) = 8 [ ( ( ( )] ki+ k i ki+ +k i iǫi t sin x i )sin x i )exp. (17) 2 2 h i It is remarked that in 2Dand3D cases some of theplasmon wavefuntions become degenerate with the same energy values. VI. NUMERICAL ANALYSIS AND DISCUSSION Figure 1 shows the variations of the number density and plasmon parameters with the chemical potential µ 0 of the arbitrary degenerate free electron Fermi gas. Figure 1(a) indicates two distinct classical (µ 0 < 0) and quantum (µ 0 > 0) regimes for the number density of electrons. The large increase in the electron temperature only significantly increases the number density of electrons for chemical potential values around µ 0 = 0. The plasmon energy variation with the chemical potential is shown in Fig. 1(b) indicating typical values of few electron Volts for metals (this value for the Aluminium is 15eV) and a sharp increase with increase in the electron number density. It is also observed that the increase in temperature of the electron gas leads to significant increase in the value of plasmon energy for regions of low chemical potential values. Variations in the plasmon wavenumber and characteristic length with chemical potential are depicted in Figs. 1(c) and 1(d). It is remarked that the plasmon wavelength decreases sharply for denser gases particularly for lower gas temperature values. Figure 2 shows the discrete energy values of plasmon excitation for different values of the normalized width of the infinite potential well. The δǫ = 0 corresponds to the energy eigenvalue of ǫ = µ 0. It is remarked that the plasmon spectrum always has an energy gap of value 2ǫ p. Therefore, the quantized plasmon energy are given by ǫ = µ 0 +2ǫ p +2ǫ e n 2 in which ǫ e is the ground state energy of an electron in an infinite well of same length and n is the plasmon quantum number. It is also clearly evident that the density of states (DoS) at an specified energy range decreases as the width of the potential well decreases. Figure 3 depicts the variations in the wavenumber eigenvalues of the confined electron gas. It is remarked that, compared to the case of single electron in a potential well, the electron gas possesses two characteristic lengths [45] corresponding to wavenumbers k ± as indicated infig. 3. Thehighervaluek + correspondstothesingleelectronoscillationwhilek isdueto 9

10 thecollectiveoscillationsofelectrongas. Itisalsorevealedthatintheclassical limit(l ) the two wavenumber values merge into single value. It is also seen that decrease of the width of the potential well leads to increase/decrease of the positive/negative wavenumber branch value. This feature is shown to become more significant for higher quantum number values. It seen that the value of k vanishes at very small vales of normalized length l. Figure 4 shows the wavefunction, probability density, electrostatic potential and electric field profiles in the infinite potential well in the ground state n = 1. Figure 4(a) differs fundamentally from the ground state wavefunction of a single electron due to the influence of collective electrostatic effects of free electron gas. The probability density of the corresponding state depicted in fig. 4(b) many definite places inside the potential well where th electron number vanishes. Note that the integral R 2 (x)dx gives the fractional number of electrons distributed within the spacial range dx. The profile in Fig. 4(b) confirms that in the ground state of the free electron gas electrons are mostly distributed in the middle of the potential well. The later feature is quite similar to the ground state probability of finding a particle confines inside the impenetrable well. The electrostatic potential and field profiles shown in Fig. 4(c) and 4(d) reveal that the electrostatic energy inside the well maximized in the middle of the bix where mos of the particles exist in the ground state. In Fig. 5 we have shown the probability distribution of electrons in ground state of electron gas in infinite potentials with different widths. As it is clearly evident, the decrease in the width of the potential leads to decrease in the number of probability density maxima at which electrons are accumulated. Further decrease of the well sized eventually can lead to the ideal quantum case of single electron confinement with a profile of a single maxima similar to the case of a particle in a 1D box. Figure 6 depicts the waqvefunction profile of the free electron Fermi gas confined in a one dimensional (1D) infinite potential well at different quantum state. The spiky profiles indicate oscillations in the electron number through the potential well. However, the n = 2 and other higher energy states have radical difference with the ground state one in that there are electron gas voids at broader spacial ranges at higher energy levels. However, the envelop of the wavefunction profiles at higher quantized energy values has a close resemblance to that of a particle in a box. The time evolution of the wavefunction for the 1D free electron gas at different quantum states is depicted in Fig. 7. Figure 8, on the other hand, shows the colorful wavefunction patterns for the a 2D free electron degenerate gas at it different 10

11 plasmon energy state. Figures 8(c) and 8(d) show two different wavefunctions with the same plasmon energy, so called the degenerate states. Current model in its present form fully accounts for the long range Coulomb collective effects in the context of single-electron wave function approximation considered also in Hartree-Fock method. The later method however also accounts for the electrostatic interactions via the self-consistent field in the mean field approximation as well as other electron gas correlations called the Fermi correlations which appear as electron exchange effect. It is however a possible extension for current model to include the exchange correlation effect as density functional form which appears in recent quantum plasma literature [47]. It should be noted that one of the greatest advantages of current model compared to the more general nonlinear models is that it provides important analytic results in the linear limit and a linear combination form for the electrostatic potential and wave function solution are presented for the electron gas system with arbitrary degree of degeneracy. VII. CONCLUSION Using the effective Schrödinger-Poisson system we have studied the quantum properties of the free electron Fermi gas with arbitrary degeneracy and have found that the energy of the free electron degenerate gas confined to an impenetrable potential well is quantized. An energy gap approximately equal to two plasmon energy exist in the quantized spectrum of the gas. The system of degenerate quantum electron gas is shown to have two different characteristic length scale corresponding to single- and multi-particle effects. Current theory of plasmon excitations can be extended to investigate the plasmon interactions and propagation in different environments. [1] M. Planck, Ann. Physik 1, 73(1900) [2] E. Schrödinger, Phys. Rev. 28, 1049(1926) [3] L. de Broglie, Ann. de la Fondation 12, 1(1987) 11

12 [4] C. Kittel, Introduction to Solid State Physics, (John Wiely and Sons, New York, 1996), 7th ed. [5] N. W. Ashcroft and N. D. Mermin, Solid state physics (Saunders College Publishing, Orlando, 1976). [6] H. Haug and S. W. Koch, Quantum Theory of the Optical and Electronic Properties of Semiconductors (World Scientific, Singapore, 2004). [7] C. Hu, Modern Semiconductor Devices for Integrated Circuits (Prentice Hall, Upper Saddle River, New Jersey, 2010) 1st ed. [8] K. Seeger, Semiconductor Physics (Springer, Berlin, 2004) 9th ed. [9] P. A. Markovich, C.A. Ringhofer, and C. Schmeister, Semiconductor Equations (Springer, Berlin, 1990). [10] G. Mangredi, Phys. Plasmas 25, (2018); [11] M. Bonitz, D. Semkat, A. Filinov, V. Golubnychyi, D. Kremp, D.O. Gericke, M.S. Murillo, V. Filinov, V. Fortov, W. Hoyer, J. Phys. A, (2003). [12] E. Fermi and E. Teller, Phys. Rev. 72, 399 (1947). [13] E. Madelung, Z. Phys., (1926). [14] F. Hoyle and W. A. Fowler, Astrophys. J. 132, 565(1960). [15] S. Chandrasekhar, An Introduction to the Study of Stellar Structure, The University of Chicago Press, Chicago (1939). [16] D. Bohm and D. Pines, Phys. Rev (1953). [17] D. Pines, Phys. Rev (1953). [18] P. Levine and O. V. Roos, Phys. Rev, (1962). [19] Y. Klimontovich and V. P. Silin, in Plasma Physics, edited by J. E. Drummond(McGraw-Hill, New York, 1961). [20] G. Manfredi, How to model quantum plasmas, Fields Inst. Commun. 46, (2005); in Proceedings of the Workshop on Kinetic Theory(The Fields Institute, Toronto, Canada 2004): ph/ [21] P. K. Shukla, B. Eliasson, Nonlinear aspects of quantum plasma physics Phys. Usp (2010). [22] G. Manfredi and F. Haas, Phys. Rev. B 64, (2001); [23] F. Haas, Quantum Plasmas: An Hydrodynamic Approach (Springer, New York, 2011). 12

13 [24] G. Brodin and M. Marklund, New J. Phys. 9, 277(2007). [25] M. Marklund and G. Brodin, Phys. Rev. Lett. 98, (2007). [26] N. Crouseilles, P. A. Hervieux, and G. Manfredi, Phys. Rev. B 78, (2008). [27] F. Asenjo, V. Munoz, J.A. Valdivia, S. H. Hahajan, Phys. Plas., (2011). [28] G. Brodin, M. Marklund, and G. Manfredi, Phys. Rev. Lett., (2008). [29] M. Marklund, G. Brodin, L. Stenflo, and C. S. Liu, Europhys. Lett., (2008). [30] F. Stern, Phys. Rev. Lett., (1967). [31] Young-Dae Jung, Phys. Plasmas 2, 332(1995); doi.org/ / [32] Young-Dae Jung, Phys. Plasmas 8, 3842(2001); doi.org/ / [33] A. E. Dubinov and I. N. Kitayev, Semiconductors 51 14(2017); doi.org/ /s [34] P. K. Shukla and L. Stenflo, J. Plasma Phys (2008); doi: /s [35] P. K. Shukla and B. Eliasson, Phys. Rev. Lett. 108, (E) (2012); ibid. 109, (E) (2012). [36] M. Marklund, P.K. Shukla, Rev. Mod. Phys (2006). [37] M. Marklund, B. Eliasson, P. K. Shukla, Phys. Rev. E 76, (2007); doi: /physreve [38] S. Das Sarma and A. Madhukar, Phys. Rev. B 23, 805(1981). [39] G. Brodin and L. Stenflo, Phys. Plasmas 24, (2017); doi.org/ / [40] Z. Moldabekov, T. Schoof, P. Ludwig, M. Bonitz and T. Ramazanov, Phys. Plasmas 22, (2015). [41] Z. Moldabekov, Tim Schoof, Patrick Ludwig, Michael Bonitz, and Tlekkabul Ramazanov, Phys. Plasmas, 22, (2015); doi.org/ / [42] L. Stanton and M. S. Murillo, Phys. Rev. E 91, (2015) [43] M. Akbari-Moghanjoughi, Phys. Plasmas 22, (2015); doi: / [44] M. Akbari-Moghanjoughi, Phys. Plasmas 17, (2017); doi: / [45] M. Akbari-Moghanjoughi, Phys. Plasmas 25, (2018); doi: / [46] M. Akbari-Moghanjoughi, and B. Eliasson, Phys. Scripta, 91 10(2016). [47] J. Hurst, K. L. Simon, P. A. Hervieux, G. Manfredi and F. Haas, Phys. Rev. B 93, (2016). [48] F. Haas, G. Manfredi, P. K. Shukla, and P.-A. Hervieux, Phys. Rev. B, 80, (2009). 13