Spin-communication channels between Ln(III) bis-phthalocyanines molecular nanomagnets and a magnetic substrate.

Learning the art of exploiting the interplay between different units at the atomic scale is a fundamental step in the realization of functional nano-architectures and interfaces. In this context, understanding and controlling the magnetic coupling between molecular centers and their environment is still a challenging task. Here we present a combined experimental-theoretical work on the prototypical case of the bis(phthalocyaninato)-lanthanide(III) (LnPc2) molecular nanomagnets magnetically coupled to a Ni substrate. By means of X-ray magnetic circular dichroism we show how the coupling strength can be tuned by changing the Ln ion. The microscopic parameters of the system are determined by ab-initio calculations and then used in a spin Hamiltonian approach to interpret the experimental data. By this combined approach we identify the features of the spin communication channel: the spin path is first realized by the mediation of the external (5d) electrons of the Ln ion, keeping the characteristic features of the inner 4 f orbitals unaffected, then through the organic ligand, acting as a bridge to the external world.

The continuous trend towards ever smaller devices and components is leading to the requirement of fine control at the atomic scale over materials and interfaces. The natural solution can be found in molecular materials, which can be produced in large amounts and whose properties can be tailored by a bottom-up approach. Single molecule magnets12, featuring the spin degree of freedom3, have been demonstrated to possess specific magnetic properties and exhibit quantum phenomena such that they have been proposed to be employed in classical as well as quantum spintronics applications45. The enabling challenge here is the possibility to design specific magnetic features directly at the synthetic level and exploit them through a smart interaction with the environment, for instance a magnetic or conducting substrate. This requires a fine tailoring of the molecular characteristics, with the contemporary preservation of the core molecular properties and the existence of a sizable magnetic interaction with the external world through the organic ligands.

While the magnetization of simple paramagnetic molecules is pinned by a magnetic substrate and can only mimic its behavior6789, an intriguing case is represented by the TbPc2 molecule. Here the characteristic properties of the isolated system, such as the huge anisotropy barrier resulting from the “double-decker” structure in which the Tb ion is sandwiched by two phthalocyanines10, are essentially preserved while a sizable magnetic interaction with the substrate has been reported111213. This optimally balanced coupling is obtained by combining a lanthanide center, with well localized f-orbitals, with suitable organic ligands, hosting delocalized π-type electrons.

The microscopic understanding of the key interactions, linking the specific molecular features with its experimentally observed behavior, is a challenging task. This, however, may prompt the development of hybrid molecular devices, such as supramolecular spin valves made of TbPc2 on carbon nanotubes14 or graphene15 or molecular spin transistors1617. Along this line, previous results include the study of the superexchange interaction in chains of CoPc molecules by spin-flip spectroscopy18 and the modelling of the interaction between a magnetic molecule and a non-magnetic metallic surface with an anisotropic Kondo impurity by means of numerical renormalization group calculations19. When magnetic substrates are considered, their interplay with the molecular layers is also at the heart of the spin-interface problem3 and it is mandatory for using molecular units as probes on surfaces with exotic properties20.

We report here a X-ray absorption spectroscopy (XAS) and magnetic circular dichroism (XMCD) investigation on LnPc2 (Ln = Tb, Dy, Er) molecules deposited on Ni (111) single crystal. Since XMCD is element-sensitive, we can follow the magnetization of the Ln ions and Ni simultaneously. We observe that the strength of the magnetic molecule-substrate interaction changes with the Ln ion, proving the dependence on the degree of hybridization between the localized magnetic 4f-orbitals of the Ln and the delocalized π-type electrons of the Pc ligands. Experimental results are rationalized with the support of a full theoretical modelling, where the microscopic details of the molecular units are firstly estimated via ab initio methods, and then inserted in spin Hamiltonian used to fit the experimental curves and derive the values for the exchange coupling. By doing this, we identify the microscopic details of the spin path linking the Ln ion with the Pc ligand and then to the substrate.

We find that there is not a direct overlap between the localized 4f states and the Pc unit, but rather the interaction occurs through the auxiliary 5d orbitals of the Ln ion, which become spin-polarized by the 4f–electrons and hybridize with the Pc ligands π-orbitals. This three-step spin-communication channel [(4f)Ln ↔ (5d)Ln ↔ (Pc) ↔ surface] is responsible for the experimentally observed coupling between the Ln ions and their environment, yet preserving unaffected the pristine magnetic properties of the molecules.

Results

System characterization and high fields XMCD

