Abstract
Properties of the âelectron gasââin which conduction electrons interact by means of Coulomb forces but ionic potentials are neglectedâchange dramatically depending on the balance between kinetic energy and Coulomb repulsion. The limits are well understood1. For very weak interactions (high density), the system behaves as a Fermi liquid, with delocalized electrons. In contrast, in the strongly interacting limit (low density), the electrons localize and order into a Wigner crystal phase. The physics at intermediate densities, however, remains a subject of fundamental research2,3,4,5,6,7,8. Here, we study the intermediate-density electron gas confined to a circular disc, where the degree of confinement can be tuned to control the density. Using accurate quantum Monte Carlo techniques9, we show that the electronâelectron correlation induced by an increase of the interaction first smoothly causes rings, and then angular modulation, without any signature of a sharp transition in this density range. This suggests that inhomogeneities in a confined system, which exist even without interactions, are significantly enhanced by correlations.
Similar content being viewed by others
Main
Quantum dots10âa nanoscale island containing a puddle of electronsâprovide a highly tunable and simple setting to study the effects of large Coulomb interaction. They introduce level quantization and quantum interference in a controlled way, and can, in principle, be made in the very-low-density regime, where correlation effects are strong11. In addition, there are natural parallels between quantum dots and other confined systems of interacting particles, such as cold atoms in traps.
Therefore, we consider a model quantum dot consisting of electrons moving in a two-dimensional (2D) plane, with kinetic energy , and interacting with each other by long-range Coulomb repulsion (
). Here electrons are labelled by i, their positions are ri, and all energies are expressed in atomic units, defined by ħ=e2/ε=m*=1 (with electronic charge e, effective mass m*, and dielectric constant ε). The electrons are confined by an external quadratic potential Vext(r)=(1/2)Ï2r2 with circular symmetry and spring constant Ï. The ratio between the strength of the Coulomb interaction and the kinetic energy is usually characterized by the interaction parameter rsâ¡(Ï n)â1/2, with n being the density of electrons. For our confined system in which n(r) varies, we define rs in the same way using the mean density
. We have studied this system up to N=20 electrons. The spring constant Ï makes the oscillator potential narrow (for large Ï) or shallow (for small Ï); it thereby tunes the average density of electrons between high and low values, thus controlling rs. For example, for N=20, varying Ï between 3 and 0.0075 changes rs from 0.4 to 17.7. The radius of the dot grows significantly as rs increases, in an approximately linear fashion (see Fig. 1).
The extrapolated quantum Monte Carlo (QMC) estimator is used9. a, High density: rsâ0.4 (Ï=3.0). b, Low density: rsâ15 (Ï=0.01). Note the dramatic change in density profile with increasing rs: the electronâelectron correlation caused by stronger interactions at low density produces sharp rings. The three-ring structure agrees with that seen in the classical limit. Note the significant increase in the radius of the dot for larger rs: it changes from approximately 2.75 to 80 (in atomic units) by increasing rs from 0.4 to 15. c, Radial cut of n(r) for rsâ10 (Ï=0.02), where the three-ring structure is about to appear. The modulation is quantitatively characterized by the fractional peak height (FPH): draw the line tangential to the two outer peaks of n(r) (), then find the vertical line
along which the distance from
to n(r) is maximum, and finally define the FPH as the ratio of the two lengths
. d, FPH as a function of rs for N=20 and 7. The curve for N=20 is linear and completely featureless for rsâ²18. The solid line is a linear fit to the data. For smaller N, radial modulation in n(r) becomes stronger, leading to FPHâ1 for large rs, and a deviation from linearity occurs (the dotted line â¼rs0.41 is the best fit for N=7). For our largest rs, the FPH typically grows with decreasing N, although not always monotonically. For example, FPH is largest for N=7, which yields a âperfect crystalâ with equidistant electrons, and thus produces a peak in the addition energy (see Fig. 3a). The Monte Carlo statistical error is less than the size of the points.
In the bulk 2D electron gas, numerical work suggests a transition from a Fermi-liquid state to a Wigner crystal around rscâ30â35 (refs 2â8). On the other hand, experiments on the 2D electron gas (which include, of course, disorder and residual effects of the ions) show more-complex behaviour, including evidence for a metalâinsulator transition5.
Circular quantum dots have been studied previously using a variety of methods, yielding a largely inconclusive scenario. Many studies12,13,14 have used density functional theory or the HartreeâFock method. These typically predict charge or spin-density-wave order even for modest rs (unless the symmetry is restored after the fact14), which are thought to be unphysical. Exact diagonalization calculations15,16 can be highly accurate but are restricted to small N and rs. Path-integral quantum Monte Carlo (PIMC) has also been applied: Egger et al. 17,18 found a crossover from Fermi liquid to âWigner moleculeâ at rsâ4âa value significantly smaller than the 2D bulk rsc. Another study19, using different criteria, found a two-stage transition for rs larger than rsc. Although PIMC treats interactions accurately, it has its own systematic and statistical problems; for instance, it generates a thermal average of states with different L and S quantum numbers, only preserving Sz symmetry. To avoid these various difficulties and to clarify the scenario, we have carried out a study using the variational Monte Carlo (VMC) and diffusion Monte Carlo (DMC) techniques9, which we used previously to study both circular20,21,22 and irregular23 dots at rsâ¼2. This method is free of the problems of PIMC but is approximate in that a âfixed-nodeâ error is made. We believe that the fixed-node error is small for the range of parameters studied here (see the Methods section, and the detailed comparison in the Supplementary Information).
Results for the electron density, n(r), are shown in Fig. 1. There is a dramatic change in n(r) on increasing interaction strength: for weak interactions (Fig. 1a), the density is rather homogeneous; the small modulation seen is caused by shell effects in the orbitals of the mean-field problem. In contrast, large rs induces strong radial modulation in n(r) (Fig. 1b), resulting in the formation of rings. Interestingly, for rs>10 the number of rings for each N is the same as that seen in the classical limit24,25 (), for example, three rings for N=20. In all of the cases we consider, the density n(r) is circularly symmetric, as is the density of spin-up and spin-down electrons separately, as required in two-dimensional systems, because we work with states of definite angular momentum L.
We find that the formation of rings, and the increase in their sharpness is completely smooth. This is shown quantitatively in Fig. 1d by using the fractional peak height (FPH, defined in Fig. 1c) of the outer ring to characterize the degree of structure. In the resulting curve for FPH as a function of rs, no deviations or special value of rs can be seen.
Having established the role of strong correlations in the formation of radial rings, which may be detectable in scanning-probe measurements, we turn to angular modulationâthe issue of correlation-induced localization of the individual electrons in each of the circular rings. We therefore consider the pair densities gÏ Ïâ²(r0;r)âthe probability of finding an electron with spin Ïâ² at location r when an electron with spin Ï is held fixed at r0âand gT=gââ+gââ. In addition to radial rings, these detect any angular structure induced by the interactions.
The most prominent feature in the pair density is a hole around the location of the fixed electron. For unlike spins, it is caused purely by Coulomb repulsion (correlation hole), whereas for like spins the antisymmetry of the wavefunction plays an important role (exchange hole). For small rs, correlation is weak, so the hole in gÏ,âÏ is much smaller than that in gÏ,Ï. As rs increases, the correlation hole grows bigger, becoming the same size as the exchange hole around rsâ4â5.
Results for the pair density in the circular quantum dot are shown in Fig. 2 for an up-electron fixed on the outer ring. For N=20 at large rs, there are clear oscillations along the angular direction near r0. This signals âincipientâ Wigner localization. However, these oscillations are weak (weaker than the radial modulation) and short ranged (damped), indicating that long-range order is not yet established. As for radial modulation, the amplitude of the angular oscillations grows continuously, without any threshold value.
a, gââ(r0;r) and b, gââ(r0;r) for the N=20 ground state (L=0, S=0) with rsâ15 (Ï=0.01, r0=(57,0)). Short-range order develops near the fixed electron, indicating âincipientâ Wigner localization but not true long-range order. c, Evolution of angular oscillations along the outer ring with rs and N. The top trace shows 0.65ÃgT for rsâ6, N=20: although strong radial modulation has already appeared, leading to âring formationâ, there is almost no angular modulation. The middle trace is gT for rsâ15, N=20: clear angular structure is present, although compared with the ring modulation it is weak and short range. Spin-resolved angular structure is also shown here; note the peculiar bump at θ=Ï. The bottom trace is gT for rsâ16, N=6 (L=0, S=0): for small N, angular modulation is clearly stronger. (The y axis is shifted and scaled for N=20 for clarity.)
The evolution of the angular oscillations as a function of rs and N is illustrated in Fig. 2c. Comparing the top two traces, for rsâ6 and 15 at N=20, we see that gT is almost featureless, even for an rs substantially larger than 1, whereas short-range oscillations have set in by our largest rs. The weakness of these oscillations suggests that electrons remain more or less delocalized along the ring for N=20, up to the largest rs studied. An intriguing feature of the spin-resolved pair densities shown is the bump at θ=Ï: gââ decreases while gââ increases. This feature is present for all rsâ¥4, and grows with increasing interaction strength; we have no explanation for it at this time. Turning now to smaller N, we find that two rings are present for N=6 at large rs: the outer ring has five electrons, whereas the remaining electron is at the centre. The lower trace in Fig. 2c shows that individual electrons are better localized for small N, a behaviour that we find holds quite generally.
Next, we turn our attention to the addition energy, Î2E(N)=EG(N+1)+EG(Nâ1)â2EG(N), where EG(N) is the ground state energy of the dot with N electrons. This is accessible experimentally as the spacing between conductance peaks in a Coulomb blockade-transport measurement, and is given by the charging energy in the simplest model of a quantum dot10. Our results for Î2E(N) (normalized by Ï) for different interaction strengths are shown in Fig. 3 (rs for fixed Ï varies slightly with N). For rsâ2, Î2E(N) is similar to that of previous studies12,20,21: non-interacting âshell effectsâ produce strong peaks for closed-shell configurations (N=6, 12, 20). At larger rs, the peaks weaken considerably, reducing mesoscopic fluctuations in Î2E. For similar rs, shell effects are more strongly affected for small N, whereas their remnant persists for large N. For comparison, we plot the addition energy in the classical limit24,25 obtained from the ground-state energies in ref. 24. The remarkable similarity to our quantum result for small N at the largest rs is strong evidence for electron localization.
a, Addition energy (normalized) as a function of N for three different Ï and for the classical limit24 . As interactions strengthen because of decreasing Ï, the mesoscopic fluctuations in Î2E become weaker. Note that this happens more readily in the small N limit. Features in the Ï=0.01 trace at small N are remarkably similar to those found in the classical limit, showing that electrons are nearly localized for small N. (The zero of the y axis is offset for clarity, and the normalization of the classical trace is arbitrary.) b,c, gT for N=9, rsâ15 (Ï=0.01), keeping an electron fixed on the outer ring (r0â(37,0)), for L=0, S=3/2 (b, the usual ground state) and L=0, S=7/2 (c). Increasing rs, particularly for small N, often brings a strongly polarized state very close in energy to the âusualâ ground state. In the case shown, the two states become essentially degenerate (E=1.464651H*) within statistical errors. The extent of Wigner localization is clearly stronger for the S=7/2 state.
Strong correlations can shuffle the energy ordering of different quantum states at fixed N. However, for Ï>0.01, the ground state remains consistent with Hundâs first rule (except for N=3). For smaller Ï, we see a tendency towards violation of this rule, primarily for small N (which, in general, are more affected by strong correlations), as in the following example. For N=9, the Hundâs rule ground state has (L,S)=(0,3/2). We find that for Ï=0.01, the highly polarized state (0,7/2) becomes degenerate with the usual ground state (within our numerical accuracy). (All other (L,S) states lie higher in energy.) Note that S=7/2 requires promotion between non-interacting shells, and so lies much higher in energy in the weakly interacting limit. At large rs, this difference is overcome by the gain in interaction energy. The pair density gT for both these N=9 states is shown in Fig. 3. The more-polarized state is clearly more localized; we find that, as expected, this is generally the case because exchange keeps the electrons apart.
The scenario that emerges here is significantly different from that for the bulk. The gradual emergence of the radial oscillations is connected to the fact that the translational symmetry is necessarily broken, and so the interactions can readily amplify existing inhomogeneities. The development of the addition-energy curves further supports this point: the structure caused by quantum interference is rapidly suppressed by the interactions leading to surprisingly smooth behaviour over a wide range of rs. Thus, strong correlations and incipient localization should be taken into account for a very broad range of interaction strength.
Methods
As a starting point, we use the KohnâSham orbitals obtained from a density functional calculation done in the local density approximation. We then carry out a VMC calculation using a trial wavefunction, ΨT, which is a linear combination of products of up- and down-spin Slater determinants of the KohnâSham orbitals multiplied by a Jastrow factor. The Jastrow factor effectively describes the dynamic correlation between the electrons coming from their mutual repulsion, whereas the near-degeneracy correlation is taken into account by having more than one determinant. We optimize the Jastrow parameters and determinant coefficients by minimizing the variance of the local energy26,27. Finally, we use fixed-node DMC calculations9,28 to project the optimized many-body wavefunction onto a better approximation of the true ground state, an approximation that has the same nodes as ΨT.
The fixed-node DMC energy is an upper bound to the true energy, and only depends on the nodes of the trial wavefunction obtained from VMC. We have calculated the energy E(N,L,S) of a circular quantum dot for each N with angular momentum L and spin Sâall good quantum numbers for our model. (Sz is also a good quantum number, and all our calculations are done for Sz=S; however, E is independent of Sz.) We investigated all possible combinations of L and S for the low-lying states, and the combination yielding the lowest DMC energy, EG, was taken as the ground state for that N. For expectation values of operators that do not commute with the hamiltonianâfor example, the density or the pair densityâwe use an extrapolated estimator9,29 (denoted FQMC for an operator F), which eliminates the inaccuracy coming from the first-order error in the trial wavefunction. FQMC is defined as 2FDMCâFVMC when FDMCâ¥FVMC, and as FDMC2/FVMC otherwise.
In the multi-determinant expansion of ΨT, we only keep Slater determinants formed from the lowest energy KohnâSham orbitals for all of the results shown here. Our study is currently limited to rsâ¤18 for technical reasons. The most serious is the failure of the VMC optimization, as many Slater determinants need to be included for stronger interactions.
For two cases, corresponding to one moderate and one large rs, we have done preliminary calculations with higher orbitals by including all determinants involving promotion of two electrons across a shell gap (10 configuration-state functions for N=20). This allows for a change in the nodes of ΨT. We find that the change in the energy, as well as the change in density and pair density, is small, although somewhat larger for greater rs. Thus, we believe that the fixed-node error in our calculations is under control.
References
Giuliani, G. & Vignale, G. Quantum Theory of the Electron Liquid (Cambridge Univ. Press, Cambridge, 2005).
Tanatar, B. & Ceperley, D. M. Ground state of the two-dimensional electron gas. Phys. Rev. B 39, 5005â5016 (1989).
Attaccalite, C., Moroni, S., Gori-Giorgi, P. & Bachelet, G. B. Correlation energy and spin polarization in the 2D electron gas. Phys. Rev. Lett. 88, 256601 (2002).
Attaccalite, C., Moroni, S., Gori-Giorgi, P. & Bachelet, G. B. Phys. Rev. Lett. 91, 109902(E) (2003).
Abrahams, E., Kravchenko, S. V. & Sarachik, M. P. Metallic behavior and related phenomena in two dimensions. Rev. Mod. Phys. 73, 251â266 (2001).
Chakravarty, S., Kivelson, S., Nayak, C. & Voelker, K. Wigner glass, spin liquids, and the metal insulator transition. Phil. Mag. B 79, 859â868 (1999).
Jamei, R., Kivelson, S. & Spivak, B. Universal aspects of Coulomb-frustrated phase separation. Phys. Rev. Lett. 94, 56805 (2005).
Waintal, X. On the quantum melting of the two-dimensional Wigner crystal. Phys. Rev. B 73, 75417 (2006).
Foulkes, W. M. C., Mitas, L., Needs, R. J. & Rajagopal, G. Quantum Monte Carlo simulations of solids. Rev. Mod. Phys. 73, 33â83 (2001).
Kouwenhoven, L. P., Austing, D. G. & Tarucha, S. Few-electron quantum dots. Rep. Prog. Phys. 64, 701â736 (2001).
Reuter, D. et al. Coulomb-interaction-induced incomplete shell filling in the hole system of InAs quantum dots. Phys. Rev. Lett. 94, 26808 (2005).
Reimann, S. M. & Manninen, M. Electronic structure of quantum dots. Rev. Mod. Phys. 74, 1283â1342 (2002).
Reusch, B. & Grabert, H. Unrestricted Hartree-Fock for quantum dots. Phys. Rev. B 68, 45309 (2003).
Yannouleas, C. & Landman, U. Unified description of floppy and rigid rotating Wigner molecules formed in quantum dots. Phys. Rev. B 69, 113306 (2004).
Reimann, S. M., Koskinen, M. & Manninen, M. Formation of Wigner molecules in small quantum dots. Phys. Rev. B 62, 8108â8113 (2000).
Rontani, M., Cavazzoni, C., Bellucci, D. & Goldoni, G. Full configuration interaction approach to the few-electron problem in artificial atoms.. J. Chem. Phys. 124, 124102 (2006).
Egger, R., Häusler, W., Mak, C. H. & Grabert, H. Crossover from Fermi liquid to Wigner molecule behavior in quantum dots. Phys. Rev. Lett. 82, 3320â3323 (1999).
Egger, R., Häusler, W., Mak, C. H. & Grabert, H. Phys. Rev. Lett. 83, 462(E) (1999).
Filinov, A. V., Bonitz, M. & Lezovik, Yu. E. Wigner crystallization in mesoscopic 2D electron systems. Phys. Rev. Lett. 86, 3851â3854 (2001).
Pederiva, F., Umrigar, C. J. & Lipparini, E. Diffusion Monte Carlo study of circular quantum dots. Phys. Rev. B 62, 8120â8125 (2000).
Pederiva, F., Umrigar, C. J. & Lipparini, E. Phys. Rev. B 68, 89901(E) (2003).
Güçlü, A. D. & Umrigar, C. J. Maximum-density droplet to lower-density droplet transition in quantum dots. Phys. Rev. B 72, 45309 (2005).
Ghosal, A., Umrigar, C. J., Jiang, H., Ullmo, D. & Baranger, H. U. Interaction effects in the mesoscopic regime: A quantum Monte Carlo study of irregular quantum dots. Phys. Rev. B 71, 241306(R) (2005).
Bedanov, V. M. & Peeters, F. M. Ordering and phase transitions of charged particles in a classical finite two-dimensional system. Phys. Rev. B 49, 2667â2676 (1994).
Koulakov, A. A. & Shklovskii, B. I. Charging spectrum and configurations of a Wigner crystal island. Phys. Rev. B 57, 2352â2367 (1998).
Umrigar, C. J., Wilson, K. G. & Wilkins, J. W. Optimized trial wave functions for quantum Monte Carlo calculations. Phys. Rev. Lett. 60, 1719â1722 (1988).
Umrigar, C. J. in Quantum Monte Carlo Methods in Physics and Chemistry (eds Nightingale, M. P. & Umrigar, C. J.) 129â160 (Kluwer, Dordrecht, 1999).
Umrigar, C. J., Nightingale, M. P. & Runge, K. J. A diffusion Monte Carlo algorithm with very small time-step errors. J. Chem. Phys. 99, 2865â2890 (1993).
Liu, K. S., Kalos, M. H. & Chester, G. V. Quantum hard spheres in a channel. Phys. Rev. A 10, 303â308 (1974).
Acknowledgements
This work was supported in part by the NSF (grants DMR-0506953 and DMR-0205328). A.G. was supported in part by the funds from the David Saxon chair at UCLA.
Author information
Authors and Affiliations
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Supplementary information
Supplementary Information
Supplementary Material (PDF 106 kb)
Rights and permissions
About this article
Cite this article
Ghosal, A., Güçlü, A., Umrigar, C. et al. Correlation-induced inhomogeneity in circular quantum dots. Nature Phys 2, 336â340 (2006). https://doi.org/10.1038/nphys293
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/nphys293
This article is cited by
-
Recent progresses of quantum confinement in graphene quantum dots
Frontiers of Physics (2022)
-
Multielectron Bubbles in Liquid Helium
Journal of Low Temperature Physics (2020)
-
Coherent cyclotron motion beyond Kohnâs theorem
Nature Physics (2016)
-
A molecular state of correlated electrons in a quantum dot
Nature Physics (2008)