Quantum dot with ferromagnetic leads: a densiti-matrix renormaliation grou study C. J. Gaa, M. E. Torio, and J. A. Riera Instituto de Física Rosario, Consejo Nacional de Investigaciones Científicas y Técnicas, Universidad Nacional de Rosario, Rosario, Argentina (Dated: February 6, 8) arxiv:cond-mat/6368v [cond-mat.str-el] 3 May 6 A quantum dot couled to ferromagnetically olaried one-dimensional leads is studied numerically using the density-matrix renormaliation grou method. Several real sace roerties and the local density of states at the dot are comuted. It is shown that this local density of states is suressed by the arallel olariation of the leads. In this case we are able to estimate the length of the Kondo cloud, and to relate its behavior to that suression. Another imortant result of our study is that the tunnel magnetoresistance as a function of the quantum dot on-site energy is minimum and negative at the symmetric oint. PACS numbers: 73.63.Kv, 7.5.Qm, 7.5.-b The relentless ursuit of smaller, faster, more efficient electronic devices has led to an increasing interest in mesoscoic and nanoscoic devices behaving as quantum dots (QD), using semiconductor technology and nanotubes as comonents among other ossibilities. Some of these devices allow the exloitation of the electron sin in addition to its charge with various otential alications in the so-called sintronics. 3 In articular, a sintronic device connecting a QD to ferromagnetic leads has already been exerimentally studied. 4 These nanoscoic devices have also been growing in comlexity. For examle a device consisting of two QDs has been develoed to show a nonlocal control of one QD over the other via a RKKY interaction. 5 Another examle is the device where an aroriate arrangement of electrodes defines a section of a nanotube as the QD. 6 This QD turns out to be then connected to essentially one-dimensional leads. In the first examle, a real sace descrition of the magnetic correlations would be imortant to understand the internal working of this device. In the second examle, the effects of finite length of the nanotube sections connected to the QD would have to be considered. Among many techniques used to study these kinds of roblems, which have been develoed in the context of the Kondo effect, 7 the numerical renormaliation grou (NRG) has rovided many imortant and essentially exact results. However, this technique cannot tackle the increasingly comlex devices such as the above mentioned ones. For this reason, we roose in this aer the alication of a real sace numerical technique, the densitymatrix renormaliation grou (DMRG) 8, which has been extensively used to study quasi-one-dimensional strongly correlated electron systems. 9 This technique is ideally suited to rovide detailed real sace information, such as sin-sin correlations or electron site occuancies, required to understand nanoscoic devices. In addition, it works on finite sie systems. This is not a disadvantage when bulk leads are involved in the device because in most cases one could resort to some kind of extraolation. On the contrary, it can rovide results for finite systems such as the examle mentioned before. In this aer, we will study the simlest device in sintronics, the sin valve. This device consists of a singlelevel QD attached to two noninteracting leads. On these leads, a ferromagnetic olariation is introduced by an alied magnetic field, or in other words, by a rigid dislacement of sin-u and sin-down electron bands. The central quantities to study are the conductance when the olariation of the leads are arallel (G P ) and antiarallel (G AP ). The main measure of the erformance of the sin valve is the tunnel magnetoresistance (TMR) defined as TMR = (G P G AP )/G AP. This device has been studied before using different aroximations,,. These analytical techniques are not exact and in fact they lead to some degree of controversy around central issues such as the behavior of the conductance when the leads have a arallel olariation. Some of these controversial issues have been artially settled using NRG. 3,4 As in revious studies on this system, the Hamiltonian model for a QD located at the center of the chain is a generaliation of the Anderson model defined as H = t (c iσ c i+σ + H.c.) h L t i=,σ i=,σ (c iσ c i+σ + H.c.) h R i= i= t (c σ c σ + c σ c σ + H.c.) σ + ǫ n + Un, n, () where conventional notation was used. We adot t as the scale of energy. The magnetic field on the leads can be set in two configurations, h L = h R = h, which corresonds to a arallel (P) olariation of the leads, and h L = h R = h which corresonds to the antiarallel (AP) one. In the following, h > favors a ositive olariation. Notice that the leads are described by a real-sace tight-binding model which is in rincile more realistic than the ones imlicit in the NRG treatment. Model equation() was studied by the DMRG on L = 63, 79 and 95 chains with oen boundary conditions. Most of the results shown below were obtained for two sets of arameters: U =, t =.4 and U = 4, S i S i
4 8 6...3.4 FIG. : (Color online) The LDOS for U = 4, t =.8, and ǫ =., for various values of the olariation, antiarallel, L = 79 (solid lines), L = 95 (dashed lines); arallel configuration, from to to bottom at = : L = 79, =.3,.3,.3,.46 (solid lines); L = 95, =.,.399 (dashed lines).the inset shows the length of the comensation cloud (see text) for the P case, L = 79 (oen circles), L = 95 (filled circles). which minimies the total energy, and hence it deends on h. The olariation of each lead α (α = L, R) is defined as α = (nα nα )/(nα + nα ). Then, the olariation of the P (AP) configuration is (L +R)/ [(L R)/]. For small clusters, and secially for the U = case, it is difficult to get L = R, i.e. the same numbers of solitons on each lead. In these cases the AP case shows surious suression of the LDOS. 7 In Fig. we show the LDOS at the QD, U =, t =.4, at the symmetric oint ǫ = U/, for several values of the olariation. In this figure and all the following similar ones, we adoted a Lorentian broadening of the eaks of δ =. In the first lace, for =.3, corresonding to h =, the LDOS is the one exected in the Kondo roblem. The small slitting of the Kondo eak is a finite sie effect. For the case of an AP configuration (Fig. ), a small value of =.88 leads to a slight modification of the LDOS, that virtually remains unmodified by further increasing the olariation of the leads, at least in the range examined (.3). The largest value considered, =.3, is obtained for h =. (L = 95). On the other hand, for the P configuration (Fig. ) it can be seen as a suression of the Kondo eak for the olariation.. For small values of the eak at = remains roughly unmodified. As we will show below, this different deendence of the LDOS with for the AP and P arrangements leads to the exected behavior of a sin valve. The inset in Fig. shows, for the P magnetiation, the length of the comensation cloud 8, a ossible measure of the Kondo cloud, defined in such a way that the sum of the correlations S(j) ( ξ/ j ξ/, j = ) is equal to.9s(). 7 The reduction of ξ with, due to the better screening of the minority sin at the QD by the majority sins on the leads, is consistent with the suression of ξ/.5.3. (L=95).88.387 (L=95).5.5..7.9.3.5 - -.5.5.5.4.3.3..48.46. ξ/ 3 5 5..5.4...3.3.. - -.5 - -.5.5.5 FIG. : (Color online) The LDOS for U =, t =.4, and ǫ =.5, for various values of the olariation, AP magnetiation for the L = 79 chain, and P, L = 95. The inset shows the length of the comensation cloud (see text) for the P case, L = 95. t =.8. For both sets of arameters the effective Kondo couling at the symmetric oint is J = 4t /U =.64, although strictly seaking this relationshi is only valid for U >> t. There are, however, many imortant differences between both sets of arameters as we will show below. The main quantity that we have studied is the local density of states (LDOS),, at the QD. In the first lace, from this quantity it is ossible to evaluate the conductance in the linear resonse regime. In the second lace, recent advances in scanning tunneling microscoy (STM) have made it ossible to directly measure this quantity. Since in our DMRG calculation the groundstate vectorismeasuredwhen the twoadded sites are at the center of the chain (symmetrical configuration), the QD then is one of these two sites which are exactly treated. Then, we adot the aroximation of alying the creation and annihilation oerators at the QD on the ground state vector and then determine following the well-known continued fraction formalism. A more accurate aroach would be, after the alication of each of those creation and annihilation oerators, to run additional swees for an enlarged density matrix. 5 In any case, the truncation of the Hilbert sace is the essential source of error in the DMRG, and to estimate the recision of our aroach we have comared results for various numbers of retained states M, from M = 5 to. We would also like to stress the fact that the conductance in linear resonse is related to the LDOS near = where the aroximation is more recise. 6 In addition, we have comuted standard roerties such as the electron occuancy of each site, < ni,σ > (σ =, ) and sin-sin correlations from the QD, S(j) =< S Sj > < S >< Sj >. All the results resented corresond to the half-filled system. For the AP case, we work in the S = / sector. For the arallel olariation, we work in the S subsace
S FIG. 3: (Color online) S i for the L = 79 chain, U = 4, t =.5 ) (ρ.5.5 ) (ρ.5 (c).5 ) (ρ.8, ǫ =., for.,.,. and.3, from bottom to to on the left lead; antiarallel, arallel configuration..5 The inset shows S at the QD as a function of the olariation - -.5.5 the Kondo resonance observed in the LDOS. Qualitatively the same behavior can be observed for FIG. 4: (Color online) The LDOS for the L = 79 chain: U = the case of the arameters U = 4, t =.8, at the symmetric oint ǫ =.. In this case, the Kondo resonance (dashed line),.5 (dotted-dashed line),.5 (dotted-dotted-, t =.4, and ǫ =-.5 (solid line), -.5 (dotted line),. dashed line) for =., =.3, AP configuration, has a smaller weight as comared with the revious set of and (c) =.3, arallel configuration. arameters. In Fig., which shows the LDOS around =, it can be seen again that for the AP case the LDOS is not much changed as is increased. In contrast, for discussed above.,3 the P case, the Kondo eak is strongly suressed in the We now study the most imortant issue in this kind of resence of olariation. This suression is considerably systems, i.e., the behavior of the conductance and hence stronger than for the revious set of arameters, U =, of the TMR as a function of the alied voltage. To understand this behavior we show in Fig. 4 the LDOS on the t =.4. We have also examined the sie deendence of these results by comuting the LDOS on the L = 95 L = 79 chain for the arameters U =, t =.4, lead olariation =.3, and various values of ǫ. In Fig. 4 we chain. Although the osition of the eaks, as exected, moves to lower frequencies (in absolute value), the overall show for comarison the unolaried case. The relevant behavior with is similar to that for L = 79. As in the quantity is the LDOS at = which is roortional to revious case, the length of the comensation cloud the conductance, G = (πt/t)ρ( = ), in units of e/h. (shown in the inset in Fig. ) for the P olariation G as a function of ǫ has the tyical form of a dome as it decreases with increasing in agreement with the suression of the Kondo resonance. It is interesting to note cluster is an artifact of the Lorentian broadening δ of the can be seen in Fig. 5. However, ρ( = ) for a finite that, consistently with this behavior, the value of (S ) eaks adoted in our calculations. Although ρ( = ) also decreases from.8 for to. for.4. could be comuted in this way9, we refer to discuss Notice that in the full Kondo regime (S ).5. the results of Fig. 4 at a qualitative level. The LDOS The mean value of S i is shown in Fig. 3 for U = 4, for =.3 in the AP olariation is shown Fig. 4. t =.8, and ǫ =., for the AP and P olariations Its deendence for various values of ǫ is very similar in anels and, resectively. The inhomogeneous to the one for the unolaried case. On the other hand, way in which the sites are olaried, due to the oen for the P configuration, the LDOS, deicted in Fig. 4(c), boundary conditions of our clusters, is clearly visible. We shows at = a nonmonotonic behavior as a function have also observed this behavior for the same arameters of the gate voltage. This nonmonotonic behavior, when on the 95-site cluster, suggesting that it is not a mere translated to the conductance, is crucial to understand finite sie effect. The different behavior of the AP and P the ǫ deendence of the TMR. configurations can be traced to the behavior of the mean Our most relevant results are shown in Fig. 5. To value of the sin at the QD. The aarently different comute the conductance we use the Friedel sum rule7, behavior of S = S QD as a function of the olariation Gσ = sin πnσ, where nσ is the occuancy at the QD of is shown in the inset. For the AP arrangement S is an electron with sin σ. The Friedel sum rule is valid for almost not affected by, while for the P arrangement arbitrary olariation in the P case, excet for negligible higher order corrections3. The total conductance is S goes to the minimum ossible value by increasing. This behavior of S in the P case is well known and G = G + G. it has been shown that a magnetic field alied to the Figure 5 shows the conductance for the unolaried QD could in fact reduce the suression of the LDOS case, AP and P olariations, for U =, t =.4, and. -. AP P -.4...3.4.5 3 S S.3.. -. -..3.. -. -. -.3-4 - i 4
/e ]h [ G.5 (c) -5-4 -3 - - FIG. 5: (Color online) Conductance as a function of the gate voltage for U =, t =.4, L = 79, and = (oen circles). Results for L = 63, =.3: P configuration (oen triangles), TMR (filled triangles), and L = 95, =.35 AP (crosses), P (lusses), and TMR (stars). Same as for U = 4, t =.8. Symbols for L = 63 and L = 95 same as. Results for L = 79: AP (squares), P (diamonds), and TMR (filled circles). The values of the olariation are =.3 (L = 63), =.3 (L = 79), and =.36 (L = 95). (c) Usin (circles) and down-sin (squares) contributions to GP for the arameters of, L = 79. Fig. 5 for U = 4, t =.8. The conductance for the AP configuration is virtually indistinguishable from the unolaried one for L = 79. On the other hand, GP has a minimum at the symmetric oint, then increases as ǫ moves aart from that oint, and finally it merges with the curve for the unolaried case, decreasing for large values of the gate voltage. Then, in the most imortant region for alication uroses, i.e. near the symmetric 3 4 5 6 7 D. Goldhaber-Gordon et al, Nature 39, 56 (998) S. J. Tans, A. R. M. Verschueren, C. Dekker, Nature 393, 49 (998); K. Tsukagoshi et al, Nature 4, 57 (999). S. A. Wolf et al, Science 94, 488 (). S. Maekawa and T. Shinjo, Sin Deendent Transort in Magnetic Nanostructures. (Taylor & Francis, London, ). A. N. Pasuathy et al, Science 36, 86 (4); J. Nygaard, et al, cond-mat/4467. N. J. Craig, et al, Science 34, 565 (4). M. J. Biercuk, et al, Nano Letters 5, 67 (5). A. J. Hewson, The Kondo roblem to heavy fermions, oint, the TMR is negative with a minimum recisely at this oint. This result is in contradiction with the one reorted in Ref. 4, showing a ositive value at the symmetric oint. This difference may be traced to the model we adoted for the leads, which imlies loosing the article-hole symmetry for the P case, which is resent in Ref. 4. Notice that for the case of Fig. 5, there is a clear trend as the chain sie is increased, that is, a deeening of GP around the symmetric oint, with a subsequent decrease in the TMR. Notice that the lead olariation is slightly different in the three clusters. This trend is not that clear for the case of Fig. 5, since here the smaller value of U requires larger lattice sies. Finally, the u- and down-sin contributions to GP are shown in Fig. 5(c) for the coulings of Fig. 5. At the symmetric oint G = G. As ǫ is moved away from the symmetric oint, the conductance of the minority sin becomes dominating, reaching erfect conductance, at ǫ U and ǫ, and it accounts for the increase of the total conductance. In summary, we have shown the ossibility of alying the DMRG technique to the roblem of transort through quantum dots. For the articular roblem of a quantum dot couled to olaried leads we have shown that the LDOS at the QD is suressed by the P olariation of the leads. In addition, we have shown that the TMR as a function of the quantum dot on-site energy is minimum and negative at the symmetric oint, a result which is at variance with those reorted in revious studies. This is a central issue in these kinds of devices. In addition to this, we have shown that this technique allowed us to estimate the length of the Kondo cloud, and to relate its behavior with the suression of the Kondo eak, correlated in turn with the conductance. This emhasies the fact that real sace roerties, which are accessible by the DMRG, and are very difficult to obtain by alternative techniques, are very imortant to understand the behavior of these devices. Acknowledgments We acknowledge useful discussions with A. E. Feiguin, S. Maekawa, J. Martinek, and G. B. Martins. This work was suorted by grant No. PICT 3-49 (ANPCYT). 8 9 3 4 (Cambridge University Press, Cambridge, UK, 993). S. R. White, Phys. Rev. Lett. 69, 863 (99); Phys. Rev. B 48, 345 (993). I. Peschel, et al, Density-Matrix Renormaliation, Lecture Notes in Physics, vol. 58 (Sringer, Berlin, 999). J. Martinek et al, Phys. Rev. Lett. 9, 73 (3). B. Dong et al, J. Phys. Condens. Matter 5, 8435 (3). I. Weynman et al, Phys. Rev. B 7, 33 (5). J. Martinek et al, Phys. Rev. Lett. 9, 47 (3); Y. Utsumi et al, cond-mat/57, (5). M. S. Choi, D. S anche, and R. L oe, Phys. Rev. Lett. G [e /h], TMR.5.5 -.5-5 -4-3 - - 4 G [e /h], TMR.5.5 -.5 - -.5 - -.5.5
5 9, 566 (4); R. Lóe and D. Sánche, Phys. Rev. Lett. 9, 66 (3). 5 K. A. Hallberg, Phys. Rev. B 5, 987 (995); E. Jeckelmann, Phys. Rev. B 66, 454 (). 6 Our calculation starts to be less recise as we move away from = since then higher excited states are involved. 7 M. E. Torio, C. J. Gaa and J. A. Riera, in rearation. 8 J. E. Gubernatis, J. E. Hirsch, and D. J. Scalaino, Phys. Rev. B 35, 8478 (987). 9 In this case δ should decrease by increasing L. See e.g., H. Benthien, F. Gebhard, E. Jeckelmann, Phys. Rev. Lett. 9, 564 (4). This has been recently shown by J. Martinek et al, Phys. Rev. B 7, 3(R) (5).