Voltage-controlled Kosterlitz–Thouless transitions and various kinds of Kondo behaviors in a triple dot device?

Yong-Chen Xiong(熊永臣),Jun Zhang(張俊),Wang-Huai Zhou(周望懷),and Amel Laref,3

1 School of Science,and Advanced Functional Material and Photoelectric Technology Research Institution, Hubei University of Automotive Technology,Shiyan 442002,China

2 Department of Physics,Faculty of Science,National University of Singapore,Singapore 117551,Singapore

3 Department of Physics and Astronomy,Science Faculty,King Saud University,Riyadh 11451,Saudi Arabia

Voltage-controlled Kosterlitz–Thouless transitions and various kinds of Kondo behaviors in a triple dot device?

Yong-Chen Xiong(熊永臣)1,2,?,Jun Zhang(張俊)1,Wang-Huai Zhou(周望懷)1,and Amel Laref1,3

1 School of Science,and Advanced Functional Material and Photoelectric Technology Research Institution, Hubei University of Automotive Technology,Shiyan 442002,China

2 Department of Physics,Faculty of Science,National University of Singapore,Singapore 117551,Singapore

3 Department of Physics and Astronomy,Science Faculty,King Saud University,Riyadh 11451,Saudi Arabia

The transport property and phase transition for a parallel triple dot device are studied by adopting Wilson’s numerical renormalization group technique,focusing on the effects of level spacings between neighboring dot sites.By keeping dot 2 at the half-filled level and tuning the level differences,it is demonstrated that the system transits from local spin quadruplet to triplet and doublet sequently,and three kinds of Kondo peaks at the Fermi surface could be found,which are separated by two Kosterlitz–Thouless type quantum phase transitions and correspond to spin-3/2,spin-1,and spin-1/2 Kondo effect, respectively.To obtain a detailed understanding of these problems,the charge occupation,the spin–spin correlation,the transmission coefficient,and the temperature-dependent magnetic moment are shown,and necessary physical arguments are given.

triple quantum dot,Kosterlitz–Thouless transition,Kondo effect,critical phenomena

1.Introduction

For the last several years,systems of coupled multi-quantum dots have gained significant interests both experimentally[1–6]and theoretically,due to rapid progress in spintronics[7–9]and quantum information.[10,11]An important example of a multi-dot system is the triple quantum dots(TQD).Such structure exhibits many interesting quantum phenomena in the strong correlation limit.For instance,the Aharanov–Bohm(AB)effect and magnetic frustration are studied in systems with triangular geometry,[12–15]and the multi-channel Kondo effect related to a number of Fermi-liquid and non-Fermi-liquid behaviors are investigated in various TQD models.[16–20]Furthermore,TQD systems are also considered as ideal models to show the Ruderman–Kittel–Kasuya–Yosida(RKKY)interaction,[21,22]the quantum interference,[23–26]and various kinds of quantum phase transition(QPT).[13,15,27–32]Basically,these behaviors depend closely on the organizations of the quantum dots and the conduction leads,as well as the coupling/interaction elements which have been taken into account.

Within the above phenomena,the level difference plays an important role as suggested in some double quantum dot (DQD)structures.For instance,a local spin triplet–singlet transition of Kosterlitz–Thouless(KT)type is found as the difference increases.[33]For a parallel DQD system,it is revealed that the triplet–doublet transition could be KT type orfirst order depends on the breaking of the spin-rotation SU(2) symmetry,[34]and the linear conductance is shown to have an asymmetric line shape of the Fano resonance when the interdot tunneling coupling is taken into account.[35]In a spinless two-level dot system,charge oscillation occurs in the presence of level spacing as the gate voltage sweeps.[36]These works are limited in the double dot structures,but less is studied in the TQD system.Compared to the double dot system,TQD structure exhibits more attractive behaviors,owing to its complex geometry and much more Feynman paths for the electron transmission.[25]Furthermore,experimentally speaking,scaling up the number of spin-1/2 qubits is an inevitable step towards the realization of quantum computation, not only for extending the qubit system,but also for performing practical quantum algorithms.[37]Very recently,plenty of experimental works have been carried out on this topic in TQD systems,[5,6,38–41]therefore,to illustrate the electronic transport,phase transition and relevant quantum phenomena in such systems also becomes an emergent task.

