Exchange Coupling of Spin-Crossover Molecules to Ferromagnetic Co

Feb 19, 2016 - *E-mail: [email protected]. .... Hervé Bulou , Fabrice Scheurer , François Bertran , Patrick Le Fèvre , Amina Taleb-Ibrahimi ...
0 downloads 0 Views 5MB Size
Subscriber access provided by GAZI UNIV

Letter

Exchange Coupling of Spin-crossover Molecules to Ferromagnetic Co Islands Saber Gueddida, Manuel Gruber, Toshio Miyamachi, Eric Beaurepaire, Wulf Wulfhekel, and Mebarek Alouani J. Phys. Chem. Lett., Just Accepted Manuscript • DOI: 10.1021/acs.jpclett.6b00172 • Publication Date (Web): 19 Feb 2016 Downloaded from http://pubs.acs.org on February 20, 2016

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry Letters is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 18

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Exchange Coupling of Spin-Crossover Molecules to Ferromagnetic Co Islands Saber Gueddida,† Manuel Gruber,∗,†,‡,§ Toshio Miyamachi,‡,∥ Eric Beaurepaire,† Wulf Wulfhekel,‡,¶ and Mebarek Alouani† Institut de Physique et de Chimie des Matériaux de Strasbourg (IPCMS), Université de Strasbourg, CNRS UMR 7504, 23 rue du Loess, BP 43, F-67034 Strasbourg Cedex 2, France, Physikalisches Institut, Karlsruhe Institute of Technology, Wolfgang-Gaede-Strasse 1, 76131 Karlsruhe, Germany, and Institute of Nanotechnology, Karlsruhe Institute of Technology, 76021 Karlsruhe, Germany



To whom correspondence should be addressed Institut de Physique et de Chimie des Matériaux de Strasbourg (IPCMS), Université de Strasbourg, CNRS UMR 7504, 23 rue du Loess, BP 43, F-67034 Strasbourg Cedex 2, France ‡ Physikalisches Institut, Karlsruhe Institute of Technology, Wolfgang-Gaede-Strasse 1, 76131 Karlsruhe, Germany ¶ Institute of Nanotechnology, Karlsruhe Institute of Technology, 76021 Karlsruhe, Germany § Present address: IEAP, Christian-Albrechts-Universität zu Kiel, 24098 Kiel, Germany ∥ Present address: The Institute for Solid State Physics (ISSP), The University of Tokyo, Kashiwa 2778581, Japan †

1 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Abstract The properties of Fe(1,10-phenanthroline)2 (NCS)2 (Fe-phen) molecules deposited on Co/Cu(111) are studied using scanning tunneling microscopy (STM) operated in ultrahigh vacuum at low temperature (4K) and ab initio calculations. Both the experimental and theoretical results are used to identify the high-spin (HS) state of Fe-phen. Additionally, the calculations reveal a strong spin-polarization of the density of states (DOS) and is validated experimentally using the spin sensitivity of spin-polarized STM. Finally, it is shown that the magnetic moment of the Fe-ion within HS Fe-phen is strongly magnetically coupled to the underlying magnetic Co through the NCS groups. These findings enable promising spintronic perspectives.

Graphical TOC Entry

Fe-phen

Co

Keywords

Cu

spin crossover, spin polarization, magnetism, density functional theory, scanning tunneling miscroscopy, scanning tunneling spectroscopy

2 ACS Paragon Plus Environment

Page 2 of 18

