Charge dynamics effects in conductance through a large semi-open quantum dot Piotr Stefan´ski,1,∗ Arturo Tagliacozzo,2 and Bogdan R. Bu lka1 1Institute of Molecular Physics, Polish Academy of Sciences, ul. M. Smoluchowskiego 17, 60-179 Poznan´, Poland 2Coherentia-INFM, Unit´a di Napoli and Dipartimento di Scienze Fisiche, 5 Universit´a di Napoli ”Federico II”, Monte S. Angelo-via Cintia, I-80126 Napoli, Italy 0 (Dated: February 2, 2008) 0 2 Fano lineshapes in resonant transmission in a quantum dot imply interference between localized n andextendedstates. Theinfluenceof thechargeaccumulated at thelocalized levels, which screens a the external gate voltage acting on the conduction channel is investigated. The modified Fano q J parameter and the resonant conduction is derived starting from a microscopic Hamiltonian. The 7 latest experiments on ”charge sensing” and “Coulomb modified Fano sensing “ compare well with 1 theresults of thepresent model. ] PACSnumbers: 73.23.Hk,73.63.Kv l l a h I. INTRODUCTION tothecontacts,t (r=L,R),theCBstructureiseventu- r - allywashedout. Weassumethatasingleedgechannelof s e Conductance peaks and Coulomb Blockade (CB) val- energy ǫ0 is centered at µ, providing a broad single par- m leys mark the charge transport across a small quantum ticle resonance peak in the transmission. The electron . dot(QD),weaklycoupledtothecontacts,whenthegate flow acrossthe dotis essentiallyuncorrelated. Neverthe- t a voltage V is varied. However, in a large QD the multi- less, the dynamics of the charge in the core levels does g m level structure cannot be ignored. The dot may display influence the resonant tunnelling, as we show in Section - a “core” level structure including levels ǫ which are lo- II of the present work. d γ calized at different places of the dot area. Hence, they The ǫγ levels hybridize with the contacts only indi- n o may not hybridize among themselves, nor directly with rectly, by weakly coupling to the ǫ0 resonance, with a c the contact leads, but they may be weakly coupled to hoppingmatrixelementtγ. Wealsoincludesomeon-site [ one or more conduction channels. The latter could be CoulombrepulsionU at the core sites, which produces a 1 of “edge” type and extend from the left (L) to the right lower Hubbard level ǫγ (with single occupancy) and an v (R)contact. Theexistenceofsuchanedgeandcorelevel upper Hubbard level (with double occupancy) ǫγ +U. 5 structure in large QDs has been proven experimentally Byincludingthechargepiledupintheǫγ stateswithin 8 by Single Electron Capacitance Spectroscopy2. In linear the Hartree-Fock (HF) approximation, we find that the 3 transport measurements, the coupling of localized and main features of the Fano landscape are not lost. There 1 delocalizedstates,called”bouncingstates”in3,hasbeen aretwomainalterations,however: a)thechargescreens 0 invoked to justify the asymmetry of the CB peaks. the gate voltage Vg by acting on the ǫ0 level with a ca- 5 On the other hand, Fano resonances are ubiquitous in pacitive interaction and it shifts the resonances; b) the 0 / the conductance of nanodevices. They arise when local- presence of the charge modifies the Fano q factor and at izedlevelsareembeddedinacontinuumofsingleparticle induces a marked asymmetry in the peaks. m extended states which contribute to transport4. Linear InSectionIIIweelaborateonmodelingtheexperimen- conductanceonasemi-openQDdisplaysmanytransmis- tal setup that has been investigatedrecently by Johnson - d sion peaks with the characteristic Fano shape, eachtime et al.7. The conductance measured through a quantum n a core level crosses the chemical potential µ5. pointcontact(QPC)”senses”theadditionofchargetoa o In a previous paper6 we have shown that a Fano-like QD, which is capacitively coupled, parallel to the QPC. c pattern also appears when the continuum of conduction Weshowthatourmodelapplieswithsmallmodifications. : v states at µ originates from electron correlations, as is Inthiscasetheedgechannelisnotinternaltothedotbut Xi the case of a broad Kondo resonance1, due to a strongly itdescribes the QPCconductance,while the ǫγ levelslo- hybridized level ǫ , lying deeply below µ. A bunch of calizedwithin the dotarefedwith chargebyanexternal r 0 a localized core levels ǫ , weakly coupled to the ǫ level, reservoir. We neglect U which implies only minor mod- γ 0 imprints the broad Kondo peak with Fano lineshapes. ification when the static screening approximation holds, The Kondo conductance tends to freeze the charge sit- butweinclude asmallhybridizationtγ ofthe corestates ting on the ǫ level, in competition with the Fano mech- to the current carrying QPC state. This allows us to 0 anism which causes particle number oscillations. Kondo monitor the full range from bare capacitive coupling or correlationsstick the chargeoccupancy of the ǫ level to “charge sensing “ (CS), to Fano resonant tunneling. 0 the value one, and the Fano factor q close to zero. It Anincreasinggatevoltageappliedtothedotshiftsthe follows that the transmission shows a dip in form of an discrete structure of energy levels across µ and the core antiresonance. levels are gradually filled by electrons. The charge ac- In opening the QD further, by increasing the coupling cumulated in the core levels capacitatively modifies the 2 gate voltage acting on the conducting channel. In the II. MODEL HAMILTONIAN OF THE “chargesensing “ case, sudden jumps appear in the con- QUANTUM DOT AND CALCULATIONS ductance, monitoring the discreteness of the charge ac- cumulation process. Additionally, when t 6= 0, a set of γ Fano resonances emerges in the conductance across the dot as a result of quantum interference. We show that the shapeofFanoresonancesis markedlyaffectedbythe The Hamiltonian of the system under consideration is charge trapping within the core levels as found in the composed of parts describing the edge state of the dot, experiment7. H , the core QD levels, H and the electrodes H : edge core el H =H +H +H , edge core el Hedge = ǫ0−Vg +(U′/2) nγσ′a†σaσ+ tα(Vg)[c†kασaσ+h.c], σ γσ′ kσα X X X H = (ǫ −V )d† d + Un n + t [d† a +h.c.], core γ g γσ γσ γ↑ γ↓ γ γσ σ γσ γ γσ X X X H = ǫ c† c . (1) el kr krσ krσ kσ,r=L,R X The edge QD state of energy ǫ is hybridized with elec- (n = d† d is the occupation number on the core 0 γσ γσ γσ trodes. The matrix element t , (r =L, R) has a given state |γσi ). r dependence on the gate voltage V , very much like what g happens in a QPC, which opens up when the gate volt- The retarded Green function for the edge state is ob- ageincreases8. The capacitivecoupling is describedby a tained by the Equation of Motion Method (EOM) with cross-interactionU′ betweentheedgeandthecorestates anadditionalHFdecouplingofthecapacitiveinteraction: −1 t2 G (ω,V ;hn i)= h0|[ω−H (V ;hn i)]−1|0i− γ . (2) 0,σ g γσ¯ edge g γσ¯ ω−ǫ +V ( γ,σ g) γ X The first term on the r.h.s. is an average on the ground r.h.s. in Eq. 4. Charge accumulation in the core states state of the system |0i (zero temperature is assumed). themselvesshifts thepositionoftheselevelswithrespect It includes hybridization with the contact leads. The to the chemical potential even further, but we neglect occupancyofthecorelevelsisselfconsistentlydetermined this capacitative modification of their energy location in from hn i = (−1/π) ℑmG (ω)dω, where the Green the calculation, since it is inessential, as core levels do σ γ,σ function of a γ-levelwith spin σ is givenin the Hubbard not participate directly in the transport across the dot. approximation9: R The linear conductance (i.e. inzero limit ofthe drain- source voltage) in units of the quantum conductance (ω−ǫ˜γ,σ)(ω−˜ǫγ,σ−U) −1 2e2/~, is10: G (ω)= +iΓ . (3) γ,σ γ,σ ω−ǫ˜ −U(1−hn i) (cid:20) γ,σ γσ¯ (cid:21) ∞ ∂f(ǫ) G = Γ (ǫ) − ρ (ω,V ;hn i)dǫ, (5) σ 0,σ g γσ¯ Thus, G has finite width ∂ǫ γ,σ σ Z−∞ (cid:18) (cid:19) X Γγ,σ =−ℑm tγ h0|[ω−Hedge(Vg;hnγσ¯i)]−1|0itγ (4) where f(ǫ) is the Fermi distribution function, Γσ(ǫ) = Γ (ǫ)Γ (ǫ)/[Γ (ǫ) + Γ (ǫ)] is the hybridization γ Lσ Rσ Lσ Rσ X of the edge state with the contact leads (Γ = Γ L R duetoindirectinteractionwithelectrodes. Thebareǫ − is further assumed) and the spectral density ρ = γ 0,σ V isalsoadditionallyshifted→ǫ˜ bytherealpartofthe −(1/π)ℑmG . g γ 0,σ 3 Itis instructiveto introduceFanoparameterq defined 0,15 as11,12: 2,0 ℜeG (ω =0,V ;hn i) 0,10 1,5 0,σ g γσ¯ qσ(Vg)=− , (6) 1,0 ℑmG (ω =0,V ;hn i) n 0,σ g γσ¯ 0,05 0,5 with gate voltage screened by the trapped charge: eV] 0,0 VgSC =Vg −(U′/2) hnγσ(Vg)i, (7) SC [Vg 0,00 -V0g,1 [e0V,0] 0,1 0,2 0,3 0,4 -0,05 Xγσ t=0.15 U Letusconsideronesinglecorelevelforthetimebeing. -0,10 t=0.3 U TheV ,VSC arethebareandscreenedvoltageactingon t=0.6 U g g -0,15 the edge level, but they move the core level as well. In -0,15 -0,10 -0,05 0,00 0,05 0,10 0,15 crossingµ,thecorelevelisgraduallyfilledwithelectrons. Vg [eV] The occupancy of the γ level is shown in the inset of FIG. 1: Gate voltage screened by the filling of the γ level Fig.(1) for various t ’s. In turn, while the occupancy γ when the latter is moved below µ = 0 ( α = 1). The inset of the core level changes, a step appears in the screened voltage, resulting from the capacitive coupling between showstheselfconsistentlycalculatedtotaloccupancynγ. The straight line is the bare Vg and is drawn for comparison. the core and the edge state. VSC is plotted in Fig. (1) vs the bare V . Two such g g kinks appear,correspondingto the single anddouble oc- 2,0 cupancy of the γ level. Their energy separation is of the order of the Coulomb repulsion U (taken as U = 0.1eV in the calculation). The dependence of the hybridiza- 1,5 tion parameter on the voltage is implemented by mul- tiplying the matrix element tmr ax = 2U by the function h] θth(Vegs)h=arp{nexesps[(osf−thVeg)k/iδn]k+,1a}n−d1,sw(=her-0e.δ2(i=n c0a.5lc)uclaotnitornosl)s 2 G [e/1,0 1,5 1,0 ihtysbproidsiitziaotnio.nTthγebkeintwkeiesnsmǫ0oaonthdedǫγwbitehcaiunsceretahseinegleocftrtohne 0,5 -Fano q000,,,055 ==01.8 hasanincreasingpreferencetotunnelresonantlythrough -1,0 =0.4 ǫ0 insteadofdwellingontheǫγ level6 whenǫ0 crossesthe -1,5-0,1 0,0Vg [eV]0,1 0,2 =0 chemical potential. Here α = U′/U is the charge sensi- 0,0 -0,1 0,0 0,1 0,2 0,3 tivity parameter. It can increase to unity when the level Vg [eV] broadening in QD becomes larger than the level spacing andthecross-interactionU′ approacheson-siteCoulomb FIG.2: ConductancevsVg for tγ =0.15U at variouscharge repulsion U. sensitivities α of the filling of the γ-level. In the inset the InFig. (2)weplotthecorrespondingconductance,for Fano q,Eq. (6), is shown. various α. The filling of the ǫ level has a marked effect γ on the shape of the Fano resonances because it modifies the Fano q factor (see inset). The non zero value of q H ≡H + ǫBb† b + t [b† d +h.c] dot core k,σ k kσ kσ kγσ B kσ γσ and its sign are responsible for the asymmetry of the and H ≡H +H . We have added an additional QPC edge el Fano dips within the broad conductance peak centered reservoir B, in conPtact with the doPt, with levels ǫB cor- k at ǫ ∼0. respondingtosingleparticleoperatorsb . Afeatureless 0 kσ and broad density of states for the reservoir is assumed. ThedotissuppliedwithelectronsfromthebathBdueto III. CHARGE AND FANO “SENSING” IN THE the hoppingparametert . The levelǫ correspondsnow B 0 QPC CONDUCTANCE to the QPC, having a gate voltage dependent hybridiza- tion with the leads. In Figs. (3, 4) we report the results Our model can also be applied straightforwardly to for a multilevel dot with γ= 13 discrete levels equally describe recent transport measurements through a QPC spaced by ∆ = 0.2 eV and symmetrically located with “sensing” the charge piled up in a QD placed at its side respect to µ = 0 for V = 0. In this hypothesis we can g andcapacitivelycoupledtoit7. Weshowthatitispossi- assume that ∆ is the addition energy of the dot, inde- ble to continuously move from a “capacitive sensing” to pendent of the electron number N. Thus, we drop the a “Fano sensing” by increasing the hybridization of the interaction term ∝ U in H and the relevant energy core single particle levels ǫ of the dot with the conduction scale for charge sensing is now ∆ itself. A static screen- γ channel ǫ of the QPC. ing in the multi-electron system is expected to be quite 0 To match with the experimental setup we use the effective, if the dot is large. The screened gate voltage Hamiltonian of Eq. (1) by reinterpreting it as follows: acting on the QPC, sensing the charge on the dot, takes 4 broadresonancecenteredatǫ = V = 0. Its widthand 0 g 2,0 characteristicasymmetryisgovernedbythe gatevoltage dependent coupling to the L,R contacts. By increasing α,themaximumoftheresonanceismovedtohigherV ’s, 1,5 =0 g =0.5 because of the voltage screening. The conductance peak =1 acquires a sawtooth pattern, which grows sharper when h] 2 G [e/1,0 αtiminecraeanseews. eTlehcetrpoantteenrnterfoslltohwesdtohte.juTmhpessohfaVpgesco,feathche overall conductance curve is similar to the one seen in 0,5 the experiment (see Fig. (2a) of7). In Fig.(4), the conductance is reported for finite hop- ping between QD andQPC (t =0.15 ∆), atα=0.0,0.5, 0,0 γ -1,5 -1,0 -0,5 0,0 0,5 1,0 1,5 1.0. Thecaseα=0isthepureFanocase. Themaximum Vg [eV] of the conductance is still centered at Vg ≈0 because tγ has been chosen to be rather small. The zero of the real FIG. 3: Conductance through the QPC with ”charge sens- part of the QPC Green’s function G defines the loca- ing” only (tγ =0) for various α parameter and tB/∆=0.15. tion of the maximum of the resonan0c,σe, and the zero of The curve for α=0 without charge sensing is also shown for the Fano q as well, according to Eq.(6). The symmetric comparison. Fanodipcorrespondingtoq =0iscenteredontopofthe conductance maximum, because a γ level sits exactly at 2,0 Vg = 0. When α increases, the condition ǫ0−VgSC ∼ 0 2 0 nowmarksthelocationoftheconductancemaximum,as 1,5 Fano q---642 ==00 .5 =1 wTehlulsa,stthheecpolnadcuecwtahnecree qen∼vel0op(seeemtohveesinasgeatinoftFoigla.(r4g)e)r. -8 =1 -1-01,5-1,0-0,50,0 0,5 1,0 1,5 voltages,together with the Fano-like structures, and the Vg [eV] h] conductance peak broadens as in Fig. (3). Appearance 2 G [e/1,0 of sharp peaks for large negative gate voltages is caused by rapid decrease of q (see inset) due to the closing of =0.5 QPC, Eq. (6). Our case with α = 1.0 compares well 0,5 withFig.(2b)of7,displayingCoulombmodifiedFanores- =0 onances. Both ref.7 and in ref.14 propose a model con- 0,0 ductance based on the grand canonical average of dot -1,5 -1,0 -0,5 0,0 0,5 1,0 1,5 configurations, which can be fitted to the experimental Vg [eV] data. Ourcalculationshowsthatthedynamicsoftheres- onant tunnelling together with some minor assumptions FIG.4: ConductancethroughtheQPCwith”chargesensing” andfinitehybridization with thedot(tγ=0.15∆) forvarious providesimplerresults,inqualitativeagreementwiththe α parameter. Inset shows corresponding Fano q, the straight experimental data. lineisforα=0andconstantcouplingtoelectrodestr =tmr ax. To conclude, we have demonstrated that a Coulomb modified Fano resonant tunneling arises from charge sensing in structures like a QPC with a dot aside, or in- now the same form as Eq. (7), but with U′ replaced side a large dot where core levels coexist with extended by the charging energy of the dot ∆. The α parameter edge states and both can be described within the same now describes the sensitivity of the capacitive coupling quantum mechanical model with minor modifications. between the QD and QPC which can be controlled by appropriate gates7. Again, the particle number on the dot n (V ) is determined by means of the Friedel γ,σ γ,σ g sum rule13. It shows the usual steps as in the inset of Acknowledgments P Fig. (1) when the gate voltage increases and the dot is filledbyelectrons. Forα=1theCSisequaltothecharg- This work was supported by the Polish State Com- ingenergy,eachtime anelectronhopsontothe dot. Fig. mittee for Scientific Research (KBN) under Grant no. (3) shows the conductance through the QPC when the PBZ/KBN/044/P02/2001 and the Centre of Excellence dot is only capacitively coupled to it (t =0), for various for Magnetic and Molecular Materials for Future Elec- γ α’s. For α=0 the dot has no effect on the QPC con- tronics within the European Commission Contract No. ductance. The transmission across the QPC displays a G5MA-CT-2002-04049. ∗ Electronic address: [email protected] 1 D. Goldhaber-Gordon et al.,Nature, 391, 156 (1998). 5 2 N.B. Zhitenev et al.,Phys. Rev.Lett. 79, 2308 (1997). 10 Y. Meir and N.S. Wingreen, Phys. Rev. Lett. 68, 2512 3 S. Lindermann et al.,Phys. Rev.B 66, 161312 (2002). (1992). 4 U. Fano, Phys.Rev. 124, 1866 (1961). 11 J.N¨ockelandA.D.Stone,Phys.Rev.B50,17415 (1994). 5 C.Fu¨hner et al.,cond-mat/0307590(2003); Phys.Rev.B 12 P. Stefan´ski, Solid St. Commun. 128, 29 (2003). 66, 161305 (2002). 13 A.C. Hewson, The Kondo Problem to Heavy Fermions, 6 P. Stefan´ski et al.,Phys.Rev. Lett.93, 186805 (2004). (CambridgeUniv.Press,Cambridge1993),D.C.Langreth, 7 A.C. Johnson et al., Phys.Rev.Lett. 93, 106803 (2004). Phys. Rev. 150, 516 (1966). 8 S.M. Cronenwett et al. PRL 88, 226805 (2002). 14 R.Berkovits et al.,cond-mat/0411478. 9 A.C. Hewson, Phys. Rev. 144, 420 (1966)