In this paper,we consider a parallel TQD system with level spacings between neighboring dot sites.With the help of the numerical renormalization group(NRG)method,and fixing dot 2 at the half-filled level,we demonstrated that the level spacings play important roles in the linear conductance and the QPT.Our main findings include the following:by tuning the level differences,the ground state of the system transitsfrom local spin quadruplet to triplet,then doublet sequentially, and three kinds of Kondo peaks at the Fermi surface could be found,corresponding to spin-3/2,spin-1,and spin-1/2 Kondo effect,respectively.Two Kosterlitz–Thouless(KT)type transitions are clarified,resulting from asymmetric Kondo couplings between conduction leads and different dots,as well as the unequal charge occupation of three dots.To handle these problems,some important physical quantities are shown,and necessary physical arguments are given.

The outline of the paper is as follows.In Section 2,we define the model Hamiltonian of the TQD system,and present the calculation algorithms and formulations.In Section 3,we show the quantum phase transitions and transport properties with respect to the increasing level spacings.Finally,a summary is given in Section 4.

2.Model and method

We illustrate the second-quantized form of the model Hamiltonian for the triple dot device in Fig.1,

where the first part is for the electrons in the dots,

Here,εiis the energy level for the i-th dot.For simplicity, we label two level spacings as Δ1=ε2?ε1and Δ2=ε3?ε2, which could be controlled by external gate voltage.is the creation operator for the electron with spin-σ(=↑or↓).U is the on-site(intra-dot)Coulomb repulsion,andni=ni↑+ni↓are the number operators.

Fig.1.Schematic view of the parallel TQD device connected to the L and R leads.

The second part is for the non-interacting electrons in the conduction lead,wherecreate a spin-σ electron of wave vector k and energy εkin lead ν(ν=R or L).By contrast,the last term is for the coupling interaction energy between the conduction lead and the dots,where Vkis the tunnel matrix element,and is assumed to be σ-independent,identical for three dots,and also symmetric to the right and left leads.

In the following discussion,we concentrate on the quantum behaviors in the strong correlation limit,hence a sophisticated theoretical technique must be adopted.Therefore,we use the celebrated NRG method[42–44]to solve Eq.(1).To simplify the problem into a more convenient form,we first assume the density of state for the conduction lead ρ0=1/(2D)and the hybridization function(dot–lead coupling)Γ=2πρ0|Vk|2are constant.Here,D is the half width of the conduction band, and is chosen for the energy unit throughout this paper.The strong correlation limit is defined as U?Γ,and the number of low-lying states kept at each iteration is about 3000.Furthermore,the discretization parameter Λ,which characterizes the logarithmic discretization of the conduction band,is set to be 1.8–2.2.

The total charge number Ntotand the local spinare defined as

where Siis the local spin operator of dot i.

The linear conductance through the device G is determined by the Landauer formula[45]

Here,f(ω)is the Fermi function and T(ω)is the transmission probability.The retarded dot Green function is written asThus

If we choose the Fermi energy at εkf=0,G at zero temperature in the limit of zero bias is given by

Finally,we define the temperature-dependent magnetic moment of the dotsμ2(T)at temperature T as

where χ(T)is the magnetic susceptibility,kBis the Boltzmann constant,g andμBare the electronic g factor and the Bohr magneton,respectively.〈...〉denotes the thermodynamic expectation value,and the subscript 0 labels the case without quantum dots.

3.Results and discussion

In this section,we study the phase transition and the electronic transport of the TQD device in the strongly correlated regime,focusing on the phenomena when the level differences sweep upwards.The on-site Coulomb repulsion,the dot-lead coupling,and the charge energy of dot 2 are fixed at U=0.1, Γ=0.01,and ε2=?U/2,respectively,throughout this paper. For convenience,we choose Δ2=Δ1/2.However,it is worth noting that our conclusions are robust for general cases,and do not really require Δ1and Δ2to satisfy the above relation.