Figure 1(a) displays a schematic view of our experimental geometry. More details on the sample preparation and preliminary characterizations are given in the Method Section and Supplementary Notes 1–2 and Supplementary Figures 1–4. Briefly, the samples were prepared in situ by evaporating ~0.3 monolayers of LnPc2 molecules on the freshly prepared Ni(111) surface. STM analysis combined with XLD measurements (see Supplementary Figure 1 and Figure 3) confirm that the molecules adsorb in a flat geometry, with the Pc plane parallel to the surface111221222324252627. XMCD measurements at the L2,3 absorption edges of Ni and the M4,5 absorption edges of Ln were performed in total electron yield mode at the ID08 beamline of the ESRF, Grenoble. The experimental configuration is such that the magnetic field is always parallel to the incident photon beam, with an angle Θ with respect to the normal of the sample surface which can be rotated (see Fig. 1(a)). The base temperature during all of our measurements was 8 K.

Figure 1

XAS and XMCD characterization of the system.

(a) Schematic view of the experimental geometry; (b–i), XAS and XMCD spectra for Ni-L2,3 (b,c), Tb-M4,5 (d,e), Dy-M4,5 (f,g), Er-M4,5 (h,i) edges in LnPc2/Ni(111) at Θ  = 0° (left panels) and Θ = 70° (right panels) incidence angles, taken at an applied external field  = 5 T and temperature T = 8 K. XAS spectra have been normalized with respect to the average L3 (M5) intensity (I+ + I−)/2, which has been set to 1.

Figure 1(b–i) shows the XAS and XMCD spectra of Ni-L2,3 edges and Ln-M4,5 edges for the LnPc2 molecules taken at normal and grazing incidence at a magnetic field  = 5 T. These spectra compare well with the ones already reported111221222324252627, demonstrating the robust character of these molecules. The XMCD intensity taken at the Ln M4,5 edges strongly depends on the field direction for TbPc2 (Fig. 1(d–e)) and DyPc2 (Fig. 1(f–g)) as previously reported111221222324252627, and also for ErPc2 (Fig. 1(h–i)): Tb and Dy show the maximum for Θ  = 0° while for Er the signal is vanishingly small and becomes larger for increasing Θ. These observations confirm that Tb and Dy have an easy axis of magnetization perpendicular to the molecular plane (Θ  = 0°) whilst Er exhibits an easy plane as observed in bulk Ln(III) bis-phathalocyanines10.

Magnetization measurements

In order to study the magnetic coupling between the molecules and the substrate, we performed isothermal magnetization cycles, which can be obtained by plotting the XMCD signals (see Supplementary Note 3 for more details) for each element as a function of the external magnetic field at a fixed temperature (8 K in our measurements). In Fig. 2(e–i) we show the element-resolved magnetization-loop for each lanthanide, with the field applied perpendicular (left panels) and at grazing angle (right panels) to the sample surface. For the ErPc2 molecules, only the curve at Θ = 70° is shown, given the vanishingly small signal observed at normal incidence. Along with the Ln magnetization of the LnPc2 molecules, the corresponding magnetization for the Ni substrate is plotted in Fig. 2(c–d).

Figure 2

Magnetization cycles.

(a,b) schematic representation of the magnetic coupling between the TbPc2 molecules and the Ni substrate in the low field region where the antiferromagnetic interaction prevails. At Θ  = 70° the magnetic moment possesses a (small) component parallel to the surface (usually negligible considering a perfectly axial Vcf), as a result of the introduction of the low symmetry crystal field. (c-i) element-resolved magnetization measurements of Ni (c,d), Tb (e,f), Dy (g,h), Er (i) for the LnPc2/Ni(111) systems, taken at Θ  = 0° (left panels) and Θ  = 70° (right panels) incidence angles. Experimental data are shown as black dots, while the theoretical fit (as described in the text) as the continuous red line. (j) Exchange coupling constant K (left axis), as determined from the fit, and magnitude of the magnetic polarization (right axis) induced on the 5d orbitals for each Ln compound calculated from DFT as described in the text. The error bars in K take into account also the uncertainty in the experimental data.

