Electron-Electron Interactions

1. Electron-Electron Interactions DragicaVasileska Professor Arizona State university

2. Classification of Scattering Mechanisms

3. Treatment of the Electron-Electron Interactions • Electron-electron interactions can be treated either in: • K-space, in which case one can separate between • Collective plasma oscillations • Binary electron-electron collisions • Real space Molecular dynamics • Bulk systems (Ewald sums) • Devices (Coulomb force correction, P3M, FMM)

4. K-space treatment of the Electron-Electron Interactions

5. Electron Gas • As already noted, the electron gas displays both collective and individual particle aspects. • The primary manifestations of the collective behavior are: • Organized oscillations of the system as a whole – plasma oscillations • Screening of the field of any individual electron within a Debye length

6. Collective excitations In the collective excitations each electron suffers a small periodic perturbation of its velocity and position due to the combined potential of all other electrons in the system. The cumulative potential may be quite large since the long-range nature of the Coulomb potential permits a very large number of electrons to contribute to the potential at a given point The collective behavior of the electron gas is decisive for phenomena that involve distances that are larger than the Debye length For smaller distances, the electron gas is best considered as a collection of particles that interact weakly by means of screened Coulomb force.

7. Collective behavior, Cont’d • For the collective description to be valid, it is necessary that the mean collision time, which tends to disrupt the collective motion, be large compared to the period of the collective oscillation. Thus: • Examples for GaAs: • ND=1017 cm-3, p =2×1013, coll >>2/p 1/coll <<3×1012 1/s • ND=1018cm-3, p=6.32×1013, coll >>2/p 1/coll<<1013 1/s • ND=1019cm-3, p=2×1014, coll >>2/p 1/coll<<3×10131/s

8. Collective Carrier Scattering Explained Consider the situation that corresponds to the mode q=0, when all electrons in the system have been displaced by the same amount u, as depicted in the figure below: u - --- - -- - - - - - -- E d + + + + + + + + + + + + + +

9. Collective Carrier Scattering Explained Because of the positive (negative) surface charge density at the bottom (top) slab, an electric field is produced inside the slab. The electric field can be calculated using a simple parallel capacitor model for which: The equation of motion of a unit volume of the electron gas of concentration n is:

10. Collective Carrier Scattering Explained • Comments: • Plasma oscillation is a collective longitudinal excitation of the conduction electron gas. • A PLASMON is a quantum of plasma oscillations. PLASMONS obey Bose-Einstein statistics. • An electron couples with the electrostatic field fluctuations due to plasma oscillations, in a similar manner as the charge of the electron couples to the electrostatic field fluctuation due to longitudinal POP.

11. Collective Carrier Scattering Explained • The process is identical to the Frӧhlich interaction if plasmon damping is neglected. Then: • Note on qmax: • Large qmax refers to short-wavelength oscillations, but one Debye length is needed to screen the interaction. Therefore, when qmax exceeds 1/LD, the scattering should be treated as binary collision. • qc=min(qmax,1/LD)

12. Collective Carrier Scattering Explained • Importance of plasmon scattering • Plasma oscillations and plasmon scattering are important for high carrier densities • When the electron density exceeds 1018 cm-3 the plasma oscillations couple to the LO phonons and one must consider scattering from the coupled modes www.engr.uvic.ca/.../Lecture%207%20-%20Inelastic%20Scattering.ppt

14. Electron-Electron Interactions(Binary Collisions) This scattering mechanism is closely related to charged impurity scattering and the interaction between the electrons can be approximated by a screened Coulomb interaction between point-like particles, namely: Then, one can obtain the scattering rate in the Born approximation as one usually does in Brooks-Herring approach.

15. Binary Collissions To write the collision term, one needs to define a pair transition rate S(k1,k2,k1’,k2’), which represents the probability per unit time that electrons in states k1 and k2 collide and scatter to states k1’ and k2’, as shown diagramatically in the figure below: k2’ k1’ r1 r2 k2 k1

16. Binary Collisions, Cont’d The pair transition rate is defined as: Since the interaction potential depends only upon the distance between the particles, it is easier to calculate M12 in a center-of-mass coordinate system, to get:

17. Binary Collisions Scattering Rate To evaluate the scattering rate due to binary carrier-carrier scattering, one weights the pair transition rate that a target carrier is present and by the probability that the final states k1’ and k2’ are empty: Note that a separate sum over k1’ is not needed because of the momentum conservation -function. For non-degenerate semiconductors, we have:

18. Binary Collisions Scattering Rate In summary:

19. Incorporation of the electron-electron interactions in EMC codes For two-particle interactions, the electron-electron (hole-hole, electron-hole) scattering rate may be treated as a screened Coulomb interaction (impurity scattering in a relative coordinate system). The total scattering rate depends on the instantaneous distribution function, and is of the form: Screening constant • There are three methods commonly used for the treatment of the electron-electron interaction: • A. Method due to Lugli and Ferry • B. Rejection algorithm • C. Real-space molecular dynamics

