Journal of
Materials Chemistry C

Materials for optical, magnetic and electronic devices

www.rsc.org/MaterialsC

<figure><img src="image_1.png" /></figure>

ISSN 2050-7526

<figure><img src="image_2.png" /></figure>

PAPER

Yanhua Lan, Svetlana Klyatskaya et al.

Magnetic interplay between two different lanthanides in a tris-phthalocyaninato complex: a viable synthetic route and detailed investigation in the bulk and on the surface
<figure><img src="image_3.png" /></figure>

CrossMark

← click for updates

Cite this: J. Mater. Chem. C, 2015, 3, 9794

# Magnetic interplay between two different lanthanides in a tris-phthalocyaninato complex: a viable synthetic route and detailed investigation in the bulk and on the surface†

Yanhua Lan,*a Svetlana Klyatskaya,*a Mario Ruben,ab Olaf Fuhr,ac Wolfgang Wernsdorfer,de Andrea Candini,f Valdis Corradini,f Alberto Lodi Rizzini,g Umberto del Pennino,g Filippo Troiani,g Loïc Joly,b David Klar,h Heiko Wende,h and Marco Affrontefg

Future applications of molecular units in quantum information technologies require a fine control at the single molecule level. This includes the choice of each functional element, the intramolecular interaction and the robustness of molecules when dispersed on a substrate. Keeping these goals in mind, we designed and synthesized a heterometallic phthalocyaninato-complex including two different lanthanides in each moiety, namely [PcDyPcTbPc*] (Pc being phthalocyanines; and Pc* being 2,3,9,10,16,17,23,24-octahexyl-substituted phthalocyanines). Full magnetic characterization was performed down to the mK temperature range on bulk microcrystals by means of AC susceptibility, DC magnetization (including microSQUID) and specific heat measurements. A weak, yet sizeable, interaction between the two lanthanides is clearly detected by different techniques, altering the magnetic behavior of the single lanthanide as observed in the parent [LnPc₂] complexes. Isolated [PcDyPcTbPc*] molecules dispersed on HOPG and the Au surface by liquid phase deposition are proven to maintain their main chemical and magnetic features by combined XPS, XAS and XMCD analysis and to lie with one Pc ligand flat to the surface. Opening of a small but sizable hysteresis loop at 1.8 K is directly observed on both Tb and Dy sites proving the retention of magnetization at the single molecule level.

Received 5th July 2015.
Accepted 17th July 2015

DOI: 10.1039/c5tc02011e

www.rsc.org/MaterialsC

## Introduction

Bis-phthalocyanines of lanthanides, LnPc₂ complexes, are attracting much interest since they have been shown to be

prototypical molecular units for the realization of hybrid nano-devices and surface studies.¹ In this simple molecular unit, the lanthanide is sandwiched between two phthalocyanines with a well-defined ligand field that in turn gives rise to axial magnetic anisotropy with a huge energy barrier – in the case of [TbPc₂]⁻ as demonstrated in the earlier work of Ishikawa *et al.*² In its neutral form [LnPc₂]⁰ one unpaired electron is delocalized in the two Pc ligands. The neutral derivative is robust enough to be sublimed or deposited by liquid phase. It is typically adsorbed laying with the Pc flat on metal surfaces due to van der Waals interactions³,⁴ or by π-π stacking on the graphitic substrate.⁵ Scanning tunneling microscopy (STM) studies on [TbPc₂] on Au(111) surfaces revealed the Kondo effect mostly localized on the organic ligand,⁶ while X-ray magnetic dichroism (XMCD) on TbPc₂ showed that coupling is possible with magnetic Ni,⁷ᵃ Co⁷ᵇ and Mn⁷ᶜ substrates while preserving the magnetic anisotropy of the single Ln ion. Because of the high anisotropy barrier, one expects to find an opening of the hysteresis loop in magnetization cycles at low temperature, with both quantum tunnelling and relaxation processes related to different features of the single lanthanide, also at a single molecule level as previously observed in Fe₄.⁸ XMCD actually detects the element

<figure><img src="image_4.png" /></figure>

<footnote>a Institute of Nanotechnology (INT), Karlsruhe Institute of Technology (KIT), 76344 Eggenstein-Leopoldshafen, Germany. E-mail: yanhua.lan@kit.edu, svetlana.klyatskaya@kit.edu</footnote>

<footnote>b Universite de Strasbourg, Institut de Physique et de Chimie des Materiaux de Strasbourg, Campus de Cronenbourg, 23 Rue du Loess, 67034 Strasbourg Cedex 2, France</footnote>

<footnote>c Karlsruhe Nano Micro Facility (KNMF), Karlsruhe Institute of Technology (KIT), 76344 Eggenstein-Leopoldshafen, Germany</footnote>

<footnote>d Univ. Grenoble Alpes, Inst NEEL, 25 rue des Martyrs, F-38000 Grenoble, France</footnote>

<footnote>e CNRS, Inst NEEL, F-38000 Grenoble, France</footnote>

<footnote>f CNR Institute of Nanoscience S3, via G. Campi 213A, 41125 Modena, Italy</footnote>

<footnote>g Università di Modena e Reggio Emilia, Dipartimento di Fisica, Informatica, Matematica, via G. Campi 213A, 41125 Modena, Italy</footnote>

<footnote>h Faculty of Physics and Center for Nanointegration Duisburg-Essen (CENIDE), University of Duisburg-Essen, Lotharstraße 1, D-47048 Duisburg, Germany</footnote>

