1 Solving Schrödinger Equation for Three-Electron Quantum Systems
Total Page:16
File Type:pdf, Size:1020Kb
Solving Schrödinger Equation for Three-Electron Quantum Systems by the Use of The Hyperspherical Function Method Lia Leon Margolin 1 , Shalva Tsiklauri 2 1 Marymount Manhattan College, New York, NY 2 New York City College of Technology, CUNY, Brooklyn, NY, Abstract A new model-independent approach for the description of three electron quantum dots in two dimensional space is developed. The Schrödinger equation for three electrons interacting by the logarithmic potential is solved by the use of the Hyperspherical Function Method (HFM). Wave functions are expanded in a complete set of three body hyperspherical functions. The center of mass of the system and relative motion of electrons are separated. Good convergence for the ground state energy in the number of included harmonics is obtained. 1. Introduction One of the outstanding achievements of nanotechnology is construction of artificial atoms- a few-electron quantum dots in semiconductor materials. Quantum dots [1], artificial electron systems realizable in modern semiconductor structures, are ideal physical objects for studying effects of electron-electron correlations. Quantum dots may contain a few two-dimensional (2D) electrons moving in the plane z=0 in a lateral confinement potential V(x, y). Detailed theoretical study of physical properties of quantum-dot atoms, including the Fermi-liquid-Wigner molecule crossover in the ground state with growing strength of intra-dot Coulomb interaction attracted increasing interest [2-4]. Three and four electron dots have been studied by different variational methods [4-6] and some important results for the energy of states have been reported. In [7] quantum-dot Beryllium (N=4) as four Coulomb-interacting two dimensional electrons in a parabolic confinement was investigated. Energy spectra, charge and spin densities, and electron-electron correlations in a harmonic oscillator potential were obtained by the use of Exact Diagonal Approximation method. However, all above mentioned methods are applicable only if confinement energy is much larger than electron-electron interaction energy. In lateral quantum dots, defined by metallic gates in a 2D electron gas confinement energy and electron-electron interaction energy are almost the same order, therefore above listed approximate methods cannot provide adequate description of the system. Theoretical study of the physical properties of quantum dots as a function of external magnetic field is a very important problem to solve since it will allows us to tune physical properties of these dots by experimentally changing external magnetic field frequency. 1 In the next two sections a new model-independent approach for the description of three electron quantum dots in 2D space is developed by the use of the Hyperspherical Function Method (HFM) presented in [8-10]. The wave functions of three electron system are expanded in a complete set of three body hyperspherical functions, and potential energy of the interaction between confined electrons is described by the logarithmic function. The HFM allows us to separate the center of mass and relative motion of electrons and obtain ground state energies as a function of an external magnetic field by solving Schrödinger’s equation. 2. Mathematical Modeling of Three Electron Quantum Dots in 2 D Space by the Use of The Hyperspherical Function Method Theoretical studies of three electron quantum dots have been carried out only for Coloumb interacting electrons. However, due to the fact that the solution of Poisson equation in 2 D space for three electron quantum system is represented by the logarithmic function, it is extremely important to describe electron-electron interactions with logarithmic potential. Solving Schrodinger equation for two dimensional electrons in a parabolic confinement with Hyperspherical Function Method (HFM) allows us to separate the center of mass movement and consider logarithmic potential of electron -electron interactions. Hamiltonian for three Electrons in parabolic confinement can be written as: 2 3 1 e 1 3 ∗ 2 2 ( (1) H = ∑ pi − Ai + m ω0 ri + ∑V ()rij i=1 2meff c 2 i≠ j Where pi is the generalized momentum of the i-th particle, Ai vector potential of the magnetic field at the point occupied by i-th particle 1 1 A = B × r = B()− y , x 0, j 2 i 2 i i ω m is an effective mass of an electron, ω is strength of confinement, r = 0 unit of eff 0 0 m* length, and V(r) is electron-electron interaction potential. If we substitute logarithmic potential in (1.1) we will obtain the following expression: 2 ∗ 2 3 m ω 3 rij H = − ∇2 + r 2 − β ln −ω L ∑ i i ∑ L z i=1 2meff 2 i≠ j r0 2 Where eB is Larmor frequency, and 2 2 1 2 ωL = ∗ ω = (ω0 + ωL ) 2m c ˆ ˆ Let’s introduce mass-scaled Jacobi coordinates X i,Y j, and R of three particle systems defined by: 2/1 m m X = j k r − r i ()j k (m j + mk )µ /1 2 m (m + m ) m r + m r Y = j j k − r + j j k k (2) i j Mµ m j + m k 1 3 R = m r ∑ i i Mµ i=1 1/ 2 m m m Where M= m+ m+ m {i, j , k ,} ={12,, 3} , and µ = 1 2 3 is a reduced mass of 1 2 3 M three particle system. For the systems of three identical particles (electrons) m m1 = m2 = m3 = m , µ = and Jacobi coordinates can be found easily. 3 Schrodinger equation with Hamiltonian (1) can be written as: 2 1 2 2R2 1 2 X 2 Y 2 U X Y E − ∇R + ω − ()∆ X + ∆Y + ω ()+ + 123 (), Ψ()()3,2,1 = Ψ 3,2,1 (3) 4 2µ 4 This equation enables us to represent Eigenfunctions of three electrons as a following product: ψ ( 3,2,1 ) = ϕ(x, y)φ(R)σ (s1,s2 ,s3 ) (4) Where σ (s ,1 s ,2 s3 ) are spin functions identifying parity of both φ(R) and ϕ(x, y) functions, ϕ(x, y) describes the relative motion of electrons , φ(R) describes the movement of the center of mass. If we substitute (4) into (3) we will receive the equation of the linear harmonic oscillator for φ(R) . − 1 ∇2 + ω 2R2 φ R = E φ R (5) ( 4 R ) ( ) R ( ) Where E = E + E and ω = eB is Larmor frequency and L = l + l . Hamiltonian R x, y L 2mc x y in the equation (5) coincides with one electron Hamiltonian in 2D parabolic confinement and gives us the following Fock-Darvin energy levels (See Fig. 1): CM E m,n = ω(2n + m +1)− ωLm (7) As for the relative movements of free electrons ϕ(x, y) if we substitute (1.5) into (1.3) we will obtain the following equation 3 2 1 2 2 2 tot (8) − ()∆ x + ∆ y + ω ()x + y +U123 ()()()y,x − ωL Lz ϕ y,x = E y,x ϕ y,x 2µ 4 In order to solve equation (1.8) for the relative motion of the electrons let’s introduce the Hyperspherical coordinates in Four dimensional Euclidean space. X Y i i − ∞ ≤ ρ ≤ +∞ ; nx = ; n y = ; 0 ≤ α ≤ π / 2 ; {ρ,,,Ωi} = ρ α ix i, y i i | X | i Y| | { } i i 2 2 2 2 2 2 2 ρ = X 1 + Y1 = X 2 + Y2 = X 3 + Y3 ; | X i |= ρcos α Y|; i |= ρ sinα (9) where xi, y i define directions of X i and Yi vectors. Relationship between Hyperspherical coordinates can be written in expanded way: X x1 = ρ cos α cos ϕ 1 ; | Y | = ρ sinα cosϕ ; 1 2 X 2 = ρ cos α sin ϕ 1 ; | Y2 | = ρ sinα sinϕ2; (10) | Y3 | = ρ sinα. X 3 = ρ cos α ; Connection between different sets of Jacobi coordin ates (2) can be represented as: X k = X i cosφik + Yi sinφik (11) yk = −X i sinφik + Yi cosφik The angle φik can be easily found using the following formula. 2/1 m M φ = arctg()−1 p j (12) ik m m k i Where p can be either odd or even depending on the parity of the {i,k, j} particle permutations. Obtaining φik from (12) isn’t difficult for the systems of three identical particles. 4 Motion of CM 12 n=1 10 8 ) ω ω ω ω (h 6 cm E n=0 4 2 0 0 0.5 1 1.5 2 2.5 3 3.5 ωωω/ωωωL CM.grf ωL Fig 1.1 a) Energy levels of the movement of the center of mass as a function of ω when n=0,1 and m=0,1,…5 Dotted lines represent Landau levels Kinetic energy operator for three particle systems in Hyperspherical coordinates can be represented as: ∂ 2 3 ∂ K 2 (Ω ) T = + − i (13) ∂ρ2 ρ ∂ρ ρ2 2 K (Ωi ) is the square of the four-dimensional space angular momentum ∂ 2 ∂ 1 1 K 2 ()Ω = + + 2ctg2α − l 2 + l 2 (14) i ∂α 2 ∂α cos2 a 1 sin 2 α 2 5 2 2 ∂ where li = 2 (i=1,2) represent squares of the impulses of the corresponding Jacobi ∂ϕ j vectors Let’s rewrite equation (6) for the relative motion of free electrons in Hyperspherical coordinates: ∂ 2 3 ∂ K 2 (Ω ) 2µ 2µ i U y,x 1 2 2 E Lrel , 2 + − 2 − 2 123 ()− ω ρ + 2 −ωL Z ψ ()ρ α = 0 (15) ∂ρ ρ ∂ρ ρ 4 2 Where U123 ( y,x ) represents potential energy of interacting particles. K (Ωi ) is an angular part of four dimensional Laplace operator moment with the eigenstates of K(K+1), and eigenfunctions that create complete set of the basic ortonormal Hyperspherical functions: Φ 1ll 2 (Ω) = N n P 1ll 2 (α )eil1ϕ1 eil2ϕ2 (16) KLM 1ll 2 n Where ν ν −m ν + l ν + l 2m+l (2 ν −m)+ l 1ll 2 2 1 1 2 Pn ()()α = ∑ −1 ()()cosα sinα m=0 m ν − m L = (l1 + l2 ) 2 K = 2n + l1 + l1 , N = (2n + l1 + l2 +1)n!(n + l1 + l2 )! [2π (n + l1 )!(n + l2 )!] Let’s expand Ψ(ρ,α) wave functions in Hyperspherical basis: −3 / 2 1ll 2 1ll 2 ψ(ρ,α) = ∑ρ χ KL ( ρ )Φ KLM ( Ωi ) (17) Kl1l2M after substitution of (17) into (15) for hyper radial wave functions we will receive infinite system of the following equations: 2 2 ∂ (K + )1 − ,0 25 ′′ ′′ 2 1ll 2 1ll 2 1ll 2 1ll 2 2 − χ + 2 ϕRL (ρ) = WKK ′LL′MM ′ (ρ)ϕK L′′ (ρ) (18) ∂ρ ρ ∑ K 1ll 2′′′ M ′ Where: 1ll 2 1ll ′′ 2 2µ 1ll 2 1ll ′′ 2 W ′ ′ ′ (ρ) = ∫ Φ KLM (Ωi )U123 ()α Φ ′ ′ ′ (Ωi d) Ωi (19) KK LL MM 2 K L M 2µ and χ2 = − E (20) 2 In order to calculate overlapping integral between have functions defined on the different sets of Jacobi coordinates we have to use Reinal-Revai unitary transformation coefficients defined by the formula; ~~ 1ll 2 i ~~ j 1ll 2 i Φ KLM (Ωi ) = 〈 1ll 2 | 1ll 2 〉 KL Φ KLM (Ωi ) (21) ∑~~ 1ll 2 6 ~~ j Where 〈 1ll 2 | 1ll 2 〉 KL represent unitary transformation coefficients first introduced by Reinal and Revai [11].