Spin Control Induced by Molecular Charging in a Transport Junction

Dec 12, 2017 - Carlos García,. § ... Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián, Spain. ...
2 downloads 0 Views 3MB Size
Letter pubs.acs.org/NanoLett

Cite This: Nano Lett. 2018, 18, 88−93

Spin Control Induced by Molecular Charging in a Transport Junction Sujoy Karan,*,†,‡ Carlos García,§ Michael Karolak,⊥ David Jacob,¶,∥ Nicolás Lorente,§,# and Richard Berndt† †

Institut für Experimentelle und Angewandte Physik, Christian-Albrechts-Universität zu Kiel, 24098 Kiel, Germany Institute of Experimental and Applied Physics, University of Regensburg, 93053 Regensburg, Germany § Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián, Spain ⊥ Institut für Theoretische Physik und Astrophysik, Universität Würzburg, Am Hubland, 97074 Würzburg, Germany ¶ Departamento de Física de Materiales, Universidad del País Vasco, UPV/EHU, Av. Tolosa 72, 20018 San Sebastián, Spain ∥ IKERBASQUE, Basque Foundation for Science, Maria Diaz de Haro 3, 48013 Bilbao, Spain # Centro de Física de Materiales CFM/MPC, CSIC-UPV/EHU, Paseo Manuel de Lardizabal 5, 20018 Donostia-San Sebastián, Spain

Nano Lett. 2018.18:88-93. Downloaded from pubs.acs.org by TULANE UNIV on 01/21/19. For personal use only.



S Supporting Information *

ABSTRACT: The ability of molecules to maintain magnetic multistability in nanoscalejunctions will determine their role in downsizing spintronic devices. While spin-injection from ferromagnetic leads gives rise to magnetoresistance in metallic nanocontacts, nonmagnetic leads probing the magnetic states of the junction itself have been considered as an alternative. Extending this experimental approach to molecular junctions, which are sensitive to chemical parameters, we demonstrate that the electron affinity of a molecule decisively influences its spin transport. We use a scanning tunneling microscope to trap a meso-substituted iron porphyrin, putting the iron center in an environment that provides control of its charge and spin states. A large electron affinity of peripheral ligands is shown to enable switching of the molecular S = 1 ground state found at low electron density to S = 1/2 at high density, while lower affinity keeps the molecule inactive to spin-state transition. These results pave the way for spin control using chemical design and electrical means. KEYWORDS: Molecular spintronics, electronic transport, spin excitations, electron affinity, single molecule, STM

T

he electron transport through molecular junctions is usually determined by the molecule-electrode contact geometry and the corresponding hybridization of molecular orbitals with the electrodes.1−5 Ferromagnetic electrodes, that provide spin-polarized current, are conducive to a stronger hybridization,6−8 which can even lead to the quenching of molecular spin in the junction.8 If the spin is not quenched, the coupling of the molecular spin to the electrode magnetization reduces its relaxation time. This limitation does not appear to hold when nonmagnetic leads are used. There have been examples demonstrating the detection and control of magnetic anisotropy and spin state in two-terminal device geometry9−15 as well as in the junction of a three-terminal settings where electrical gating by a third terminal was used for stabilizing a particular state.16−19 Structural changes may also affect molecular spin states or their conductance fingerprints.20−22 Here we show that the polarizability of molecular ligands can provide control of the magnetic state in different transport regimes with no electrical gating required. We built junctions with R4-type iron porphyrins (FeTPP, Figure 1), involving different meso sub-substituents, namely, sulfonyl hydroxide (−SO3H) and methoxy (−OCH3), in a scanning tunneling microscope (STM). Figure 2 summarizes the results of tetraphenylporphine sulfonate (denoted FeTPP-S © 2017 American Chemical Society

Figure 1. Schematic of meso-substituted porphyrins. Molecules host an Fe ion at the center of tetraphenylporphyrin with four meso-subgroups occupying para positions of the phenyl rings. H atoms are not shown. Phenyl rings are substituted with sulfonyl hydroxide (−SO3H) or methoxy (−OCH3) groups. Pink and gray colors are used in STM images below for these molecules. Directions connecting opposite pyrroles correspond to c2 axes.

below). When deposited (see Methods) on Au(111), the molecules appear to lie flat featuring a quatrefoil shape in STM Received: August 9, 2017 Revised: November 17, 2017 Published: December 12, 2017 88