<footnote>† Electronic supplementary information (ESI) available: Full experimental details, crystallographic data, magnetic data, low temperature magnetization measurements, heat capacity measurements, and further spectroscopic data and analysis. CCDC 1054639 (1). For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/c5tc02011e</footnote>
specific magnetization of individual magnetic centers even when isolated on surfaces, although it turns out that the measurement time is relatively long⁹ and the open hysteresis loop can only be observed at very low temperature.¹⁰ [TbPc₂] molecules have also been used to fabricate hybrid spintronic devices such as molecular spin valves with carbon nanotubes¹¹ or graphene¹² or single molecule transistors.¹³ Here single electrons and even nuclear spin states have been addressed. The read out of an individual spin and coherent manipulation has been clearly demonstrated a new method to encode quantum bits with molecular spin.¹⁴

The next step on this line is the realization of two-qubit gates. To this end, molecules comprising two distinguishable magnetic centers weakly coupled among them have been suggested as possible candidates.¹⁵ Two different magnetic centers are required if one wants to address them independently so as to engineer a universal CNOT quantum logic gate. Recently compound [CeEr] prepared by Aromí and co-workers¹⁶ is deemed to meet the quagre requirements, in which both ions have a doubly degenerate magnetic ground state and can be addressed individually. Their isotopes have mainly zero nuclear spin, which enhances the electronic spin coherence. So far there are several reports concerning the preparation of heterodinuclear lanthanide complexes.¹⁷ Of them, two examples¹⁷,ᵇ are related to complexes of phthalocyanines or its substituted derivatives. The synthesis and the magnetic properties of homometallic dinuclear [PcLnPcLnPc] complexes have been reported by Ishikawa¹⁸ first, and then by Yamashita and co-workers,¹⁹ along with preliminary XMCD studies on homometallic [PcLnPcLnPc] on the surface;²⁰ however, the physical properties of heterometallic complexes of this class towards its deposition onto substrates have not been reported yet. With the aim to realize molecular spintronic devices with two-qubit gates, we have designed and synthesized an asymmetric heterometallic phthalocyaninato-Dy(III)-Tb(III) complex with one of the Pc rings decorated with flexible long hydrocarbon chains.

By exploiting a divergent stepwise protocol, we report here the synthesis, single crystal structure and detailed magnetic studies on [PcLn₁PcLn₂Pc*] with Ln₁ ≠ Ln₂ (1) and Ln₁ = Ln₂ (2). We start presenting thermodynamic properties measured on bulk micro-crystals by AC susceptibility and DC magnetization microSQUID and heat capacity to demonstrate that a weak ferromagnetic coupling is established when Ln₁ = Tb and Ln₂ = Dy. Secondly, we report X-ray photoelectron spectroscopy (XPS), X-ray absorption spectroscopy (XAS) and XMCD measurements of [PcDyPcTbPc*] molecules, hereafter named [Dy,Tb], dispersed on Au(111) and HOPG surfaces to prove the robustness of this molecule and to show the opening of a hysteresis loop at 1.8 K. The results demonstrate that the main magnetic features of the molecule are preserved also when they are dispersed on the substrate.

## Experimental

### Synthesis

The complexes PcLn₁PcLn₂Pc* (Scheme 1) were synthesized in a divergent stepwise protocol. First, the reaction of 1,2-dicyano-4,5-di(hexyl)benzene (3) with lithium in dry methanol led to the

di-lithium salt of the Pc ligand with 2,3,9,10,16,17,23,24-octahexyl substituents, Pc* (4). In the presence of Ln(acac)₃, the reaction mixture was worked up by the *in situ* fusing with excess of [Ln₁Pc₂] (5 or 6) affording the targeted crude product. Pure complex 1 as a dark green solid was separated from the crude mixture by column chromatography (basic alumina oxide) with an eluent of CH₂Cl₂/MeOH (10:1). The complexes (1) and (2) are composed of three Pc²⁻ ligands and two Ln³⁺ ions, resulting in a neutral complex with a closed shell π-electron system. 1% of (1) diluted sample was prepared by mixing 1 and 2 in a ratio of 1:100 in CH₂Cl₂. With slow evaporation of the resulting solution, dark green needles were grown for micro-SQUID measurements. Single crystals of 1 suitable for X-ray diffraction analysis were obtained by slow diffusion of EtOH into a solution of complex (1) in hexane as fine dark green needles. In order to steer a preferential absorption side of the heteronuclear complex, only one of the Pc-rings, Pc*, is decorated with flexible long hydrocarbon chains. Further details on synthesis and chemical characterization are reported in the ESI.†

### Physical characterization

**Spectroscopic characterization and single crystal X-ray diffraction.** Complexes are fully characterized by means of UV-vis and near-IR spectroscopies, and mass spectrometry (MALDI ToF) (Fig. S1, ESI†). The experimental and calculated values of the relative abundance of the isotopic ions of the main peak in HR MALDI ToF are absolutely identical demonstrating the absence of any homometallic Dy₂ or Tb₂ complexes. Finally, the structure was confirmed by X-ray crystal structure analysis (Fig. S2, ESI†). When Ln₁ ≠ Ln₂, the X-ray chemical selectivity offers the unique opportunity to study separately the properties of the two lanthanide ions and hence to investigate the intra-molecular Ln–Ln coupling.

**Magnetic measurements.** Magnetic susceptibility measurements were obtained using a Quantum Design SQUID magnetometer MPMS-XL. DC susceptibility measurements were obtained over the temperature range 1.8–300 K under an applied field of 1000 Oe. Magnetization was measured with fields of 0–70 kOe. AC susceptibility measurements were measured with an oscillating AC field of 3 Oe and AC frequencies ranging from 1 to 1500 Hz. Magnetic data were corrected for the sample holder contribution. Low temperature magnetization measurements were carried out on