20. (A) Method Due to Lugli and Ferry • This method starts form the assumption that the sum over the distribution function is simply an ensemble average of a given quantity. • In other words, the scattering rate is defined to be of the form: • The advantages of this method are: • 1.The scattering rate does not require any assumption on the form of the distribution function • 2. The method is not limited to steady-state situations, but it is also applicable for transient phenomena, such as femtosecond laser excitations • The main limitation of the method is the computational cost, since it involves 3D sums over all carriers and the rate depends on k rather on its magnitude.

21. (B) Rejection Algorithm • Within this algorithm, a self-scattering mechanism, internal to the interparticle scattering is introduced by the following substitution: • When carrier-carrier collision is selected, a counterpart electron is chosen at random from the ensemble. • Internal rejection is performed by comparing the random number with:

22. If the collision is accepted, then the final state is calculated using: • where: • The azimuthal angle is then taken at random between 0 and 2. • The final states of the two particles are then calculated using:

25. Real-Space Treatment of the Electron-Electron Interactions Bulk Systems Semiconductor Device modeling

26. Bulk Systems

27. (C) Real-space molecular dynamics • An alternative to the previously described methods is the real-space treatment proposed by Jacoboni. • According to this method, at the observation time instant ti=it, the total force on the electron equals the sum of the interparticle coulomb interaction between a particular electron and the other (N-1) electrons in the ensemble. • When implementing this method, several things need to be taken into account: • 1. The fact that N electrons are used to represent a carrier density n = N/V means that a simulation volume equals V = N/n. • 2. Periodic boundary conditions are imposed on this volume, and because of that, care must be taken that the simulated volume and the number of particles are sufficiently large that artificial application from periodic replication of this volume do not appear in the calculation results.

28. Using Newtonian kinematics, the real-space trajectories of each particle are represented as: • and: • Here, F(t) is the force arising from the applied field as well as that of the Coulomb interaction: • The contributions due to the periodic replication of the particles inside V in cells outside is represented with the Ewald sum:

29. Simulation example of the role of the electron-electron interaction: • The effect of the e-e scattering allows equilibrium distribution function to approach Fermi-Dirac or Maxwell Boltzmann distribution. • Without e-e, there is a phonon ‘kink’ due to the finite energy of the phonon

30. Semiconductor Device Modeling

31. Ways of accounting for the short-range Coulomb interactions Long-range Coulomb interactions are accounted for via the solution of the Poisson equation which gives the so-called Hartree term If the mesh is infinitely small, the full Coulomb interaction is accounted for However this is not practical as infinite systems of algebraic equations need to be solved To avoid this difficulty, a mesh size that satisfies the Debye criterion is used and the proper correction to the force used to move the carriers during the free-flight is added

32. Earlier Work – k-space treatment of the Coulomb interaction Good for 2D device simulations Requires calculation of the distribution function to recalculate the scattering rate at each time step and the screening which is time consuming Implemented in the Damocles device simulator

33. K-space Approach

34. Present trends – Real-space treatment • Requires 3D device simulator, otherwise the method fails • There are several variants of this method • Corrected Coulomb approach developed by Vasileska and Gross • Particle-particle-particle-mesh (p3m) method by Hockney and Eastwood • Fast Multipole method

35. Real Space Treatment Cont’d • Corrected Coulomb approach and p3m method are almost equivalent in philosophy • FMM is very different • Treatment of the short-range Coulomb interactions using any of these three methods accounts for: • Binary collisions + plasma (collective) excitations • Screening of the Coulomb interactions • Scattering from multiple impurities at the same time which is very important at high substrate doping densities

36. 1. Corrected Coulomb approach A resistor is first simulated to calculate the difference between the mesh force and the true Coulomb force Cut-off radius is defined to account for the ions (inner cut-off radius) Outer cut-off radius is defined where the mesh force coincides with the Coulomb force Correction to the force is made if an electron falls between the inner and the outer radius The methodology has been tested on the example of resistor simulations and experimental data are extracted

37. Corrected Coulomb Approach Explained

38. Resistor Simulations

39. MOSFET: Drift Velocity and Average Energy

40. 2. p3m Approach

41. Details of the p3m Approach

42. Impurity located at the very source-end, due to the availability of Increasing number of electrons screening the impurity ion, has reduced impact on the overall drain current.

43. 3. Fast Multipole Method Different strategy is employed here in a sense that Laplace equation (Poisson equation without the charges) is solved. This gives the ‘Hartree’ potential. The electron-electron and electron-ion interactions are treated using FMM The two contributions are added together Must treat image charges properly. Good news is that the surfaces are planar and the method of images is a good choice

44. Idea

45. The philosophy of FMM:Approximate Evaluation

46. Ideology behind FMM

47. Simulation Methodology

48. Method of Images

49. Resistor Simulations

50. More on the Electron-Electron Interactions for Q2D Systems Exchange-correlation effects Screening of the Coulomb interaction potential