Page 3 of 18

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Organic spintronics is a fertile research field aiming at using organic materials to host spin transport while benefiting of new advantages such as material flexibility, light weight and potential low production cost. 1 Furthermore, various experiments have recently shown that organic molecules on a ferromagnetic substrate can strongly enhance the spin polarization of the molecule-metal interface or “spinterface”. 2–5 Alternatively, spin-crossover (SCO) complexes have gathered considerable attention for their potentialities as intrinsic switch to build nanoscale electronic components. 6–10 Indeed, SCO complexes can undergo a spin transition upon application of various stimuli such as temperature, light, pressure and magnetic field. 11 In this work, we experimentally and theoretically explore the synergy between the two research fields by studying the (spin) properties of single SCO molecules deposited on a ferromagnetic substrate. The molecular system of interest is Fe-phen (see Fig. 1(a)). This prototypical SCO complex 12 undergoes a thermal SCO from the low-spin (LS) state S=0 to the HS state S=2 when heated above 175 K. 13 The transition is accompanied by changes in the Fe-N bond lengths and angles 14 that lead to changes in molecular conformation. 6,15 We deposited Fe-phen on Co bilayer islands grown over a Cu(111) surface (see methods). The Co islands form triangular-like shapes with opposite orientations corresponding to fcc and hcp stacking structures with a ratio of 3:2. 16,17 In Fig. 1(b), we show a scanning tunneling microscopy (STM) topography recorded at 4K with a home built instrument 18 where the Co bilayer islands are clearly visible. We additionally observe protrusions over the Co islands corresponding to molecular adsorbates upon Fe-phen deposition. In Fig. 1(c) we show the three different types of adsorbates that were observed. Molecules of type I and II are of similar sizes and both present an elongated shape. While type I molecules have their maximum apparent height in the center, type II molecules exhibit two protrusions on their sides (see height profiles along the long axis of the molecules in Fig. 1(d)). We note that the overall shape of type I and II molecules present similarities to that of Fe-phen deposited on Cu(100), Cu(111) and CuN/Cu(100). 6,15 We therefore suppose that type I and II molecules are Fe-phen molecules, possibly in different spin states. In addition, we believe that, similarly

3 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

to Fe-phen on Cu(100), Cu(111) and CuN substrates, the NCS groups of Fe-phen bind to the Co/Cu(111) substrate while the phenanthroline (phen) groups face upwards and are seen with the STM as two lobes. The third kind of molecular adsorbate is much smaller than type I and II adsorbates and is compatible with a single phen group. 19 We note that we observed single phen groups in a non-negligible proportion. We foresee two possible scenarios accounting for the presence of the isolated phen groups: i) phen groups were co-evaporated with Fe-phen and/or ii) some Fe-phen molecules reacted with the Co/Cu(111) substrate and were decomposed into phen and the remaining part of Fe-phen. The precise origin of the single phen groups is out of the scope of this paper, and in the following, we will focus on type I and type II molecules. In Fig. 1(e), we present the dI/dV (V ) spectra recorded on the center of type I and type II molecules and on the Co surface. The Co reference spectrum resembles that of previous studies. 17,20,21 The spectrum over the center of type I molecule does not reveal clear spectroscopic features while that on the center of type II molecule reveals peaks at the following energies: -1.18 V, -0.57 V, -0.26 V, +0.25 V. Since we do not yet know the energy position of the molecular orbitals of HS and LS Fe-phen, we are unable to identify type I and type II molecules to HS and LS Fe-phen. To circumvent this limitation, we performed ab initio density functional theory (DFT) calculations of HS and LS Fe-phen on Co/Cu(111). In Fig. 2 we show the calculated density of states (DOS) within the GGA+U+vdW method for LS and HS Fe-phen (see methods and SI). Within the measured energy interval, the LS Fe-phen DOS has a single wide peak centered at 0.28 eV. These states mainly arise from the phen groups with small contributions from the Fe ion (see Fig. 2(a)). The calculated DOS for HS Fe-phen presented in Fig. 2(b) also exhibits a wide peak for the unoccupied states (centered at 0.28 eV). However, the peak is much wider than that of LS Fe-phen. We observe an additional peak, for HS Fe-phen, located at -0.90 eV. As expected, the difference between HS and LS Fe-phen mainly arises from the Fe DOS. We now compare the calculated DOS to the experimental dI/dV spectra presented in Fig. 1(e). The dI/dV for type II

4 ACS Paragon Plus Environment

Page 4 of 18

Page 5 of 18

(a)

(b)

Fe N C S H (e)

(c) I

phen II

dI/ dV ( arb. unit s)

6

Height (pm)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

4

I

II

2 Co

120 (d) 80 40 0.0

I II

0 -1 .0

1.0 Length (nm)

0 .0

1 .0

Sample volt age ( V)

