当前位置:首页 >> >>

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


...Nuclear and Neutron Matter In a Relativistic Hartree-Fock ....pdf
Properties of Hot Nuclear and Neutron Matter In a Relativistic Hartree-Fock Theory_电子/电路_工程科技_专业资料。办声 A to , e e l a r P h y s ...
第44 卷第 11 期 2015 年 11 月 当代化工 Contemporary Chemical Industry Vol.44,No.11 November,2015 用 Hartree-Fock 法预测取代苯亚硒酸的 pKa 值张梦阳...
昌吉 学 院学报 2014 年第 4期 浅析《 计算材料 学》 教学 中的 Hartree-Fock 方法 以H 2 为例探讨 H a rt r e e-Fo C k 方法 朱应 涛 ...
Hartree_Fock__Method_自然科学_专业资料。摘自:http...ground-state energy of quantum many-body system.... Exact broken-symmetry ... 8页 免费 Point ...
...Density in Hartree, Hartree-Fock, and Density Fu....pdf
Hartree-Fock energy of a... 暂无评价 7页 免费 Broken Symmetry in Densi...ned angular momentum quantum numbers L and Lz , the exact density n(r) ...
Relativistic Shock Acceleration A Hartree-Fock Approach.pdf
We compare our results with accurate numerical solutions obtained using the ...equation of state via the QJ method, 2 3 Hartree-Fock approximation and ...
...in the Integer Quantum Hall Effect within Hartree-Fock ....pdf
Received 21 August 2005, revised 27 August 2005 Published online Revision : 1.13, compiled 2/2/2008 Key words Quantum-Hall effect, Hartree-Fock, ...
...D Accuracy of the Hartree-Fock method for Wigner....pdf
quantum dots at high magnetic fields are studied by the Hartree-Fock (HF)...obtained with the broken-symmetry HF wave functions become exact in the in...
...Hartree-Fock calculations of conductance for a quantum ....pdf
of conductance for a quantum point co_专业资料。...state that can be assessed by the Hartree-Fock ...Exact broken-symmetry ... 8页 免费 Tensor ...
Tensor interaction in Hartree--Fock calculations.pdf
The uni?ed description of ground and excited nuclear states within a Hartree-Fock plus Random Phase Approximation theoretical scheme, rules out the ...
...by the Symmetry-Unrestricted Cranked Hartree-Fock Calcula_....pdf
Symmetry-Unrestricted Cranked Hartree-Fock Calcula_...asymptotic quantum numbers (?Λ = 1, ?nz = 2...Exact broken-symmetry ... 8页 免费 ...
Dirac Hartree-Fock for Finite Nuclei Employing realistic ....pdf
Dirac Hartree-Fock for Finite Nuclei Employing realistic Forces_专业资料。We ...many-body theory within the framework of non-relativistic quantum mechanics. ...
Bimodal fission in the Skyrme-Hartree-Fock approach.pdf
Bimodal fission in the Skyrme-Hartree-Fock approach_专业资料。Spontaneous-fission properties of 256Fm, 258Fm, and 260Fm isotopes are studied within the ...
An Hartree-Fock study of Potassium adsorption on TiO 2 ....pdf
An Hartree-Fock study of Potassium adsorption on TiO 2 surfaces_专业资料。...Non-irreducible atoms (those symmetry related to atoms in the assymetric ...
Vibrational Analysis of TeCl. II. A 4 HartreeFock, MP2, and ....pdf
quantum chemical methods w Hartree Fock, second-...without the inclusion of the exact 1997 John ...d The computed geometry has tetrahedral symmetry. ...
Slow long-range decay of bound Hartree-Fock orbitals and ....pdf
The existence of extra nodes in solutions of Hartree-Fock equations was ...motivated by spintronics and solid state quantum computers based on spin ...
...positronium affinities in the work-formalism Hartree-Fock ....pdf
work-formalism Hartree-Fock approximation_专业资料。...using variational and di?usion quantum Monte Carlo...Exact broken-symmetry ... 8页 免费 ...
...localized Hartree-Fock density-functional approach for the....pdf
Spin-dependent localized Hartree-Fock density-...dependent of the orbital symmetry of the state. ...iσ (r) of an electron with quantum numbers n...
Particle-drip lines from the Hartree-Fock-Bogoliubov theory ....pdf
Particle-drip lines from the Hartree-Fock-... the usual HFB solutions describe even nuclei ...the exact positions of drip lines from the ...

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

copyright ©right 2010-2021。