We investigate the spin- and energy dependent tunneling through a single
organic molecule (CoPc) adsorbed on a ferromagnetic Fe thin film, spatially
resolved by low-temperature spin-polarized scanning tunneling microscopy.
Interestingly, the metal ion as well as the organic ligand show a significant
spin-dependence of tunneling current flow. State-of-the-art ab initio
calculations including also van-der-Waals interactions reveal a strong
hybridization of molecular orbitals and surface 3d states. The molecule is
anionic due to a transfer of one electron, resulting in a non-magnetic (S= 0)
state. Nevertheless, tunneling through the molecule exhibits a pronounced
spin-dependence due to spin-split molecule-surface hybrid states.Comment: Version of Submission, 18-03-201

Atodiresei, Nicolae

Bluegel, Stefan

Brede, Jens

Caciuc, Vasile

Hoffmann, Germar

Kuck, Stefan

Lazic, Predrag

Morikawa, Yoshitada

Wiesendanger, Roland

English

arXiv

We investigate the spin- and energy-dependent tunneling through a single organic molecule (CoPc) adsorbed on a ferromagnetic Fe thin film, spatially resolved by low-temperature spin-polarized scanning tunneling microscopy. Interestingly, the metal ion as well as the organic ligand show a significant spin dependence of tunneling current flow. State-of-the-art ab initio calculations including also van der Waals interactions reveal a strong hybridization of molecular orbitals and substrate 3d states. The molecule is anionic due to a transfer of one electron, resulting in a nonmagnetic (S = 0) state. Nevertheless, tunneling through the molecule exhibits a pronounced spin dependence due to spin-split molecule-surface hybrid states

Brede, J.

Atodiresei, N.

Kuck, S.

Lazic, P.

Caciuc, V.

Morikawa, Y.

Hoffmann, G.

Blügel, S.

Wiesendanger, R.

Juelich Shared Electronic Resources

