当前位置:首页 >> >>

Exact broken-symmetry states and Hartree-Fock solutions for quantum dots at high magnetic f

arXiv:cond-mat/0409346v1 [cond-mat.mes-hall] 14 Sep 2004

Exact broken-symmetry states and Hartree-Fock solutions for quantum dots at high magnetic ?elds
B. Szafran a,b, F.M. Peeters a, S. Bednarek b, and J. Adamowski b
a Departement b Faculty

Natuurkunde, Universiteit Antwerpen (Campus drie Eiken), B-2610 Antwerpen, Belgium

of Physics and Nuclear Techniques, AGH University of Science and Technology, Krak?w, Poland o

Abstract Wigner molecules formed at high magnetic ?elds in circular and elliptic quantum dots are studied by exact diagonalization (ED) and unrestricted Hartree-Fock (UHF) methods with multicenter basis of displaced lowest Landau level wave functions. The broken symmetry states with semi-classical charge density constructed from superpositions of the ED solutions are compared to the UHF results. UHF overlooks the dependence of the few-electron wave function on the actual relative positions of electrons localized in di?erent charge puddles and partially compensates for this neglect by an exaggerated separation of charge islands which are more strongly localized than in the exact broken-symmetry states. Key words: quantum dots, Wigner molecules, PACS: 73.21.La, 73.20.Qt

Preprint submitted to Elsevier Science

2 February 2008

At high magnetic ?elds the electron systems in circular quantum dots form Wigner molecules [1,2] in the internal structure of the system. Deformation [3,4] of the circular symmetry allows the molecules to appear in the laboratory frame only [4] when the classical counterpart [5] of the quantum system possesses a single lowest-energy con?guration. Otherwise the charge density at high ?eld is a superposition of equivalent semi-classical densities and the quantum system undergoes symmetry transformations when the magnetic ?eld is increased [4]. These transformations are associated with level crossings at which the ground state is two-fold degenerate. Superposition of the states of the degenerate levels allows [4] to extract the semi-classical brokensymmetry charge density into the laboratory frame. On the other hand the unrestricted Hartree-Fock (UHF) produces [1] broken-symmetry states for Wigner molecules. In the in?nite magnetic ?eld limit UHF gives exact [6] results for the total energy. At ?nite magnetic ?elds for which exact brokensymmetry eigenstates exist the artifactal symmetry breaking cannot be blamed to the inaccuracy of the UHF [4]. In this paper we look for the e?ects neglected by UHF comparing the ED and UHF solutions for elliptical and circular dots. We assume a spin-polarization of electrons at high magnetic ?eld (0, 0, B) oriented perpendicular to the quantum dot plane and use the Landau gauge. In the ED, described in detail in Ref. [4], the single electron wave functions used for construction of the Slater determinants are obtained via diagonalization of the single electron Hamiltonian in the multicenter basis [7,8,9] of M displaced lowest Landau level wave functions

ψ(r) =

ck exp{?α[(x ? xk )2 + (y ? yk )2 ]/4 + +ieB(x ? xk )(y + yk )/2? }, h (1)

where α is treated as variational parameter. In the present UHF approach the one-electron orbitals (1) are optimized self-consistently. We study up to N = 4 electrons, use the material data of GaAs and a basis of 12 centers (xk , yk ) put on an ellipse with a size determined variationally. Basis (1) of displaced lowest Landau level wave functions reproduces [4] also higher Fock-Darwin bands. Contrary to previous multicenter HF calculations [6,9] using a single center per electron, the present HF approach produces results which are exact in the UHF [10] sense. Classical system of three electrons in an elliptical con?nement potential with hωx = 3 meV and hωy = 4 meV possesses two equivalent lowest-energy con?g? ? urations [cf. inset of Fig. 1] and the quantum system undergoes parity transformations [4] with the magnetic ?eld [cf. Fig. 1]. Superposition [4] of the two lowest-energy eigenstates √ ΨBS = (Ψeven + eiφ Ψodd )/ 2 (2) 2