Figure 1: (a) Three dimensional model of Fe-phen molecule. (b) STM topography of Cu(111)/Co/Fe-phen. The black background corresponds to the Cu(111) surface while orangish triangular-like structures are Co bilayer islands. The yellow spots on the Co islands are molecular adsorbates. (c) STM topography of type I, type II and phen molecules on a Co island. (d) Height profiles along type I and II molecules. (e) dI/dV spectra recorded on the center of type I and II molecules and on the Co island. The image sizes are (b) 92×60 nm2 (V =1.0 V, I=100 pA) and (c) 3.2×5.4 nm2 (V =1.0 V, I=100 pA).

5 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

molecule exhibits a very wide resonance for positive sample voltage (unoccupied states) as well as a resonance located at -1.18 V. The position and the width of these resonances are in qualitative agreement with the calculated DOS for HS Fe-phen. Furthermore, the dI/dV for type II molecule has two additional peaks located at -0.57 V and -0.26 V. These peaks possibly come from states localized at the molecule/metal interface induced by the Co Shockley surface states. 22 Our calculated model has infinite lateral dimensions while the Co islands are of finite size. Therefore, the calculated Shockley-induced states can differ from the experimental one. Considering the qualitative agreement between the dI/dV spectrum of type II molecule and calculated DOS for HS Fe-phen, type II molecule is identified as HS Fe-phen. Note that further evidences for this identification are given below. Concerning type I molecule one may tentatively ascribe it to LS Fe-phen by deduction. However, there is less agreement with the calculations. Below the Fermi energy, the calculated DOS for LS Fe-phen has no strong peak but a small feature around -0.7 eV in agreement with experimental dI/dV at negative voltage. The peak at +0.2 eV could correspond to the shoulder observed in the dI/dV . We also recall the presence of isolated phen (see Fig. 1(c)). Type I molecules could therefore also be a residue of Fe-phen dissociation leading to the isolated phen. Considering the doubts associated to type I molecule, in the following, we only discuss type II molecules, i.e. HS Fe-phen molecules. Confirming the qualitative agreement between the experimental dI/dV and calculated DOS of HS Fe-phen, we stress the strong spin-polarization of the calculated HS Fe-phen DOS. Indeed, the spin-up DOS (solid blue line in Fig. 2(b)) differs from the spin-down DOS (dashed red line in Fig. 2(b)). The spin polarization reaches 84% at the Fermi level, i.e. higher than that of the Co alone in agreement with previous spinterface studies. 4,5 Actually, a spin-polarized DOS is expected for FeII complexes in the HS state since the FeII ion carries a magnetic moment. 11,23 To be of practical spintronic usage, the Fe-ion magnetic moment must be coupled to that of the underlying Co at finite temperature. We estimate the magnetic coupling by performing calculations of HS Fe-phen on Co/Cu(111) with the

6 ACS Paragon Plus Environment

Page 6 of 18

(a) Density of states (states/eV/Cell)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

6

Fe-phen

LS

3 0 1

Fe

0 2

phen

0 1

NCS

0 -1.5

Spin up Spin down -1.0

-0.5 0.0 E-EF (eV)

0.5

1.0

(b) Density of states (states/eV/Cell)

Page 7 of 18

6

Fe-phen

HS

3 0 2

Fe

0 phen

2 0 1 0 -1.5

NCS

Spin up Spin down -1.0

-0.5 0.0 E-EF (eV)

0.5

1.0

Figure 2: Calculated spin-dependent DOS of (a) LS Fe-phen and (b) HS Fe-phen adsorbed on Co/Cu(111) using the GGA+U+vdW method. For each sub-figure, the top panel indicates the spin-dependent DOS integrated over the entire molecule. The three next panels report the integrated spin-dependent DOS over the Fe ion, the phen and NCS groups. The spin-up (spin-down) DOS are represented in solid blue (dashed red).