DOI: 10.1021/acs.nanolett.7b03411 Nano Lett. 2018, 18, 88−93

Letter

Nano Letters

Figure 2. (a) Constant−current image (I = 100 pA, sample voltage V = 2.56 V, 3.4 × 4.6 nm2) of FeTPP-S on Au(111). (b) Conductance measured while the STM tip is brought closer to the center of a FeTPP-S molecule. Positive displacements indicate a reduced tip-molecule distance. (c) Spectra of the differential conductance dI/dV recorded at displacements indicated with stars in panel b. The inset schematically shows the experimental setting. Zero displacement (Δz = 0) corresponds to I = 100 pA and V = 0.1 V. The lower spectra (black curves) have been multiplied by 2000, 170, and 17, respectively. (d) dI/dV data recorded at a molecular lobe. The lower three spectra have been scaled by factors of 5000, 300, and 20, respectively.

a Kondo temperature about 96 K (this value41 may differ from the actual Kondo temperature42) in good agreement with previous reports on similar molecules43. This large value translates into a field of exceeding 15 Tesla to split the Kondo peak, unattainable in STM setups. The appearance of the Kondo feature would require a change of either the magnetic anisotropy or the spin state of the molecule to have a degenerate ground state. The data tentatively rule out the former scenario, since the transition from two excitation steps to a Kondo-dip is abrupt, and no intermediate regime with a splitting of the Kondo feature was observed, in contrast to earlier results from a Co-terpyridine derivative.44 Our calculations suggest that it is indeed a S = 1/2 ground state (details below), ruling out inelastic spin flips by virtue of Kramers theorem that forces a degenerate S = 1/2 ground state. When lobes are contacted, the spectrum remains virtually unchanged (Figure 2d). On the lobes, where the tip mainly couples to the molecular π system, the inelastic steps appear more symmetrical around zero-bias than at the center of the molecule, where the d orbitals of Fe, particularly dz2, prevail. We therefore moved to less polarizable methoxy (−OCH3) groups as meso-substituents to investigate the influence of the peripheral ligands on the Kondo physics, while keeping the environment of the Fe atom unchanged. The molecules (denoted FeTPP-Me below) agglomerate into small domains (Figure 3a) along ⟨112̅⟩ and ⟨11̅0⟩ directions starting from fcc regions of the herringbone reconstruction. Individual molecules are imaged as “crabs” with “claws” representing the methoxy groups (Figure 3b). Both tunneling and contact regimes were explored under experimental conditions similar to those used before. Data from a molecule with protruding opposite pyrrole

topographs (Figure 2a) with one of their c2 axes preferentially oriented along ⟨11̅0⟩ direction of the substrate. However, deviations up to ≈ ± 15° from this orientation are observed. Figure 2b shows a conductance-displacement curve acquired while bringing the STM tip closer to the center of a molecule until contact was reached. The decreasing width of the tunneling barrier is reflected by an exponential increase of the conductance followed by an abrupt rise indicating the transition from tunneling to contact. Next, differential conductance (dI/ dV) spectra were recorded while holding the tip at predetermined heights. Spectra acquired in the tunneling regime (Figure 2c, black curves) exhibit steps at about ±7 meV. We attribute these features to inelastic spin-flip excitations, similar to those of iron porphyrin and octaethylporphyrin,23,24 which are consistent with our calculations revealing a S = 1 state of the molecule exhibiting magnetic anisotropy (details further). The possibility that these features are due to vibrational excitation may be excluded because the changes of conductance exceed several tens of percent that are typical of magnetic excitations as shown for similar molecular systems.23,25−28 In contrast, the largest conductance changes caused by vibration excitation rarely reach 20% (in the cases of Ce−H29,30 and CO-molecular cascades31) because of partial cancellation of the changes of the elastic and inelastic conductances.32 Once the tip is brought into contact with the iron center, the spectrum (light blue) develops a dip as typical of the Kondo effect.33,34 Further approach of the tip does not change the shape of the spectrum until destruction of the junction occurs. The spectrum may be fit by a Frota line35,36 (red curve) as expected for the spin-1/2 Kondo effect.37−40 Our fit points out 89