<figure><img src="image_5.png" /><figcaption>**Scheme 1** Synthesis of [Pc–Ln₁–Pc–Ln₂–Pc*] complexes, Ln₁ = Dy, Ln₂ = Tb (1), Ln₁ = Ln₂ = Y (2), (a) Li/MeOH; (b) Ln(acac)₃; [Ln₃Pc₂]⁰ (5 – Ln = Tb, 6 – Ln = Y); 1-chloconaphthalene, reflux.</figcaption></figure>
single crystals with an array of micro-SQUIDs between 0.04–1.6 K using an applied field with sweep rates of 0.001–0.280 T s⁻¹.²¹

**Specific heat measurements.** Heat capacity measurements were performed by means of the QD-PPMS7T system using the two-tau relaxation method. 1 mg of microcrystalline sample was mixed with Apiezon N grease and glued on the calorimeter.

**Deposition and surface analysis.** Sub-monolayers (MLs) were obtained by immersing HOPG, Au(111) single crystals or the Au/mica flamed-annealed surface in a 10⁻⁵ M solution of (1) using dichloromethane (DCM) as the solvent, then rinsed in DCM and blow dried in a flux of nitrogen gas. Thick films (TF) were obtained by drop casting the saturated solution on the same substrates. STM and XPS²² were used to check that the desired two dimensional distribution of nanometric entities was actually obtained. Room temperature STM image acquisition was carried out in constant current mode under typical conditions of 2.0 V and the lowest achievable current (30 pA) in order to minimize dragging and damaging of the soft organic materials by the scanning tip. XPS measurements were performed using an Omicron hemispherical analyzer (EA125) and a non-monochromatized Al-Kα X-ray source (hν = 1486.6 eV).

**XAS and XMCD measurements.** Soft XAS and XMCD measurements were performed at the SIM-X11MA beamline of Swiss Light Source (SLS), Paul Scherrer Institut (PSI), Villigen (CH) proposal number 20140289. The lowest sample temperature reached was ~1.8 K²³ and the base pressure of the experimental chamber was 1.0 × 10⁻¹⁰ mbar. We paid attention to avoid any sample degradation induced by radiation exposure, working at a very low flux (below 10¹⁰ photons per s) and by strictly monitoring the XAS spectra throughout all the experiments for detecting even the smallest indication of sample damaging. XMCD measurements at the Tb-M₄.₅ and Dy-M₄.₅, edges were performed in total electron yield mode using circularly and linearly polarized light with ~100% polarization rate and with external magnetic fields µ₀H up to 6 T, applied parallel to the incident photon beam (see Scheme 2).

The dichroic XMCD signal (expressed in %) is evaluated by taking the difference between the two XAS spectra obtained with different X-ray circular polarizations (σ↑↓-σ↑↑) and

dividing by the height of the average of the two polarizations (edge jump). Hence, the XMCD value expressed in percentage (%) does not depend on the amount of the material analyzed and thus the signal of the ML can be directly compared with that obtained on the thick films, whilst the quantification of the dichroic signal in terms of absolute magnetic moments is less straightforward. According to that reported in the literature for Tb³⁺ and Dy³⁺ ions, we simply assume here that the area of the XMCD curve is proportional to the magnetization. We also checked that the shape of the dichroic signal does not change with the magnetic field, therefore the height of the M₅ peak can be considered proportional to the area of the dichroic signal. We can therefore obtain fast detection of the dichroic signal by simply measuring the edge (E) and pre-edge (P) intensities, when the magnetic field is swept under isothermal conditions:

$$ \frac{[(E\sigma^{\uparrow\downarrow} - P\sigma^{\uparrow\downarrow}) - (E\sigma^{\uparrow\uparrow} - P\sigma^{\uparrow\uparrow})]}{1/2[(E\sigma^{\uparrow\downarrow} - P\sigma^{\uparrow\downarrow}) + (E\sigma^{\uparrow\uparrow} - P\sigma^{\uparrow\uparrow})]} $$

here E is the intensity of the XAS at the energy which corresponds to the maximum of the XMCD signal and P is the level of the background preceding the absorption edge (see Fig. S16a, ESI†).

## Results and discussion

### Structural description

Complex (1) crystallizes in the monoclinic space group P2₁/c with four molecules per unit cell (Table S1, ESI†). The terbium and dysprosium ions occupy central positions in the complex and are eightfold coordinated by the isoindole nitrogen atoms (Niso) of the phthalocyanine ligands: Tb–N2, N4, N6, N8 (Pc* in gray) and N10, N12, N14, N16 (Pc in orange); Dy–N10, N12, N14, N16 (Pc in orange) and N18, N20, N22, N24 (Pc in red). Both outer Pc and Pc* ligands of 1 are equally distorted from planarity and, therefore, adopt a biconcave shape. The central terbium ion lies 1.782 Å from the Niso mean plane of the Pc ligand and 1.262 Å from the Niso mean plane of the Pc* ligand bearing peripheral substituents, whereas dysprosium ion lies 1.262 Å from the Niso mean plane of the outer Pc ligand (in red, Fig. 1) and 1.727 Å from the Niso mean plane of the middle Pc ligand (in orange, Fig. 1). The calculated distances between the three Niso mean planes are 2.984 and 3.150 Å respectively. The intramolecular distance between Dy and Tb ions in the same unit cell is 3.508 Å, which is slightly smaller than for the homonuclear [PcLnPcLnPc] complex (3.52 Å).¹⁸ᵃ The twist angles between the outer rings Pc/Pc* and the center Pc ring were determined to be 46.4° and 32.2° respectively, (Fig. S3a, ESI†), causing a pseudo fourfold axis (the direction of the uniaxial magnetic anisotropy) perpendicular to the Pc rings.