yields a broken-symmetry (BS) charge density with a distinct electron separation. Fig. 1 shows that in contrast to the exact ground-state energy the UHF energy estimate is a smooth function of the magnetic ?eld. The charge densities of considered states are shown in Figs. 2(a) and 2(b) for two magnetic ?eld values corresponding to the even-odd energy crossing presented in Fig. 1. The phase φ in the BS state is chosen such that the electrons are localized at the classical Wigner molecule positions. Notice that in the UHF the separation of electrons is more pronounced than in the exact BS states. Fig. 2(c) shows the pair correlation function (PCF) [2] for the UHF and the exact BS state corresponding to the charge density of Fig. 2(a) with the position of one of the electrons ?xed at two di?erent locations: in the center and on the edge of the central charge puddle. In contrast to the exact BS state in the UHF wave function the two electrons are insensitive to the actual position of the third electron in its charge puddle. This is a consequence of the single-determinantal form of the UHF wave function, and can be easily explained for two electrons. In the spin-polarized two electron Wigner molecule the UHF spatial wave function is given by Ψα (r1 )Ψβ (r2 )?Ψβ (r1 )Ψα (r2 ), where |Ψα |2 and |Ψβ |2 are the charge densities of separate charge puddles α and β. Wave functions Ψα and Ψβ are orthogonal due to the vanishing overlap between the puddles. The calculation of the PCF gives (up to a constant) PCF(ra , rb ) = |Ψα (ra )|2 |Ψβ (rb )|2 + |Ψα (rb )|2 |Ψβ (ra )|2 . For rb inside puddle β the second term of the sum vanishes and the remaining one signi?es that the probability of ?nding an electron in point ra inside puddle α is independent of the actual position of the second electron in puddle β. The UHF self-consistency is reached only in one of the two classical orientations [cf. inset of Fig. 1] in which the UHF energy is minimal. On the other hand, the exact BS states can be oriented under an arbitrary angle [cf. Fig. 2(d)] by modifying the phase φ in Eq. (2). Moreover, since the BS state is constructed with states of opposite parities all the plots in Fig. 2(d) correspond to the same value of the kinetic, potential and electron-electron interaction energies equal to the arithmetic average of the expectation values for Ψodd and Ψeven states. Fig. 3 shows the two lowest energy levels and the UHF total energy calculated with respect to the lowest Landau level for the elliptical dot with hωx = 3 ? meV and hωy = 4 meV. For these values the classical counterpart of the four? electron system is unique and conform with the symmetry of the con?nement potential, so that the Wigner crystallization is visible in the exact quantum ground-state for an arbitrary magnetic ?eld after the MDD decay. In this case the MDD decay is a continuous process and appears at the anticrossing around 6 T. The inset to Fig. 3 shows the charge densities calculated with the UHF and 3

ED methods. In ED the charge density in between the charge maxima takes on larger values than in UHF in which the separation of electrons is more distinct like for N = 3 [cf. Fig. 2(a) and (b)]. The average radius of the charge puddle as obtained for the two-electron ground state in the circular quantum dot (? ωx = hωy = 3 meV) is displayed h ? in Fig. 4. The ED and UHF values are similar below the MDD breakdown (B < 5.6 T). The exact value present discontinuities at the angular momentum transitions. After the MDD decay the UHF value is close to the average around which ED and BS results oscillate, but at higher ?elds it becomes an upper bound for these oscillations. The inset of Fig. 4 shows the exact BS and UHF charge density for B = 11.8 and 28.9 T. BS charge densities for the two values of the magnetic ?eld have been obtained from superpositions of the degenerate states with angular momenta -5,-7 and -9,-11 (in h units), ? respectively. The charge maxima in the exact BS solutions are less strongly separated. For B = 28.9 T the charge density islands of the exact BS have crescent shape while the UHF charge puddles are more oval and the distance between them is larger. The di?erence in shape of the separated charge islands in circular dots is largest for N = 2. For larger N the charge puddles in the exact solutions are less spectacularly spread. Fig. 5 shows the comparison of the exact BS and UHF charge densities for N = 3. The BS plots for N = 3 correspond to degeneracy of states with angular momentum ?9? , ?12? for B = 7.5 T and h h ?12? , ?15? for B = 15 T. h h In summary, we have presented a comparison of the UHF and the ED results for the charge density of Wigner molecules in circular and elliptical dots. For the comparison we have used the broken-symmetry states obtained from the superposition of the exact eigenstates. We have found that the UHF exaggerates the separation of the electron charge densities in the laboratory frame. In this way the UHF method partially compensates for the overlooked correlations related to the reaction of electrons on their actual position in the separated charge density islands. This reaction is of smaller importance for larger magnetic ?elds, for which the charge density islands shrink to points, which explains the vanishing of the UHF energy overestimation in the in?nite magnetic ?eld [6]. Due to the exaggerated electron separation the charge islands forming Wigner molecules calculated in the UHF shrink with the magnetic ?eld faster than in the exact broken-symmetry states. Acknowledgments This paper was supported in part by the Polish Ministry of Scienti?c Research and Information Technology in the framework of the solicited grant PBZ-MIN-008/P03/2003, the Flemish Science Foundation (FWO-Vl), the Belgian Science Policy, and the University of Antwerpen (VIS and GOA). One of us (BS) is supported by the Foundation for Polish Science 4