Fig.2.(color online)(a)Charge occupation〈ni〉on each dot,(b)spin correlation〈SiS j〉between dotsiand j,and(c)local spinof the dots at zero temperature as functions of Δ1.Here,Γ=0.01,U=0.1, ε2=?U/2 and Δ2=Δ1/2.

Let us start our discussion from the quantum phase transition for the TQD system.In Figs.2(a)–2(c),we depict the charge occupation〈ni〉of each dot,the spin correlation〈SiSj〉 between dotsiand j,and the localspinon the dots as functions of Δ1.We can see when the level differences are absent, three dots are singly occupied with〈ni〉=1.0,since the intradot repulsion U is strong enough.In this case,〈SiSj〉≈0.20 and≈3.20,implying the local spins on three dots are arranged parallelly and a local spin quadruplet is generated,owing to the ferromagnetic RKKY interaction mediated by the Kondo exchange coupling between electrons on the leads and those on the dots.It is noted that〈SiSj〉andcould not reach the ideal values of 1/4 and 15/4 for spin quadruplet,as a finite U/Γ is adopted.With increasing Δ1,the level ε1decreases,while ε3increases but with a lower speed.Therefore,〈n1〉grows to 2.0 gradually,while〈n2〉and〈n3〉retain nearly at 1.0.In this process,changes continuously to about 1.75,describing that the ground state is dominated by the spin triplet.In this triplet,dot 1 is almost doubly occupied,hence〈S1S2〉and〈S1S3〉decrease to about 0,while〈S2S3〉has nearly no change.As Δ1increases continuously,〈n3〉decreases to the regime n3≈0.In this case,the system maps to the single impurity case and the ground state is a local spin doublet with〈SiSj〉≈0 and≈0.75.

Fig.3.(color online)(a)–(c)Transmission coefficient T(ω)at zero temperature for various Δ1.The remaining parameters are the same as in Fig.2.

Figures 3(a)–3(c)show the transmission coefficient T(ω) at zero temperature for various Δ1.It is seen that the spectral weight is symmetric to the Fermi level ω=0 when Δ1=0, since Eq.(2)satisfies particle–hole(p–h)symmetry.Two Coulomb peaks located at ω=±U/2 are found,referring to the process of annihilating(creating)an additional electron (hole)on the bonding orbitalBesides,a Kondo peak could also be seen at the Fermi level, which reaches the unitary limit of 2e2/h corresponding to full conductance.It could be attributed to the process where the itinerant electron screens a local spin-1/2 degree of freedom and results in a partially screened spin-3/2 Kondo effect.[15,21]As Δ1increases,the spectral weight moves to the left and is not symmetric to ω=0,for the p–h symmetry is broken and Ntotis away from a triply occupied state.In this process,the amplitude of the Kondo resonance peak decreases, implying that the linear conductance is reduced slightly(see Δ1=0.043 in panel(a)).As Δ1exceeds the first critical point=0.0462,a sharp Kondo peak appears at ω=0 (see Δ1=0.05 in Fig.3(b)),which is related to the partially screened spin-1 Kondo effect,and originates from the screen-ing of the local spin triplet between dot 2 and 3 by the conduction electrons.With increasing Δ1,the spectral weight moves away from ω=0 and the Kondo peak is broadened,indicating a decrease of the conductance(e.g.,Δ1=0.07).When Δ1grows larger than the second critical=0.0863,another sharp Kondo peak is observed at the Fermi surface(see Δ1=0.10 in panel(c)),which results from the antiferromag-netic Kondo coupling between the conduction leads and dot 2, since dot 1 is almost doubly occupied and dot 3 is empty.As Δ1grows continuously,the Kondo peak keeps at ω=0,for the TQD is stable in the spin doublet(see Δ1=0.16).