Complex (1) crystallizes with molecules of hexane in the crystal lattice. The shortest intermolecular Dy–Tb distance along the *a* axis was determined to be 14.40 Å, while being 17.95 Å within two neighboring columns, Fig. S2 (ESI†). Each molecule of (1) is rather well separated from the neighboring molecule due to the *n*-hexyl chains with all Dy–Tb-axes pointing collinearly. The sizes of (1) were estimated to be ~15 Å (non-substituted Pc) and 27 Å (the Pc-ring with hexyl chains), on the

<figure><img src="image_6.png" /><figcaption>**Scheme 2** Schematic view of the system under investigation. [Dy,Tb] (1) molecule deposited on Au(111) or HOPG surfaces. The beam incidence angle theta (θ) can be varied from normal incidence (θ = 0°) to grazing incidence (θ = 60°).</figcaption></figure>

<figure><img src="image_7.png" /></figure>
<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /><figcaption>**Fig. 1** ORTEP plot of one molecule of [Dy,Tb] (1) with ellipsoids drawn at 50% level of probability for all non-hydrogen atoms, indicating the numbering scheme; top view (left), atoms labelled 'n' represent the carbon atoms Cn. Side view (right), the chains at C12/21/28 have been omitted for clarity. There is disorder in the hexylchains at C5, C12 and C21. The disorder model and the positions of the hydrogen atoms are excluded. Selected bond lengths [Å]: Tb(1)–N(7) 2.334(4), Tb(1)–N(1) 2.343(4), Tb(1)–N(5) 2.344(4), Tb(1)–N(3) 2.355(4), Tb(1)–N(11) 2.614(4), Tb(1)–N(9) 2.623(4), Tb(1)–N(13) 2.635(3), Tb(1)–N(15) 2.636(4), Dy(1)–N(19) 2.339(4), Dy(1)–N(17) 2.341(5), Dy(1)–N(21) 2.344(4), Dy(1)–N(23) 2.352(4), Dy(1)–N(13) 2.582(4), Dy(1)–N(11) 2.591(4), Dy(1)–N(9) 2.592(3), Dy(1)–N(15) 2.596(4).</figcaption></figure>

basis of the distance between two hydrogen atoms of *n*-hexyl chains at both ends and between the hydrogen atoms of the upper non-substituted Pc rings, respectively. The height of ~7 Å is determined as the distance between overlapping hydrogen atoms at meta-positions of upper and lower rings (Fig. S3b, ESI†). In the crystal structure, (1) is arranged in columns along the *a* axis in a π-π stack arrangement (3.54 Å, Fig. S2, ESI†).

## Magnetic susceptibility

The static magnetic susceptibility of (1) was measured on a poly-
crystalline sample in the temperature range 1.8–300 K in an applied
magnetic field of 1000 Oe. The χT value of 25.77 cm³ K mol⁻¹ at
300 K is very close to the expected value of 25.99 cm³ K mol⁻¹ for the
system containing one isolated Dy³⁺ (*J* = 15/2, *g* = 4/3, ⁶H₁₅/₂, *C* =
14.17 cm³ K mol⁻¹) and one Tb³⁺ ion (*J* = 6, *g* = 3/2, ⁷F₆, *C* =
11.82 cm³ K mol⁻¹) (Fig. 2). Upon cooling, the χT product slightly
decreases to reach a minimum value of 23.38 cm³ K mol⁻¹ at
22 K as a consequence of thermal depopulation of the Stark

sublevels of the anisotropic Ln³⁺ ions. Below 22 K, the χT product sharply increases to reach 32.37 cm³ K mol⁻¹ at 1.8 K, indicative of ferromagnetic coupling between Dy³⁺ and Tb³⁺ ions. The presence of ferromagnetic interactions between the two Ln³⁺ ions is also evident in the in-phase AC magnetic susceptibility at low frequency (Fig. 2 inset). This observation is consistent with the results reported for the homometallic dinuclear Tb-phthalocyaninato-complexes.¹⁹

The dynamic properties of (1) were also investigated using
temperature and frequency dependent AC susceptibility mea-
surements. The strong frequency dependence of both in-phase
and out-of-phase components observed in zero DC field below
30 K manifests slow magnetization dynamics (Fig. 3a and b).
Taking the maximum of the χ''(T) curve as a criterion to char-
acterize the relaxation, the blocking temperature at 1000 Hz is

<figure><img src="image_10.png" /><figcaption>**Fig. 2** Temperature dependence of the χT product of complex (1) at 1000 Oe (with χ being the molar susceptibility defined as M/H). Inset: frequency and temperature dependence of the in-phase AC magnetic susceptibility.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Fig. 3** Temperature dependence of the in-phase (a) and out-of-phase (b) components of the AC magnetic susceptibility of complex (1) in zero dc field at different frequencies; the solid lines are a guide to the eye. Argand plots in zero dc field between 2 and 6 K (c) and 7 and 13 K (d); the solid lines represent the least-squares fit obtained using a generalized Debye model. The parameters are discussed in the text.</figcaption></figure>

<figure><img src="image_12.png" /></figure>
<figure><img src="image_13.png" /><figcaption>Fig. 4 The relaxation time $\tau$ as a function of $1/T$ of complex (1) extracted from AC susceptibilities between 1.8 and 15 K under zero dc field. The two solid lines represent two thermally activated regimes with $\Delta E_1 = 6.9$ K and $\tau_1 = 1.5 \times 10^{-3}$ s in the temperature range of 1.8–5.0 K and $\Delta E_2 = 43.2$ K and $\tau_2 = 6.9 \times 10^{-6}$ s between 10 and 15 K.</figcaption></figure>