7 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Fe-ion magnetic moment either parallel or anti-parallel to that of Co. We find an energy difference between the two magnetic configurations of ∆E = E AF M − E F M = 22.6 meV, which implies that i) the magnetic moment of the Fe-ion is preferentially ferromagnetically coupled to that of the underlying Co and ii) the coupling is relatively strong and can be observed at the measurement temperature of ≈ 4 K. Since the magnetic referential of HS Fephen is linked to that of the underlying Co, we can use the spin sensitivity of spin-polarized STS to experimentally confirm the spin-polarization of HS Fe-phen. To do so, we used a Co coated W tip, which possesses a dominant out-of-plane magnetization that is parallel or anti-parallel to the magnetization of the Co islands. The exact magnetic orientation of the tip is unknown, but as in previous spin-polarized STS studies on Co/Cu(111), 21,24–26 Co islands presenting the strongest (lowest) dI/dV signal -0.38 eV are referred to as parallel (anti-parallel). We then acquired dI/dV spectra over HS Fe-phen adsorbed on Co islands. In Fig. 3(a), we show two HS Fe-phen dI/dV spectra for HS Fe-phen lying on parallel and antiparallel spin polarized Co islands, respectively represented in dashed red and solid blue. The amplitudes and the widths of the peaks of the HS Fe-phen dI/dV spectra clearly depend on the magnetic orientation of the underlying Co islands. As pointed out in other spin-polarized STS studies, 25,27,28 the difference in the dI/dV spectra arises from the asymmetry between the spin up and spin down DOS, weighted by the spin-polarization of the tip. Interestingly, the experimental dI/dV presented in Fig. 3(a) are very similar to the calculated DOS of HS Fe-phen presented in Fig. 2(b). For the unoccupied states, we see that the spin-up DOS peak is higher but narrower than that of the spin-down DOS. In addition, the spin-down DOS has a peak around -0.90 eV that is absent for the spin-up DOS. If one assumes that the tunneling current within the STM junction for a parallel (anti-parallel) tip/sample spin polarizations is dominated by spin-up electrons, we have a qualitative agreement between the ab initio calculations and the experiments. We now stress several points: i) as for the HS Fe-phen dI/dV spectrum acquired with a W tip (see Fig. 1(e)), the calculations fail to reproduce the peak at around -0.5 eV; ii) for the experimental HS Fe-phen dI/dV , the

8 ACS Paragon Plus Environment

Page 8 of 18

Page 9 of 18

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

“parallel” (“anti-parallel”) Co islands had a hcp (hcp) while the calculations were carried out using a fcc stacking. Concerning i), we believe that the peak(s) located around -0.5 eV originate from surface states as described above for the comparison of Fig. 1(e) and Fig. 2(b). ii) The different stacking structures is known to induce an energy shift of the Co occupied states by about 0.03 eV. 20,24 It is also known that the same states can shift by more than 0.1 eV with the size of the islands. 20 Therefore, for the acquired dI/dV spectra, we took care to select islands of similar sizes but ignored the stacking structure of the islands. This refrains us to appreciate shifts lower than 0.03 V. Nevertheless, the experimental data support the calculations concerning two points: i) HS Fe-phen has a spin-polarized DOS close to the Fermi level and ii) the magnetic moment of HS Fe-phen is coupled to that of the underlying Co. For completeness, we also theoretically estimated the conductance of the system composed of a HS Fe-phen molecule on Co/Cu(111) and an ideal Co tip (see Fig. 3(b)). The conductance is estimated by the Jullière model using the calculated DOS of HS Fe-phen and the independently calculated DOS of an ideal Co STM tip (see SI). For the unoccupied states, we observe that the conductance of an anti-parallel spin polarization has a narrower but higher peak than that of a parallel spin polarization. This remains in qualitative agreement with the experimental data and any deviation can be easily understood since i) we calculated an idealized STM tip while the exact shape of the experimental tip apex is unknown, ii) the DOS are calculated for zero voltage and iii) the tunneling probability is not the same for all electrons but depends on their in-plane k vector. Various studies already demonstrated the magnetic coupling between a molecular magnetic moment and a ferromagnetic substrate. 30–34 However, in the present case, the Fe-ion that carries the magnetic moment, is not in direct contact with the substrate but at a distance of 5.1 Å. This distance prohibits direct exchange interaction. In order to have more insights about the origin of the magnetic coupling, we show calculated magnetization density of HS Fe-phen on Co/Cu(111) in Fig. 4. The Fe-ion carries a magnetic moment of 3.5 µB in

9 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

(a)

Normalized dI/ dV ( arb. unit s)

Parallel Ant i-parallel

-1 .5

-1 .0

-0 .5 0 .0 0 .5 Sample volt age ( V)

1 .0

(b)

Est imat ed conduct ance ( arb. unit s)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Parallel Ant i-parallel

-1 .5

-1 .0

-0 .5 0 .0 E-EF ( eV)

0 .5