To get more information about the phase transition,we show T(ω)on the triplet side near=0.0462 in Fig.4(a). One finds that the width of the spin-1 Kondo peak W1enlarges with increasing Δ1.Here,the width of the spin-1 Kondo peak is defined as the half width at half maximum of the Kondo peak.Figure 4(b)suggests that W1depends exponentially on the distance to the critical pointand can be adequately described using an exponential functionwhere the fitting parameters are given by P1=?0.0192,P2=3.9409,and P3=0.0552. It is seen that the fitting function(solid line)agrees very well with our NRG results(scatter dots).This behavior describes that the quadruplet–triplet transition atis a KT transition.Figure 4(c)gives T(ω)for various Δ1on the doublet side near.One may also find that the width of the spin-1/2 Kondo peak W2approaches zero exponentially as Δ1is close to,and can be described by an exponential functionHere,the fitting parameters are P4=?0.0019,P5=6.4015,P6=0.0431.Figure 4(d)indicates that the fitting function(solid line)agrees very well with the NRG results(scatter dots),illustrating the triplet–doublet transition is also of the KT type.

To exhibit more information about the low-temperature scenario,we presentμ2(T)for different Δ1in Fig.5.The upmost curve is for the case of Δ1=0.It is seen that TQD goes through four different regimes as T decreases.For example,when T is high enough(e.g.,T＞U),the electrons on three dots are independent,then each dot contributes 1/8 toμ2(T).As T decreases to T＜U(e.g.,T～0.01),the electrons are then in the local moment regime,thus the full-and zero-occupied states are strongly suppressed.Therefore, each of them contributes 1/4 toμ2(T).When T decreases to the order of the RKKY interaction(e.g.,T～0.001),the electrons on three dots form a spin quadruplet and contribute Sdot(Sdot+1)/3=5/4 toμ2(T).As a result,μ2(T)grows to a higher spin stage.Here,μ2(T)in the local moment regime (μ2(T)～0.62)and the RKKY regime(μ2(T)～0.80)are smaller than the ideal values,for the TQD is not in a pure ground state in these temperature regimes.When T reduces to the Kondo temperature scale(e.g.,T～10?7),the itinerant electron screens a local spin-1/2 degree of freedom,thus one finds an obvious drop in μ2(T)and μ2(T)～0.67 at low temperature.As Δ1increases,the height ofμ2(T)in the RKKY regime is reduced,since the possibilities of generating the RKKY interaction decreases due to increasing double occupancy on dot 1.On the other hand,μ2(T)at low temperature μ2(T=0)decreases exponentially.When Δ1is large enough (e.g.,Δ1=0.052),μ2(T=0)≈0.25,since in this case the spin triplet is partially screened by the conduction leads,and the remaining spin-1/2 degree of freedom contribute toμ2(T). As Δ1grows continuously,μ2(T=0)reduces exponentially again.When Δ＞,μ2(T)→0 at low temperature,for in this regime the local spin singlet is totally screened by the conduction electrons.

Fig.4.(color online)(a)T(ω)on the triplet side near the first critical points=0.0462.The curves from top to bottom on the left side are for Δ1=0.051 to 0.060 in steps of 0.001.(b)Width of the spin-1 Kondo peak W1 (scatter dots)and its fitting exponential functions(solid lines).(c)T(ω)on the doublet side near the second critical point=0.0863.The curves from top to bottom on the right side are for Δ1=0.091 to 0.099 in steps of 0.001.(d) Width of the spin-1/2 Kondo peak W2(scatter dots)and its stimulant functions (solid lines).The remaining parameters are the same as in Fig.2.

Fig.5.Total magnetic moment of the TQDμ2(T)as a function of temperature T according to different Δ1,the curves from top to bottom are for Δ1=0,0.04,0.046 to 0.050 in steps of 0.001,0.052(dash line),and Δ1=0.08,0.087 to 0.09 in steps of 0.001,0.094,0.1,respectively.The remaining parameters are the same as in Fig.2.

To explore the physical origination of the KT transitions, we change Eq.(1)to a three impurities spin-1/2 Kondo model by adopting the Schrieffer–Wolff transformation.[?,21]The effective Hamiltonian then could be written as