detected at 3.9 K. Besides the well-defined maxima in the out-of-phase components in the low temperature range, a broad maximum can also be observed when the frequency was increased up to 800 Hz at higher temperature (Fig. S4a, ESI†).

The relaxation time ($\tau$) against $1/T$ was extracted from the frequency sweeping AC data between 1.8 and 15 K and is plotted in Fig. 4 and Fig. S7 (ESI†). The nonlinear relaxation plot suggests the presence of multiple relaxation pathways. Data analysis²⁴ revealed that the thermal variation of $\tau$ can be approximated by two exponential laws, characteristics of thermally activated regimes, with $\Delta E_1 = 6.9$ K and $\tau_1 = 1.5 \times 10^{-3}$ s in the temperature range of 1.8–5.0 K and $\Delta E_2 = 43.2$ K and $\tau_2 = 6.9 \times 10^{-6}$ s between 10 and 15 K. Two activated regimes have also been observed in the case of homometallic [PcTbPcTbPc] complexes.¹⁹

For both [TbPc₂] and [DyPc₂] complexes only one relaxation time is observed with one characteristic activation energy and saturation at low temperature due to the presence of quantum tunneling. The observation of the thermally activated regime with a low energy barrier ($\Delta E_1 = 6.9$ K) denotes the partial suppression of quantum tunneling – with respect to the single [LnPc₂] complexes (see below) and the presence of low energy barrier that, interestingly, is comparable to the splitting of the Ising doublet due to the intramolecular interaction.

To see if it is possible to further characterize these relaxation processes, an Argand plot of the in-phase versus out-of-phase susceptibilities was constructed. If only one single relaxation process is active, the plot would have a semicircular shape with a vanishing $\alpha$ value. Fitting the data using a generalized Debye model leads instead to a $\alpha$ value of 0.26–0.30 between 2 and 5 K (Fig. 3c and Table S2, ESI†) and 0.13–0.24 above 10 K (Fig. 3d and Table S2, ESI†), suggesting that relaxation mechanism operating in the two temperature regimes is more complex than single activated processes.

## Low temperature magnetization

Isothermal magnetization was measured at very low temperature on small crystals of both (1) and 1% of (1) diluted in [Y,Y]

(2) by micro-SQUID as a function of magnetic field applied along the principal anisotropy axis. Opening of hysteresis loops is evident due to blocking of magnetization at the typical sweeping rates used in these experiments (Fig. 5a and b and Fig. S9, ESI†). The general behavior of (1) is similar to those reported for homometallic [PcLnPcLnPc] complexes.¹⁸,¹⁹ The intramolecular interactions between Dy³⁺ and Tb³⁺ have some clear influence on the step height and width of the hysteresis loops. For comparison, the behavior (i.e. the sum of the two contributions) of isolated [TcPc₂] and [DyPc₂] complexes is also reported in Fig. S10 (ESI†), showing clear differences.

An approximate Zeeman diagram of the ground states of [Dy,Tb] (1) is calculated²⁵ (Fig. 5c) based on the energy splitting of 4.8 K obtained from point-dipole calculations (see below). The red levels correspond to the ferromagnetic state, with Tb and Dy pointing in the same direction, whereas the green levels correspond to the antiferromagnetic alignment. The dotted lines represent the tunnel transitions at avoided level crossings and the big arrows correspond to direct relaxation between energy levels, which becomes faster and faster at higher fields. The nuclear spin of Dy is neglected for clarity, because the hyperfine coupling of Dy is much smaller than that of Tb and some isotopes of Dy have no nuclear spin.

The Zeeman diagram reveals that the large step at zero field corresponds to a co-tunneling process where both Tb and Dy flip their magnetic moments. Four steps are visible at this low field transition (Fig. 5d), which are due to the four nuclear spin states of the Tb ion. Such “hyperfine steps” were also observed for the [TbPc₂] complexes (Fig. S11b, ESI†).²b However, because of intramolecular coupling of Tb–Dy, the field separation between these steps are about two times smaller than for the [TbPc₂] complexes.

The crossing between the ferromagnetic ground state and the antiferromagnetic excited state induces a fast relaxation at about 0.4 T (Fig. 5c), which is confirmed by the hysteresis loops of the 100% sample, showing clearly the step at about the same field (Fig. 5a). This is not seen in the diluted sample because the direct relaxation processes are enhanced by small disorder, reversing all spins before reaching 0.4 T (Fig. 5d).

## Specific heat

In order to get quantitative evaluation on the strength of the coupling we have measured specific heat at low temperature. The specific heat of non-magnetic [Y,Y] (2) can be taken as reference since it shows only pure lattice contribution down to the lowest temperatures (Fig. 6). An anomaly with a maximum around 1.2 K is well visible for [Dy,Tb] (1) instead. This shifts with the magnetic field (Fig. S12, ESI†). A fit with a simple two-level Schottky $C(T)$ gives an effective energy gap of 3.8 K. In the point-dipole approximation, an energy splitting of ≈4.8 K can be estimated between two configurations of the two Ising spins of Dy and Tb separated by 3.5 Å. The agreement between the two values is reasonably good taking into account that we neglect the degeneracy of the levels. Heat capacity measurements on 10% [Dy,Tb] (1) diluted in the non-magnetic [Y,Y] (2) sample (Fig. S13, ESI†) shows a small, but sizable Schottky anomaly with maximum at essentially the same position (1.2 K)

<figure><img src="image_14.png" /></figure>
<figure><img src="image_15.png" /></figure>