DOI: 10.1021/acs.nanolett.7b03411 Nano Lett. 2018, 18, 88−93

Letter

Nano Letters

out-of-plane excited state. This interpretation coincides with our magnetic-anisotropy energy (MAE) calculations (Supporting Information, section 2), which show the easy-plane character, that is, Sz = 0 ground state, of the S = 1 system. The contact regime is modeled using the above system plus an FCC pyramid of gold atoms grown along the crystallographic [111] direction (see Methods). The apex atom is placed at 2.35 Å from the Fe atom for both molecules. During the relaxation of forces on the tip apex and the molecule, the Fe atom is displaced toward the surface, until the vertical distance with the surface plane is 2.40 Å. These values correspond to typical Fe−Au covalent-bond distances and represent a strong, nondestructive interaction with the metallic tip. This simulation setup corresponds to the large measured contact conductance. Despite the similarity in molecular structure, particularly in the frontier orbitals and spin distribution calculated for the tunneling range, the contact regime shows qualitative differences. The presence of the tip brings in enough electrons and a change in the electrostatic balance for the FeTPP-S molecule, while it basically leaves FeTPP-Me unaltered. Indeed, FeTPP-S takes a full electron, while FeTPP-Me does not. The electron is largely placed in the C and N atoms next to Fe. In a first approximation, FeTPP-S looks much like its anion, Figure 4e and f, where the Fe core holds a magnetic moment of 2 Bohr magnetons (S = 1), while the C and N rings polarize their extra electron to create an antiferromagnetic coupling. Consequently, the anion is a doublet (S = 1/2). For the FeTPP-S trapped between tip and substrate, the Fe atom captures charge (∼0.2e, Table 1), reducing its magnetic moment. The charge on C and N carries significant magnetic moment, oriented opposite to the total moment, such that the molecule also becomes a doublet in contact. However, FeTPP-Me captures less charge with a consequently smaller reduction of its magnetic moment, Table 1. These results show that the valence state of the Fe ion remains unchanged regardless of the electrode configuration, despite the different charge transfer to the ligands. The ligands behave as “innocent” ligands.47 Calculations for a tip-apex centered 2.45 Å above a benzene ring of FeTPP-S do not show significant changes of the molecular states except for those caused by the new constraint to the geometry of the molecule (Supporting Information, section 4). The magnetic moment (Table 1, last row) is hardly changed from the tunneling case, although there is a small charge transfer comparable to the one of the FeTPP-Me molecule in the contact range. These results fully agree with unaltered excitation steps observed experimentally. The main effect of the tip contacting the Fe atom is to force an extra electron into the FeTPP-S molecule, reducing its spin to S = 1/2, while FeTPP-Me maintains its S = 1 configuration.

Figure 3. (a) STM image (11.5 × 11.4 nm2, I = 100 pA, V = 0.05 V) of a square-shaped island of FeTPP-Me molecules on Au(111). (b) Image of a FeTPP-Me chain. (c) Conductance measured during tip approach to the molecular center. (d) dI/dV spectra recorded at displacements indicated with colored stars in panel c. The upper spectra (black) have been multiplied by 4500, 275, and 65, respectively.

units in the ⟨110̅ ⟩ direction are shown in Figure 3d. The tunneling characteristics of the molecule are virtually identical to those of FeTPP-S. However, in the contact regime, the spectrum retains the inelastic steps. These results are remarkable raising the question how the proximity of metallic leads influences the magnetic state of the molecule and which chemical parameter becomes predominant at contact. Density functional theory (DFT) calculations give a clear picture of the different conductance regimes. In the tunneling regime, FeTPP-Me and FeTPP-S have similar electronic structures and bind in the same way to Au(111). The interaction is basically driven by dispersion forces and negligible charge transfer takes place (see ΔQ in Table 1).45 As a consequence, the spin distribution over the molecule is hardly perturbed by the substrate. Most of the spin is localized at the Fe atom, Table 1 and red isocontour in Figures 4a−c, with some small redistribution of spin in the adjacent C and N atoms. The molecule is in a S = 1 ground state, in good agreement with a Fe(II) configuration of the molecule.45,46 The experimental conductance shows the signature of an excitation from the ground state spin lying in the molecular plane to an

