ElectronElectron Interactions DRAGICA VASILESKA PROFESSOR ARIZONA STATE UNIVERSITY

Electron-Electron Interactions DRAGICA VASILESKA PROFESSOR ARIZONA STATE UNIVERSITY

Classification of Scattering Mechanisms

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, P 3 M, FMM)

K-space treatment of the Electron-Electron Interactions

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

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 longrange 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.

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 Ga. As: ND=1017 cm-3, p =2× 1013, coll >>2 / p 1/ coll <<3× 1012 1/s ND=1018 cm-3, p =6. 32× 1013, coll >>2 / p 1/ coll <<1013 1/s ND=1019 cm-3, p =2× 1014, coll >>2 / p 1/ coll <<3× 1013 1/s

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 d - - - - E + + + +

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:

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.

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)

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 -%20 Inelastic%20 Scattering. ppt

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.

Binary Collissions To write the collision term, one needs to define a pair transition rate S(k 1, k 2, k 1’, k 2’), which represents the probability per unit time that electrons in states k 1 and k 2 collide and scatter to states k 1’ and k 2’, as shown diagramatically in the figure below: k 2’ k 1’ r 1 k 1 r 2 k 2

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 M 12 in a center-of-mass coordinate system, to get:

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 k 1’ and k 2’ are empty: Note that a separate sum over k 1’ is not needed because of the momentum conservation -function. For nondegenerate semiconductors, we have:

Binary Collisions Scattering Rate In summary:

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 interaction: A. Method due to Lugli and Ferry B. Rejection algorithm C. Real-space molecular dynamics

(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 3 D sums over all carriers and the rate depends on k rather on its magnitude.

(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:

• 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:

Real-Space Treatment of the Electron-Electron Interactions BULK SYSTEMS SEMICONDUCTOR DEVICE MODELING

Bulk Systems

(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.

• 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:

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

Semiconductor Device Modeling

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 freeflight is added

Earlier Work – k-space treatment of the Coulomb interaction Good for 2 D 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

K-space Approach

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

Real Space Treatment Cont’d Corrected Coulomb approach and p 3 m 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

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

Corrected Coulomb Approach Explained

Resistor Simulations

MOSFET: Drift Velocity and Average Energy

2. p 3 m Approach

Details of the p 3 m Approach

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.

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

Idea

The philosophy of FMM: Approximate Evaluation

Ideology behind FMM

Simulation Methodology

Method of Images

Resistor Simulations

More on the Electron-Electron Interactions for Q 2 D Systems EXCHANGE-CORRELATION EFFECTS SCREENING OF THE COULOMB INTERACTION POTENTIAL

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

Space Quantization • Poisson equation: e 3 e 2 e 1 e 0 d 2 V (z) e 2 H = n(z) + N (z)] [ 2 a dz e sc • Hohenberg-Kohn-Sham Equation: (Density Functional Formalism) VG>0 e 1’ e 0’ EF z-axis [100] (depth) D 4 -band 2 h ¶ 2 * 2 + Veff (z) y n (z) = e n y n (z) 2 mz ¶z [ ] D 2 -band V (z) = V (z) + V (z) eff H xc im [100]-orientation: Finite temperature generalization of the LDA (Das Sarma and Vinter) D 2 -band : mz=ml=0. 916 m 0, mxy=mt=0. 196 m 0 D 4 -band: mz=mt, mxy= (ml mt)1/2

Exchange-Correlation Effects HF E = EHF + E corr = EHF E + kin exchange + E corr Total Ground State Energy of the System Accounts for the error made with the Hartree-Fock Approximation for the Ground State Energy Accounts for the reduction of the Ground State Enery due to the inclusion of the Pauli Exclusion Principle Ways of Incorporating the Exchange-Correlation Effects: ï Density-Functional Formalism (Hohenberg, Kohn and Sham) ï Perturbation Method (Vinter)

Subband Structure Importance of Exchange-Correlation Effects Vasileska et al. , J. Vac. Sci. Technol. B 13, 1841 (1995) (Na=2. 8 x 1015 cm-3, Ns=4 x 1012 cm-2, T=0 K) Exchange-Correlation Correction: ï Lower subband energies ï Increase in the subband separation ï Increase in the carrier concentration at which the Fermi level crosses into the second subband ï Contracted wavefunctions Thick (thin) lines correspond to the case when the exchange-correlation corrections are included (omitted) in the simulations.

Subband Structure Comparison with Experiments Kneschaurek et al. , Phys. Rev. B 14, 1610 (1976) Infrared Optical Absorption Experiment: Si. O 2 Al-Gate LED far-ir Si-Sample radiation Vg Transmission-Line Arrangement

Subband Structure Comparison with Experiments Experimental data from: F. Schäffler and F. Koch (Solid State Communications 37, 365, 1981) Unprimed ladder Primed ladder

Screening of the Coulomb Interaction

What is Screening? l. D - Debye screening length - + r 3 D: - 1 r screening cloud æ rö 1 expç - ÷ r è l. D ø - - Example: - - Ways of treating screening: • Thomas-Fermi Method static potentials + slowly varying in space • Mean-Field Approximation (Random Phase Approximation) time-dependent and not slowly varying in space

Diagramatic Description of RPA Polarization Diagrams Bare interaction Effective interaction (or ‘dressed’ or ‘renormalized’) Proper (‘irreducible’) polarization parts = + bare pair-bubble + +. . .

Screening: Simulation results are for: Na=1015 cm-3, Ns=1012 cm-2 Re lative Polarization Function: 0) (q, 0) /P (0) P (00 00(0, 0) Screening Wavevectors : 2 s qnm (q, 0) = - e P (0) nm(q, 0) 2 k T=300 K 2 D-Plasma Frequency: (q) = pl e 2 N q s xy 2 km*

Screening: Form-Factors: Na=1015 cm-3, Ns=1012 cm-2 ¥ ¥ 0 0 Fij, nm(q) = ò dz'y j*(z)y i(z)G˜(q, z, z')y *m(z')y n(z' ) Diagonal form-factors => j z z' i Off-diagonal form-factors m n