<figure><img src="image_16.png" /><figcaption>**Fig. 5** Sweeping rate dependence of magnetization measured as a function of applied magnetic field for (a) 100% of (1) and (b) 1% of (1) diluted in [Y,Y] (2). (c) Zeeman diagram of the ground states of [Dy,Tb] (1). The dotted lines indicate resonant tunnel transitions and the large arrows indicate the possibility of direct transitions with phonon emission. (d) Enlargement at low-fields of (b) showing four steps corresponding to ground state tunnel transition of the Tb³⁺ ion, which depends on the four nuclear spin states.</figcaption></figure>

<figure><img src="image_17.png" /><figcaption>**Fig. 6** Specific heat in R units ($R = 8.314$ J mol⁻¹ K⁻¹ is the gas constant) measured on 0.5 mg of microcrystals of [Dy,Tb] (1, black open circle) and non-magnetic [Y,Y] (2, blue open squares). The red line is the calculated Schottky anomaly contribution.</figcaption></figure>

as in the case of the non-diluted (1) derivative, indicating that intermolecular interactions are not particularly relevant here. DFT calculations³ suggest that through-bond exchange interaction may well occur in Ln-phthalocyaninato-complexes. Yet, in this case the neutral molecule (1) has no electrons delocalized over the Pc planes and this pathway is also expected to give small contribution. Thus, we conclude that dipolar Dy–Tb intra-molecular interaction is the dominant coupling mechanism.

## Molecules on the surface

While (micro-)SQUID measurements detect the magnetization of the whole molecule, XMCD measurements are element-sensitive

and allow one to discern the magnetization of the different lanthanides. X-ray absorption and magnetic dichroism have been measured on [Dy,Tb] (1) molecules dispersed on HOPG and Au surfaces. Preliminary XPS analysis (Fig. S15, ESI†) demonstrated that molecules are intact after the deposition by liquid phase and the chemical stoichiometry is preserved as reported in Table S3 (ESI†). STM images (Fig. S14b, ESI†) show isolated spots with a size consistent with what determined by X-ray diffraction indicating that the molecules do not aggregate. Note that in such a situation intermolecular interaction can be neglected, while the interaction with different substrates can be significant and needs to be checked. Thus, to further test the integrity of the lanthanide core of the molecule, we compared the Tb-M₄.₅ and Dy-M₄.₅ XAS collected on the MLs on Au(111) and on HOPG with the corresponding spectra obtained on the thick films and on the powder (Fig. S16, ESI†). The XAS and XMCD spectra of the powder perfectly match those obtained on the thick films (Fig. S16a, ESI†). This confirms that the (1) core is stable in solution. Secondly, the XAS/XMCD spectra of the ML on gold perfectly match those obtained on MLs on HOPG (Fig. S17, ESI†) for both normal and grazing incidence (θ = 0° and θ = 60°). This shows that the interaction with the gold surface does not affect the valence electronic structure of the core (i.e. Tb and Dy ions remain trivalent). The different intensities of the XMCD signals at 0° and 60° (Fig. S17, ESI†) indicate a strong magnetic anisotropy of the molecules, consistently with the fact that they are flat on the surface, with the Pc plane parallel to the substrate, as actually depicted in Scheme 2. This point is further confirmed by X-ray linear dichroism (XLD) at the Tb M₅, Dy M₅, and N K edges at θ = 60° (Fig. S18, ESI†).
<figure><img src="image_18.png" /></figure>

<figure><img src="image_19.png" /><figcaption>Fig. 7 XMCD signal of complex (1) measured at the Tb and Dy $M_5$ edges respectively at 1.8 K as a function of applied magnetic field. Molecules of (1) have been dispersed on the HOPG substrate and lie with one Pc flat on the substrate. The external magnetic field was applied at an angle of 0° (left panel) and 60° (right panel) with respect to the HOPG surface. Continuous and dashed lines are the calculations with and without dipolar interaction, respectively.</figcaption></figure>

The magnetization loop was derived by measuring the XMCD signal as a function of an applied magnetic field at the Tb and Dy $M_5$ edges for the ML on HOPG (Fig. 7), the ML on gold (Fig. S19, ESI†) and a thick film of derivative (1) (Fig. S20, ESI†). Two different angles are compared: normal incidence ($\theta = 0°$) and grazing incidence ($\theta = 60°$). To evidence the effect of intramolecular Dy–Tb interaction, we evaluated the magnetization of the single ion with the crystal field parameters reported in ref. 18. For comparison we also plotted in Fig. 7 the estimated value of magnetization in the presence of dipolar Dy–Tb interaction, using the Dy–Tb distance value as estimated by X-ray diffraction. The effect of dipolar interaction is small but quite visible. The agreement between the experimental data

<figure><img src="image_20.png" /><figcaption>Fig. 8 XMCD signal of complex (1) measured at the Tb and Dy $M_5$ edges respectively at 1.8 K as a function of applied magnetic field. This zoom of Fig. 7 evidences the opening of hysteresis loops in the *M*(*H*) cycle.</figcaption></figure>

and the calculated magnetization is remarkable, which confirms that through-space (dipolar) intramolecular interaction is active and that the crystal field parameters are close to those deduced for bulk samples.