Table 1. Magnetic Moment for Full System and Fe Atom Using Bader Charges, Fe Charge Difference with Respect to Gas-Phase Molecule, and Affinity Energy as the Difference of Total Energy between Anion and Neutral Gas-Phase Moleculea tunneling −OCH3 −SO3H −SO3H lobe

contact

μTotal (μB)

μFe (μB)

ΔQ (e)

μTotal (μB)

μFe (μB)

ΔQ (e)

affinity (eV)

2.05 1.93

2.28 2.15

−0.013 0.012

1.67 0.99 1.82

2.00 1.54 2.08

0.11 0.19 0.12

−2.13 −3.00 −3.00

Magnetic moment in Bohr magnetons, μB; full system, μTotal; Fe atom, μFe; ΔQ of Fe, in electrons, positive values mean the adsorbed molecule is negatively charged with respect to the gas-phase one. For the two first cases, the contact regime is simulated by locating the pyramidal tip apex at 2.35 Å above the Fe atom. On the lobes of the FeTPP-S molecule, the tip is centered 2.45 Å above a benzene ring (the tunneling regime is not considered since the calculation does not include any tip in the tunneling configuration). a

90

DOI: 10.1021/acs.nanolett.7b03411 Nano Lett. 2018, 18, 88−93

Letter

Nano Letters

Figure 4. (a−c) Calculated spin-density isosurfaces (0.01 e/Å3) for (a) FeTPP-Me, (b) FeTPP-S on Au(111), and (c) gas-phase FeTPP-S. Red and yellow indicate spin directions. These calculations represent the system for a tunneling junction where the STM tip does not perturb the surface. The three systems show qualitative resemblance. (d, e) Corresponding results for the “contact” situation, calculated with a [111] pyramid of 10 Au atoms placed above the Fe atom. The pyramid is depicted with smaller gray balls in order not to obscure the molecule. (f) Negatively charged gas-phase FeTPP-S molecule, which resembles panel e, indicating a larger charge transfer for FeTPP-S than for FeTPP-Me at contact.

This is clearly seen in the experiment, where the differential conductance displays a zero-bias anomaly for FeTPP-S characteristic of a S = 1/2 Kondo effect. However, FeTPP-Me maintains unaltered spin transitions at the same energy, reflecting an unchanged S = 1 configuration. The qualitative differences between the molecules despite their electronic similarity can be traced back to their different electron affinities, which result from the different polarizabilities of the −OCH3 and −SO3H substituents. We evaluate the affinity as the energy difference between the ground states of the anion and the neutral molecule. As expected, the −SO3H group is highly polarizable, more than the −OCH3 group, and the affinity of FeTPP-S is 0.87 eV larger. Consequently, FeTPPS easily captures and stabilizes one extra electron in an environment of increased electron density. In conclusion, the tunneling conductance spectra of FeTPPMe and FeTPP-S show magnetic excitations corresponding to a S = 1 system with easy-plane anisotropy. When the tip contacts the Fe core of the molecules, FeTPP-Me shows unchanged magnetic structure (S = 1), but FeTPP-S displays a zero-bias anomaly typical of a S = 1/2 system. The change from S = 1 to 1 /2 in the case of FeTPP-S is due to more polarizable end groups resulting in larger electron affinity, which stabilizes an electron from the electron-rich metallic environment at contact. The lower electron affinity of FeTPP-Me however preserves its spin even at contact. The main observation is that the observed Kondo features of FeTPP-S are only compatible with a S < 1 state of the molecule in good agreement with the computed value of its total spin. This work shows that a degree of control on the molecular magnetism can be achieved by judiciously choosing their ligands and their electronic environment. Methods. Sample Preparation. Purified iron meso-tetra(4sulfonatophenyl)porphine chloride and meso-tetra(4methoxyphenyl)porphine chloride were dissolved separately in aqueous methanol in the presence of 1 vol % acetic acid. Singly charged molecular cations originated from the