1 .0

Figure 3: (a) Experimental normalized dI/dV acquired over Fe-phen lying on “parallel” (dashed red) and “anti-parallel” (solid blue) Co islands. We performed the normalization using the procedure described in Ref. 29. (b) Conductance estimated by the Jullière model of HS Fe-phen and an ideal Co tip for parallel (dashed red) and anti-parallel (solid blue) spin polarizations.

10 ACS Paragon Plus Environment

Page 10 of 18

Page 11 of 18

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

comparison to 3.6 µB for the gaz-phase molecule. In addition, we observe that the nitrogen atoms of the NCS groups as well as the sulfur atoms have an induced magnetization. To a limited extend, the N atoms of the phen groups are also magnetized. The (induced) magnetization on the Fe, N and S within NCS is of the same sign than that of the underlying Co. The ferromagnetic coupling between the Fe-ion and the Co surface may be explained as follow: the S atoms are magnetized by and coupled to the underlying Co. S and N within NCS are magnetically coupled by double-exchange interaction through the C atoms. Finally the N of the NCS groups couples to Fe by direct-exchange interaction. The remarkable point for this system is the apparent large magnetic-coupling energy of 22.6 meV despite the large distance between the Fe and the Co. We also calculated the magnetization density of LS Fe-phen (see SI). Interestingly, the Fe ion has an induced magnetic moment.

Figure 4: Calculated magnetization density of HS Fe-phen using the GGA+U+vdW method. The positive magnetization density is red, and the negative is blue. The values of magnetic moments of iron, nitrogen, carbon and sulfur are also indicated. In conclusion, we have experimentally and theoretically shown that the magnetic moment of a spin-crossover complex (here Fe-phen) can be strongly magnetically coupled to the underlying ferromagnetic substrate (here Co). We have additionally shown that the spinpolarization of the DOS is relatively important close to the Fermi level revealing promising 11 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

spintronic perspectives.

Methods Sample preparation Clean and flat Cu(111) surfaces were prepared by several cycles of Ar sputtering (1.5 keV) and subsequent annealing to 500◦ C. The Co was deposited from a e-beam source to the Cu(111) sample held at room temperature. Fe-phen molecules were subsequently sublimed at 180◦ C while maintaining the substrate at a temperature of 50◦ C. Tips STM tips were prepared by chemical etching of a tungsten wire and subsequent flashing in ultrahigh vacuum. The tip used for non spin-polarized measurements were repeatedly indented into the Cu(111) substrate in order to sharpen the tip. For spin-polarized contrast, the flashed W tip was covered by about 12-13 ML of Co. Differential conductance acquisition Differential conductance dI/dV spectra were obtained with a lock-in technique using a modulation voltage ranging from 15 mV to 30 mV and a frequency of 5.03 kHz. When ramping the voltage, the feedback loop was opened. STM topographies All the STM topographies were acquired using the constant-current mode. Ab initio calculations The first-principle calculations were carried out using the Vienna Ab initio Simulation Package (VASP) 35 using the plane-wave pseudo-potential method. 36 For the exchange-correlation functional, we adopted the generalized gradient approximation (GGA) with the Perdew, 12 ACS Paragon Plus Environment

Page 12 of 18

Page 13 of 18

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Burke and Ernzerhof parametrization. 37 In order to account for the strong correlated electronelectron interactions within the central iron atom of the investigated molecule, we used the GGA+U method 38 as implemented in the VASP code, 39 in which an additional on-site Hubbard-U term is included on the iron atom. We also included van der Waals (vdW) interactions in the calculations 40,41 (GGA+U+vdW), which is known to give more reliable structural, electronic and magnetic properties of metallo-organic molecules on a metallic substrate. 42 In this work, we used an effective U of 3 eV and an exchange parameter J=0.9 eV for the 3d electrons of the Fe atom within the molecule. These values have been used in a previous study 23 and correctly describe the properties of the free and adsorbed Fe-phen molecule on a metallic substrate. The kinetic energy cutoff for the expansion of the wave functions is 450 eV at the Gamma point of the Brillouin zone.