For high magnetic fields, both the molecules and the Ni moments are aligned parallel to the external field and follow a conventional magnetization curve. However, for small fields, the magnetization curves measured on the lanthanides are generally opposed to ; in addition, they present an abrupt change of the slope in correspondence of the field Bs at which the Ni magnetization saturates (Bs = 0.5 T for Θ = 0° and 0.2 T for Θ = 70°). These features reveal a coupling of the lanthanide’s magnetic moment to the Ni substrate of antiferromagnetic character. Close to  = 0, the coupling with the underlying Ni crystal forces the lanthanide moment to be antiparallel to the substrate (see Fig. 2(a,b)), while at high fields the Zeeman interaction prevails and the magnetizations are all parallel to the external field. When these two interactions compensate each other, the magnetization of the molecule is zero, and the actual magnetic field value at which this condition is satisfied (which we call Bexc) can be considered as a qualitative estimation of the coupling strength. We note that the magnetic coupling can be tuned by changing the Ln ion as the interaction strength displays a clear trend in going from TbPc2 to ErPc2. This can be better visualized considering for instance the values of Bexc at Θ = 70° for the three compounds, namely Bexc = 0.90 T, 0.35 T and 0.25 T for Tb, Dy and Er, respectively. We shall address this observation in the following sections.

Theoretical modeling of the magnetic interaction

In order to proceed with a quantitative analysis of the Ln magnetization curves reported in Fig. 2 to extract information about the strength and character of the molecule-surface coupling, we modeled the system with the following spin Hamiltonian11:

where L and S are the orbital and spin angular momentum operators of the 4f electrons, the second term describes the Zeeman interaction between the Ln total magnetic moment and the magnetic field, and the last term describes the exchange-like magnetic coupling between the lanthanide 4f-electronic spin S and the Ni magnetization M via the coupling constant K. The Ni magnetization M is normalized by taking the ratio with its saturation value, and it is taken from the experimental curves presented in Fig. 2(a–b)). Strong spin-orbit coupling and strong on-site electron repulsion within the 4f-shell are taken into account by projecting the Hamiltonian onto the Ln ground atomic spin-orbit multiplet , where JMAX = L + S (i.e. |7F6 > for TbIII, |6H15/2 > for DyIII, and |4I15/2 > for ErIII). Within such subspace, it is always true that S/S = L/L = J/J, i.e. the three operators are all equivalent, up to a multiplicative constant. Finally, VCF is the microscopic crystal field potential exerted by the charge distribution of the ligands on the 4f electrons. Given the pseudo D4 symmetry of LnPc2 molecules, it is customary to approximate VCF to an axially symmetric electrostatic potential containing only three fitting parameters, which were first optimized by Ishikawa et al.10. However, such an approximation (as the X-ray structure has no formal local point-group symmetry, although not too far from an idealized D4 structure) can be inaccurate especially when the angular dependence of the magnetization is considered28. Therefore, we perform a full ab initio calculation of all 27 crystal field parameters for the 4f electrons of the LnPc2 molecules, using the experimental (X-ray) structure measured for the TbPc2 compound in bulk29. The ab initio methodology to describe the crystal field spectrum in rare earth complexes is well established303132 and it consists of the Complete Active Space Self-Consistent Field (CASSCF) and State-Interacting (RASSI) methodologies32, accounting for strong electron correlation and strong spin-orbit coupling, respectively (more details about the calculations are in the Methods section and in the Supplementary Note 4).

Inserting the optimized ab initio low-symmetry crystal field potential VCF in the Hamiltonian of Equation (1), we obtain the fit reported as solid lines in Fig. 2, panels (e–i). The quality of the fit is remarkably good, despite the use of just a single fit parameter (the molecule-surface coupling constant K) plus a multiplicative factor which scales the computed value for the Ln magnetization with the measured XMCD intensity. The best fits are obtained for KTb ~ 1.6 K, KDy ~ 1.1 K and KEr ~ 0.9 K, finding a quantitative estimation for the strength of the molecule-surface interaction and confirming the experimental observation that it depends on the specific Ln ion and decreases by going from TbPc2 to ErPc2. The visibility of this effect at the moderate temperature of 8 K should not be surprising, as the effect of the temperature is only to broad the thermal population of the eigenstates of Equation (1). At small external fields already sufficient to saturate the Ni substrate magnetization, the two lowest lying levels (which would be almost degenerate in the absence of coupling) will be separated by an energy roughly proportional to SLnKLn. As a consequence, the magnetization will deviate from a conventional curve up to temperatures even much higher than K.

Our fit procedure of the TbPc2 magnetization curve on Ni employs realistic crystal field parameters obtained from ab initio calculations, and without the introduction of a strongly anisotropic coupling interaction, which was recently proposed in Reference 11 to account for the experimental observation that Bexc at Θ = 70° for TbPc2 is (slightly) greater than Bexc at Θ = 0°. While the effect of low symmetry harmonics in the crystal field potential is indeed to mix the M = ±6 ground doublet with all the excited states (see Supplementary Note 4 for a description of the ground state wave functions) and also introduce J-mixing, only MJ-mixing was considered in the current approach, and this did not introduce significant improvements in the fit, although we expect that overall state mixing can be even more pronounced for the actual situation with molecules on substrate that are probably more distorted than what we estimated here.