With increasing Δ1,Jkiare asymmetric.On the other hand,the charge occupation in three dots also become different.The symmetric broken by these features are the origination of the KT transitions,similar to a two-level dot system with level difference[33]and inter-dot Coulomb repulsion.[34]Furthermore,it is worth noting that one can estimate the critical points by considering the energy levels of an isolated TQD model.For instance,the energy level for the spin quadruplet can be written as EQ=ε1+ε2+ε3=?3U/2.While that for the triplet is ET=2ε1+ε2+ε3+U=?U?3Δ1/2, and the level for the doublet is given by ED=2ε1+ε2+U=?U/2?2Δ1.As Δ1increases,the energy levels for the triplet and the singlet are pulled down,and eventually degenerate with the quadruplet at Δ1=U/3 and the triplet at U respectively.As a result,two QPTs of the KT type can be found atand

4.Summary

In conclusion,we have studied the phase transition and Kondo behavior in a triple quantum dot device with parallel organization in the strongly correlated regime.We concentrate on the effect of the level differences Δ1=ε2?ε1and Δ2=ε3?ε2.By keeping dot 2 at the half-filled level and tuning the level differences,it is shown that the ground state of the system transits from local spin quadruplet to triplet,then doublet sequentially,three kinds of Kondo peaks at the Fermi surface could be found in the transmission coefficient,which correspond to spin-3/2,spin-1,and spin-1/2 Kondo effect,respectively.Two KT type phase transitions are clarified,resulting from asymmetric Kondo couplings between the conduction leads and different dots,as well as unequal electron occupation of three dots.We believe our work not only clarifies the effect of the level difference on the phase transition and the electronic transport in a triple dot structure,but it may also afford useful guidance for spintronics and molecular electronics devices.

[1]Vidan A,Westervelt R M,Stopa M,Hanson M and Gossard A C 2004 Appl.Phys.Lett.85 3602

[2]Schr?er D,Greentree A D,Gaudreau L,Eberl K,Hollenberg L C L, Kotthaus J P and Ludwig S 2007 Phys.Rev.B 76 075306

[3]Amaha S,Hatano T,Teraoka S,Tarucha S,Tokura Y,Miyazaki T,Oshima T,Usuki t,Yokoyama N 2008 Appl.Phys.Lett.92 202109

[4]Gaudreau L,Kam A,Granger G,Studenikin S A,Zawadzki P and Sachrajda A S 2009 Appl.Phys.Lett.95 193101

[5]Takakura T,Pioro-Ladrière M,Obata T,Shin Y S,Brunner R,Yoshida K,Taniyama T and Tarucha S 2010 Appl.Phys.Lett.97 212104

[6]Gaudreau L,Granger G,Kam A,Aers G C,Studenikin S A,Zawadzki P,Pioro-Ladrière M,Wasilewski Z R and Sachrajda A S 2012 Nat. Phys.8 54

[7]?uti? I,Fabian J and Sarma S D 2004 Rev.Mod.Phys.76 323

[8]Prinz G A 1998 Science 282 1660

[9]Wolf S A,Awschalom D D,Buhrman R A,Daughton J M,von Molnar S,Roukes K L,Chtchelkanova A Y and Treger D M 2001 Science 294 1488

[10]Loss D and DiVincenzo D P 1998 Phys.Rev.A 57 120

[11]DiVincenzo D P 2005 Science 309 2173

[12]Delgado F,Shim Y P,Korkusinski M,Gaudreau L,Studenikin S A, Sachrajda A S and Hawrylak P 2008 Phys.Rev.Lett.101 226810

[13]Wang W Z 2008 Phys.Rev.B 78 235316

[14]Anderson P W 1973 Mater.Res.Bull.8 153

[15]Wang W Z 2007 Phys.Rev.B 76 115114

[16]Kuzmenko T,Kikoin K and Avishai Y 2003 Europhys.Lett.64 218

[17]?itko R and Bon?a J 2008 Phys.Rev.B 77 245112

[18]Mitchell A K,Logan D E and Krishnamurthy H R 2011 Phys.Rev.B 84 035119

[19]Mitchell A K,Galpin M R,Wilson-Fletcher S,Logan D E and Bulla R 2014 Phys.Rev.B 89 121105(R)