Acknowledgement We thank G. Rogez for synthesizing the molecule. We acknowledge funding from the FrenchGerman University and the Baden-Württemberg Stiftung in the framework of the Kompetenznetz für Funktionale Nanostrukturen (KFN), from the Institut Carnot MICA’s “Spinterface” grant, from the Agence Nationale de la Recherche ANR-09-JCJC-0137 and ANR-11-LABX0058 NIE and from the International Center for Frontier Research in Chemistry. This work was performed using HPC resources from GENCI-CINES Grant 2014-gem1100 and Strasbourg HPC Mesocenter.

Supporting Information Available Additional figures and information. This material is available free of charge via the Internet at http://pubs.acs.org/.

Note: The authors declare no competing financial interest. 13 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

References (1) Dediu, V. A.; Hueso, L. E.; Bergenti, I.; Taliani, C. Spin Routes in Organic Semiconductors. Nat. Mater. 2009, 8, 707–716. (2) Barraud, C.; Seneor, P.; Mattana, R.; Fusil, S.; Bouzehouane, K.; Deranlot, C.; Graziosi, P.; Hueso, L.; Bergenti, I.; Dediu, V. et al. Unravelling the Role of the Interface for Spin Injection into Organic Semiconductors. Nature Phys. 2010, 6, 615–620. (3) Sanvito, S. Molecular spintronics: The Rise of Spinterface Science. Nature Phys. 2010, 6, 562–564. (4) Djeghloul, F.; Ibrahim, F.; Cantoni, M.; Bowen, M.; Joly, L.; Boukari, S.; Ohresser, P.; Bertran, F.; Le Fèvre, P.; Thakur, P. et al. Direct Observation of a Highly SpinPolarized Organic Spinterface at Room Temperature. Sci. Rep. 2013, 3, 1272. (5) Djeghloul, F.; Garreau, G.; Gruber, M.; Joly, L.; Boukari, S.; Arabski, J.; Bulou, H.; Scheurer, F.; Hallal, A.; Bertran, F. et al. Highly Spin-Polarized Carbon-Based Spinterfaces. Carbon 2015, 87, 269–274. (6) Miyamachi, T.; Gruber, M.; Davesne, V.; Bowen, M.; Boukari, S.; Joly, L.; Scheurer, F.; Rogez, G.; Yamada, T. K.; Ohresser, P. et al. Robust Spin Crossover and Memristance Across a Single Molecule. Nat. Commun. 2012, 3, 938. (7) Gopakumar, T. G.; Matino, F.; Naggert, H.; Bannwarth, A.; Tuczek, F.; Berndt, R. Electron-Induced Spin Crossover of Single Molecules in a Bilayer on Gold. Angew. Chem., Int. Ed. 2012, 51, 6262–6266. (8) Meded, V.; Bagrets, A.; Fink, K.; Chandrasekar, R.; Ruben, M.; Evers, F.; BernandMantel, A.; Seldenthuis, J. S.; Beukman, A.; Van der Zant, H. S. J. Electrical Control over the Fe(II) Spin Crossover in a Single Molecule: Theory and Experiment. Phys. Rev. B 2011, 83, 245415. 14 ACS Paragon Plus Environment

Page 14 of 18

Page 15 of 18

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