It is possible that also other mechanisms contribute to the observed angular dependence of Bexc. For instance, the role of the spin s = 1/2 π-vacancy delocalized over the Pc ligands of the molecules could be taken into account since the spin radical may be involved in linking the magnetic interaction between the Ln center and the magnetic substrate3334.

DFT calculations of the intramolecular Ln-Pc interaction

In our experiments we tune the LnPc2 – Ni interaction by changing the Ln ion inside the molecule. Since the organic ligand is the same for the three molecules, the character of the Ln orbital wave-function should play a key role in determining the magnetic anisotropy and interactions. We are therefore addressing the intramolecular Ln-Pc interaction only, corroborating the intuitive picture that the Ln-Ni coupling is indeed mediated by the organic part of the molecule.

To get a deeper insight on the intramolecular coupling, we have performed DFT calculations on TbPc2, DyPc2, ErPc2 molecules in gas phase with the Perdew-Burke-Ernzerhof (PBE)35 exchange correlation functional as implemented in Quantum Espresso36 including static correlation effects using the LDA + U formalism37. Further computational details are given in the Methods section.

The calculated projected local density of states (LDOS) of the LnPc2 complexes are shown in Fig. 3, where we plot the states with f and d orbital character of the Ln ions, as well as the p-states of the N and C ions closer to the Ln ion, where is left most of the spin-polarization of the s = 1/2 spin radical delocalized over the two phthalocyanine parts.

Figure 3

DFT calculation of spin-resolved local density of states (LDOS) of the LnPc2 molecules.

(a) spin-resolved LDOS projected on the 4f orbitals of Ln = Tb, Dy and Er. (b) spin-resolved LDOS projected on the 5d orbital of Tb (the one for Dy and Er is very similar) and on the 2p orbitals of the nearest-neighbours (to the Ln ions) N and C atoms (which is only weakly modified by the choice of the Ln) in the phthalocyanine planes. In each panel, the upper (lower) side is relative to the majority (minority) spin channel. d-character in (b) has been enhanced for better comparison.

Already for Tb, the f orbitals are well below the Fermi level, due to the inclusion of U, and, moving from to Dy to Er, progressive filling of the minority f states pushes the electronic levels further down in energy. We observe a partial occupation of the 5d orbitals together with a small spin polarization in all the three cases with the same sign of the f spin moments. This implies that there is a hybridization between the f and d states in the Ln ions, which in our LDOS becomes evident for Tb, in the region of −4 eV from EF in the minority spin channel. The values of this induced polarization are plotted in Fig. 2(j), being 0.07 μB, 0.06 μB and 0.03 μB, respectively for the Tb, Dy and Er ions. We note that the trend going from Tb to Er is the same as observed experimentally for the coupling strengths. Since the Ln 5d orbitals occupy the same energy region of (and hybridize with) the p-states of the C and N atoms in the phthalocyanines, we conclude that such d-mediated interaction is the mechanism through which the spin information in the f orbitals propagates to the delocalized states of the two phthalocyanines, which, in turn, communicate to the external world by direct contact (in this case with the spin-polarized surface states). Our picture is in agreement with the one set forward in dinuclear complexes, where the role of 5d orbitals has been advocated in order to explain the exchange coupling through Oxygen between Gd and Cu ions38. We stress that this result has implications also for different studies focusing on LnPc2 and analogue compounds, where the magnetic coupling between the Ln ions to the external world is relevant1314151617 and, more generally, could help in the realization of novel Ln-based molecular magnets39.

Conclusions

We systematically studied how the interaction between LnPc2 single ion magnets and a ferromagnetic substrate is tuned by changing the Ln ion. With the support of composite ab initio, DFT, and spin Hamiltonian theoretical and computational approaches, we have identified the microscopic origin of the spin communication channels between the central Ln ion and the magnetic substrate, and characterized it as a three-step process: the f-electrons hybridize with the 5d orbitals of the metal ion, which in turn act as a bridge to the delocalized states in the organic part of the molecules, ultimately mediating the communication with the magnetic surface. We employed an ab initio calculated low-symmetry crystal field potential in our spin Hamiltonian to reproduce the experimental results evaluating the strength K of the magnetic coupling between lanthanide ion and surface, which is found to decrease in going from Tb to Dy to Er. This is explained in terms of the modulation of the hybridization between (4f) ad (5d) orbitals on these metal ions, as reproduced by the DFT calculations. These systems constitute a test bed for more general models which aim at a microscopic and quantitative description of the magnetic interactions between individual molecular compounds and magnetic substrates.

