arXiv:1102.2798v1 [cond-mat.mes-hall] 14 Feb 2011 ipysprse yculn ftesi omolecu- to spin the of One be coupling to by [8]. transport suppressed arise Kondo simply can QST-assisted QST- such SMMs a expect to might electrodes specific attached resonance to Kondo with tunneling the Combined exchange on anisotropy. the relies transverse weak term of anisotropy presence dom- magnetic a from uni-axial arising barrier inant energy the through (QST). QST spin-tunneling cou- This quantum a by induced such effect probing Kondo SMM sensitively the for of stud- pling candidate dynamic been not their One have far SMM ied. so an through 7], crystal transport [6, on SMM effect addressed distor- of been molecular has measurements such samples magnetization of effect on transport the tions a Although me- in terms demonstrated [3]. anisotropy recently setup these the of is tuning aspect chanical Hamiltonian. char- new the is further in which A terms SMM, spin-quadrupole sensitively the by of depends acterized anisotropy transport is magnetic the result the key regime on A either in [5]. that control [2], time-dependent possible is for moment ac- spin allowing the virtually of only excitation proper- [4] inelastic are arises and interaction states magnetic exchange charge spin-spin effective enhanced such the cessible, When of exhibit states [2]. can ties charge the several which instance, to For SMM and access study allows magnetism. to tunneling molecular possibilities their new control provides quan- single- with transport as interplay tum The to spin- spin-anisotropy. large referred with [4]. single, moment a setup all constitute STM (SMM), simplicity magnets an molecule for in systems, surface molecular insulating a Such trans- in an for embedded on arises atoms network situation magnetic similar through very to- A break port atoms a 3]. in such contacted [2, molecule holding junction magnetic single groups a ligand in gether provided, by is environment instance, an Such for [1]. few current to electron a respond an environment how anisotropic an of in problem atoms magnetic fundamental the address systems rnpr esrmnso aoee-ie magnetic nanometer-sized on measurements Transport neato fsi n irtosi rnpr hog sin through transport in vibrations and spin of Interaction 1 ntttfu hoeicePyi,Jhn ofagGoe Wolfgang Johann Physik, f¨ur Theoretische Institut (1) ASnmes 56.h 36.v 85.35.-p 73.63.Kv, 85.65.+h, numbers: PACS side-bands. a Kondo effect vibrational Kondo of the int dependence enhance The strongly effect. can QST-Kondo vibrations a molecular and anis (QST) transverse spin-tunneling small quantum and barrier distortio anisotropy vibrational longitudinal quantized with spin anisotropic its 3 ntttfu hoi e ttsice hsk WHAa RWTH Physik, Statistischen der f¨ur Theorie Institut (3) esuylna lcrntasottruhasingle-molecu a through transport electron linear study We spin 4 ee rubr ntttadJR udmnaso Info of Fundamentals - JARA and Gr¨unberg Institut Peter (4) ovbain saseictp of type specific a is vibrations to 2 a lnkIsiuefrPlmrRsac,518Mainz, 55128 Research, Polymer for Institute Planck Max (2) .May F. oshnsetu ¨lc,545Juih Germany J¨ulich, 52425 J¨ulich, Forschungszentrum (1 , 2) .R Wegewijs R. M. , (3 oedsotn nSMwt afitgrsi.Strik- spin. half-integer vibrational with quantized even SMM ingly, a an by distorting SMM mode an of anisotropy netic spin-isotropic for considered [9–11]. been molecules in has vibrations regime to Kondo coupling though the even not far, have so coupling studied been spin-vibration transverse longi- between and anisotropy vibrational competition tudinal the of possible the effect increase and dynamical to fluctuations the However, tends [7]. this barrier as vibrations lar oalclvbainlmd ihfeunyΩ h total The Ω. frequency with reads couples mode Hamiltonian spin vibrational the local the that a assume to in We SMM temperature virtually. the on changes low charge only the at where regime, leads blockade Coulomb electronic to coupled ln hc rnpr n unu ants a be may effects. avenue nano-mechanical magnetism with quantum possible combined and a transport indicates which vibrations and along scattering discrete en- spin by This Kondo QST of interplay SMMs. the through of transport hancement higher temperature investigation low a experimental of for generally, relevant is More spin-vibration which from coupling result fluctuate. may temperature to thereby vibra- QST-Kondo and spin by SMM, the the generated allow distort dynamically which fluctuations is conductance. tional effect linear the Kondo in arise This can peak Kondo induced ecnie eeol afitgrvle ftesi mag- spin the of nitude values half-integer only here consider We Here ei ed[] n ae omn nteitgrspin integer the on comment later and [8], field netic xs hc ecos ob the be to choose we which axis, , 4) nti ae ecnie h ouaino h mag- the of modulation the consider we paper this In oe n method. and Model so h oeue eso ht ept the despite that, show We molecule. the of ns h-nvri¨t 03 rnfr/an Germany Frankfurt/Main, the-Universit¨at, 60438 n .Hofstetter W. and , toy irtoa utain a induce can fluctuations vibrational otropy, S H z S SMM didc naoaosmgei field magnetic anomalous an induce nd H stepoeto ftemlcl’ pno t easy its on spin molecule’s the of projection the is emge SM n h nepa of interplay the and (SMM) magnet le o hc hr saKnoeeta eomag- zero at effect Kondo a is there which for , hn 25 ahn emn and Germany Aachen, 52056 chen, K ihu ttctases anisotropy transverse static without rlyo pnsatrn,QTand QST scattering, spin of erplay + = = − X kσ − DS D mto Technology, rmation ǫ H kσ z 2 ′ S + a = z 2 kσ † l-oeuemagnets gle-molecule + Germany (1) 1 2 H ecnie M strongly SMM a consider We a E SMM kσ 2 1 E S + ′ + 2 z + J S + ai and -axis S + 2 H S + · K − 2 s S with . − Ω + 2 S b ± † b b † = + QST- a , S b x ± iS (1) (2) y . 2 case. Starting from an isolated molecule, the longitudinal anisotropy D splits the eigenstates of Sz into the inverted parabolic magnetic spectrum, which is sketched in Fig. 1 for the representative case of S = 3/2 used throughout this work. The zero-field splitting (ZFS) corresponds to the energy difference between the ground-state and the first magnetic excitation, and equals δ = (2S 1)D for E = 0 (and E′ = D′ = 0). The transverse anisotropy− E breaks the continuous rotational symmetry about the easy-axis of the SMM, thereby causing spin-tunneling through the barrier. As shown in [7], vibrational modes modulate the magnetic anisotropy and can significantly FIG. 1. (Color online). Magneto-mechanical excitation contribute to the observed magnetic splittings. Here spectrum of the SMM with S = 3/2. Magnetic states, ly- we additionally consider the dynamical effects of such ing on an inverted parabola due to the static longitudinal coupling by allowing the dominant anisotropy parame- anisotropy (D), each have vibrational excitations (n) on the ters to depend linearly on the vibrational mode coor- upright parabolas. Magnetic transitions induced by Kondo † spin-scattering (J), static anisotropy (E) and longitudinal dinate Q = (b + b )/√2 through coupling coefficients ′ ′ ′ ′ † (D ) and transverse (E ) spin-vibration coupling are indi- D and E . Here the operator b (b ) relaxes (excites) cated. Magnetic states in the two different Kramer’s sub- the vibration by one quantum. Thus, when the SMM spaces are marked blue and orange. vibrates it lowers its symmetry and QST is enhanced. Importantly, this also holds for virtual quantum vibra- tions. We note that recently such a linear dependence Static anisotropy. Even without the vibration (E′ = of the D parameter on the pitch angle coordinate [12] D′ = 0) or transport (J = 0) the SMM eigenstates are in the tetra-iron (III) “propeller”-SMM used in [2] has not spin-eigenstates when E = 0. However, for the typi- been measured. The conduction electron states, repre- cal case of moderate transverse6 anisotropy E
FIG. 2. (Color online). Logarithmic color plot of the SMM FIG. 3. (Color online). Effect of the spin-vibration coupling 2 spectral function A(ω) normalized to the value A0 = π /4 on the QST-Kondo peak: SMM spectral functions shown for achieved for a symmetrically coupled quantum channel. This various combinations of zero and non-zero values of the pa- ′ ′ corresponds to the linear conductance at ω = eV and T = 0. rameters. Non-zero values used: D = 0.04D, E = 0.16D − ′ Parameters S = 3/2, D = 5 · 10 4W , E = 0.1D, E = 0, and E = 0.02D. Remaining parameters as in Fig. 2. ′ Ω = 0.5D, J = 0.2W and D is varied. The dashed line represents the renormalized ZFS from Eq. (3). E′ = 0 (black), or, E′ = 0 and D′ = 0 (blue), does not 6 both the height and width of the peak (without splitting change this result. Only in the latter case, vibrational side-peaks appear in the spectrum at ω = Ω (renormal- it). ± Dynamic anisotropy. The anisotropic couplings of the ized to smaller value due to strong J). Strikingly, when large spin to the vibrational mode D′ and E′ are how- both types of couplings are non-zero, a pronounced QST ever also of importance [7], especially if the vibrational Kondo peak appears, even though there is no transverse mode frequency energy Ω is low. The simplest effect magnetic anisotropy E = 0 (red). This vibration-induced of the longitudinal vibrational coupling D′ is a pola- QST Kondo effect is the central result of this work. We ronic shift which is different for each magnetic level. For now first explain why it requires the presence of both lon- E = E′ = 0 one can shift the vibrational coordinate Q gitudinal and transverse couplings referring to processes 2 sketched in Fig. 1. The Kondo effect is related to fluc- (or the operator b) in HSMM by an Sz -dependent amount, resulting in an effective Hamiltonian with eigenvalues tuations between degenerate states of the SMM which 2 ′2 4 EM,n = DM (D /Ω)M +Ωn, where M = S, .., S are in opposite Kramer’s subspaces [8]. Since E = 0, to and n =0−, 1, .. are− the quantum numbers of the spin− and reach states on the other side of the anisotropy barrier, a ′ the polaron, respectively. Each SMM eigenstate is thus vibration-induced spin-tunneling of type E is required, dressed by vibrational excitations as sketched in Fig. 1. which however, involves a virtual vibrational excitation. The effective static energy barrier opposing spin inver- An exchange scattering process (J), which changes the sion is changed in shape and increased in height. This Kramer’s subspace, cannot change the vibrational num- always results in an effectively enhanced ZFS δ ∆, ber at low energy, c.f. [9]. This is why the two processes → E′ and J result only in a QST Kondo side-peak in Fig. 3 3 D′2 ∆= δ + (4S(S2 S + 1) 1) . (3) which is split at Ω and suppressed due to the inability − 2 − Ω to reach the vibrational ground-state. Only when a lon- ′ Based on this simple picture one may expect that the gitudinal spin-coupling D is present as well, the virtual QST Kondo effect is suppressed by coupling to vibra- vibrational excitation can coherently reach the ground tions. For fixed E′ = 0 this is indeed the case as can be state and a full zero-bias QST Kondo anomaly can de- seen in Fig. 2: Due to the moderate but non-vanishing velop, as the red curve in Fig. 3 shows. One may say that static E = 0.1D a zero-bias QST Kondo peak (red) oc- due to the quantum-fluctuations of the vibrational mode curs which is increasingly suppressed with the coupling of the SMM, the magnetic symmetry is broken in virtual D′ due to an increased barrier. In addition two side-peaks intermediate states, allowing for Kondo exchange scatter- are found at the effective ZFS ω ∆ given by Eq. (3) ing, c.f. [10]. Therefore even for SMMs which have van- as indicated by the dashed lines≈± (and renormalized to ishing static E due to symmetry, spin-fluctuations may slightly smaller value due to strong J). These peaks are result in pronounced transport features due to the in- the remnants of the S =3/2 high-spin Kondo effect. terplay of exchange scattering and intra-molecular spin- In Fig. 3 we now focus on the dynamical effect of the vibration coupling. Clearly, the vibrational fluctuations spin-vibration coupling by first considering zero trans- can further assist the QST Kondo effect when it already verse anisotropy (E = 0). Without coupling to the vi- is present due to static E: this results in a higher Kondo brations there is no QST Kondo peak. Switching on spin- temperature as shown in the green curve of Fig. 3. vibration coupling of only one type, either D′ = 0 and In Fig. 4 we show the dependence of the QST 6 4
FIG. 4. (Color online). QST-Kondo temperature TK in units of bandwidth 2W (log color scale) determined from the NRG level flow for S = 3/2 and J = 0.2W . (a) Static QST FIG. 5. (Color online). Magnetic field evolution of the spec- Kondo effect: TK as function of the static anisotropy pa- tral function (log color scale) for the parameters of the red ′ rameters D, E without spin-vibration coupling. Dashed con- curve of Fig. 3: S = 3/2, E = 0.02D, D = 0.04D and ′ tour lines indicate that TK increases with the ratio E/D [8]. E = 0.16D. Dashed lines mark the QST-Kondo (red) and (b) vibration-induced QST Kondo effect: TK as function of ′ ′ vibrational side-peak (yellow) both evolving with an anoma- spin-vibration couplings D and E with Ω = 0.5D for static lous g-factor and zero-field split high-spin Kondo peak (blue). anisotropy, D = 5·10−4W and E = 0. The dashed line marks ′ ′ the regime D > 0.4D where D starts to suppress TK , note the constant D here and also the offset: TK = 0 for either ′ ′ E = 0 or D = 0. Results for small finite E are similar. instead a zero-bias conductance dip (for 0 Kondo temperature TK on the longitudinal and trans- as studied here could occur. verse static anisotropies and their respective vibration- Conclusion. We have studied the interplay of spin and ′ ′ couplings D and E . TK grows as both the transverse vibration on the conductance through a single-molecule anisotropies are increased since they both enhance QST. magnet. Whereas longitudinal coupling to the vibra- Whereas D always suppresses the QST Kondo effect, tion increases the zero-field splitting, suppressing the Fig. 4(a), its fluctuations first enhance TK by allow- quantum spin-tunneling Kondo peak, a vibrationally in- ing QST between the Kramers degenerate ground-states. duced quantum spin-tunneling Kondo effect can occur at Eventually they will suppress TK if the vibrational con- zero bias if transverse coupling is present as well. The tribution to the anisotropy barrier S4D′2/Ω, protecting transition to virtual vibrational excited states and the the SMM against these processes, increases too much. transverse spin-mixing in these virtual states results in Finally we show in Fig. 5 how the experimentally ac- a Kondo effect, even in the absence of static transverse cessible evolution of the conductance with a magnetic anisotropy. The interplay with vibrations thus can in- field reveals the different origin of the various peaks. For crease the quantum spin-tunneling Kondo temperature simplicity, we consider the field to be along the easy for a given static anisotropy, which may motivate further axis as in some experiments [2, 3] and added the Zee- experimental investigation of low temperature transport though single-molecule magnets. The measurable mag- man term SzHz to HSMM in Eq. (1) absorbing the g- factor into− the magnetic field. As the magnetic field is netic field evolution of the conductance reveals that vi- increased, the QST Kondo peak is weakened and splits brational side-bands acquire an anomalous g-factor. We acknowledge A. Cornia, J. Kortus and J. Paaske for stim- with the anomalous g-factor, ωKondo =2SHz [18]. This clearly indicates the origin of the| QST| Kondo effect, since ulating discussions and support from NanoSci-ERA. the ground-state Kramer’s doublet M = S is split by ∆M = 2S. Strikingly, the vibrational side-peaks± have the same,± strong field dependence, as they correspond to a similar transition offset in energy by Ω: ω = | vib| [1] L. Bogani and W. Wernsdorfer, Nature Materials 7, 179 Ω+2SHz. In contrast to this, the high-spin Kondo peak (2008). (ZFS) evolves much slower in the magnetic field, inde- [2] A. S. Zyazin, J. W. van den Berg, E. A. Osorio, H. S. pendent of the spin magnitude S: ωZFS = ∆+ Hz, van der Zant, N. P. Konstantinidis, F. May, M. Leijnse, signaling that it corresponds only to| a transition| with W. Hofstetter, M. R. Wegewijs, C. Danieli, et al., Nano ∆M = 1. Comparing with the above formulas for the Lett. 10, 3307 (2010). peak evolution,± indicated by dashed lines in Fig. 5, we [3] J. J. Parks, A. R. Champagne, T. A. Costi, W. W. Shum, A. N. Pasupathy, E. Neuscamman, S. Flores-Torres, P. S. conclude that the QST Kondo and high-spin Kondo ef- Cornaglia, A. A. Aligia, C. A. Balseiro, et al., Science fects are distinguishable, especially for SMM with large 328, 1370 (2010). spin. Finally, we remark that for an SMM with integer [4] A. F. Otte, M. Ternes, K. von Bergmann, S. Loth, spin S there is no QST Kondo effect at zero field, but H. Brune, C. P. Lutz, C. F. Hirjibehedin, and A. J. Hein- 5 rich, Nature Physics 4, 847 (2008). Rev. B 74, 235409 (2006). A. F. Otte, M. Ternes, S. Loth, C. P. Lutz, C. F. Hir- P. S. Cornaglia, G. Usaj, and C. A. Balseiro, Phys. Rev. jibehedin, and A. J. Heinrich, Phys. Rev. Lett. 103, B 76, 241403 (2007). 107203 (2009). [12] L. Gregoli, C. Danieli, A.-L. Barra, P. Neugebauer, [5] S. Loth, K. von Bergmann, M. Ternes, A. F. Otte, G. Pellegrino, G. Poneti, R. Sessoli, and A. Cornia, C. Lutz, and A. Heinrich, Nature Physics 6 (2010). Chem. Eur. J. 15, 6456 (2009). [6] M. N. Leuenberger and D. Loss, Phys. Rev. B 61, 1286 [13] G. Gonz´alez, M. N. Leuenberger, and E. R. Mucciolo, (2000). Phys. Rev. B 78, 054445 (2008). T. Pohjola and H. Schoeller, Phys. Rev. B 62, 15026 [14] A. A. Aligia, C. A. Balseiro, and C. R. Proetto, Phys. (2000). Rev. B 33, 6476 (1986). [7] J. Kortus, T. Baruah, N. Bernstein, and M. R. Pederson, H. Lustfeld, Physica B+C 106, 213 (1981). Phys. Rev. B 66, 092403 (2002). [15] T. A. Costi, Phys. Rev. Lett. 85, 1504 (2000). [8] C. Romeike, M. R. Wegewijs, W. Hofstetter, and [16] R. Bulla, T. A. Costi, and T. Pruschke, Rev. Mod. Phys. H. Schoeller, Phys. Rev. Lett. 96, 196601 (2006). 80, 395 (2008). [9] J. Paaske and K. Flensberg, Phys. Rev. Lett. 94, 176801 [17] W. Koller, A. C. Hewson, and D. Meyer, Phys. Rev. B (2005). 72, 045117 (2005). [10] K. Kikoin, M. N. Kiselev, and M. R. Wegewijs, Phys. [18] C. Romeike, M. R. Wegewijs, W. Hofstetter, and Rev. Lett. 96, 176801 (2006). H. Schoeller, Phys. Rev. Lett. 106, 019902(E) (2011). [11] C. A. Balseiro, P. S. Cornaglia, and D. R. Grempel, Phys. [19] M. N. Leuenberger and E. R. Mucciolo, Phys. Rev. Lett. 97, 126601 (2006).