[1] S.M. Reimann and M. Manninen, Rev. Mod. Phys. 74 (2002) 1283. [2] P.A. Maksym, H. Immamura, G.P. Mallon, and H. Aoki, J. Phys. Condens. Matter 12 (2000) R299. [3] M. Manninen, M. Koskinen, S.M. Reimann, and B. Mottelson, Eur. Phys. J. D. 16 (2001) 381. [4] B. Szafran, F.M. Peeters, S. Bednarek, and J. Adamowski, Phys. Rev. B 69 (2004) 125344. [5] V.M. Bedanov and F.M. Peeters, Phys. Rev. B 49 (1994) 2667. [6] B. Szafran, S. Bednarek, J. Adamowski, M. Tavernier, E. Anisimovas, and F.M. Peeters, Eur. Phys. J. D 28 (2004) 373. [7] J. Kainz, S.A. Mikhailov, A. Wensauer, and U. R¨ssler, Phys. Rev. B 65 (2002) o 115305. [8] C. Yannouleas and U. Landman, Phys. Rev. B 66 (2002) 115315. [9] B. Szafran, S. Bednarek, and J. Adamowski, Phys. Rev. B 67 (2003) 045311; J. Phys.: Cond. Matter 15 (2003) 4189. [10] B. Reusch and H. Grabert, Phys. Rev. B 68 (2003) 045309.


y [nm]


E [meV] - 3 hωc / 2


19 even 18 broken odd 17
-20 -20 0 x [nm] 20





B [T]

Fig. 1. Energy of the lowest even (dashed curve) and odd parity (solid curve) levels, the broken symmetry state (dotted curve) and UHF energy calculated for N = 3, h ? ωx = 3 meV and hωy = 4 meV. Inset shows the two equivalent classical con?gu? rations.


100 nm


B=8.75 T



even odd B=17.25 T




B=8.75 T

-20 -40 20 0 40 0 40


0 0 40 0 40


-20 -40

π/8 3/8π (d)


B=17.25 T


Fig. 2. (a) The charge densities for the even, the odd, the broken-symmetry (BS) state and the UHF for N = 3, hωx = 3 meV, ? ωy = 4 meV at B = 8.75 T. (b) Same ? h as (a) at B = 17.25 T. (c) Pair correlation function plots for UHF and BS state for an electron ?xed at (0,13.4 nm) [left panel] and (-3 nm,10.4 nm) [right panel] marked by crosses. (d) BS states for given shifts of phase φ [Eq. (2)] with respect to the BS plot in (b). Plots (a), (b) and (d) have the same scale given by the length bar in (a).




E [meV] - 2 hωc

100 nm





32 4

8 12 16 20

B [T]

Fig. 3. Lowest energy levels for N = 4, hωx = 3 meV and hωy = 4 meV (solid ? ? lines) and total UHF energy (dotted) line. Inset: Charge density calculated in the UHF (two upper plots) and the ED (two lower plots) for B = 8 T (left side) and B = 20 T (right side).



80 nm

11.8 T


28.9 T

2 1/2


<r >






B [T]

Fig. 4. Average radius of the charge puddle (< r 2 >1/2 ) for two electrons in a circular quantum dot (? ωx = hωy = 3 meV). Dotted, dashed and solid lines correspond h ? to UHF, BS, and ED results. Inset: charge density obtained in the UHF (two upper plots) and in the exact broken-symmetry state (two lower plots) for B=11.8 T (left side) and 28.9 T (right side).


UHF 100 nm

7.5 T BS

15 T

Fig. 5. Charge density for N = 3 electrons in a circular quantum dot (? ωx = hωy = 3 meV) calculated with UHF (upper plots) and for exact BS states h ? (lower plots) for B = 7.5 T (left plots) and B = 15 T (right plots).



量子化学 非限制Hartree-Fock方程
量子化学 非限制Hartree-Fock方程_物理_自然科学_专业资料。非限制 Hartree-Fock 方程 1. 波函数 波函数是描述电子运动的函数,一般也称为轨道,因为描述电子运动...
百万核 Hartree-Fock 自洽场迭代(SCF) 来源: 作者: 发布时间:2015-07-28 阅读:389 次 研究背景 研究 HIV 与硫酸茚地那韦相互作用机制,对 HIV 药物设计具有...

All rights reserved Powered by 甜梦文库 9512.net

copyright ©right 2010-2021。