[20]Tooski S B,Ram?ak A and Bu?ka B R 2016 Physica E 82 366

[21]?itko R and Bon?a J 2006 Phys.Rev.B 74 045312

[22]Xiong Y C,Wang W Z,Luo S J,Yang J T and Huang H M 2016 J. Magn.Magn.Mater.399 5

[23]Ladrón de Guevara M L and Orellana P A 2006 Phys.Rev.B 73 205303

[24]Webb R A,Washburn S,Umbach C P and Laibowitz R B 1985 Phys. Rev.Lett.54 2696

[25]Fu H H and Yao K L 2010 J.Appl.Phys.108 084510

[26]Huang R,Ming S W and Wang Y 2012 Chin.Phys.Lett.29 47201

[27]?itko R and Bon?a J 2007 Phys.Rev.B 76 241305(R)

[28]Mitchell A K,Jarrold T F and Logan D E 2009 Phys.Rev.B 79 085124

[29]Mitchell A K and Logan D E 2010 Phys.Rev.B 81 075126

[30]Mitchell A K,Jarrold T F,Galpin M R and Logan D E 2013 J.Phys. Chem.B 117 12777

[31]Xiong Y C,Huang J and Wang W Z 2012 J.Phys.:Condens.Matter. 24 455604

[32]Xiong Y C,Wang W Z,Yang J T and Huang H M 2015 Chin.Phys.B 24 027501

[33]Hofstetter W and Schoeller H 2001 Phys.Rev.Lett.88 016803

[34]Wang W Z 2011 Nanotechnology 22 205203

[35]Ding G H,Kim C K and Nahm K 2005 Phys.Rev.B 71 205313

[36]SindelM,Silva A,Oreg Y and von Delft J 2005 Phys.Rev.B 72 125316

[37]Steffen L,Salathe Y,Oppliger M,Kurpiers P,Baur M,Lang C,Eichler C,Puebla-Hellmann G,Fedorov A and Wallraff A 2013 Nature 500 319

[38]Medford J,Beil J,Taylor J M,Rashba E I,Lu H,Gossard A C and Marcus C M 2013 Phys.Rev.Lett.111 050501

[39]Medford J,Beil J,Taylor J M,Bartlett S D,Doherty A C,Rashba E I,DiVincenzo D P,Lu H,Gossard A C and Marcus C M 2013 Nat. Nanotechnol.8 654

[40]Takakura T,A Noiri,Obata T,Otsuka T,Yoneda J,Yoshida K and Tarucha S 2014 Appl.Phys.Lett.104 113109

[41]Noiri A,Yoneda J,Nakajima T,Otsuka T,Delbecq M R,Takeda K, Amaha S,Allison G,Ludwig A,Wieck A D and Tarucha S 2016 Appl. Phys.Lett.108 153101

[42]Krishna-Murthy H R,Wilkins J W and Wilson K G 1980 Phys.Rev.B 21 1003

[43]Krishna-Murthy H R,Wilkins J W and Wilson K G 1980 Phys.Rev.B 21 1044

[44]Bulla R,Costi T A and Pruschke T 2008 Rev.Mod.Phys.80 395

[45]Meir Y,Wingreen N S and Lee P A 1993 Phys.Rev.Lett.70 2601

[46]Schrieffer J R and Wolff P A 1966 Phys.Rev.149 491

2 May 2017;published online 31 July 2017)

10.1088/1674-1056/26/9/097102

?Project supported by the National Natural Science Foundation of China(Grant No.11504102),the Scientific Research Items Foundation of Hubei Educational Committee,China(Grant Nos.Q20161803 and B2016091),the Doctoral Scientific Research Foundation(Grant No.BK201407),and the Major Scientific Research Project Pre-funds of Hubei University of Automotive Technology,China(Grant No.2014XY06).

?Corresponding author.E-mail:xiongyc lx@huat.edu.cn

?2017 Chinese Physical Society and IOP Publishing Ltd http://iopscience.iop.org/cpb http://cpb.iphy.ac.cn