dissolution of Cl ions were electrosprayed onto clean Au(111) surfaces in ultrahigh vacuum.48 STM. Imaging and contact experiments were performed using a scanning tunneling microscope operated at ∼5 K. An electrochemically etched tungsten wire, indented in situ into the substrate to get it coated with gold, was used as the tip. The differential conductance spectra were recorded adding a sinusoidal modulation (2 mVrms, 1.2 kHz) to the bias. DFT. Ab initio calculations have been performed in the framework of density functional theory (DFT) as implemented in the VASP code (details in the Supporting Information, section 1). The molecules were relaxed on a three-layer slab with a 10 × 10 Au (111) unit cell. To model the tip, a [111] pyramid of 10 atoms was chosen. Because of the lack of coordination, the DFT electronic structure predicts magnetic moments on the edge atoms. However, when the tip contacts the molecule, the magnetism appears only in the apex atom on top of the Fe atom. This shows that, in our calculations, the pyramid is interacts with the substrate and does not develop magnetism contrary to the free-standing pyramid. The calculation of molecular spins within DFT can be problematic. Typically spin contamination is a minor problem,49 but the total spin is not a good quantum number in the monoconfigurational implementations of DFT. However, Sz is a good quantum number. Fixing different Sz solutions in the absence of spin−orbit coupling (hence keeping spin degeneracies) leads to different total energies, showing that the only possible ground state compatible with Sz = 1/2 corresponds to the antiferromagnetic coupling of the ligand and Fe spins.



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.nanolett.7b03411. Density functional theory calculations; magnetic anisotropy; total spin-dependent charge on adsorbed mole91

DOI: 10.1021/acs.nanolett.7b03411 Nano Lett. 2018, 18, 88−93

Letter

Nano Letters



