Volume 10, Issue 2 Articles In and Out Trott's Corner Beyond Sudoku New Products New Publications Calendar News Bulletins New Resources Classifieds Download This Issue Editorial Policy Staff and Contributors Submissions Subscriptions Advertising Back Issues Contact Information 
Electronic Structure of MultiElectron Quantum Dots
IntroductionQuantum dots are artificially fabricated atoms, in which charge carriers are confined in all three dimensions just like electrons in real atoms. Consequently, they exhibit properties normally associated with real atoms such as quantised energy levels and shell structures. These properties are described by the electron wavefunctions whose evolution is governed by the Schrödinger equation and the Pauli exclusion principle. There are many methods available to solve the Schrödinger equation for multiple electrons. They roughly fall into the categories of the diagonalisation method, meanfield densityfunctional theory, and the selfconsistent field approach. One of the first theoretical studies of quantum dots was by Pfannkuche et al. [1], who compared the results of HartreeFock selfconsistent calculations and exact diagonalisation of the Hamiltonian for two electrons in a circularly symmetric parabolic potential. They found good agreement between the two methods for the triplet state but marked differences for the singlet state, indicating important spin correlations were not included properly in their HartreeFock model. This suggests that the proper treatment of electron spins is crucial for correctly obtaining the electronic structures in quantum dots. Examples of selfconsistent field approaches in the literature include Yannouleas and Landman [2, 3], who studied circularly symmetric quantum dots using an unrestricted spinspace HartreeFock approach, and McCarthy et al. [4], who developed a HartreeFock Mathematica package. Macucci et al. [5] studied quantum dots with up to 24 electrons using a meanfield localdensityfunctional approach, in which the spin exchangecorrelation potential was approximated by an empirical polynomial expression. Lee et al. [6] also studied an Nelectron quantum dot using densityfunctional theory, where the generalisedgradient approximation was used for exchangecorrelation potentials. Exchange interaction comes directly from the antisymmetrisation of wavefunctions as required by Pauli's exclusion principle. In densityfunctional theory this is a major problem since the mathematical object is the electron density rather than the electron wavefunction, making evaluation of the exchange interaction intrinsically difficult. Diagonalisation approaches in the literature include Ezaki et al. [7, 8], Eto [9, 10], Reimann et al. [11], and Reimann and Manninen [12], who each applied a brute force approach by numerically diagonalising the Nelectron Hamiltonian using Slater determinants composed of singleelectron eigenstates as the basis functions. This approach, namely the configuration interaction (CI) method, takes into account the full interaction and correlation of the electrons in the system as long as the numerical results converge with an increasing number of basis functions. However, such an approach involves the calculation of a very large number of interaction integrals and the inversion of large matrices, which can be prohibitively expensive in terms of computer resources. Reimann et al. [11] employed matrices of dimensions up to 108,375 with 67,521,121 nonzero elements for a sixelectron quantum dot. Calculations for any higher number of electrons were not considered numerically viable using the conventional CI formalism, even with stateoftheart computing facilities. We have recently developed a spinadapted configuration interaction (SACI) method to study the electronic structure of multielectron quantum dots [13]. This method is based on earlier work by quantum chemist R. Pauncz [14], which expands the multielectron wavefunctions as linear combinations of antisymmetrised products of spatial wavefunctions and spin eigenfunctions. The SACI method has an advantage over using Slater determinants in that a smaller basis is used. This reduces the computational resources required, allowing calculations for dots with more than six electrons on a desktop computer. After some theoretical background, we present results from a Mathematica package which employs the SACI method to calculate the energy levels and wavefunctions of multiple electrons confined by a 2D harmonic well and under the influence of a perpendicular magnetic field. Mathematica enables the exact calculation of interaction integrals which greatly improves the speed and accuracy of calculations.


About Mathematica  Download Mathematica Player © Wolfram Media, Inc. All rights reserved. 