(9) Bousseksou, A.; Molnár, G.; Salmon, L.; Nicolazzi, W. Molecular Spin Crossover Phenomenon: Recent Achievements and Prospects. Chem. Soc. Rev. 2011, 40, 3313–3335. (10) Molnár, G.; Salmon, L.; Nicolazzi, W.; Terki, F.; Bousseksou, A. Emerging Properties and Applications of Spin Crossover Nanomaterials. J. Mater. Chem. C 2014, 2, 1360. (11) Gütlich, P.; Garcia, Y.; Goodwin, H. A. Spin Crossover Phenomena in Fe(II) Complexes. Chem. Soc. Rev. 2000, 29, 419–427. (12) König, E.; Madeja, K. Unusual Magnetic Behaviour of Some Iron(II)–bis-(1,10phenanthroline) Complexes. Chem. Commun. 1966, 61–62. (13) Davesne, V.; Gruber, M.; Miyamachi, T.; Da Costa, V.; Boukari, S.; Scheurer, F.; Joly, L.; Ohresser, P.; Otero, E.; Choueikani, F. et al. First Glimpse of the Soft X-Ray Induced Excited Spin-State Trapping Effect Dynamics on Spin Cross-Over Molecules. J. Chem. Phys 2013, 139, 074708. (14) Gallois, B.; Real, J. A.; Hauw, C.; Zarembowitch, J. Structural Changes Associated with the Spin Transition in bis(isothiocyanato)bis(1,10-phenanthroline)iron: A SingleCrystal X-Ray Investigation. Inorg. Chem. 1990, 29, 1152–1158. (15) Gruber, M.; Davesne, V.; Bowen, M.; Boukari, S.; Beaurepaire, E.; Wulfhekel, W.; Miyamachi, T. Spin State of Spin-Crossover Complexes: From Single Molecules to Ultrathin Films. Phys. Rev. B 2014, 89, 195415. (16) Vázquez de Parga, A. L.; García-Vidal, F. J.; Miranda, R. Detecting Electronic States at Stacking Faults in Magnetic Thin Films by Tunneling Spectroscopy. Phys. Rev. Lett. 2000, 85, 4365–4368. (17) Heinrich, B. W.; Iacovita, C.; Rastei, M. V.; Limot, L.; Ignatiev, P. A.; Stepanyuk, V. S.; Bucher, J. P. A Spin-Selective Approach for Surface States at Co Nanoislands. EPJ B 2010, 75, 49–56. 15 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(18) Zhang, L.; Miyamachi, T.; Tomanić, T.; Dehm, R.; Wulfhekel, W. A Compact subKelvin Ultrahigh Vacuum Scanning Tunneling Microscope with High Energy Resolution and High Stability. Rev. Sci. Instrum. 2011, 82, 103702. (19) Gopakumar, T. G.; Bernien, M.; Naggert, H.; Matino, F.; Hermanns, C. F.; Bannwarth, A.; Mühlenberend, S.; Krüger, A.; Krüger, D.; Nickel, F. et al. Spin-Crossover Complex on Au(111): Structural and Electronic Differences Between Mono- and Multilayers. Chem. Eur. J. 2013, 19, 15702–15709. (20) Rastei, M. V.; Heinrich, B.; Limot, L.; Ignatiev, P. A.; Stepanyuk, V. S.; Bruno, P.; Bucher, J. P. Size-Dependent Surface States of Strained Cobalt Nanoislands on Cu(111). Phys. Rev. Lett. 2007, 99, 246102. (21) Schmaus, S.; Bagrets, A.; Nahas, Y.; Yamada, T. K.; Bork, A.; Bowen, M.; Beaurepaire, E.; Evers, F.; Wulfhekel, W. Giant Magnetoresistance Through a Single Molecule. Nature Nanotech. 2011, 6, 185–189. (22) Heinrich, B. W.; Limot, L.; Rastei, M. V.; Iacovita, C.; Bucher, J. P.; Djimbi, D. M.; Massobrio, C.; Boero, M. Dispersion and Localization of Electronic States at a Ferrocene/Cu(111) Interface. Phys. Rev. Lett. 2011, 107, 216801. (23) Gueddida, S.; Alouani, M. Spin Crossover in a Single Fe(phen)2 (NCS)2 Molecule Adsorbed onto Metallic Substrates: An Ab Initio Calculation. Phys. Rev. B 2013, 87, 144413. (24) Pietzsch, O.; Kubetzka, A.; Bode, M.; Wiesendanger, R. Spin-Polarized Scanning Tunneling Spectroscopy of Nanoscale Cobalt Islands on Cu(111). Phys. Rev. Lett. 2004, 92, 057202. (25) Iacovita, C.; Rastei, M. V.; Heinrich, B. W.; Brumme, T.; Kortus, J.; Limot, L.; Bucher, J. P. Visualizing the Spin of Individual Cobalt-Phthalocyanine Molecules. Phys. Rev. Lett. 2008, 101, 116602. 16 ACS Paragon Plus Environment

Page 16 of 18

Page 17 of 18

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