Spin- and Energy-Dependent Tunneling through a Single Moleculewith Intramolecular Spatial ResolutionJens Brede,1 Nicolae Atodiresei,2,3,* Stefan Kuck,1 Predrag Lazic´,2,4 Vasile Caciuc,2 Yoshitada Morikawa,3Germar Hoffmann,1,† Stefan Blu¨gel,2 and Roland Wiesendanger11Institute of Applied Physics, University of Hamburg, 20355 Hamburg, Germany2Institut fu¨r Festko¨rperforschung and Institute for Advanced Simulation, Forschungszentrum Ju¨lich, 52425 Ju¨lich, Germany3The Institute of Scientific and Industrial Research, Osaka University, Osaka, Japan4Rudjer Boskovic Institute, Zagreb, Croatia(Received 18 March 2010; published 21 July 2010)We investigate the spin- and energy-dependent tunneling through a single organic molecule (CoPc)adsorbed on a ferromagnetic Fe thin film, spatially resolved by low-temperature spin-polarized scanningtunneling microscopy. Interestingly, the metal ion as well as the organic ligand show a significant spindependence of tunneling current flow. State-of-the-art ab initio calculations including also van der Waalsinteractions reveal a strong hybridization of molecular orbitals and substrate 3d states. The molecule isanionic due to a transfer of one electron, resulting in a nonmagnetic (S ¼ 0) state. Nevertheless, tunnelingthrough the molecule exhibits a pronounced spin dependence due to spin-split molecule-surface hybridstates.DOI: 10.1103/PhysRevLett.105.047204 PACS numbers: 85.75.d, 68.37.Ef, 73.22.f, 75.50.XxMolecular based systems are fascinating, yet promisingcandidates for nanoscale spintronic devices and open via-ble routes toward quantum computing [1]. Previous experi-ments on the spin transport through break junctions [2] andspin valves [3] unveil exciting new frontiers of molecularmagnetism. Much effort is dedicated to understand theproperties of organic-magnetic interfaces. To this end,spatially averaging techniques show substantial spin injec-tion into, as well as long spin-coherent transport through-out films of metal-phthalocyanines [4,5]. However,detailed and quantitative access to different constituentsof a single molecule is desirable, though challenging.Scanning tunneling microscopy (STM) is well establishedas a probe of a local spin [6–13] in an atomically well-defined environment.Iacovita et al. recently performed a spin-polarized STM(SP-STM) study of a Co-Phthalocyanine (CoPc) in contactwith a ferromagnetic cobalt nanoisland [14]. Stackingcontrast, spin-dependent scattering, edge states, meso-scopic relaxations as well as the adsorbate induced modi-fication create a complex environment [15] toward under-standing the influence of the substrate on molecular mag-netism. After careful selection of electronically equivalentCo nanoislands a ferromagnetic exchange interaction be-tween the molecular spin and the cobalt lead was success-fully deduced, both theoretically and experimentally.In this Letter we demonstrate a significant spin polar-ization for a CoPc molecule in contact with a ferromag-netic Fe thin film due to molecule-substrate hybridizationeven though the molecule loses its net spin. As confirmedby SP-STM, an energy- and site-dependent spin polariza-tion from inversion to amplification is resolved on thesubmolecular scale. State-of-the-art density functional the-ory (DFT), which includes the decisive role of van derWaals (vdW) interactions, reveals both the magnetic andelectronic nature of the molecule coupled to the ferromag-netic substrate. Even though the net spin of the molecule islost due to a transfer of one electron, spin splitting isrecovered through the local bonding of molecular orbitalswith Fe 3d bands.Simulations were carried out in the DFT [16] formalismwith a plane wave implementation as provided by the VASPcode [17]. Pseudopotentials used were generated with theprojector augmented wave method [18] by using the PBEgeneralized-gradient exchange-correlation energy func-tional [19] (GGA). A slab consisting of two Fe and threeW atomic layers, with a (5 7) in-plane surface unit cell(22:20 A 22:42 A) modeled the molecule-surface sys-tem. The kinetic energy cutoff of the plane waves was set to500 eV while the Brillouin zone was sampled by the point. In our calculations we employed a Gaussian smear-ing with a broadening of 0.11 eV. Optimized molecule-surface geometries were obtained by relaxing all moleculardegrees of freedom and those of the Fe overlayers byincluding long-range vdW interactions in a semiempiricalway [20,21]. The threshold of the calculated forces wasfixed to 0:001 eV= A.Experiments were performed in an ultrahigh vacuumsystemwith an STM operated at 6 K [22]. In situ chromiumcoated tungsten tips with out-of-plane spin sensitivity [23]were utilized for all measurements. 1.8 ML Fe was in situdeposited on a W(110) single crystal kept at 400 K.Afterward, the magnetic sample was cooled down to mea-surement temperature and spin sensitivity was confirmedby imaging the alternating spin-up () and spin-down ()domains of the 2nd layer portion of the Fe film [23]. CoPcPRL 105, 047204 (2010) P HY S I CA L R EV I EW LE T T E R Sweek ending23 JULY 20100031-9007=10=105(4)=047204(4) 047204-1  2010 The American Physical Societymolecules were in situ thermally sublimated from a home-built crucible onto the precooled substrate to ensure spa-tially well separated molecules.All STM images were recorded in constant-currentmode at the set current I and sample bias voltage U. Thelocal spin polarization of the atomic-scale tunnel junctioncan directly be determined from height variations s re-sulting from different magnetization directions of elec-tronically identical areas as introduced by Wiesendangeret al. [24]:PtPs cos ¼ eaffiffiffips  1eaffiffiffips þ 1¼ tanha2ffiffiffiffips a2ffiffiffiffipszfflfflfflfflfflffl}|fflfflfflfflfflffl{s<1 A:The polarization of tip (Pt) and sample (Ps) are defined asthe normalized difference of the energy integrated spin-upand spin-down states (PtðsÞ :¼ tðsÞ  tðsÞ; tðsÞ þ tðsÞ ¼1).  is the angle between the respective magnetizationdirections,  the mean local tunneling barrier height, anda ffi 1 eVð1=2Þ A1.Figure 1(a) shows a representative 2nd ML area of Feafter deposition of CoPc molecules. The two dislocationlines indicate the crystallographic axes. Isolated CoPcsadsorb in three well-defined orientations (denoted as A,B, and C). The line profile [Fig. 1(b)] along the indicatedpositions in (a) quantitatively illustrates the height varia-tion within the atomically flat Fe film (s ¼ 0:3 A) andthe domain-dependent molecular appearance across theorientation A CoPc on an up and down domain. Thisvariation is due to the relative alignment of tip and samplemagnetization being either parallel or antiparallel. Whilethe center of the molecule (Co site) has the same apparentheight on both domains (s ¼ 0 A), we see a profoundcontrast at the ligand site (s ¼ 0:3 A). For an up do-main the ligand exhibits an increased apparent heightrelative to the surrounding Fe film, which is in contrastto the apparent height for a molecule on the down domain.Thereby, orientation A [25] represents a highly symmet-ric configuration (C2v symmetry) with one molecular axisoriented exactly along the h001i direction. While there aredifferences in intramolecular contrast amongst the threemolecular orientations [Figs. 1(d) to 1(i)], all three con-figurations show similar domain-dependent contrast varia-tions. In the following, we focus on molecules inconfiguration A which are not in the vicinity of dislocationlines or domain walls. As the molecular appearance con-trollably switches under an applied external magnetic field[26], the magnetic contrast can unambiguously be attrib-uted to the spin-polarized local density of states (SP-LDOS) of the molecule-surface system.First-principles calculations within DFT (GGA) clarifythe origin of the observed spin dependence. Calculationsfor the A geometry show the preferred top adsorption, i.e.,the Co atom on top of a surface Fe atom. Previous resultsfor a Co island revealed an adsorption at the bridge site[14], while the hollow site was deduced for a Au(111)substrate [12]. Furthermore, calculations cover both spin-orbit coupling (SOC) as well as the considerable role ofvdW interactions. Results without SOC are discussed interms of molecular orbitals (MOs) to clearly illustrate theeffect of vdW interactions alone.Figure 2 depicts the change in electronic structure of aCoPc with increasing molecule-substrate interaction.Figure 2(a) shows the SP-LDOS for the free molecule. Inagreement with previous work [9,12,27,28] the origin ofspin splitting is an unpaired electron in a MO with dz2contribution situated at the Co site and the total molecularspin S is 12 . Figure 2(b) presents the effect of the surface,when vdW forces are neglected during the relaxation pro-cess. The molecule adsorbs 3.1 A˚ above the surface re-maining flat. The SP-LDOS shows a hybridization of MOsand substrate 3d states. According to the spatial extensionFIG. 1 (color online). (a) SP-STM image of CoPcs on mag-netic iron domains (;  ). (b) Line profile as indicated in (a).(d)–(i) High-resolution images of CoPcs in all three orientations(A, B, C) on oppositely magnetized domains. (c) Adsorptionconfiguration of a CoPc in orientation A. (a) and (d)–(i) 0.05 V,200 pA; (a) 240 A 280 A; (d)–(i) 25 A 25 A. The colorscale for (a) is nonlinear to illustrate the domain structure, whilethe color scale for (d)–(i) is linear.PRL 105, 047204 (2010) P HY S I CA L R EV I EW LE T T E R Sweek ending23 JULY 2010047204-2of the MOs perpendicular to the molecular plane (i.e., with character), the hybridization of MOs containing the Codz2 atomic orbital is the strongest, followed by those in-cluding d and pz atomic contributions. The d and MOs are only slightly broadened, compared to the freemolecule case, as they are localized in the molecular plane.The spin splitting of the molecule-surface hybrid states isreduced due to a transfer of an electron from the substrateto the CoPc. As a result, the formerly unoccupied dz2 typeMO becomes occupied and the total molecular spin isquenched (0 b).The role of vdW forces is crucial as it brings the mole-cule 0.5 A˚ closer to the surface [Fig. 2(b)] and distorts themolecular geometry [Fig. 2(d)]. This new adsorption ge-ometry has a drastically different electronic structure dueto the overall changes of hybridization between moleculeand substrate. The SP-LDOS [Fig. 2(c)] shows not onlyMOs with a  character strongly hybridize with spin-polarized Fe 3d states of the same symmetry to form broadspin-split bands but also  type MOs are significantlyaffected by the interaction with the surface. Again, a trans-fer of one electron from the surface to the molecule anni-hilates the molecular spin, but spin splitting is recovereddue to the local molecule-surface bonding at different partsof the molecule. The newly formed molecule-surface hy-brid states have, within a given energy interval [E, E0], anunbalanced, locally varying electronic charge in the up anddown channels which is mapped in SP-STM (E ¼ EF,E0 ¼ EF þ eU).For a direct comparison of the first-principles calcula-tions with constant-current SP-STM images, isochargesurfaces above the CoPc are extracted from the spatialvariation of the energy integrated SP-LDOS. This ap-proach mimics the experimental situation of a local andspin-sensitive tip probing the charge density above thesurface at constant current and thereby, accounts for thevariation of decay lengths of the different states into vac-uum. Figure 3 gathers calculated spin-dependent isochargesurfaces and experimental SP-STM images for U ¼0:05 V. The isocharge surfaces reproduce all features ofthe experimental SP-STM images, only if vdW forces aretaken into account [Fig. 3(c)]. Corrections due to the in-clusion of spin-orbit coupling effects are minor [Fig. 3(d)].Figure 4 illustrates quantitatively the local spin polar-ization. To compensate for the tip polarization, we normal-ized the experimental polarization above the molecule(PexpCoPc) relative to the spin polarization of the free surface(PexpFe ),P ¼ PexpCoPcPexpFe¼ PtPCoPc cosPtPFe cos¼ PsimCoPcPsimFe: (1)The normalized polarization (P) is therefore, to a first ap-proximation, independent of the tip polarization and candirectly be compared to our simulated polarization (Psim).For comparison, both experimental and theoretical spa-tial maps of the polarization at three characteristic energiesare given in Fig. 4. We see an excellent agreement ofFIG. 2 (color online). The geometry and electronic structurefor a free CoPc (a) and adsorbed on an Fe surface (b) without and(c) with vdW forces included during the relaxation.(d) Molecular deformation due to vdW forces. The C6 ringsand two of the outer N are twisted by up to 0.3 A˚ toward thesurface and away from the plane defined by the Co ion and fourinner N. All H are pointing away from the surface.FIG. 3. Experimental and simulated SP-STM images at U ¼0:05 V and 22 A 22 A for both spin directions. (a) Averagedexperimental images. (b)–(d) Isocharge surfaces; (b) conven-tional DFT (GGA) (without vdW), (c) DFTþ vdW,(d) DFTþ vdW and including SOC.PRL 105, 047204 (2010) P HY S I CA L R EV I EW LE T T E R Sweek ending23 JULY 2010047204-3experimental findings and DFT calculations, when thevdW relaxed adsorption geometry is used. Remarkably,not only reduction (1>P> 0) but also inversion (P< 0)and amplification (P> 1) of spin polarization can equallybe observed. The P above the C6 ring shows an inversionfor all energies, while the spin polarization near the Co siteranges from inversion (at 0:2 V) to amplification (at0.75 V).In conclusion we investigated single CoPc moleculesadsorbed on a ferromagnetic Fe surface and locally ob-served the energy-dependent spin transport with intramo-lecular spatial resolution. The tunneling current above theentire molecule—metal ion, as well as the organic ligand—shows a high, locally varying spin polarization rangingfrom inversion up to amplification with respect to theferromagnetic Fe film. Optimized molecule-surface ge-ometries for the first-principles calculations were obtainedby including long-range van der Waals interactions duringthe relaxation, leading to an excellent agreement of theo-retical and experimental data. Therefore, the observed spinpolarization is identified as a unique property of themolecule-substrate hybrid states, by our combined DFTand SP-STM approach. Our observations imply that furtherprogress in the field of molecular spintronics does notexclusively rely on tailoring magnetic molecular proper-ties; the clever design of spin-active hybrid interfaces isequally important.This work is funded by the DFG (SFB 668-A5, SPP1243and GrK 611), the JSPS, Alexander von Humboldt foun-dation, the EU project SpiDME. The computations wereperformed on JUROPA and JUGENE supercomputers atthe Ju¨lich Supercomputing Centre, ForschungszentrumJu¨lich (Germany).*n.atodiresei@fz-juelich.de†ghoffman@physnet.uni-hamburg.de[1] L. Bogani and W. Wernsdorfer, Nature Mater. 7, 179(2008); M.N. Leuenberger and D. Loss, Nature(London) 410, 789 (2001).[2] J. Park et al., Nature (London) 417, 722 (2002).[3] A. J. Drew et al., Nature Mater. 8, 109 (2009).[4] T. Suzuki, M. Kurahashi, X. Ju, and Y. Yamauchi, J. Phys.Chem. B 106, 11 553 (2002); T. Suzuki, M. Kurahashi,and Y. Yamauchi, ibid. 106, 7643 (2002).[5] M. Cinchetti et al., Nature Mater. 8, 115 (2009).[6] F. Meier, L. Zhou, J. Wiebe, and R. Wiesendanger,Science 320, 82 (2008).[7] C. F. Hirjibehedin et al., Science 317, 1199 (2007).[8] Y. Yayon et al., Phys. Rev. Lett. 99, 067202 (2007).[9] X. Chen et al., Phys. Rev. Lett. 101, 197208 (2008); Y.-S.Fu et al., ibid. 103, 257202 (2009).[10] N. Tsukahara et al., Phys. Rev. Lett. 102, 167203(2009).[11] D. Wegner et al., Phys. Rev. Lett. 103, 087205 (2009).[12] A. Zhao et al., Science 309, 1542 (2005).[13] T. Balashov et al., Phys. Rev. Lett. 102, 257203 (2009).[14] C. Iacovita et al., Phys. Rev. Lett. 101, 116602 (2008).[15] O. Pietzsch, A. Kubetzka, M. Bode, and R. Wiesendanger,Phys. Rev. Lett. 92, 057202 (2004); O. Pietzsch et al., ibid.96, 237203 (2006); M.V. Rastei et al., ibid. 99, 246102(2007); M. Sicot et al., Phys. Rev. B 77, 035417 (2008); H.Oka et al., Science 327, 843 (2010).[16] P. Hohenberg and W. Kohn, Phys. Rev. 136, B864(1964).[17] G. Kresse and J. Hafner, Phys. Rev. B 49, 14251 (1994);G. Kresse and J. Furthmu¨ller, ibid. 54, 11 169 (1996).[18] P. E. Blo¨chl, Phys. Rev. B 50, 17 953 (1994).[19] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett.77, 3865 (1996).[20] S. Grimme, J. Comput. Chem. 27, 1787 (2006).[21] N. Atodiresei, V. Caciuc, P. Lazic´, and S. Blu¨gel, Phys.Rev. Lett. 102, 136809 (2009).[22] C. Wittneven, R. Dombrowski, S. H. Pan, and R.Wiesendanger, Rev. Sci. Instrum. 68, 3806 (1997).[23] A. Kubetzka, M. Bode, O. Pietzsch, and R. Wiesendanger,Phys. Rev. Lett. 88, 057201 (2002).[24] R. Wiesendanger et al., Phys. Rev. Lett. 65, 247 (1990).[25] Orientation A is only observed when prepared on a pre-cooled surface, indicating a metastable configuration.[26] See supplementary material at http://link.aps.org/supplemental/10.1103/PhysRevLett.105.047204 for addi-tional experimental data.[27] X. Ge et al., J. Am. Chem. Soc. 131, 6096 (2009).[28] P. A. Reynolds and B.N. Figgis, Inorg. Chem. 30, 2294(1991).FIG. 4 (color online). Local spin polarization at (a)–(c) threerepresentative energies. Raw and average (over multiple imagesand according to symmetry [26]) experimental data are com-pared with DFT simulations incl. vdW and incl. spin-orbitCoupling. Line profiles follow high-symmetry directions asindicated in the sphere model inset in (a). Simulated dataincluding SOC are only given in (b) as the SOC correctionsfor (a) and (c) are negligible. Insets: 22 A 22 A.PRL 105, 047204 (2010) P HY S I CA L R EV I EW LE T T E R Sweek ending23 JULY 2010047204-4

Spin- and Energy-Dependent Tunneling through a Single Molecule with Intramolecular Spatial Resolution

Jens Brede

Nicolae Atodiresei

Stefan Kuck

Predrag Lazić

Vasile Caciuc

Yoshitada Morikawa

Germar Hoffmann

Stefan Blügel

Roland Wiesendanger

Crossref

Physical Review Letters

Spin- and energy-dependent tunneling through a single molecule with
  intramolecular spatial resolution

Spin- and energy-dependent tunneling through a single molecule with intramolecular spatial resolution

Abstract

Similar works

Full text

Available Versions

Juelich Shared Electronic Resources

Crossref