magnetic - doped quantum dots@xcite have attracted considerable theoretical and experimental interests over the last two decades .
diluted magnetic ii - vi and iii - v semiconductor ( dms ) quantum dots were fabricated with single - electron control@xcite .
a rich variety of different magnetic and optical properties were discovered@xcite .
in such structures one can explore the physical properties coming from inter - carrier interactions and the interaction of the carriers with the magnetic ion .
this system promises to be relevant for future quantum computing devices , where for instance the spin of the magnetic ion is used as a quantum bit .
more recently , electrically active devices were fabricated in which a single manganese ion is inserted into a single quantum dot @xcite with control of the amount of charge in the quantum dot and consequently the possibility of the control of magnetism of single mn - doped quantum dots .
investigation of the exact electronic structure of a two - dimensional quantum dot confined by a parabolic potential containing several electrons and a single magnetic impurity ( in this paper @xmath1 ) in the presence of an applied magnetic field is a new topic . in a recent investigation@xcite
, a three - dimensional ( 3d ) cd(mn)se quantum dot containing several electrons , where only the low - energy levels of the single - electron problem were taken into account , was investigated in the presence of a magnetic field . here
we will extend this work and include all relevant energy levels in order to obtain a convergent solution for the ground state ( and also excited spectrum ) of the system .
it is known that in the absence of a magnetic ion , an external magnetic field is able to change the spin polarized state of weakly interacting electrons @xmath2 in a quantum dot in such a way that in the ground state it maximizes the total spin of the system : i.e. @xmath3 .
if the inter - particle interaction is strong , even without an applied magnetic field the electrons may already be polarized . however , with increasing magnetic field and in the case the inter - particle interaction is strong , the total spin of the system can be unusually reduced by the magnetic field@xcite .    in the present study ,
we investigate theoretically the few - electron two - dimensional confined quantum dot system that contains a single magnetic ion in the presence of an external magnetic field taking into account a sufficient large number of single - particle orbitals such that numerical
exact " results are obtained .
we explore how sensitive the whole system is to the position of the magnetic ion in the quantum dot and to the presence of a magnetic field and investigate the competition between the following three energies : i ) the interaction of the magnetic ion with the electrons , ii ) the interaction of the magnetic ion with the magnetic field , and iii ) the interaction of the external field with the electrons .
these terms affect the spin polarization of the @xmath2 electrons in the quantum dot .
explicit studies of a @xmath2-correlated - electron system interacting with a single magnetic ion in nonzero magnetic field are very rare in the literature .
recent theoretical@xcite and experimental work@xcite has focused either on a small number of electrons using the exact diagonalization approach at zero field ( refs.@xcite ) for a @xmath4d quantum dot or at nonzero field@xcite including only the lowest single - particle states for a 3d system or the exciton states relevant for optical spectroscopy of self - assembled magnetic - doped quantum dots@xcite .    here
, we will examine thoroughly the exact properties of the system containing several correlated electrons and a single magnetic impurity in the presence of a magnetic field . in our numerical
exact " diagonalization approach we include an arbitrary number of single - particle states to guarantee the accuracy of our results .
we investigate the influence of the strength of the inter - particle interaction and the position of the magnetic ion on the ground state of the system .
we predict the interesting phenomenon that the magnetic ion ferromagnetically couples with the electrons in a region below a critical magnetic field and antiferromagnetically with the electrons above this critical field .
thermodynamic properties as magnetization , susceptibility , and the heat capacity are investigated as function of magnetic field and temperature .    ) , and the coulomb interaction strength @xmath5 .
the inset is the first derivative of @xmath6 with respect to @xmath7 that highlights the non - linearity of @xmath6 in certain regions of the magnetic field.,width=347 ]    this paper is organized as follows .
section ii introduces the model and the numerical method . in section iii
, we present our numerical results for the many - particle ground state and investigate correlations through the appearance of vortices in the many - electron wave function .
iv addresses the many - particle spectrum and in sect .
v we present results for different thermodynamic quantities .
our discussion and conclusions are presented in sect .
a quantum dot containing @xmath2 electrons with spins @xmath8 confined by a parabolic potential and interacting with a single magnetic ion ( @xmath1 ) with spin @xmath9 and a magnetic field is described by the following hamiltonian :    @xmath10 + \sum_{i <   j=1}^{n_e}\frac{e^2}{4\pi\epsilon_0\epsilon|\overrightarrow{r_i}-\overrightarrow{r_j}| }   \nonumber\\+   \frac{1}{2}\hbar\omega_c\left(g_e m^ *   { s}_z+g_{mn } m^ *   m_z\right )   -j_c\sum_{i=1}^{n_e}\overrightarrow{s_i}\cdot   \overrightarrow{m } \delta(\overrightarrow{r_i}-\overrightarrow{r}).\end{aligned}\ ] ]    the vector potential @xmath11 is taken in the symmetric gauge : @xmath12 where the magnetic field @xmath13 points perpendicular to the plane of the interface .
the confinement frequency @xmath14 is related to the confinement length by : @xmath15 . @xmath16 and @xmath17 are the land g - factor of the host semiconductor and the magnetic ion , respectively . the dimensionless coulomb strength is defined as : @xmath18 with @xmath19 the effective bohr radius .
the cyclotron frequency is : @xmath20 .
@xmath21 , @xmath22 , and @xmath23 are the projections of the total angular momentum of the electrons , their spins , and the magnetic ion in the direction of the magnetic field .
the electrons and the magnetic impurity in the quantum dot interact via the contact exchange interaction with strength @xmath24 .
we use the set of parameters@xcite that is applicable to cd(mn)te which is a ii(mn)vi quantum dot with typical lateral size of about tens of nanometers .
the dielectric constant @xmath25 , effective mass @xmath26 , @xmath27 @xmath28 , @xmath29 , @xmath30 , @xmath31 , and @xmath32 about tens of nanometers ( @xmath33 corresponding to tens of @xmath34 ) .
for example , with @xmath35 @xmath34 gives @xmath36  .
we rewrite the hamiltonian in second - quantized form :    @xmath37,\end{aligned}\ ] ]    where the first term is the single - particle energies @xmath38 for an electron in state @xmath39 with spin @xmath40 and the second term is the coulomb interaction .
the third term is the electron and magnetic ion zeeman energy .
the last sum is the electron-@xmath41 interaction in which the first term describes the difference between the number of spin up and spin down electrons , and the last two terms describe the energy gained by flipping the electron spin along side with flipping the spin of the magnetic ion .
@xmath42 , and @xmath43 are the z - projection , raising and lowering operators , respectively , of the magnetic ion spin ( we consider @xmath44 ions which have a spin of size @xmath45=5/2 ) .
the single - particle states in a parabolic confinement potential define a complete basis of fock - darwin orbitals @xmath46 and spin functions @xmath47 : @xmath48 with the fock - darwin orbitals : @xmath49 in hamiltonian ( [ e : secondquantized ] ) , @xmath39 denotes a set of quantum numbers @xmath50 with @xmath51 the radial and azimuthal quantum numbers , respectively .
the effective length @xmath52 in the presence of a magnetic field is defined through the hybrid frequency @xmath53 where @xmath54{1+(\omega_c/2\omega_0)^2}$ ] .
the single - particle orbital energy is given by : @xmath55 the interaction parameters between the electrons and the magnetic ion in the quantum dot is expressed : @xmath56 as a product of two fock - darwin orbitals calculated at the position of the magnetic ion .
we construct the many - particle wave function following the configuration interaction ( ci ) method .
@xmath57 where @xmath58 is the @xmath59-th state of the non - interacting many electron wave function determined by @xmath2 electrons with @xmath2 different sets of quantum numbers @xmath60 and the single scatterer with one of the six states of the @xmath44 .
@xmath61=@xmath62 stands for the coordinates and spin of a single electron . in second quantization representation , the state @xmath58 , which is a slater determinant composed of single - electron states ,
can be translated into a ket vector @xmath63 grouping a total of @xmath2 electrons into @xmath64 electrons with z - component of spin up and @xmath65 electrons with z - component of spin down : @xmath66 here @xmath67 and @xmath68 are the indices of the single - electron states for which each index is a set of two quantum numbers ( radial and azimuthal quantum numbers ) , as mentioned above .    ) and @xmath69 .
the upper inset is the first derivative of the main plot with respect to @xmath7 to visualize more clearly the cusps .
the lower inset is a zoom of the rectangular regions of the main plot .
( b ) the z - projection of the total spin of the magnetic impurity ( solid black curve ) , of the electrons ( dashed red curve ) , and the total angular momentum ( dash - dotted green curve).,width=340 ]    the ci method , which is in principle exact if a sufficient number of states are included , is limited to a small number of electrons due to computational limitations . for a larger number of electrons and/or magnetic ions , other approaches that e.g.
are based on spin - density functional theory ( sdft ) using e.g. the local spin density approximation ( lsda ) as was used in ref.@xcite is able to handle a large number of electrons and/or magnetic ions .
the lsda is exact only in the case of the homogeneous electron gas , and in practice , works well also in most inhomogeneous systems .
however , in really highly correlated few - particle systems as discussed in this paper , the lsda might fail or be at least less accurate .     and the coulomb strength @xmath69 . the upper inset is a zoom of the small magnetic field region and the lower one presents the average of @xmath23 , @xmath22 , and @xmath21 for the ground state at very small magnetic fields.,width=347 ]
we first explore the magnetic field dependence of the total zeeman energy : @xmath70 which is the difference in energy in the presence and without a magnetic ion .
it consists of three terms : the difference of the zeeman energy of the electrons @xmath71 , the zeeman energy describing the interaction of the magnetic ion having spin @xmath72 with the magnetic field , @xmath73 , and the exchange interaction of the ion with the electrons , @xmath74 .
@xmath75 is just the so called  local zeeman splitting term " as discussed in ref.@xcite .
this sum is basically the difference in the zeeman energy of the electrons between the cases with and without a magnetic ion plus the energy contribution of the magnetic ion .    , @xmath76 , and @xmath77 for the coulomb strength @xmath69 .
the horizontal line at @xmath78 separates the @xmath75 plane into ferromagnetic ( fm ) and antiferromagnetic ( afm ) regions .
, width=340 ]    for @xmath79 , we find a total zeeman energy that appears linear in magnetic field . a similar linear behavior is also found in the cases with @xmath80 but with different slopes ( see fig .  [ localzeeman1234e ] ) .
let us suppose , in hamiltonian ( [ e : hamiltonian ] ) , that the contribution from the last term ( the local zeeman energy or the exchange interaction term ) is zero .
for instance this is the case when a magnetic ion is located at the center of the quantum dot having three - electrons in the partially filled @xmath81 shell , in which the first two - electrons fully fill the @xmath82 shell and the remaining one is in either of the orbitals of the @xmath81 shell
. then a perfect linear behavior of the total zeeman energy is found .
a closer look to @xmath6 gives us a slightly different picture , as is provided by taking the derivative ( see the inset of fig .
[ localzeeman1234e ] ) .
notice that the total zeeman term @xmath6 has pronounced cusps and the number and positions of these cusps is different for different number of electrons @xmath2 .
there exists one at @xmath83 for a two - electron quantum dot , one at @xmath84 for a three - electron quantum dot , and two at @xmath85 and @xmath86 for the four - electron quantum dot , with the magnetic ion located at @xmath87 .
the three - electron system exhibits a much richer behavior when we increase the coulomb interaction strength to @xmath69 as seen in fig .
[ zeeman3e0.5l0andlzmzsz](a ) where we placed the magnetic ion at @xmath76 .
cusps , which are highlighted in the two insets in fig .
[ zeeman3e0.5l0andlzmzsz](a ) , appear when the total angular momentum and/or the total z - projection of the spin of the electrons change abruptly with magnetic field .
notice that the total zeeman energy of a two - electron quantum dot in the presence of the magnetic ion does not produce a similar behavior due to the fact that the z - projection of the total spin is zero making the main contribution ( from the zeeman spin term of the @xmath44-impurity ) negligible .    .
the inset in the upper plot are the averages of the three quantities : @xmath23 , @xmath22 , and @xmath21 with the ion located at @xmath76 at small magnetic fields .
the horizontal line separates the plane into ferromagnetic ( fm ) and antiferromagnetic ( afm ) regions in both two plots .
the schematic diagram in ( b ) explains why the exchange energy is almost zero in the case @xmath88 .
, width=340 ]    the coulomb strength and the position of the magnetic impurity affects the total zeeman energy and influences the number and the position of the cusps .
the first pronounced cusp appears at lower magnetic field for larger coulomb interaction strength .
this is a consequence of the competition between the coulomb energy and the energy gap of the single - particle problem .
larger coulomb strength ( smaller energy gap ) leads to stronger electron - electron correlation and consequently the electrons are more clearly separated from each other .
it results into a high probability for finding the electrons to occupy higher energy states .
that also means that the system transfers to a configuration with larger @xmath22 and @xmath21 at smaller applied field .    .
the horizontal line separates the plane into ferromagnetic ( fm ) and antiferromagnetic ( afm ) regions .
the inset shows the total @xmath89-projection of the @xmath23 , @xmath22 , and @xmath21 of the electron system with the magnetic ion located at the center of the dot.,width=340 ]    .
the horizontal line separates the plane into ferromagnetic ( fm ) and antiferromagnetic ( afm ) regions .
the inset shows @xmath22 alone vs. magnetic field.,width=340 ]    we also found that the ground - state energy is sensitive to the presence of the magnetic field . in zero magnetic field
, the ground state receives contributions from many different configurations having z - projection of the magnetic ion , @xmath23 , from @xmath90 to @xmath91 .
when a magnetic field is applied , the ground state favors states with projection of the spin of the magnetic ion down and the states with @xmath92 give the main contribution to the ground state .    , @xmath93 , and @xmath4 at @xmath94 .
the horizontal dotted line at @xmath78 separates the plane into ferromagnetic ( fm ) and antiferromagnetic ( afm ) regions.,width=347 ]    fig .
[ zeeman3e0.5l0andlzmzsz](b ) shows the average of the three quantities @xmath23 , @xmath22 , and @xmath21 of the three - electron quantum dot as a function of the magnetic field for a coulomb interaction strength @xmath69 .
we realize that with increasing @xmath95 the @xmath96 ( @xmath97 ) exhibits jumps at smaller critical @xmath7 ( compare green dash - dot - dot curve in fig .
[ localzeeman1234e ] ) .
now we direct our attention to the very small magnetic field behavior .
there exists a very small region of the magnetic field where the total spin of electrons and the total spin of the magnetic ion are oriented parallel , we found this earlier in ref.@xcite for the zero magnetic field case .
these results are now extended to nonzero magnetic field .
this is made more clear in the upper inset of fig .
[ twoterms_3e_0.5l0 ] where the crossing point of the two terms : the zeeman energy of the magnetic ion ( @xmath98 ) and the exchange interaction ( @xmath99 ) , occurs at @xmath100 ( converted to @xmath101 for the considered system ) .
this ferromagnetic coupling extends further , up to @xmath102 ( see the lower inset of fig .
[ twoterms_3e_0.5l0 ] ) .
( black squares referred to ( 1 ) ) and ( @xmath103 ) ( red circles referred to ( 2)).,width=347 ]    it is worth noticing that this ferromagnetic coupling is extended to a much larger magnetic field range ( up to @xmath104 ) if we move the magnetic ion to the center of the quantum dot ( see fig .  [ exchange_3e_pos ] ) . this can be understood as follows .
when the magnetic ion is located at the center of the dot and the magnetic field is very small the absolute value of @xmath98 always dominates over @xmath99 .
this is opposite to the case when the magnetic ion is located at @xmath76 .
recall that in ref.@xcite we found for zero field that the exchange zeeman energy is minimum when the magnetic ion is at the center of the quantum dot and approximately zero at positions very far from the center of the quantum dot .
( black squares referred to ( 1 ) ) and @xmath105 ( red circles referred to ( 2 ) ) . ,
width=347 ]    fig .
[ exchange_3e_pos ] tells us that the magnetic field where the antiferromagnetic coupling between electron and the magnetic ion starts depends on the position of the magnetic ion in the quantum dot .
the system with the magnetic impurity located at @xmath77 exhibits an antiferromagnetic coupling for @xmath106 , that is larger than @xmath107 in case of @xmath76 .
we have discussed the appearance of antiferromagnetism in a three - electron quantum dot .
now we go back to the two simpler cases with the number of electrons @xmath108 ( see fig .  [ exc1e2e ] ) .
let us first discuss the results for @xmath79 as given in fig .
[ exc1e2e](a ) .
the antiferromagnetic coupling between the electron and the magnetic ion starts at smaller magnetic field as the magnetic ion is moved .
this is different from the previous results for @xmath109 .
the reason is as follows : for the quantum dot with a single electron , the electron tends to accommodate permanently the @xmath82-shell with @xmath110 ( see the inset of fig .
[ exc1e2e](a ) ) in the ground state while the exchange parameter in the @xmath82 shell ( @xmath111 ) is found to be maximum right at the center of the quantum dot . moving the magnetic ion away from the center of the dot
, this @xmath111 is found to be smaller and as a consequence the exchange electron-@xmath44 interaction becomes smaller than the electron zeeman energy , leading to an antiferromagnetic coupling at smaller magnetic field .
the story for @xmath88 electrons ( see fig .
[ exc1e2e](b ) ) is now interesting since the two electrons accommodate the @xmath82 shell with spins antiparallel making the total spin of the electron zero in the ground state with almost unit probability .
this leads to zero contribution to the first term written in the last line in [ [ e : secondquantized ] ] for diagonal elements .
therefore , the main contribution ( even very small ) to the exchange energy is now expected to come from coupling with configurations where one of the electrons ( spin down ) stays in the @xmath82 level and the other occupies higher level ( see the schematic diagram in fig .  [ exc1e2e](b ) ) . in this diagram , the magnetic ion is assumed to be located at the center of the dot with spin down ( @xmath90 ) .
the coupling of the electron ( spin up ) in the @xmath82 orbital with an electron from either of the @xmath81 shell is zero .
the only non - zero coupling is with an electron with the quantum numbers ( 1,0 ) of the @xmath112 shell ( as shown in the diagram and this quantum state would change if the ion is located away from the center of the dot ) with the amount of about @xmath113 .
this picture remains valid until the magnetic field is high enough to excite one electron from the @xmath82 shell to a higher quantum state forming the ground state with two up spins antiferromagnetically coupling with the magnetic ion .
for smaller coulomb interaction strength , the antiferromagnetic behavior occurs at larger magnetic field since the two electrons repel each other less and consequently they stay longer antiparallel in the @xmath82 shell .     calculated for a three - electron quantum dot with the magnetic ion ( white cross ) located at @xmath76 in ( a ) , ( b ) and @xmath114 in ( c ) when the coulomb strength @xmath69 in a magnetic field @xmath115 ( see ( a ) ) and @xmath116 ( see ( b ) and ( c ) ) .
( d ) is the magnetic ion position dependence of @xmath74 with the magenta circles indicating the respective calculated positions of the magnetic ion in ( a ) , ( b ) , and ( c ) .
red color region is for highest value of radial density.,width=340 ]    from the @xmath88 result we may ask a question : whether a four - electron quantum dot has similar properties since the number of electrons in both cases are even and one may have a situation where the total spin of the electrons is zero . indeed , if the magnetic ion is located at the center of the quantum dot even though for @xmath117 the outer shell is half - filled , this is possible as illustrated in fig .
[ exc4e ] where the antiferromagnetic coupling occurs at @xmath118 at which the total spin of the electrons reaches the maximum value @xmath119 . in this case
, the first two electrons will occupy the @xmath82 shell and the remaining two will occupy two of the five orbitals of the @xmath81 and @xmath112 shells .
this picture holds at small magnetic field . however , there is a big difference in the exchange energy as compared to the previous case of @xmath88 when the ion is shifted away from the center of the dot , e.g. in this plot at @xmath76 .
the exchange energy is much larger than the result obtained for @xmath120 because when the ion is out of the center of the quantum dot , the two remaining electrons at higher orbitals have a non - zero contribution in the diagonal exchange elements dominating the exchange energy of the ground state .
this is the reason why the antiferromagnetic transition occurs at smaller magnetic field ( @xmath85 ) as compared to the case when the ion is at the center of the dot ( @xmath118 ) although the pictures of the @xmath23 and @xmath21 transition in these cases are similar .    to complete the picture for few - electron quantum dot system
, we will discuss the afm behavior for the system with the highest number of electrons , @xmath121 , where we were able to obtain accurate numerical results .
we focus on the small magnetic field region ( see fig .  [ exc5e ] ) .
for the magnetic ion in the center of the dot , the fm coupling is dominant in the shown magnetic field region because the diagonal exchange matrix elements dominate over the zeeman energies of the electrons and of the magnetic ion .
this is different for the cases with the magnetic ion displaced a bit from the center of the quantum dot .
the fm - afm occurs at @xmath122 and @xmath123 for @xmath124 and @xmath125 , respectively .
it is similar to the cases for the system with @xmath126 due to the zero coupling between the orbitals from the @xmath127 shell with the @xmath82 orbital . to observe the afm behavior for the system with the magnetic ion located at the center of the dot where the diagonal exchange elements are almost zero , it is crucial to include enough quantum orbitals ( that rapidly increases the size of the hamiltonian matrix resulting in very time consuming calculations ) so that one allows the electrons to jump to higher energy levels and having parallel spins as previously shown for the case @xmath117 ( see fig .
[ exc4e ] ) . in that case
, the four - electron system exhibits an anti - ferromagnetic coupling with the magnetic ion at the magnetic field where the total @xmath89-projection of the spin is maximum @xmath128 .
the system is strongly polarized .
for the case @xmath121 , up to @xmath129 , the total @xmath130 and the total @xmath131 .
the inset in fig .
[ exc5e ] supports the afm behavior for the out of center @xmath44 as obtained in the main plot .
now we change the coulomb interaction strength and explore the magnetic behavior as function of the position of the magnetic ion . from fig .
[ exchangepos_3e_severallamb ] , it is clear that when reducing the coulomb interaction the system undergoes a ferromagnetic to antiferromagnetic transition at gradually larger magnetic fields for @xmath1 positions that are closer to the center of the dot .
we see that @xmath132 has a peak structure with a maximum at some specific position of the magnetic ion , e.g. see the peak for the case @xmath5 ( the blue full circles ) .
however , it is certain that at high magnetic field , the system is always antiferromagnetic .    the fm - afm phase diagram for a three - electron quantum dot in ( @xmath133 ) space
is shown in fig .
[ pdlamb - omegac ] for two different positions of the magnetic ion .
when the magnetic ion is in the center of the quantum dot ( black curve with squares in fig .
[ pdlamb - omegac ] ) the critical magnetic field increases as the coulomb interaction strength decreases .
the reason is that increasing the coulomb interaction helps the electrons to approach closer the magnetic ion and therefore the critical magnetic field for the system to transit to the antiferromagnetic phase decreases .
calculated for the quantum dot of fig .
[ denlambda1 ] where one electron ( orange arrow ) is pinned at one of the three most - probable positions of the electrons as obtained in fig .
[ denlambda1 ] .
it is at @xmath134 in ( a ) for fig .
[ denlambda1](a ) , @xmath135 in ( b ) for fig .
[ denlambda1](b ) , and @xmath136 in ( c ) for fig .
[ denlambda1](c ) , respectively .
( b ) the same correlation function for the case the position of the fixed electron is closer to the ion ( at @xmath137 ) as compared to ( b ) .
the position of the magnetic ion is indicated by the white cross.,width=294 ]    now we move the magnetic ion away from the center of the quantum dot and we obtain the phase diagram as shown by the red curve ( with solid circles ) in fig .
[ pdlamb - omegac ] . for @xmath138 , the stability of the fm phase with respect to an applied magnetic field
is strongly reduced and a small magnetic field turns the three - electron system into the afm phase . notice
that for sufficient strong electron - electron interaction ( i.e. @xmath139 ) we obtain practically the same fm - afm phase diagram as for the case the @xmath44-ion is located in the center of the quantum dot .
a remarkable re - entrant behavior is found in the region @xmath140 and @xmath141 where with increasing @xmath95 we go from an antiferromagnetic to ferromagnetic and back to antiferromagnetic phase .
this unusual behavior is understood as follows .
as the impurity is moved away from the center of the quantum dot the exchange matrix will have many nonzero off - diagonal terms that leads to a smaller fm - afm critical transition magnetic field .
now let us turn our attention to the region @xmath138 .
for very small coulomb interaction strength the electrons will repel each other only weakly and are therefore pulled towards the magnetic ion ( the nonzero exchange matrix elements increase strongly ) resulting in a very small fm - afm magnetic field . for @xmath139 ,
the electrons become more strongly correlated and the critical field stays about @xmath142 from the result for @xmath143 .
if one moves the ion further and further away from the center , the @xmath144 transition line moves to larger @xmath95 values .
for example , for @xmath69 and the magnetic ion located at @xmath76 the fm - afm critical transition occurs at @xmath145 , which is much smaller than @xmath146 found for @xmath147 .     and the coulomb interaction strength is @xmath5.,width=359 ]
the dependence of the ferromagnetic - antiferromagnetic transition of a three - electron quantum dot system on the position of the magnetic ion is summarized in the phase diagram shown in fig .
[ pdlamb - rmn ] for two different magnetic fields @xmath148 .
we can predict that with slightly larger ( smaller ) magnetic field the slope of the curve will be larger ( smaller ) . from fig .
[ pdlamb - omegac ] , we already learned that the fm - afm transition magnetic field is largest for the ion in the center of the quantum dot as is also seen in fig .
[ pdlamb - rmn ] .
the re - entrant behavior of the afm phase as function of @xmath95 is found for small @xmath149 values , i.e. when the @xmath44-ion is not too far from the center of the quantum dot , in case the magnetic field is not too small . the critical point ( @xmath150)=(@xmath151 ) for @xmath105 moves down ( up ) with increasing ( decreasing ) magnetic field .      in high magnetic field
, the magnetic ion tends to attract electrons because they are oppositely polarized . because the exchange interaction is small as compared to the coulomb interaction ,
the electrons and magnetic ion are arranged in such a way that the electrons repel each other and also try to be as close to the magnetic ion as possible .
this picture holds above the fm - afm critical magnetic field .    to show this behavior explicitly , we studied the radial density and the radial pair correlation functions .
their respective operators are defined as : @xmath152 and @xmath153 we plot in fig .
[ denlambda1 ] the radial density of a three electron quantum dot that is polarized in high magnetic field for the case that the coulomb strength is @xmath69 and the magnetic ion is located at two different positions for two magnetic fields .
the electrons and the magnetic ion are antiferromagnetically coupled . the strength of that coupling can be seen from fig .  [ denlambda1](d ) in which we plot the magnetic ion s position dependence of the exchange energy at two magnetic fields @xmath154 and @xmath155 .
those magnetic fields are typical in the sense that the exchange term is found to be very large ( @xmath115 ) or the correlation between the electrons very high ( @xmath116 ) .
density plots are shown for @xmath156 at @xmath76 and ( @xmath157 ) .
we observe three distinct peaks of maximum probability .
they are found at : @xmath134 , @xmath158 , and @xmath159 in fig .
[ denlambda1](a ) ; @xmath135 , @xmath160 , and @xmath161 in fig .
[ denlambda1](b ) ; @xmath136 , @xmath162 , and @xmath163 in fig .
[ denlambda1](c ) .
these figures show clearly the interplay effect where the three electrons on the one hand try to be close to the magnetic ion and on the other hand repel each other via the coulomb potential energy .
it results in the merging of the radial density such that the higher the exchange energy the larger the merging of the local maxima in the electron density and the smaller the correlations .
[ denlambda1](d ) gives an idea about the variation of @xmath74 with the position of the magnetic ion and reaches a maximum at @xmath76 for @xmath116 . in fig .
[ denlambda1](c ) , the three electrons are less attracted to the magnetic ion via the antiferromagnetic coupling as compared to that in fig .  [ denlambda1](b ) .
this is due to the fact that the @xmath74 for the case shown in fig .
[ denlambda1](b ) is larger than that in fig .
[ denlambda1](c ) .
the electrons are therefore found more correlated in the latter case presented by the extended red region in fig .
[ denlambda1](c ) . thereby , correlation between electrons in fig .
[ denlambda1](c ) is expected to be the highest and in fig .
[ denlambda1](a ) the smallest .
the position of the magnetic ion affects the ground - state property as is made clear in fig .
[ paircorrpos ] .
we fix the spin state and the position of one electron ( indicated by the orange arrow ) and the position of the magnetic ion ( the white cross ) .
the magnetic field is such that @xmath115 in fig .
[ paircorrpos](a ) and @xmath116 in the others .
it also reflects the fact that the system in fig .  [
denlambda1](a ) exhibits the smallest correlation as compared to the other two .
this illustrates the point raised above about the density . at the magnetic field @xmath116 ,
the electrons are strongly polarized resulting in the red regions of the up - up spin pair correlation function that tends to surround the magnetic ion .
we see that the three electrons are most likely to localize around some specific positions defining a triangle with the three electrons at the three vertices , while they are attracted to the magnetic ion .
when we locate one electron at a position closer to the magnetic ion , see fig .
[ paircorrpos](b ) , the two peaks decrease in amplitude as compared to those in fig .
[ paircorrpos](b ) .
the addition energy ( often called the chemical potential ) is defined as the increase of the energy of the quantum dot system when an electron is added : @xmath164 .
this quantity can be measured experimentally and is plotted in fig .
[ chimicalpot ] as function of the magnetic field .     with the magnetic ion ( the white cross ) located at @xmath76 in ( a),(b ) and at @xmath114 in ( c ) for coulomb strength @xmath69 .
two fixed electrons ( indicated by the two magenta circles ) are located at the two peaks appearing in the radial density : @xmath135 and @xmath160 in ( a ) , ( c ) ; and @xmath135 and @xmath161 in ( b ) .
, width=347 ]    there are several cusps appearing in the addition energy curves as a consequence of changes in the ground state .
these changes are due to variations in the z - projection of the total spin of the electrons and/or the z - projection of the total angular momentum of the system when the magnetic field increases beyond some specific values .
the presence of the magnetic ion leads to more cusps and the position of these cusps is also influenced by the number of electrons and the position of the magnetic ion .
the cusps are from either of the two systems in the study .
for instance , the green triangles in fig .
[ chimicalpot ] are for @xmath165 has two cusps at @xmath83 and @xmath84 .
the cusp at the point @xmath83 comes from the change in the configuration of the average of the total z - projection spin and the total z - projection of angular momentum @xmath166 of the two - electrons in the quantum dot from ( 0,0 ) to ( 1,1 ) . while the other cusp @xmath84 comes from the change of the phase of the three - electron quantum dot from @xmath167 to @xmath168 .
it is similar to the case for @xmath169 ( the blue left pointing triangles ) where the cusp appears at @xmath84 . at this point
, we observed a change from configuration ( @xmath170 ) to ( @xmath171 ) .
the remaining one , @xmath172 , is from the four - electron case when its configuration changes from ( @xmath171 ) to ( @xmath173 ) .
another way to obtain information on the correlations that are present in the many - particle wave function is to investigate the vortex structure . at a vortex
the many - body wave function is zero and is characterized by a change of phase of @xmath174 when we go around this point .
levels of the three ( the top ) , four ( the middle ) , and five ( the bottom ) electron quantum dot with the magnetic impurity located at @xmath76 .
the coulomb strength is @xmath69 with magnetic fields @xmath175 ( the blue squares ) and @xmath94 ( the magenta circles ) . the inset in the top figure shows the average of @xmath23 ( the dark - blue triangles ) and @xmath22 ( the red triangles ) as a function of @xmath7 of the sixth energy level.,width=321 ]    the zeros of the wave function are similar to flux quanta when e.g. the wave function corresponds to the order parameter in a superconductor .
the fixed electrons and the zero of the wave function follows closely the displaced electron and one may say that the electron plus its zero form a composite - fermion object .
the composite - fermion@xcite ( and references therein ) is a collective quasi - particle that consists of one electron bound to an even number of vortices ( flux quanta ) .
the composite - fermion concept introduces a new type of quasi - particle that is used to understand the fractional quantum hall effect in terms of the integer quantum hall effect of these composite fermions .    to obtain the zeros of the wave function of the system with @xmath2 electrons , we fix @xmath176 electrons at some positions inside the quantum dot and
leave the remaining one free .
the resulting reduced wave function gives the probability to find the remaining electrons at different positions in the quantum dot and the zeros of this function are those points where the phase of the wave function changes by @xmath174 .
as an example , we investigate the situation of a three - electron quantum dot .
levels are shown ) of the one - electron quantum dot with the magnetic impurity located at @xmath76 and @xmath69 ( the inset is a zoom of the first @xmath177 levels at low magnetic fields and the inset in the inset shows the first @xmath178 levels , it is a zoom of the rectangular region.,width=347 ]    figs .  [ vortice_lamb1 ] shows the vortex pictures of a three - electron quantum dot containing a magnetic ion located at positions ( see the white cross ) that are identical to its positions in fig .  [
denlambda1](a ) .
two , among three , electrons are fixed at the respective peaks in the electron density .
red and black regions are referring to the highest ( 2@xmath179 ) and lowest phase @xmath180 , respectively .
those plots show that there are always two vortices near the pinned electrons positions .
for example , the number of vortices pinned to each electron in the case @xmath109 at @xmath181 is @xmath4 describing the system at filling factor @xmath182 .
notice also that one of the vortices appears to be pinned at a position very close to the @xmath44 ion .
we realize that moving the magnetic impurity to a different position changes the relative positions of the vortices that are pinned to the electrons with respect to one another , as shown in figs .
[ vortice_lamb1](a ) and ( c ) .
as the electrons are antiferromagnetically coupling to the magnetic ion this kind of movement consequently depends on the position of the magnetic ion .
but now for the two - electron quantum dot .
the inset in the inset is a zoom of levels @xmath155th to @xmath177th for magnetic field close to zero.,width=347 ]    in the case @xmath69 , we found that the average of the maximum z - projection of the total angular momentum is @xmath183 and the two vortices appear at the external field @xmath184 .
apparently , the larger @xmath95 the smaller @xmath7 for which the first two vortices appear at the pinned electrons .
in the presence of an external magnetic field , the many - fold degeneracy of the energy spectrum of the system is lifted . fig .
[ spectrum3eimp0.5l0](a ) illustrates that point for the case of three electrons . in the absence of the interaction between the electrons and the magnetic ion and in the absence of a magnetic field ( blue squares )
, the energy spectrum is @xmath155-fold degenerate for the first @xmath155 lowest energy levels , the next level is then @xmath185-fold degenerate , and the next @xmath155-fold degenerate , and so on .
the origin of this was explained in ref.@xcite and is due to the coupling of the electrons and the magnetic ion .
when the magnetic field is different from zero , see red circles in fig .
[ spectrum3eimp0.5l0 ] , the degeneracy is lifted . in the inset of fig .
[ spectrum3eimp0.5l0 ] we plot @xmath186 ( the magenta triangles ) and @xmath97 ( the dark - blue ones ) as a function of magnetic field for the sixth level .
the average of @xmath186 and @xmath97 change abruptly as compared to those found for the ground - state energy , e.g. see fig .  [ zeeman3e0.5l0andlzmzsz ] .
@xmath23 and @xmath97 of the sixth state jump between two different values , e.g. @xmath187 and @xmath188 for @xmath186 and @xmath93 and @xmath189 for @xmath97 as function of the field .
this is a consequence of anti - crossings of energy levels as will be apparent later .
the result for four- and five - electron quantum dots are also shown in fig .
[ spectrum3eimp0.5l0 ] .
we see the degeneracy of @xmath190 , @xmath191 , @xmath192 , and @xmath178 for the first @xmath193 levels in the case @xmath117 and of @xmath155 , @xmath185 , @xmath155 , @xmath185 , @xmath185 , and level @xmath193 has the same degeneracy with the next energy level beyond the first @xmath193 for the case @xmath121 in @xmath194 tesla .
but now for the three - electron quantum dot.,width=347 ]    to have a clearer picture of the energy spectrum of the quantum dot system we plot in figs .
[ spectrum - b1eimp0.5l0 ] ,  [ spectrum - b2eimp0.5l0 ] , and  [ spectrum - b3eimp0.5l0 ] the magnetic field dependence of the first @xmath195 energy levels for @xmath196 , respectively .
the spectra at small magnetic fields is enlarged ( see insets ) to see the zeeman splitting and the nearly - linear behavior of the energy levels .
remember that this is due to the coupling of the electron spins with the magnetic ion spin .
for @xmath79 the first two levels for @xmath194 are @xmath155- and @xmath185-fold degenerate ( @xmath155-fold degenerate is due to the ferromagnetic coupling of the @xmath82-shell electron spin @xmath197 with the magnetic ion spin @xmath91 and @xmath185-fold degenerate of that electron now with spin @xmath198 to the magnetic ion with spin @xmath91 ) , respectively . a closer inspection ( see fig .
[ spectrum - b1eimp0.5l0 ] ) tells us that these @xmath178-levels are exchange split into two bundles of @xmath155- and @xmath185-levels ( the inset in the inset in fig .
[ spectrum - b1eimp0.5l0 ] ) .
notice that there is a first large energy gap at very small fields between the first @xmath178 levels and the next @xmath177 ones as seen in fig .
[ spectrum - b1eimp0.5l0 ] while that kind of gap appears between the first @xmath191 and the next @xmath199 for @xmath88 ( fig .  [ spectrum - b2eimp0.5l0 ] ) .
for @xmath109 ( fig .
[ spectrum - b3eimp0.5l0 ] ) , this kind of gap appears after the first @xmath177 energy levels .
the origin is the coupling of the third electron , which can reside at either two states of the @xmath81 shell while the @xmath82 shell is already fully filled , with the magnetic ion with @xmath191 @xmath89-components of the spin at very small fields , i.e. the intra - shell ( @xmath81 ) exchange interaction . for @xmath88
the electron ground state corresponds with a filled @xmath82-shell , i.e. @xmath200 , and therefor for @xmath194 only a @xmath191-fold degeneracy , as shown in fig .
[ spectrum - b2eimp0.5l0 ] , is found due to the @xmath89-component of the @xmath44-spin .
the next level is @xmath190-fold degenerate at @xmath194 ( @xmath190 come from the ferromagnetic coupling of the two - electron system with total spin @xmath201 to the magnetic ion with spin @xmath91 ) ( see the inset in the inset in fig .  [ spectrum - b2eimp0.5l0 ] ) , etc .     for three values of the coulomb strength @xmath5 , @xmath93 , and @xmath202.,width=328 ]    with increasing magnetic field
, we see that for @xmath79 there is periodically an opening of energy gaps in the spectrum .
similar energy gaps have been found earlier ( as an example see e.g ref.@xcite ) for a quantum dot without a magnetic impurity and are a consequence of the electron with two - fold spin degeneracy filling the equally - gaped - energy single - particle quantum states with different sets of the radial and angular quantum numbers .
notice that for @xmath203 , these gaps have disappeared .     and
@xmath5.,width=328 ]    the spectra exhibit a lot of crossings and anti - crossings , the number of them has increased as compared to the quantum dot case without a magnetic ion because of the zeeman splitting of the @xmath44-spin . when the applied field increases the gaps in the spectrum of @xmath79 are still open and appear more often than in the cases of @xmath203 .
once again , we see a lot of cusps in the energy levels and that reminds us to abrupt changes in the configuration of the system with magnetic field as discussed before for the ground state .
we first calculate the magnetization and susceptibility of the system : @xmath204 and @xmath205 at zero temperature .
the magnetization of a quantum dot with the magnetic ion located at @xmath76 with @xmath203 electrons is plotted in fig .
[ mag2and3e ] .
we see several jumps that are a consequence of changes in the ground state , e.g. changes in @xmath96 ( see previous section ) . for example , the magnetization of the three - electron quantum dot as plotted in fig .
[ mag2and3e](b ) for the case @xmath5 and the magnetic ion at @xmath76 has a step at @xmath206 .
consequently , the susceptibility also has a peak at @xmath206 .
the same thing happens at @xmath207 and @xmath208 for @xmath209 in the magnetization and the susceptibility .    for non - zero temperature ,
the temperature dependence of the magnetization and susceptibility is defined by : @xmath210 , @xmath211 , respectively .
the statistical average @xmath212 is calculated as : @xmath213 where the sum is over the energy levels as displayed in e.g. fig .
[ spectrum3eimp0.5l0 ] .
these quantities are explored in figs .
[ magtempo3e ] for @xmath109 and a few different temperatures ( including the zero temperature case ) .
with increasing temperature the jumps become smoother .
a very low magnetic field peak shows up because for @xmath214 we have @xmath215 at @xmath216 .
an important quantity that is related to the storage of energy is the heat capacity : @xmath217 the heat capacity is investigated as a function of the coulomb strength @xmath218 , temperature @xmath219 , magnetic field , and the position of the impurity @xmath156 .     and the magnetic field @xmath220 ( the three lower curves ) and @xmath221 ( the three upper curves).,width=328 ]    we plot in fig .  [ cvseveralpos_lamb1 ] the specific heat for two values of the magnetic field , i.e. @xmath220 and @xmath222 , and three typical positions of the magnetic ion .
for weak fields , the three electrons start to polarize and we see that the position of the main peak moves towards higher temperature as the magnetic ion is moved away but not too far from the center of the dot .
for the high magnetic field case the three electrons are strongly polarized and we see a different behavior in the shift of the main peak .
this results from the change of the statistical average of the energy levels at different fields .    , @xmath5 , at several temperatures : @xmath223 , @xmath224 , @xmath225 and @xmath226 k.,width=302 ]
now we examine the behavior of the heat capacity at a specific temperature as function of magnetic field . figs .
[ cv - b_lamb0.4_severalt ] ,  [ cv - btemp=5.8and11.6 ] , and  [ cv - bseveralt_4e ] are the plots of the magnetic field dependence of the heat capacity of three ( at two coulomb interaction regimes ) and four - electron quantum dots at some specific temperatures and two different @xmath227 .
the peak at small magnetic fields broadens and moves to higher fields with increasing temperature .
the heat capacity exhibits a number of peaks and a clear minimum around e.g. @xmath228 as shown in figs .
[ cv - b_lamb0.4_severalt ] .
remember that this field corresponds to a cusp in the energy versus magnetic field behavior as discussed previously in subsection iii.a .
at very low temperatures , this cusp still affects the heat capacity through the sharpness of the minimum as shown in the figure and this gradually becomes small at high temperatures . in fig .
[ cv - btemp=5.8and11.6 ] , we see a very interesting behavior of the heat capacity at @xmath229 : the single peak becomes a double peak as temperature increases from @xmath230 to @xmath231 . this is due to the cusps now occurring around this field in the low - energy levels of the spectrum of the three - electron quantum dot system as observed in fig .
[ spectrum - b3eimp0.5l0 ] .
besides , the structure of the heat capacity is more complex ( more peaks ) with increasing @xmath95 .
this is made clear if one looks back to the previous discussion related to figs .
[ localzeeman1234e ] and  [ zeeman3e0.5l0andlzmzsz](a ) .    ,
@xmath69 , and the temperatures are @xmath232 and @xmath233.,width=321 ]    for the case @xmath117 , the heat capacity exhibits more peaks as compared to the case @xmath109 and the behavior of the peaks with increasing temperature is also very different .
temperature affects the heat capacity of the system in the sense that it increases the peak values and separates them in magnetic field .
the coulomb interaction strength changes the structure of the magnetic field dependence of the heat capacity and is illusive through figs .
[ cv - b_lamb0.4_severalt ] and  [ cv - btemp=5.8and11.6 ] .
the peak of the heat capacity for the case with smaller coulomb interaction strength appears at higher magnetic field as compared to the case with larger one .
due to the presence of the magnetic ion ( and electron - electron interaction ) , electrons in the ground state do not always completely polarize in the presence of an external magnetic field .
the configurations are mixed consisting of electrons having spins up and down .
but for very large magnetic field , the magnetic ion tends to pull the electrons closer to the ion forming a ring - like electron density profile .
these are the consequences of the interplay of several effects such as the zeeman effect ( on the electrons and the magnetic ion s spins ) , the coulomb repulsion , the spin exchange interaction , etc .
this competition results in a crossover from ferromagnetic- to antiferromagnetic coupling between the electrons and the magnetic ion at some specific magnetic field .
interestingly , this effect is observed to appear at higher magnetic field when we move the magnetic ion further from the origin of the quantum dot
. a re - entrant behavior of the fm - afm transition is found as function of the coulomb interaction strength when the magnetic ion is moved out ( but not too far ) from the center of the quantum dot .
, @xmath5 plotted for several temperatures.,width=302 ]    the energy levels exhibit cusps as function of the magnetic field which correspond to changes of the configuration of the system as expressed by the values of ( @xmath234 ) .
these cusps move to lower magnetic field with increasing coulomb interaction strength .
the number of cusps increases with increasing number of electrons .
these cusps show up in the addition energy .    the transformation of the electron system to those of composite fermions is studied . in high magnetic fields ,
the electrons attach an even number of quantized vortices which we made clear by examining the many - body ground - state wave function in the presence of a magnetic ion . unlike the case without a magnetic ion where
all the vortices are tightly bound to the electrons , when we fix the electrons at different positions the system of vortices stays pinned to the electrons and moves with the electrons but the relative positions of the vortices are modified .
the contribution of the local zeeman splitting energy to the total energy of the system in large external fields is very small as compared to the contributions from the other parts .
however , a slight movement of the position of the magnetic ion inside the quantum dot affects the result , slightly .
with increasing applied magnetic field , each time the system jumps to a different @xmath235 configuration leads to the appearance of a peak in the thermodynamic quantities as e.g. the susceptibility and the heat capacity . in the presence of the magnetic ion ,
the structure of peaks in the heat capacity changes with the position of the magnetic ion . as temperature increases ,
these peaks split into two peaks and become smoother .
this work was supported by fwo - vl ( flemish science foundation ) , the eu network of excellence : sandie , and the belgian science policy ( iap ) .
99 j. k. furdyna , j. appl . phys . * 64 * , r29 ( 1988 ) .
d. l. klein , r. roth , a. k. l. lim , a. p. alivisatos , and p. l. mceuen , nature ( london ) * 389 * , 699 ( 1997 ) .
y. f. chen , j. h. huang , w. n. lee , t. s. chin , r. t. huang , f. r. chen , j. j. kai , and h. c. ku , appl .
. lett . * 90 * , 022505 ( 2007 ) .
k. chang , j. b. xia , and f. m. peeters , appl .
. lett . * 82 * , 2661 ( 2003 ) .
a. o. govorov , phys .
b * 72 * , 075358 ( 2005 ) .
paul i. archer , steven a. santangelo , and daniel r. gamelin , nano lett .
* 7 * , 1037 ( 2007 ) .
j. fernndez - rossier and l. brey , phys .
* 93 * , 117201 ( 2004 ) ; l. besombes , y. lger , l. maingault , d. ferrand , h. mariette , and j. cibert , phys .
93 * , 207403 ( 2004 ) ; y. lger , l. besombes , l. maingault , d. ferrand , and h. mariette , phys .
* 95 * , 047403 ( 2005 ) ; @xmath236 phys .
b * 72 * , 241309(r ) ( 2005 ) ; l. maingault , l. besombes , y. lger , c. bougerol , and h. mariette , appl .
89 * , 193109 ( 2006 ) .
t. schmidt , m. scheibner , l. worschech , a. forchel , t. slobodskyy , and l. w. molenkamp , j. appl . phys . * 100 * , 123109 ( 2006 ) .
p. wojnar , j. suffczyski , k. kowalik , a. golnik , g. karczewski , and j. kossut , phys .
b * 75 * , 155301 ( 2007 ) .
j. fernndez - rossier and ramn aguado , phys .
* 98 * , 106805 ( 2007 ) .
y. lger , l. besombes , j. fernndez - rossier , l. maingault , and h. mariette , phys .
lett . * 97 * , 107401 ( 2006 ) ; l. maingault , l. besombes , y. lger , h. mariette , and c. bougerol , phys .
( c ) * 3 * , 3992 ( 2006 ) ; m. m. glazov , e. l. ivchenko , l. besombes , y. lger , l. maingault , and h. mariette , phys .
b * 75 * 205313 ( 2007 ) .
m. b. tavernier , e. anisimovas , f. m. peeters , b. szafran , j. adamowski , and s. bednarek , phys .
b * 68 * , 205305 ( 2003 ) ; m. b. tavernier , e. anisimovas , and f. m. peeters , phys . rev .
b * 70 * , 155321 ( 2004 ) ; @xmath236 phys . rev .
b * 74 * , 125305 ( 2006 ) ; t. stopa , b. szafran , m. b. tavernier , and f. m. peeters , phys . rev .
b * 73 * , 075315 ( 2006 ) .
shun - jen cheng , phys .
b * 72 * , 235332 ( 2005 ) .
fanyao qu and pawel hawrylak , phys .
lett . * 95 * , 217206 ( 2005 ) ; _ ibid_. * 96 * , 157201 ( 2006 ) .
nga t. t. nguyen and f. m. peeters , phys .
b * 76 * , 045315 ( 2007 ) .
ramin m. abolfath , pawel hawrylak , and igor zuti , phys .
lett . * 98 * , 207203 ( 2007 ) .
j. k. jain , phys .
* 63 * , 199 ( 1989 ) ; @xmath236 phys .
b * 41 * , 7653 ( 1990 ) ; @xmath236 * 42 * , 9193(e)(1990 ) .
h. saarikoski , a. harju , m. j. puska , and r. m. nieminen , phys .
lett . * 93 * , 116802 ( 2004 ) .
s. tarucha , d. g. austing , t. honda , r. j. van der hage , and l. p. kouwenhoven , phys .
* 77 * , 3613 ( 1996 ) .