(26) Fei, X.; Wu, G.; Lopez, V.; Lu, G.; Gao, H.-J.; Gao, L. Spin-Dependent Conductance in Co/C60/Co/Ni Single-Molecule Junctions in the Contact Regime. J. Phys. Chem. C 2015, 119, 11975–11981. (27) Wiesendanger, R. Spin Mapping at the Nanoscale and Atomic Scale. Rev. Mod. Phys. 2009, 81, 1495–1550. (28) Schwöbel, J.; Fu, Y.; Brede, J.; Dilullo, A.; Hoffmann, G.; Klyatskaya, S.; Ruben, M.; Wiesendanger, R. Real-Space Observation of Spin-Split Molecular Orbitals of Adsorbed Single-Molecule Magnets. Nat. Commun. 2012, 3, 953. (29) Yamagishi, Y.; Nakashima, S.; Oiso, K.; Yamada, T. K. Recovery of Nanomolecular Electronic States From Tunneling Spectroscopy: LDOS of Low-Dimensional Phthalocyanine Molecular Structures on Cu(111). Nanotechnology 2013, 24, 395704. (30) Scheybal, A.; Ramsvik, T.; Bertschinger, R.; Putero, M.; Nolting, F.; Jung, T. A. Induced Magnetic Ordering in a Molecular Monolayer. Chem. Phys. Lett. 2005, 411, 214–220. (31) Wende, H.; Bernien, M.; Luo, J.; Sorg, C.; Ponpandian, N.; Kurde, J.; Miguel, J.; Piantek, M.; Xu, X.; Eckhold, P. et al. Substrate-Induced Magnetic Ordering and Switching of Iron Porphyrin Molecules. Nat. Mater. 2007, 6, 516–520. (32) Bernien, M.; Miguel, J.; Weis, C.; Ali, M. E.; Kurde, J.; Krumme, B.; Panchmatia, P. M.; Sanyal, B.; Piantek, M.; Srivastava, P. et al. Tailoring the Nature of Magnetic Coupling of Fe-Porphyrin Molecules to Ferromagnetic Substrates. Phys. Rev. Lett. 2009, 102, 047202. (33) Gruber, M.; Ibrahim, F.; Boukari, S.; Isshiki, H.; Joly, L.; Peter, M.; Studniarek, M.; Da Costa, V.; Jabbar, H.; Davesne, V. et al. Exchange Bias and Room-Temperature Magnetic Order in Molecular Layers. Nat. Mater. 2015, 14, 981–984.

17 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(34) Gruber, M.; Ibrahim, F.; Boukari, S.; Joly, L.; Da Costa, V.; Studniarek, M.; Peter, M.; Isshiki, H.; Jabbar, H.; Davesne, V. et al. Spin-Dependent Hybridization between Molecule and Metal at Room Temperature through Interlayer Exchange Coupling. Nano Lett. 2015, 15, 7921–7926. (35) Kresse, G.; Furthmüller, J. Efficient Iterative Schemes for Ab Initio Total-Energy Calculations Using a Plane-Wave Basis Set. Phys. Rev. B 1996, 54, 11169–11186. (36) Blöchl, P. E. Projector Augmented-Wave Method. Phys. Rev. B 1994, 50, 17953–17979. (37) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. (38) Anisimov, V. I.; Aryasetiawan, F.; Lichtenstein, A. I. First-Principles Calculations of the Electronic Structure and Spectra of Strongly Correlated Systems: the LDA + U Method. J. Phys.: Condens. Matter 1997, 9, 767. (39) Bengone, O.; Alouani, M.; Blöchl, P.; Hugel, J. Implementation of the Projector Augmented-Wave LDA+U Method: Application to the Electronic Structure of NiO. Phys. Rev. B 2000, 62, 16392–16401. (40) Grimme, S. Semiempirical GGA-type Density Functional Constructed with a LongRange Dispersion Correction. J. Comput. Chem. 2006, 27, 1787–1799. (41) Bučko, T.; Hafner, J.; Lebègue, S.; Ángyán, J. G. Improved Description of the Structure of Molecular and Layered Crystals: Ab Initio DFT Calculations with van der Waals Corrections. J. Phys. Chem. A 2010, 114, 11814–11824. (42) Javaid, S.; Lebègue, S.; Detlefs, B.; Ibrahim, F.; Djeghloul, F.; Bowen, M.; Boukari, S.; Miyamachi, T.; Arabski, J.; Spor, D. et al. Chemisorption of Manganese Phthalocyanine on Cu(001) Surface Promoted by van der Waals Interactions. Phys. Rev. B 2013, 87, 155418. 18 ACS Paragon Plus Environment

Page 18 of 18