当前位置:首页 >> >>

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


Exact broken-symmetry states and Hartree-Fock solutions for ....pdf
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 ?...
RDX中14N核四极耦合常数的Hartree-Fock计算_电子/电路_工程科技_专业资料。用... Exact broken-symmetry ... 59人阅读 8页 免费 化学理论计算中Hartree .....
Exact broken-symmetry states and Hartree-Fock solutions for ....unkown
Exact broken-symmetry states and Hartree-Fock solutions for quantum dots at high magnetic fields B. Szafran a,b, F.M. Peeters a, S. Bednarek b, ...
Brueckner-Hartree-Fock study of circular quantum dots_图文_....pdf
Brueckner-Hartree-Fock study of circular quantum dots - We calculate ground state energies in the...... Hartree-Fock calculations are close to the exact ....
Hartree_Fock__Method_自然科学_专业资料。摘自:http...ground-state energy of quantum many-body system.... Exact broken-symmetry ... 8页 免费 Point ...
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 ...
...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 ...
...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, ...
There are No Unfilled Shells in Hartree-Fock Theory.pdf
lled in an exact Hartree-Fock calculation. Our theorem generalizes to any...states and ground state energies of quantum mechanical systems such as atoms...
Hartree-Fock ground state of the composite fermion metal.pdf
Hartree-Fock ground state of the composite fermion metal_专业资料。Within the Hartree-Fock approximation the ground state of the composite fermion metal is ...
Analytical Hartree-Fock gradients with respect to the cell ....pdf
Analytical Hartree-Fock gradients with respect to ...The traditional quantum chemical methods are dif?... without symmetry, in such a way that three ...
第44 卷第 11 期 2015 年 11 月 当代化工 Contemporary Chemical Industry Vol.44,No.11 November,2015 用 Hartree-Fock 法预测取代苯亚硒酸的 pKa 值张梦阳...
3320第六讲(2005217) Hartree, Fock, Hohenberg, Kohn, and Sham_....pdf
3320第六讲(2005217) Hartree, Fock, Hohenberg, ...(Quantum Monte Carlo) 2005.2.17 3.320 材料...Exact broken-symmetry ... 8页 免费 Hohenberg-...
A Hartree-Fock Study of Persistent Currents in Disordered ....pdf
where the Hartree-Fock ground state energy is ...remarkable agreement with the exact results of [4...Quantum Monte Carlo St... 8页 免费 ...
...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页 免费 ...
...fluid. VIII. Time-dependent Hartree-Fock approximation ....pdf
VIII. Time-dependent Hartree-Fock approximation ...The first-order solutions Comment on ”Dynamics ...Broglie - Bohm theory depend on quantum states. ...
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 ...
...on Hartree Fock Roothaan Calculations for Ground States of....pdf
Comment on Hartree Fock Roothaan Calculations for ...Truhlat, Int.J.Quantum Chem., 62(1997)1. 2...Exact broken-symmetry ... 8页 免费 Hartree-...

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

copyright ©right 2010-2021。