In Fig. 8 we plot a zoom of magnetization curves reported in Fig. 7 in a smaller magnetic field range. Sizable opening of hysteresis is visible at both Dy and Tb sites for $\theta = 0°$, although the coercive fields are smaller than what observed in micro-SQUID measurements. That is due in part to the higher temperature and also to the slower sweeping rate in XMCD experiments (minutes for one cycle). Note that the small but finite opening of the hysteresis is visible for $\theta = 0$ for 1 deposited on HOPG but it is vanishingly small for $\theta = 60°$ and for MLs of 1 deposited on Au(111) (data reported in Fig. S19, ESI†) and for a TF (Fig. S20, ESI†). The vanishingly small hysteresis observed in 1 on HOPG at $\theta = 60°$ and TF (Fig. S20, ESI†) is not surprising due to the presence of a transverse field. In the case of 1 on gold (Fig. S19, ESI†) it is clear that interaction with the metal surface affects somehow the magnetic properties of 1. Thus, we conclude that the magnetic anisotropy of both Dy and Tb is essentially preserved when isolated molecular units are deposited on HOPG while some change can occur when 1 is deposited on gold. It is worth stressing however that it is the first time – to our knowledge – that memory effects are independently observed on two distinct magnetic centers in the same molecular units on the timescale of experiments at synchrotron radiation facility (minutes).

## Conclusions

By using a stepwise synthesis, we were able to obtain a hetero-metallic [Dy,Tb] complex (1) and determine its molecular structure by single crystal X-ray diffraction. The control of the nature of each moiety opens up the possibility to distinguish lanthanides with different magnetic features in the same molecular units. The weak but finite ferromagnetic coupling between the two magnetic centers affects the magnetic properties of this derivative. To show this, we performed a systematic study of DC magnetization, AC susceptibility and specific heat at low temperature in microcrystalline bulk samples. Data analysis shows the presence of two thermally activated regimes and quantum tunneling particularly active at very low temperatures. Due to the intramolecular interaction, mostly of dipolar origin, the low temperature magnetization dynamics is affected by a co-tunneling process.

Complex (1) is soluble in solution and we deposited it on HOPG and Au surfaces in order to study how the magnetic features are preserved when molecules are dispersed on the surface. XAS and XMCD experiments, performed at low temperatures, show that these molecules preserve their stoichiometry and the magnetic coupling between the two $Ln^{3+}$ is still active within isolated molecular units. Opening of a small but sizable hysteresis loop was detected at 1.8 K on isolated molecules dispersed on HOPG. This is a remarkable result since, except for the case of $Fe_4^8$ and $LnPc_2$ on HOPG,¹⁰ few
<figure><img src="image_21.png" /></figure>

experiments have been reported so far to directly prove the retention of magnetization in isolated molecular objects. We believe that these results may contribute to define strategies to process quantum information at the single molecule level.¹⁴

## Acknowledgements

This work has been partially supported by the EU by the FP7 FET-Proactive project MoQuaS No. 610449, the SFB TR 88 “3MET” and the Agence Nationale de la Recherche project MolQuSpin, No. ANR-13-BS10. We acknowledge PSI for provision of synchrotron facility and Armin Kleibert for technical assistance. We also thank Prof. Annie K. Powell and Dr Valeriu Mereacre for the allocation of technical equipment.

## Notes and references

1. Y. Lan, S. Klyatskaya and M. Ruben, in *Lanthanides and Actinides in Molecular Magnetism, "Bis-Phthalocyaninato Lanthanide(m) Complexes - From Molecular Magnetism to Spintronic Devices"*, ed. R. Layfield and M. Murugesu, Wiley-VCH, Weinheim, 2015.

2. (a) N. Ishikawa, M. Sugita, T. Okubo, N. Tanaka, T. Iino and Y. Kaizu, *Inorg. Chem.*, 2003, **42**, 2440–2446; (b) N. Ishikawa, M. Sugita and W. Wernsdorfer, *Angew. Chem., Int. Ed.*, 2005, **44**, 2931–2935; (c) F. Branzoli, P. Carretta, M. Filibian, G. Zoppellaro, M. J. Graf, J. R. Galan-Mascaros, O. Fuhr, S. Brink and M. Ruben, *J. Am. Chem. Soc.*, 2009, **131**, 4387–4396.

3. L. Vitali, S. Fabris, A. M. Conte, S. Brink, M. Ruben, S. Baroni and K. Kern, *Nano Lett.*, 2008, **8**, 3364–3368.

4. T. Komeda, H. Isshiki, J. Liu, K. Katoh, M. Shirakata, B. K. Breedlove and M. Yamashita, *ACS Nano*, 2013, **7**, 1092–1099.

5. (a) M. Gonidec, R. Biagi, V. Corradini, F. Moro, V. De Renzi, U. del Pennino, D. Summa, L. Muccioli, C. Zannoni, D. B. Amabilino and J. Veciana, *J. Am. Chem. Soc.*, 2011, **133**, 6603–6612; (b) R. Biagi, J. Fernandez-Rodriguez, M. Gonidec, A. Mirone, V. Corradini, F. Moro, V. De Renzi, U. del Pennino, J. C. Cezar, D. B. Amabilino and J. Veciana, *Phys. Rev. B: Condens. Matter Mater. Phys.*, 2010, **82**, 224406.

6. T. Komeda, H. Isshiki, J. Liu, Y.-F. Zhang, N. Lorente, K. Katoh, B. K. Breedlove and M. Yamashita, *Nat. Commun.*, 2011, **2**, 217.

7. (a) A. Lodi Rizzini, C. Krull, T. Balashov, J. J. Kavich, A. Mugarza, P. S. Miedema, P. K. Thakur, V. Sessi, S. Klyatskaya, M. Ruben, S. Stepanow and P. Gambardella, *Phys. Rev. Lett.*, 2011, **107**, 177205; (b) D. Klar, S. Klyatskaya, A. Candini, B. Krumme, K. Kummer, P. Ohresser, V. Corradini, V. de Renzi, R. Biagi, L. Joly, J.-P. Kappler, U. Del Pennino, M. Affronte, H. Wende and M. Ruben, *Beilstein J. Nanotechnol.*, 2013, **4**, 320–324; (c) A. Lodi Rizzini, C. Krull, T. Balashov, A. Mugarza, C. Nistor, F. Yakhou, V. Sessi, S. Klyatskaya, M. Ruben, S. Stepanow and P. Gambardella, *Nano Lett.*, 2012, **12**, 5703–5707.