(18) Osorio, E. A.; Moth-Poulsen, K.; van der Zant, H. S. J.; Paaske, J.; Hedegård, P.; Flensberg, K.; Bendix, J.; Bjørnholm, T. Nano Lett. 2010, 10, 105−110. (19) Vela, S.; Verot, M.; Fromager, E.; Robert, V. J. Chem. Phys. 2017, 146, 064112. (20) Karan, S.; Jacob, D.; Karolak, M.; Hamann, C.; Wang, Y.; Weismann, A.; Lichtenstein, A. I.; Berndt, R. Phys. Rev. Lett. 2015, 115, 016802. (21) Hiraoka, R.; Minamitani, E.; Arafune, R.; Tsukahara, N.; Watanabe, S.; Kawai, M.; Takagi, N. Nat. Commun. 2017, 8, 16012. (22) Kuang, G.; Zhang, Q.; Lin, T.; Pang, R.; Shi, X.; Xu, H.; Lin, N. ACS Nano 2017, 11, 6295−6300. (23) Heinrich, B. W.; Ahmadi, G.; Müller, V. L.; Braun, L.; Pascual, J. I.; Franke, K. J. Nano Lett. 2013, 13, 4840−4843. (24) Rubio-Verdú, C.; Sarasola, A.; Choi, D.-J.; Majzik, Z.; Ebeling, R.; Calvo, M. R.; Ugeda, M. M.; Garcia-Lekue, A.; Sánchez-Portal, D.; Pascual, J. I. arXiv:1708.01268v1 2017. (25) Chen, X.; Fu, Y.-S.; Ji, S.-H.; Zhang, T.; Cheng, P.; Ma, X.-C.; Zou, X.-L.; Duan, W.-H.; Jia, J.-F.; Xue, Q.-K. Phys. Rev. Lett. 2008, 101, 197208. (26) Tsukahara, N.; Noto, K.; Ohara, M.; Shiraki, S.; Takagi, N.; Takata, Y.; Miyawaki, J.; Taguchi, M.; Chainani, A.; Shin, S.; Kawai, M. Phys. Rev. Lett. 2009, 102, 167203. (27) Gauyacq, J.-P.; Novaes, F. D.; Lorente, N. Phys. Rev. B: Condens. Matter Mater. Phys. 2010, 81, 165423. (28) Gauyacq, J.-P.; Lorente, N.; Novaes, F. D. Prog. Surf. Sci. 2012, 87, 63. (29) Pivetta, M.; Ternes, M.; Patthey, F.; Schneider, W. D. Phys. Rev. Lett. 2007, 99, 126104. (30) Hofer, W. A.; Teobaldi, G.; Lorente, N. Nanotechnology 2008, 19, 305701. (31) Heinrich, A. J.; Lutz, C. P.; Gupta, J. A.; Eigler, D. M. Science 2002, 298, 1381. (32) Morgenstern, K.; Lorente, N.; Rieder, K.-H. Phys. Status Solidi B 2013, 250, 1671. (33) Schneider, W.-D.; Berndt, R. J. Electron Spectrosc. Relat. Phenom. 2000, 109, 19−31. (34) Scott, G. D.; Natelson, D. ACS Nano 2010, 4, 3560−3579. (35) Frota, H. O. Phys. Rev. B: Condens. Matter Mater. Phys. 1992, 45, 1096. (36) Prüser, H.; Wenderoth, M.; Weismann, A.; Ulbrich, R. G. Phys. Rev. Lett. 2012, 108, 166604. (37) Prüser, H.; Wenderoth, M.; Dargel, P. E.; Weismann, A.; Peters, R.; Pruschke, T.; Ulbrich, R. G. Nat. Phys. 2011, 7, 203. (38) Zhang, Y.; Kahle, S.; Herden, T.; Mayor, C. S. M.; Schlickum, U.; Ternes, M.; Wahl, P.; Kern, K.; et al. Nat. Commun. 2013, 4, 2110. (39) Karan, S.; Li, N.; Zhang, Y.; He, Y.; Hong, I.-P.; Song, H.; Lü, J.T.; Wang, Y.; Peng, L.; Wu, K.; Michelitsch, G. S.; Maurer, R. J.; Diller, K.; Reuter, K.; Weismann, A.; Berndt, R. Phys. Rev. Lett. 2016, 116, 027201. (40) Karan, S.; Berndt, R. Phys. Chem. Chem. Phys. 2016, 18, 9334. (41) Different definitions of the Kondo temperature exist. For the comparability to the Kondo temperatures (TK) extracted from previous STM experiments we use the definition relating TK to the half width at half maximum of the symmetric resonance. This gives kBTK = 2.54Γ, where kB is the Boltzmann constant and the width Γ of the Frota line is defined in ref 36. (42) Zhang, Y.-h.; Kahle, S.; Herden, T.; Stroh, C.; Mayor, M.; Schlickum, U.; Ternes, M.; Wahl, P.; Kern, K. Nat. Commun. 2013, 4, 2110. (43) Iancu, V.; Deshpande, A.; Hla, S.-W. Phys. Rev. Lett. 2006, 97, 266603. (44) Parks, J. J.; Champagne, A. R.; Costi, T. A.; Shum, W. W.; Pasupathy, A. N.; Neuscamman, E.; Flores-Torres, S.; Cornaglia, P. S.; Aligia, A. A.; Balseiro, C. A.; Chan, G. K.-L.; Abruna, H. D.; Ralph, D. C. Science 2010, 328, 1370−1373. (45) Gottfried, J. M. Surf. Sci. Rep. 2015, 70, 259. (46) Borgogno, A.; Rastrelli, F.; Bagno, A. Dalton Trans. 2014, 43, 9486.

cule; tip contacting FeTPP-S lobe; spatial evolution of dI/dV spectra; dI/dV spectra of molecules over wider range of sample voltages (PDF)

AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. ORCID

Sujoy Karan: 0000-0002-1319-7733 Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS Financial support of the Deutsche Forschungsgemeinschaft (DFG) via Sonderforschungsbereich 677 and Forschergruppe 1162, MINECO (Grant No. MAT2015-66888-C3-2-R) and FEDER funds is acknowledged.



REFERENCES

