ALKAN PHYSICS LETTERS
. L, 3 (3), pp. 135 - 142 (1995) 1 July 1995
G�OUND STATE CORRELATIONS IN QUASI-ONE-DIMENSIONAL
ELECTRON AND CHARGED BOSE SYSTEMS
E. H. TURECi and B. TANATAR Department of Physics, Bilkent University,
Bilkent, 06533 Ankara, TURKEY.
(received 17 April 1995 ; accepted 23 May 1995)
Abstract. - Ground-state properties of quantum confined quasi-a d'
(Q1D) systems in which the charge carriers are electrons or charged b ne- 1!'1e!ls1onal
Th · b · th' k osons 1s ·
gated. e charge earners are assumed to em a zero !c ness xy-plane w'th 1nvesti-confinement potential. We assume that the charge earners are in their I 1 a harmonic and employ many-body techniques such as the Hartree-Fock approximation owest subband, approximation (RPA), and the self-consistent scheme of Singwi et al. t� ra:idorn-phase ground-state properties of this model system. Similar calculations are also c culate the thin the Hubbard approximation for a double-layer Q1D electron system in th Performed wi
limit. e no-tunneling
,'·INTRODUCTION
111.ecent advances in the fabrication technology have made possible the realizat · [ \�imensional
(Q1D)
semiconductor electron systems. Adding a lateral con�on
l]
of qua.sir.
electron gas as occurs in GaAs heterostructures, results in the quantizatio nernent to the l� except for one space dimension. Because of the reduced density of stat n of electronic�ts have been detected_ with promi�ing possibilities for device applications.esQ many novel
(�ms have been the subject of extensive theoretical study because of the enh lD electron ('ttelation effects in lower dimensions
[2, 3, 4].
Most of these investigations anced exchange i�collective modes, dielectric response, and quasi-particle properties includinc
�;centrated on
·r\cts
but relatively little attention is paid to the ground state correlation � e rnany-body!� '
S In Ql0
1�·
syste111s.. 1'he ground state properties of charged-particle systems obeying quantum 8
!� a major problem in many-body theory. Although a vast amount of lit tat1stics have
i . ( I . 1· I eratur .
, arged-ferm1ons
e.g.,
the e ectron gas), comparatively 1tt e work [6] has b e exists for,�arged bosons, mainly because they have no simple experimental realizatio een devoted to :�arged Bose gas is of interest because of its intrinsic simplicity and its relation�- However, the
��l>orta.nt ele�tron gas. Recently, the importance of disorder for the experirnen
°
the Physically. 'Perconductmg films led to the suggestion that a 20 charged Bose fluid · ts on ultra-tl ·
136
BALKAN PHYSICS LETTERS
There has bee
na lot of i
nterest in the one-
dimensional electron systems from the �a-called
Lutti
nger liqui
dpoint of view [8] which ·uses a linear energy-momentum relatior.. Her� we take
the Fermi liquid approach for the electro
ns as argue
dby Hu an
dDas Sarma (4]. Experimental
evide
nce from light scatteri
ng measureme
nts [9] support this point of view for quantum wires
as occur in semiconductor structures.
In this paper we stu
dy the grou
ndstate properties of a QID system in which the charge
carriers are either electrons or charge
dboso
ns, with the aim of applying stan<lard many-body
tech
niqu
esi
na comparative ma
nner. For the case of electro
ns, we assume the partides are all
confined to the low
est subba
nd. Begi
nni
ng from the Hartree-Fock (for electrons) and random
phase approximatio
n(RPA), we employ the self-consistent appro
ach of Singwi d al. [10)
(STLS) to inv
estigate the ground state e
nergy a
ndvarious correlation functions. STLS scheme
was applie
dto QID electro
ns usi
ng differe
nt interaction mo
dels [2, 11] an
dto charged bosons
with a contact i
nteraction (12]. As an i
nteresti
ng application of the STLS method, we also
co
nsider a double-layer QI D system, a
ndcalculate the groun
dstate energy as a function of the
I
,,
layer separatio
ndista
nce with the
no tu
nneling assumption within the Hubbar
dapproximation.
�
Correlations i
ndouble-layer electron systems have been calculated mai
nly"for 20 systems and
li
qua
ntum-wells [13, 14] with i
nteresting results. The ·aPA is k
nown to
describe the dielectric
!
properti
esof i
nteracting particles quite well i
nthe high-
de
nsity limit. It is believe
dthat RPA
\J
becomes less accurate as the dime
nsionality of the system is reduce
d. As a computatienally
ij·.,·, ..
i'•.1.
ma
nageable formulatio
n, the STLS method gives remarkably good results for the correlation
�
energy in 20. Our aim is to check the extent to which STLS approach
describes the ground
�
state properti
esof Q 1 D systems.
The r
est of this paper is orga
nized as follows. In the
next sectio
n, we
describe the model
use
dQID system, an
doutline the calculation of excha
nge-correl
ation energies of electron and
charged Bose systems. I
nsectio
n3, we present an
ddiscuss our results.
2. MODEL AND THEORY
Semico
nducti
ng QlD systems are realize
dexperimentally by several groups i11
different sha
pes. In our model system, the charge .carriers, bosons or fermions, are assume
dto occupy a zero
thick
ness zy.. plane with a harmo
nic co
nfi
neme
nt pote
ntial [3, 14] in they-direction so that the
subband e
nergies are
fn=
O(n+l/2) where n
=
0, 1, 2, ··.,
andO
describes the !:3tre
ngth of the
confini
ng pote
ntial. Furthermore, we make the assumptio
nthat the partic.1 es occupy the lowest
subban
d, i
nthe case of electrons. For a system of particles co
nfi
ned to one dimeusio
nwe define
the average i
nterparticle dista
nce a , such that the ID de
nsity of particles is given by n
=
1/2a.
As the scali
ng unit for our system, we choose the effective Bohr radius aB -= foh
2/(me
2) where
to is the backgrou
nd dielectric constant an
dm is the effective ban
dma..'-.s. The
dimension!�
de
nsity parameter is
defi
ne
dby r
11 :::;a/as
=
1r/(4kFas), where kF is the Fermi wave vector
for the corresponding noni
nteracti
ng Fermi system which is relate
dto the particle density by
n
=
2kF/'lr_· We shall use as to scale le
ngths an
deffective ry
dbergs (R
y
=
1/(2mat)) to
scale e
nergies. The effective {2oulomb i
nteraction is give
nby (14] V{q) = (e
2/to)
F(q),
where
F(q) = exp(q
2/20)/(o(q
2/20), 'in which /(o(z) is the. mo
difie
dBessel function of the first kind
an
d.q
and O are expresse
din units of as a
ndR
y
, respettively. The interaction e
nergy (p.;
particle) of the charge
dparticle system is
.. fl. TURECI et al. : GROUND STATE CORRECTIONS 137
E
1(
'°N
= ;
lo
dq F(q)[S(q) - 1],
(1)
,here S( �ptibility
q)
is the static structure factor, which in turn is expressed through the dynamicS(q)
= -
r3f
00d
w Imx(q,
w).
1r
lo
(2)
fbe dynamic susceptibility is expressed in a general way as
(
Xo(q,
w)
X q,
w)
=
1 - V(q)[l - G(q)]xo(q,w)
(3)
IJ.ere
xo(q,w)
is the susceptibility of a non-interacting system [2]. The above formalism to �ct1late the interaction energy holds irres�ecti�e of the particle st�tistics. The local-fieldfac-G(q) describes the short-range correlations m the system. Takmg_G == O amounts to the :dom-phase approximation (RPA), whereas G
=
1 results in the Hartree-Fock approxima tio . In general, a wave vector dependent G(q)
accounts for the correlations neglected by the R;A, and has origins in the m�ny-body ii:iteractions. �he abo�e relations are applicable to both fermions and bosons, provided one uses the respective non-mteracting susceptibility XAn approximation scheme to calculate
G(q),
based on a physically appealing pict ois O"ided by Singwi et al. [10] (STLS). The choice of G(q) in the STLS depends on an a:re tz !:tiich decouples tl_1e two-part�cle _dis�ribt�tion fun�tion in a p�o?uc� of two one-particle d�:iri� t,ution functions tunes the pau distribution function [5] and 1t 1s given byr
$1
00 k · q V ( k)G(q)=-
-
1r -oodk-
q V(q)
2-(S(q-k)-1].
(4)
S'nce
• G(q) dep· ends on r lS(q),
· t t Eqs. (2-4) need to be solved self-consistently. Using th H t fS
( ) . ( e ar r ee-Fi ck expression ,or t 1e static s rue ure actor, HF q , m Eq. 4) results in the H bb d 0 --imation to the local-field correction, GH(q). Different approachesto co t u ar
ap-prv.,. ns ruct G(q) in
electron and charged Bose systems are �lso pos�ible [15J. .
For charged boson systems, the non-mteractmg response function is given by ( w) _
2n(q
2/2m)
xo
q'- ( . )2 ( 2/2 )2 .
w + i7J - q mThe density response function x( q,
w)
has the same form as that of electrons. Th B" .1espression
for the excitation energy q2 /[2mS(q
)] identified with the poles 0� t�e-Feynm�nsusceptibility yields [6 ,7], · dynamic
S(q)
=
[1
+2F(q)[l -
G(q)]]-1/2raq
(5)
Similar to the electron system, G
=
0 retains the RPA for the structure factor. In tlof STLS S(q) and
G(q)
have to be solved self-consistently. le framework As a' commonly used experimental system, we now consider a double-lasystem.·
We assume that the particles are confined in two zero-thickness playerQID
electron nes, separated by1 38 BA LKAN PHYSICS UTTERS
a distance
d,
and there is no tunneling between different layers. Additional harn1onic confinement in the lateral direction renders a double-layered quantum wire system. Denoting by 1 and 2 the different layers, the interaction energy is given by
E 2 /00
N
= ;Jo
dq
{ F1 1(q)[S1 1 (q)
- 1 ] + F1 2(q)S12(q)} ,(6 )
where Fu(q) is the same as F(q) for a single-layer system, and inter-Jaye� Couiomb interaction ·
lS
The static structure factors Sij(q) arr evaluated from the dynamic response funct ions [13]
_ l ( [XO ( q , W )] - I - ef> I 1 (
q)
:
-q,
1 2 ( q) ) X (q,w) = 1 ' -<1>1 2( q) [xo( q, w )]- - </>22( q) (7 ) (8)in which ,/Ji;(q) = Vij(q) [1 - Gi; (q)]. The RPA is recovered if the local-fidd factors are set
to zero,
viz.,
Gi; (q)=
O. Formally, the above set of equations to calculate the mteraction.energy of a double-layer system is similar to those of two-component plasma, first considered
by Vashishta
et al.
(16].3� RESULTS AND DISCUSSION
Experimentally realized semiconducting quantum wires typical1y hav<' rarrin densities n
-105 - 101 cm- •. We choose a GaAs material (dielectric constant c0 � 1 3) ! o present our
numerical calculations, so that the density parameter r3 r.J 0.05 - 5. We -1 lso n .... sume that the
electrons occupy only the lowest subba11d, although generalization of our forn1,1.lism to higher sub bands is possible [5]. We use the same parameters for Ql D chargeJ b, ,�ons !or comparison. Figure 1 shows the interaction energy (per particle) for Q 1 D electro1 1 and charged boson
systems with harmonic potential parameter
n
=
5 Ry, in various appro, ima1 .. lllS. In Fig. la,we consider the electron system. Shown by the dashed line is the Hartree Fock <>r the exchange energy. Dotted and solid lines indicate the correlation energies in RPA :wd llubbard appro
ximation, respectively. Similar correlation energies for the charged hos• 1 1 1 syst ,·m is illustrated
in Fig. 1 b. We observe that for the same Coulomb potential model, ti,,· ek, tron correlation energy is lower than that for charged bosons. although the many-body , · ,rrert ions to the RPA are roughly the same for both systems.
In Fig. 2, we display the static structure fact.or S(
q)
that enters the c:, :culation of interactionenergies. RPA ( dotted lines) and self-consistent STLS approximation : solid lines) to S(q) in
a Q lD elec_tron system is depicted in Fig. 2a for 7'8 = 0.5 (lower cun · s) and 7•3
=
1 (uppercurves). Similar results of
S(q)
for a slightly different model of a qu .. ntum wire, within theSTLS scheme were reported'··by Borges el
al.
[5] . Fig. 2b shows the r . rrespondingS(q)
for oQ 1
D
charged Bose system with same parameters.In the calculation of ground-state energy of a double-layer Q I D . :(·ctron systems, rat.her than employing the full STLS scheme, we resort to the Hubbard appr. xirnat.ion for the former
E. H. TU RECI et al. : G RO U N D STATE CO RRECTIONS
is
computationally more demanding. In this case, the local-filed factors Gij are given byB
G
r3
loo
k . q �;(k) HF
]
Gii (q)
= - -
1r - � dk --y- u .. ( ) [Sii (q - k) - 6ii . q Vj3 q0 ...
.--..---r--r-,--l'""""T--,---r--r-.,...r--r-r--r--r--r-,--,- 1
(a)
0=5 Ry
� �-
m-2
�-3
-4
0.0
0.2
0.40.6
0.8
1.0
0
- 1
(b}
0=5 Ry
� �-
m-2
�-3
-4
0.0
0.2
0.40.6
0.8
1.0
raFig. I . - ( a) I nteract.ion energy per particle for a QJ D electron system in a confining potential. Solid , dotted, and dashed lines are calculated using Ga (q) ( Hubbard ) , G
=
0 ( RPA), and G=
1 ( H artree-Fock) , respectively. (b} Interaction energy for charged bo sons. Solid and dotted lines indicate Huh bard appro ximation and RPA , respectively.1 39
(9)
The Hubbard approximation takes only the exchange into account where the Pauli hole is used for electrons in the same layer,. and Coulomb correlations are not included. Since
sr{
= O, for electrons in different layers we get G12=
0, in this approximation. We displaythe
interaction energy as a function of the density parameter 1·3 , and separation distance dfor a double-layer Q 1 D electron system in Fig. 3. The density dependence of
E
/N
is showni
i�
Fig. 3a for RPA (dotted line) and Hubbard approximation (solid line). We assume the density ..
of electrons on each layer is the same, and inter-layer distance is taken as d = q8 . Similar
to the case of dou ble-layer 20 electron systems [13], magnitude of the RPA energy remains
larger than that of H ubbard approximation. This is a consequence of RPA dielectric func tion's not incorporating the short-range correlations of electrons. Since we have not attempted
140
BALKAN PHYSICS LETTERS
to calculate the full STLS energy, the local-field factors reduce the effective interaction onlyslightly, and the difference is small suggesting the need of better approxim:itions. In
Fig. 3b,
the i nteraction energy per particle obtained from the Hubbard approximatian (sohd line) and the RPA (dotted line) are shown as a function of the layer separation. Th'..· energy attains 8 constant value for large separation distances which is about twice that for two independent layers. This implies that correlations between layers are important only for smal l distances.
1 .2 ... ---�---1 .0
0.8
_.. 0.6 rn0.4
0.2
0=5 Ry
(a)
0.0 i...ll�-L....L..-L...L-L..-L-.L...1"--1...L...I-..L...1...L....L-.L...I�0
24
6 810
qaB
1 .2 ,...,__,.._.-.--.,..., __________ _ 1 .00.8
'?i.f
o .. a
0=5 Ry
0.4 0.2 i (b) 0.0 -...__._--.. ... _.__.__,__,_-'-...._"--'__.__._.&....I__._-.J0
24
6 8 1 0qaB
Fig.2. - (a) The static structure factor S(q) for elec
trons. Solid and dotted lines indicate Hubbard appro
ximation and
RPA,
respectively.(b)
S{q) for chargedbosons.
We have applied the self-consistent scheme of Singwi
et al.
[2] t.,, calr ulate the groundstate energies of Q l D electron and charged Bose systems. Corrections imroduced by local
field factors to the interaction energy are approximately the same for
:Joth
:-;ystems. Hubbardapproximation is used to investigate the density dependence of ground ::tatf' ,·nergy of a-�ouble layer Ql D electron system. We find that, to obtain reliable result�· beyc ,id the RPA, it is neces�ary t.o develop better approximations such as the full STLS scl: · · me. .\lthough we have
presented results of ground-stat.e energy calculations in Q I D system� for � :
=
5 Ry, it shouldt.
H. TURECI et al. : GRO U N D STATE CORRECTIONSz
... (a) 0=5Ry
d=aB···
I:!- 1 0
- 12 ___.__.__...._.__.__..___.__._.._.__._._.._.._.__,_._.__.__.
0.0
0.2
0.40.6
0.8
1.C
rs-6 ,_ ______ __,...,. ____ __,._. __________ __
(b) 0=5Ry
r11=1-8
- 10
---
---·
·
- 12 ._.._._...,.,..._ ... _,_...._ ... _,_..._�_.__._..__.
0.0
0.2
0.4
0.6
0.8
1 .0
Fig.3. - (a) Interaction energy per particle as a func tion of electron density for a double-layer QI D elec tron system. Solid and dotted lines are for Hubbard approximation and RPA, respectively. The layer se
paration is d
=
as. (b) Interaction energy as a function of layer separation, at r$
=
1. Solid and dottedlines are for Hubbard approximation and RPA, re spectively.
Acknowledgements. - E. H . T. gratefully acknowledges a TUBITAK scholar h'
s lp.
REFERENCES
141
[1] Petroff, P. M., Gossard, A. C. , Lagan, R. A. and Wiegmann, W. W. A I p
41 635, 1982; Goiii, A. R., Pinczuk, A., Weiner, J. S., Calleja J. 'M
PP ])h
!s. Lett.,
Pf�iffer, L. N. and West, K. W. , Phys. Rev. Lett., 67, 3298, 1991; Dern�l T
ennis., B. S.,
D. P., Grambow, P. and Ploog, K., Phys. Rev., B38, 12, 732, 1988. .. ' ·, Heitmann,.
[2] Williams, P. F. and Bloch, A. N., Phys. Rev., B 10, 1097, 1974; Friesen W
1 42 BALKAN PH YSICS LETTERs
768, 1 99 1 .
(3] Gold, A . and Ghazali, A . , Phys. Rev. , B41, 7626, 1 990; Das Sarma, S and Lai,
W . Y.,
Phys. Rev. , B32, 1401 , 1985.
[4] Hu, B. Y.-K. and Das Sarma, S., ·Phys. Rev. Lett., 68, 1750, 1 992; Phy� Rev . B48
,
546 9,1 993.
(5] Knorr, W. and Godby, R. W. , Phys. Rev. , B50, 1 779, 1994; Borges,
:\ .N . : Degani, M.
H. and Hipolito, 0., Superlatt. Microstruc., 13, 375, 1993.
[6] Fo)dy, L. L., Phys. Rev., 124, 649, 1 96 1 ; Hines, D. F. and Frankel, :\ . E. : Phys
. Rev.,
B20, 972 , 1979; Um, C. I . , Kahng, W. H . , Yim, E. S. and George, · 1'.· F . . Ph
ys. Rev.,
B41 , 259, 1 990.
[7] Gold, A . , Z. Phys. , B8·3, 429, 1 99 1 ; ibid., 87, 1 69, 1 992.
[8] Schulz, H . J . , Phys. Rev . Lett. , 71, 1 864, 1 993; Glazman , L. I., Ruzin . I. l'vl and
Shklov
skii, B. I., Phys. Rev., B45, 8454, 1 992.
[9] Strenz, R., Bqckelmann, U . , Hirler, F., Abstreiter, G., Bohm, G. and Weiman
n, G., Phys.
Rev . Lett. ,