Methods

Sample preparation and synchrotron measurements

Experiments were carried out at the ID08 beamline of the European Synchrotron Radiation Facility (ESRF) in Grenoble, France. The Ni(111) single crystal (Surface Preparation Laboratory) was used as the substrate. Before molecule deposition, the surface was cleaned by repeated cycles of Ar+ sputtering (Energy = 2 KeV for 20 minutes and E = 0.8 KeV for 10 minutes) and annealing (Temperature = 800 °C for 5 minutes). The quality of the surface was checked by Low Energy Electron Diffraction (LEED). A ~0.3 monolayer of LnPc2 molecules was evaporated after long degassing of the powders, keeping the evaporator temperature at 420 °C at a base pressure of 1.0 × 10−9 mbar and monitoring the thickness with an in situ quartz microbalance. XMCD measurements at the L2,3 absorption edges of Ni and the M4,5 absorption edges of Ln were performed in total electron yield mode. The magnetic field was applied parallel to the incident photon beam, at an angle Θ with respect to the normal of the sample surface (see Fig. 1(a) for a schematic picture).

The STM investigation as well as further characterization of our sample (linear dichroism absorption measurements, ex situ XPS and Raman spectroscopies) are presented in Supplementary Notes 1–2 and Supplementary Figures 1–4.

DFT calculations

We performed our calculations using the projector augmented wave method (PAW)40 as recently implemented in the pseudopotentials made by Andrea Dal Corso41. We have taken into account the influence of the on-site Coulomb interaction imposing static correlation effects on the f electrons of the Ln ions36 (values of U and J have been taken by ref. 42, we used the same value of U = 9 eV and J = 1 eV for the Tb, Dy and Er ions).

No spin-orbit coupling is included in the calculation. We modeled our systems with a tetragonal supercell of 20 × 20 X 14 Å. Regarding the kinetic and charge density cutoffs, we have chosen 121 and 620 Ry, 125 and 499 Ry, 135 and 578 Ry, respectively for TbPc2, DyPc2 and ErPc2. A smearing parameter of 8 × 10−4 Ry ensured a correct description of the spin radical in the two Pcs. The version of Quantum Espresso used in all the calculations was the 5.0.3.

Ab-initio strategy for the evaluation of the crystal field parameters of the LnPc2 molecules

We perform an explicitly correlated full ab initio calculation of the crystal field states in TbPc2 using the software MOLCAS 8.032. As we did not impose any symmetry in the crystal field, we employed the experimental (X-ray) structure measured for the TbPc2 compound in bulk29. It is important to remark that the main finding of our approach (i.e. the activation of the new tunneling mechanisms) is a consequence of the introduction of the low symmetry harmonics and it is therefore very robust with respect to the microscopic details of the molecular structures used in the calculations. First we separately optimize state-averaged Complete Active Space Self-Consistent Field (SA-CASSCF) wavefunctions for all relevant spin symmetries obtained by the Ln-ion f electrons in the active space composed of the seven 4 f atomic orbitals. The optimization takes into account scalar-relativistic Douglas-Kroll Hamiltonian as implemented in the software MOLCAS 8.0. The resulting multiconfigurational wavefunctions, having a well-defined spin symmetry, are then non-perturbatively mixed by suitable spin-orbit Hamiltonian accounting for spin-other-orbit interactions in a mean-field manner (AMFI).

The split J-multiplet thus obtained is then projected onto a full crystal field Hamiltonian, with Stevens operators quantized along the principal magnetic axis obtained from direct calculation of the ground state g-tensor, which provides all 27 crystal-field parameters evaluated ab initio31. Such approach is of course sensitive to the choice of basis set, but previous experience has shown that using the relativistically-corrected ANO-RCC-DZP basis set on the lanthanide ion, and an ANO-RCC-DZ basis set on the lighter elements, leads to quite accurate results, especially when the ground doublet is well-known from experiments to be very close to a pure atomic state30. This is well known to be the case for the three molecules studied here. More details on the ab-initio calculation of the crystal field parameters are given in Supplementary Note 4.

Additional Information

How to cite this article: Candini, A. et al. Spin-communication channels between Ln(III) bis-phthalocyanines molecular nanomagnets and a magnetic substrate. Sci. Rep. 6, 21740; doi: 10.1038/srep21740 (2016).