Dragica Vasileska Professor - PowerPoint Presentation

1. Dragica VasileskaProfessorArizona State universityElectron-Electron Interactions

2. Classification of Scattering Mechanisms

3. Treatment of the Electron-Electron InteractionsElectron-electron interactions can be treated either in:K-space, in which case one can separate betweenCollective plasma oscillationsBinary electron-electron collisionsReal space Molecular dynamicsBulk systems (Ewald sums)Devices (Coulomb force correction, P3M, FMM)

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

5. Electron GasAs 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 oscillationsScreening of the field of any individual electron within a Debye length

6. Collective excitationsIn 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 pointThe collective behavior of the electron gas is decisive for phenomena that involve distances that are larger than the Debye lengthFor 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’dFor 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/sND=1018 cm-3, p =6.32×1013, coll >>2/p  1/coll <<1013 1/sND=1019 cm-3, p =2×1014, coll >>2/p  1/coll <<3×1013 1/s

8. Collective Carrier Scattering ExplainedConsider 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:+ + + + + + + + + + + + + +- - - - - - - - - - - - - -udE

9. Collective Carrier Scattering ExplainedBecause 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 ExplainedComments: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 ExplainedThe 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 ExplainedImportance of plasmon scatteringPlasma oscillations and plasmon scattering are important for high carrier densitiesWhen the electron density exceeds 1018 cm-3 the plasma oscillations couple to the LO phonons and one must consider scattering from the coupled modeswww.engr.uvic.ca/.../Lecture%207%20-%20Inelastic%20Scattering.ppt

13.

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 CollissionsTo 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:k1k2k2’k1’r1r2

16. Binary Collisions, Cont’dThe 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 RateTo 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 RateIn summary:

19. Incorporation of the electron-electron interactions in EMC codesFor 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 constantThere are three methods commonly used for the treatment of the electron-electron interaction:A. Method due to Lugli and FerryB. Rejection algorithmC. Real-space molecular dynamics

20. 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.(A) Method Due to Lugli and Ferry

21. 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:(B) Rejection Algorithm

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:

23.

24.

25. Bulk SystemsSemiconductor Device modelingReal-Space Treatment of the Electron-Electron Interactions

26. Bulk Systems

27. (C) Real-space molecular dynamicsAn 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 interactionsLong-range Coulomb interactions are accounted for via the solution of the Poisson equation which gives the so-called Hartree termIf 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 interactionGood for 2D device simulationsRequires calculation of the distribution function to recalculate the scattering rate at each time step and the screening which is time consumingImplemented in the Damocles device simulator

33. K-space Approach

34. Present trends – Real-space treatmentRequires 3D device simulator, otherwise the method failsThere are several variants of this methodCorrected Coulomb approach developed by Vasileska and GrossParticle-particle-particle-mesh (p3m) method by Hockney and EastwoodFast Multipole method

35. Real Space Treatment Cont’dCorrected Coulomb approach and p3m method are almost equivalent in philosophyFMM is very differentTreatment of the short-range Coulomb interactions using any of these three methods accounts for:Binary collisions + plasma (collective) excitationsScreening of the Coulomb interactionsScattering from multiple impurities at the same time which is very important at high substrate doping densities

36. 1. Corrected Coulomb approachA resistor is first simulated to calculate the difference between the mesh force and the true Coulomb forceCut-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 forceCorrection to the force is made if an electron falls between the inner and the outer radiusThe 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 MethodDifferent 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 togetherMust 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. Exchange-correlation effectsScreening of the Coulomb interaction potentialMore on the Electron-Electron Interactions for Q2D Systems

51. Exchange-Correlation Correction to the Ground State Energy if the System

52. Space Quantization Poisson equation:d2VH(z)dz2=e2escn(z)+Na(z)[] Hohenberg-Kohn-Sham Equation: (Density Functional Formalism)Finite temperature generalizationof the LDA (Das Sarma and Vinter)Veff(z)=VH(z)+Vxc(z)+Vim(z)-h22mz*¶2¶z2+Veff(z)yn(z)=enyn(z)[]EFVG>0e0e1e2e3e0’e1’z-axis [100](depth)[100]-orientation:D2-band : mz=ml=0.916m0, mxy=mt=0.196m0D4-band: mz=mt, mxy= (ml mt)1/2D2-bandD4-band

53. Exchange-Correlation EffectsE=EHF+Ecorr=EkinHF+EexchangeHF+EcorrTotal Ground StateEnergy of the SystemHartree-Fock Approximationfor the Ground State EnergyAccounts for the error madewith the Hartree-Fock ApproximationAccounts for the reductionof the Ground State Enerydue to the inclusion of thePauli Exclusion PrincipleWays of Incorporating the Exchange-Correlation Effects:Density-Functional Formalism (Hohenberg, Kohn and Sham)Perturbation Method (Vinter)

54. Subband StructureImportance of Exchange-Correlation EffectsExchange-Correlation Correction:Lower subband energiesIncrease in the subband separationIncrease in the carrier concentration at which the Fermi level crosses into the second subbandContracted wavefunctionsVasileska et al., J. Vac. Sci. Technol. B 13, 1841 (1995)(Na=2.8x1015 cm-3, Ns=4x1012 cm-2, T=0 K)Thick (thin) lines correspond to thecase when the exchange-correlationcorrections are included (omitted) inthe simulations.

55. Subband StructureComparison with Experiments Kneschaurek et al., Phys. Rev. B 14, 1610 (1976)Infrared Optical AbsorptionExperiment:far-irradiationLEDSiO2Al-GateSi-SampleVgTransmission-Line Arrangement

56. Subband StructureComparison with Experiments Unprimed ladderPrimed ladderExperimental data from: F. Schäffler and F. Koch (Solid State Communications 37, 365, 1981)

57. Screening of the Coulomb Interaction

58. What is Screening?+---------lD - Debye screening lengthWays of treating screening:Thomas-Fermi Method static potentials + slowly varying in spaceMean-Field Approximation (Random Phase Approximation) time-dependent and not slowly varying in spacer3D: 1r 1rexp-rlDæ è ç ö ø ÷ screeningcloudExample:

59. Diagramatic Description of RPAPolarization DiagramsEffective interaction (or ‘dressed’ or ‘renormalized’)Bare interaction =+++ . . .bare pair-bubbleProper (‘irreducible’)polarization parts

60. Screening:Simulation results are for: Na=1015 cm-3, Ns=1012 cm-2Relative Polarization Function:P00(0)(q,0)/P00(0)(0,0)Screening Wavevectors:qnms(q,0)=-e22kPnm(0)(q,0)T=300 K2D-Plasma Frequency: wpl(q)=e2Nsq2kmxy*

61. Screening:Form-Factors: Na=1015 cm-3, Ns=1012 cm-2Diagonal form-factorsOff-diagonal form-factorsFij,nm(q)=dzdz'0¥ò0¥òyj*(z)yi(z)˜ G (q,z,z')ym*(z')yn(z') =>inmjz'z