8. (a) M. Mannini, F. Pineider, C. Danieli, F. Totti, L. Sorace, Ph. Sainctavit, M.-A. Arrio, E. Otero, L. Joly, J.-C. Cezar, A. Cornia and R. Sessoli, *Nature*, 2010, **468**, 417–421; (b) L. Malavolti, V. Lanzilotto, S. Ninova, L. Poggini, L. Cimatti, B. Cortigiani, L. Margheriti, D. Chiappe, E. Otero, P. Sainctavit, F. Totti, A. Cornia, M. Mannini and R. Sessoli, *Nano Lett.*, 2015, **15**, 535–541.

9. S. Stepanow, J. Honolka, P. Gambardella, L. Vitali, N. Abdurakhmanova, T.-C. Tseng, S. Rauschenbach, S. Tait, L. V. Sessi, S. Klyatskaya, M. Ruben and K. Kern, *J. Am. Chem. Soc.*, 2010, **132**, 11900–11901.

10. D. Klar, A. Candini, L. Joly, S. Klyatskaya, B. Krumme, P. Ohresser, J. P. Kappler, M. Ruben and H. Wende, *Dalton Trans.*, 2014, **43**, 10686–10689.

11. M. Urdampilleta, S. Klyatskaya, J.-P. Cleuziou, M. Ruben and W. Wernsdorfer, *Nat. Mater.*, 2011, **10**, 502–506.

12. A. Candini, S. Klyatskaya, M. Ruben, W. Wernsdorfer and M. Affronte, *Nano Lett.*, 2011, **11**, 2634–2639.

13. R. Vincent, S. Klyatskaya, M. Ruben, W. Wernsdorfer and F. Balestro, *Nature*, 2012, **488**, 357–360.

14. S. Thiele, F. Balestro, R. Ballou, S. Klyatskaya, M. Ruben and W. Wernsdorfer, *Science*, 2014, **344**, 1135–1138.

15. G. Aromì, D. Aguila, P. Gamez, F. Luis and O. Roubeau, *Chem. Soc. Rev.*, 2012, **41**, 537–546.

16. D. Aguilà, L. A. Barrios, V. Velasco, O. Roubeau, A. Repollés, P. J. Alonso, J. Sesé, S. J. Teat, F. Luis and G. Aromí, *J. Am. Chem. Soc.*, 2014, **136**, 14215–14222.

17. (a) K. Katoh, R. Asano, A. Miura, Y. Horii, T. Morita, B. K. Breedlove and M. Yamashita, *Dalton Trans.*, 2014, **43**, 7716–7725; (b) P. Zhu, N. Pan, R. Li, J. Dou, Y. Zhang, D. Y. Y. Cheng, D. Wang, D. K. P. Ng and J. Jiang, *Chem. – Eur. J.*, 2005, **11**, 1425–1432; (c) J.-P. Costes and F. Nicodème, *Chem. – Eur. J.*, 2002, **8**, 3442–3447; (d) R. Sato, K. Suzuki, M. Sugawa and N. Mizuno, *Chem. – Eur. J.*, 2013, **19**, 12982–12990.

18. (a) N. Ishikawa, T. Iino and Y. Kaizu, *J. Phys. Chem. A*, 2002, **106**, 9543–9550; (b) N. Ishikawa, T. Iino and Y. Kaizu, *J. Am. Chem. Soc.*, 2002, **124**, 11440–11447; (c) N. Ishikawa, S. Otsuka and Y. Kaizu, *Angew. Chem., Int. Ed.*, 2005, **44**, 731–733.

19. (a) K. Katoh, T. Kajiwara, N. Nakano, Y. Nakazawa, W. Wernsdorfer, N. Ishikawa, B. K. Breedlove and M. Yamashita, *Chem. – Eur. J.*, 2011, **17**, 117–122; (b) K. Katoh, H. Isshiki, T. Komeda and M. Yamashita, *Coord. Chem. Rev.*, 2011, **255**, 2124–2148; (c) K. Katoh, Y. Horii, N. Yasuda, W. Wernsdorfer, K. Toriumi, B. K. Breedlove and M. Yamashita, *Dalton Trans.*, 2012, **41**, 13582–13600.

20. A. Lodi Rizzini, C. Krull, A. Mugarza, T. Balashov, C. Nistor, R. Piquerel, S. Klyatskaya, M. Ruben, P. M. Sheverdyaeva, P. Moras, C. Carbone, C. Stamm, P. S. Miedema, P. K. Thakur, V. Sessi, M. Soares, F. Yakhou-Harris, J. C. Cezar, S. Stepanow and P. Gambardella, *Surf. Sci.*, 2014, **630**, 361–374.

21. W. Wernsdorfer, *Adv. Chem. Phys.*, 2001, **118**, 99–190.

22. *Handbook of X-ray Photoelectron Spectroscopy*, Perkin-Elemer Corporation, 1979.

23. I. Letard, P. Sainctavit, J.-P. Kappler, P. Ghigna, D. Gatteschi and B. Doddi, *J. Appl. Phys.*, 2007, **101**, 113920.

24. R. J. Glauber, *J. Math. Phys.*, 1963, **4**, 294–307.

25. M. Urdampilleta, S. Klyatskaya, M. Ruben and W. Wernsdorfer, *ACS Nano*, 2015, **9**, 4458–4464.