(1) Kiguchi, M.; Tal, O.; Wohlthat, S.; Pauly, F.; Krieger, M.; Djukic, D.; Cuevas, J. C.; van Ruitenbeek, J. M. Phys. Rev. Lett. 2008, 101, 046801. (2) Néel, N.; Kröger, J.; Limot, L.; Frederiksen, T.; Brandbyge, M.; Berndt, R. Phys. Rev. Lett. 2007, 98, 065502. (3) Quek, S. Y.; Kamenetska, M.; Steigerwald, M. L.; Choi, H. J.; Louie, S. G.; Hybertsen, M. S.; Neaton, J. B.; Venkataraman, L. Nat. Nanotechnol. 2009, 4, 230−234. (4) Schull, G.; Frederiksen, T.; Arnau, A.; Sánchez-Portal, D.; Berndt, R. Nat. Nanotechnol. 2011, 6, 23−27. (5) Requist, R.; Baruselli, P. P.; Smogunov, A.; Fabrizio, M.; Modesti, S.; Tosatti, E. Nat. Nanotechnol. 2016, 11, 499. (6) Kawahara, S. L.; Lagoute, J.; Repain, V.; Chacon, C.; Girard, Y.; Rousset, S.; Smogunov, A.; Barreteau, C. Nano Lett. 2012, 12, 4558− 4563. (7) Schmaus, S.; Bagrets, A.; Nahas, Y.; Yamada, T. K.; Bork, A.; Bowen, M.; Beaurepaire, E.; Evers, F.; Wulfhekel, W. Nat. Nanotechnol. 2011, 6, 185−189. (8) Brede, J.; Atodiresei, N.; Kuck, S.; Lazić, P.; Caciuc, V.; Morikawa, Y.; Hoffmann, G.; Blügel, S.; Wiesendanger, R. Phys. Rev. Lett. 2010, 105, 047204. (9) Gopakumar, T. G.; Matino, F.; Naggert, H.; Bannwarth, A.; Tuczek, F.; Berndt, R. Angew. Chem., Int. Ed. 2012, 51, 6262. (10) Miyamachi, T.; Gruber, M.; Davesne, V.; Bowen, M.; Boukari, S.; Joly, L.; Scheurer, F.; Rogez, G.; Yamada, T. K.; Ohresser, P.; Beaurepaire, E.; Wulfhekel, W. Nat. Commun. 2012, 3, 938. (11) Warner, B.; Oberg, J. C.; Gill, T. G.; El Hallak, F.; Hirjibehedin, C. F.; Serri, M.; Heutz, S.; Arrio, M.-A.; Sainctavit, P.; Mannini, M.; Poneti, G.; Sessoli, R.; Rosa, P. J. Phys. Chem. Lett. 2013, 4, 1546. (12) Ruiz, E. Phys. Chem. Chem. Phys. 2014, 16, 14. (13) Bairagi, K.; Bellec, O. I. A.; Kartsev, A.; Li, D.; Lagoute, J.; Chacon, C.; Girard, Y.; Rousset, S.; nad Yannick J Dappe, F. M.; Smogunov, A.; Barreteau, C.; Boillot, M.-L.; Mallah, T.; Repain, V.; et al. Nat. Commun. 2016, 7, 12212. (14) Jasper-Tönnies, T.; Gruber, M.; Karan, S.; Jacob, H.; Tuczek, F.; Berndt, R. J. Phys. Chem. Lett. 2017, 8, 1569. (15) Jasper-Tönnies, T.; Gruber, M.; Karan, S.; Jacob, H.; Tuczek, F.; Berndt, R. Nano Lett. 2017, 17, 6613. (16) Jo, M.-H.; Grose, J. E.; Baheti, K.; Deshmukh, M. M.; Sokol, J. J.; Rumberger, E. M.; Hendrickson, D. N.; Long, J. R.; Park, H.; Ralph, D. C. Nano Lett. 2006, 6, 2014−2020. (17) Zyazin, A. S.; van den Berg, J. W. G.; Osorio, E. A.; van der Zant, H. S. J.; Konstantinidis, N. P.; Leijnse, M.; Wegewijs, M. R.; May, F.; Hofstetter, W.; Danieli, C.; Cornia, A. Nano Lett. 2010, 10, 3307−3311. 92

DOI: 10.1021/acs.nanolett.7b03411 Nano Lett. 2018, 18, 88−93

Letter

Nano Letters (47) Jørgensen, C. K. Coord. Chem. Rev. 1966, 1, 164. (48) Hamann, C.; Woltmann, R.; Hong, I.; Hauptmann, N.; Karan, S.; Berndt, R. Rev. Sci. Instrum. 2011, 82, 033903. (49) Cohen, A. J.; Tozer, D. J.; Handy, N. C. J. Chem. Phys. 2007, 126, 214104.

93

DOI: 10.1021/acs.nanolett.7b03411 Nano Lett. 2018, 18, 88−93