Tailored homo- and hetero- lanthanide porphyrin dimers: a synthetic strategy for integrating multiple spintronic functionalities into a single molecule

Abstract

We present the design, synthesis and magnetic properties of molecular magnetic systems that contain all elements necessary for spin-valve control in molecular spintronic devices in a single molecule. We investigate the static and dynamic magnetic properties and quantum spin properties of butadiyne-linked homo- and hetero-nuclear lanthanide-porphyrin dimers. A heterometallated porphyrin dimer containing both Tb^III and Dy^III centres is created rationally by the stepwise oxidative homocoupling of distinct lanthanide-porphyrin monomers. Tb^III and Dy^III mononuclear porphyrin complexes, homodimers and heterodimers all exhibit slow magnetic relaxation below 10 kelvin under a static magnetic field. The coherence times for Gd^III porphyrin monomers and dimers are found to be in excess of 3.0 μs at 2 K, allowing distinct magnetic manipulations in low temperature transport experiments.

---

Introduction

Covalently linked porphyrin oligomers constitute a thriving area of research, particularly as model systems for light harvesting, molecular electronics and catalysis.¹ Hetero-metallic porphyrin systems have been studied as donor–acceptor dyads.² They would also be extremely appealing for the creation of molecular spintronic devices,³ which are typically composed of a spin polariser and a spin analyser, setting and reading out the current, respectively. The ultimate miniaturisation limit of these technology lies in the atomic and molecular level, and there is considerable research effort aimed at reducing spintronic elements down to the single-molecule level.⁴ The magnetic complexes integrated into molecular spintronic devices have mainly been single-molecule magnets (SMMs), with a particular role played by bis(phthalocyaninato) rare earth sandwich complexes.⁵ At sufficiently low temperatures, SMMs exhibit molecular hysteresis and behave as magnetic switches that can assume two classical states with a long relaxation time from one to the other.⁶ The observed slow relaxation could thus allow integrating the spintronic analyser and polariser inside a single molecule.

The development of a molecular spin valve requires the synthesis of rationally designed ad hoc systems that contain two different lanthanide-based SMMs. Recent advances in the synthesis of lanthanide hetero-dinuclear complexes have yielded prototypes for two-qubit molecular spin processors, which require coupled lanthanide centres.⁷ In contrast, for the aims of a molecular spin-valve, the spin centres must be magnetically shielded from one another, so as to operate independently. Moreover, each of the magnetic centres needs to be in contact with an electronic quantum dot, in order to optimise the spintronic transport signal. For these purposes the highly developed chemistry of porphyrin systems⁸ including homo-lanthanide dimers,^8f–h may offer the degree of flexibility and reliability that is necessary for the rational tailoring of the magnetic molecules.

However, the use of metalloporphyrin oligomers in molecular magnetism⁹ is only scarcely investigated. Moreover, the creation of reliably hetero-dinuclear lanthanide complexes, requires the layout and implementation of a new methodology, in order to produce an original and interesting architecture.¹⁰ In particular, it is necessary to build a synthetic route that includes the introduction of protecting groups of very different polarity onto the mononuclear building blocks. Only in this way the subsequent separation by chromatography of the hetero-dinuclear assemblies becomes possible, in a reliable procedure.

Here we present the first synthetic approach to the rational creation of such molecular systems, concentrating on hetero-metallic rare-earth based porphyrin dimers. We analyse the magnetic properties of the independent building blocks and then concentrate on Dy–Tb dimers, where each lanthanide component acts as an individual magnetic source. The synthetic strategy is optimised to retain the presence of two electronic quantum dots in series, in the overall linear arrangement necessary to bridge electronic nanogaps. We investigate the quantum coherent properties of these molecular systems, paving the way to the creation of coherent spintronic devices. By pulsed electron paramagnetic resonance (EPR) we quantify the spin coherence times of Gd-based monomers and dimers, the magnetic interactions and the quantum features relevant for coherent manipulation.

Results and discussion

Molecular design and synthesis

Molecules for molecular spin valves need to possess two magnetic elements that can be tuned by an external magnetic field and can retain their magnetisation state for times as long as the measurement. Different conductivities can be obtained when the spins of the metal centres adopt a parallel or antiparallel configuration. This would constitute a fundamental, conceptual step forward in miniaturising devices to a truly molecular level. In addition, this strategy has the advantage of reducing the number of costly nanofabrication steps, by removing the need of magnetic electrodes.

The main requirements that new systems need to fulfil are rather stringent: the magnetic centres need to be reliably slow-relaxing, i.e. they must both have spin-flip times (T₁) comparable to the experimental timescale (i.e.>10 s) in the temperature range used for single-molecule transport measurements (sub-K, for our purposes); high magnetic anisotropy barriers; the molecule must allow electrical conductance via single-electron processes, so that the molecular electronic quantum dots must be easy to reduce or oxidise; on the contrary, the spin centres must not easily undergo oxidation–reduction reactions, in order to preserve the magnetic state while in operation; the spin centres must be different, and show different coercive fields; the quantum dots must be connected to the spin centres; the system must be stable enough to permit deposition on different surfaces; it is desirable that a coherent counterpart can also be created, with dephasing times (T₂) as long as possible.

We use porphyrin ligands in order to exploit their rich chemistry,⁸ for their high chemical stability and tendency to crystallise (Fig. 1). Lanthanide insertion can be achieved, as previously reported,⁹ by heating the free-base porphyrin in the presence of an excess Ln^IIICl₃ (Ln = Y, Gd, Tb, Dy) in a melt of diphenylether and imidazole. While many reports in the literature⁸ use molten imidazole as a solvent, diphenylether was an essential co-solvent for dissolving our porphyrins. Lanthanides are most commonly found in the +3 oxidation state and tend to adopt high coordination numbers both in the solid state and in solution due to their large ionic radii. Since porphyrins form di-anionic ligands which coordinate to metal ions through four nitrogen atoms, the Ln^III ions are not satisfied by a single porphyrin ligand, and one or more supplementary ligands are required in the design of the molecule to balance the charge and account for the remainder of the coordination sphere.

Fig. 1 General structure of P1_Ln and the synthetic route of the heterometallated dimer P2_Dy-Tb (Ln = Tb). Reaction conditions: (i) TBAF, CHCl₃/CH₂Cl₂, 50%; (ii) TBAF, CHCl₃/CH₂Cl₂, 44%; (iii) Pd(PPh₃)₂Cl₂, CuI, 1,4-benzoquinone, toluene/ⁱPr₂NH, 48%. Ln = Y, Gd, Tb, Dy; TIPS = triisopropylsilane; CPDIPS = cyanopropyldiisopropylsilane; Ar = 3,5-bis(tert-butyl)phenyl.

For this work we used Kläui's tripodal ligand [CpCo{P(=O) (OEt)₂]₃]⁻ (Cp = cyclopentadienyl) (Fig. 1) as a redox-stable capping group for the porphyrins.^9c,d The tridentate ligand has three coordinating oxygen atoms and a single negative charge, yielding a neutral 7-coordinate lanthanide complex. This choice of ligands is fundamental to the creation of a reliable system that can operate in the non-crystalline state and in molecular electronic devices on surfaces. The tripodal shape of the Kläui ligand effectively shields the Ln^III from the environment. Moreover it confers solubility in most organic solvents, and its redox state and binding to Co(III) are resilient to a broad range of reaction and purification conditions. All these characteristics are required if this complex were to be grafted into a devices. Moreover the ligand has an encumbered geometry that is likely to favour grafting onto surfaces via the porphyrin system, thus facilitating the presence of a single geometry in the devices.

Statistical de-protection of the porphyrin monomers (P1_Ln), followed by oxidative homo-coupling, was used for the formation of the homoleptic dimers (P2_Ln2).¹¹ We prepared homoleptic dimers containing two Gd, Tb and Dy centres. The hetero-metallated dimer containing both terbium and dysprosium centres was prepared via an alternative route. An analogue of the dysprosium monomer (P1*_Dy) was prepared containing polar CPDIPS protected acetylenes. Oxidative coupling of a mixture of mono de-protected P1_Tb and P1*_Dy, resulted in a statistical mixture of three dimers which could be separated on silica based on the different polarity of the protecting groups. Successful separation was confirmed using MALDI (Fig. 2) where the hetero-dinuclear peak is clearly present with the expected isotopic pattern, and homo-dinuclear species are not present. This indicates that no metal switching occurs during the oxidative coupling of monomer units.

Fig. 2 MALDI-TOF mass spectrum of P2_Dy-Tb, acquired with a dithranol matrix. The absence of homodinuclear species (which would fall at 3188 and 3198 amu for homo-nuclear Tb^III and Dy^III complexes) is clearly observable (top), and the agreement of the isotopic pattern with the one expected for P2_Dy-Tb is shown in the bottom.

The presence of the paramagnetic lanthanide ions renders the ¹H-NMR spectra of the complexes rather uninformative. A diamagnetic yttrium analogue (P1_Y) was prepared in order to confirm the general structure by NMR spectroscopy (ESI†).

Structural characterisation

Single-crystals of the porphyrin monomers were grown by slow liquid diffusion of methanol into a CH₂Cl₂ solution of P1_Ln. The homoleptic dimers (P2_Ln2) yielded single crystals by slow evaporation of a solution of P2_Ln2 in a mixture of ethanol/CH₂Cl₂.

Single-crystals of each P1_Y, P1_Gd, P1_Tb, P1_Dy, P2_Gd2, P2_Dy2 and P2_Tb2 were obtained and select crystallographic information is provided in Table 1. P1_Gd, P1_Tb, P1_Dy and P1_Y monomers display the same overall molecular connectivity but are not isostructural, with differing packing between the monomers resulting in the adoption of different space groups. P1_Gd crystallizes in monoclinic Cc, P1_Tb and P1_Y crystallise in monoclinic P2₁/c and P1_Dy crystallises in a triclinic P1̄ space group. Of these complexes, only the packing structure of P1_Tb contains solvent molecules: one dichloromethane molecule per Tb complex. It is surprising that these monomer complexes crystallise differently because they were all crystallised under identical conditions. Additionally, for P1_Dy there are two molecules in the asymmetric unit whereas all other monomers have one molecule per asymmetric unit. In related lanthanide structures^9e the Gd^III, Tb^III and Dy^III complexes all crystallise in a monoclinic P2₁/c space group. The discrepancy in packing arrangement is a manifestation of crystal polymorphism.¹²

Table 1 Selected crystallographic angles (°) and distances (Å)

	Space group	Nearest intramolecular Ln^III–Ln^III (Å)	Nearest intermolecular Ln^III–Ln^III (Å)	PorphyrinCentroid–Ln^III–Co^III (°)	PorphyrinCentroid–Ln^III (Å)	Ln^III–P or 1–P or 2–Ln^III torsion (°)	Co^III–Ln^III (Å)
a Data from ref. 9e: LnP LOEt^− = [(LOEt)Ln-(TPP)]·0.25H2O where TPP = 5,10,15,20-tetraphenylporphyrinate, LOEt^− = [(η^5-C5H5)Co{P(=O) (OEt)2}3]^−.
P1Y	P21/c	—	11.1407(3)	177.83(4)	1.1958(9)	—	4.2699(4)
P1Gd	Cc	—	12.10900(9)	176.97(5)	1.2360(11)	—	4.3042(4)
P1Tb	P21/c	—	9.87989(19)	174.50(13)	1.233(3)	—	4.2467(14)
P1Dy(A)	P1̄	—	9.8418(3)	170.48(9)	1.2084(19)	—	4.2457(8)
P1Dy(B)	P1̄		9.8418(3)	175.49(9)	1.212(2)		4.2410(8)
P2Gd2	P1̄	13.9202(7)	10.3106(2)	174.15(10)	1.266(2)	180.0	4.272(3)
P2Dy2	P1̄	13.8726(6)	10.25340(18)	174.63(7)	1.2239(16)	180.0	4.248(2)
P2Tb2	P1̄	13.9299(5)	10.3082(1)	174.47(9)	1.2519(18)	180.0	4.270(3)
GdP LOEt^−^a^8e	C2/c	—	10.14	177	1.23	—	4.30
TbP LOEt^−^a^8e	C2/c	—	10.11	177	1.23	—	4.29
DyP LOEt^−^a^8e	C2/c	—	10.10	177	1.22	—	4.28
HoP LOEt^−^a^8e	C2/c	—	10.14	177	1.23	—	4.30

---

For simplicity only one monomer of the Dy-analogue will be described in detail, complete crystallographic data for all the compounds are available in the ESI.† For P1_Dy the Dy^III ion is sandwiched between four N atoms of the porphyrin ligand and three O atoms of the capping ligand resulting in a 7-coordinate complex (Fig. 3a). The two aryl pendants on the porphyrin are not parallel to the plane of the porphyrin, but are tilted at an angle of 67.42(16)° and 87.74(14)°. The Dy^III ion sits 1.212(2) Å out of the plane of the porphyrin and 4.2410(8) Å from the Co^III ion. The Cp ring and porphyrin appendage are near parallel to each other with a tilt angle of 4.4(3)°. The porphyrin, Dy^III and Co^III ions also have a near-linear arrangement with a porphyrin-centroid (C_Por)–Dy–Co angle of 175.49(9)°, illustrated in Fig. 3a and b.

Fig. 3 Crystal structures of representative examples of the investigated compounds, as acquired at room temperature. (a) Side view of P1_Dy; (b) top view of P1_Dy, along the Co–Dy vector; (c) side view of P2_Dy2; (d) top view of P2_Dy2. H atoms have been omitted for clarity. Triisopropylsilylethynyl and aryl groups are faded for clarity and the colour-scale is: carbon = light grey, oxygen = red, nitrogen = dark blue, phosphorous = purple, cobalt = light blue and dysprosium = pink.

Examination of the four lanthanide monomers reveals that the C_Por–Ln distance for P1_Gd, P1_Tb, P1_Dy and P1_Y decreased roughly in proportion to the atomic radius from P1_Gd through P1_Y. The different packing arrangement may explain why we do not observe any periodic trends related to intermolecular Ln–Ln (Å) distance or linearity of the complexes measured through the C_Por–Ln–Co tilt angle (Table 1).

We next examined the solid state structures of three homoleptic lanthanide dimers. In contrast to the monomers, each of P2_Gd2, P2_Tb2 and P2_Dy2 complexes crystallise in a triclinic P1̄ space group with no solvent molecules in the crystal lattice and all structures contain a centre of inversion. The nearest intramolecular Ln–Ln distance decreases moving from left to right across the lanthanide series P1_Gd −13.9202(7) Å, P1_Tb −13.9299(5) Å and P1_Dy −13.8726(6) Å. There is a centre of inversion at the centre of each porphyrin dimer. This is reflected in the Ln–C_Por–C_Por–Ln torsion angle (180° in each dimer). We expect that the centre of inversion would not be maintained in solution or on a surface as the barrier to torsional rotation in a butadiyne-linked porphyrin dimer is quite low (ΔH = 5.27 kJ mol⁻¹).¹³ There is very little deviation of the coordination environment surrounding the lanthanide ions in the monomer and dimer structures. SHAPE¹⁴ software was used to compare the polyhedral of P1_Dy, P2_Dy2, P1_Tb, P2_Tb2, P1_Gd, P2_Gd2 with other possible 7-vertex polyhedral. All 6 complexes exhibited predominately either capped octahedron or a capped trigonal prism (see full discussion in the ESI†).

The main objective of creating the P2_Dy-Tb is to obtain hetero-dinuclear systems where every molecule in the crystal contains one Dy^III centre and one Tb^III centre. We demonstrated the purity of this hetero-dimer using Matrix Assisted Laser Desorption Ionization-Time of Flight (MALDI-TOF) mass spectrometry, which is an excellent way of identifying perfectly hetero-dinuclear compounds, as they provide a starkly different mass value and isotopic pattern for hetero and homo-dinuclear molecules (Fig. 2).

Magnetic properties

Each of P1_Gd, P1_Tb, P1_Dy, P2_Gd2, P2_Dy2, P2_Tb2 and P2_Dy-Tb have an overall charge of zero. The Co^III ion is diamagnetic and the magnetic properties discussed below originate solely from the localized spin of the Ln^III ion(s). We shall investigate both the static and dynamic properties of the systems, relevant to their suitability as spin-valve components.

Static magnetic properties. The variable temperature magnetic properties of all complexes were determined using a MPMS-XL SQUID magnetometer. Fig. 4 provides the static magnetic properties for each crystalline complex, with the magnetic susceptibility χ calculated as the ratio between the magnetisation M and the applied magnetic field H. The χT products for the P1_Gd, P1_Tb and P1_Dy monomers are 8.40, 11.19 and 14.16 cm³ kmol⁻¹ respectively at room temperature. These values are in good agreement with the theoretical values of 7.88, 11.82 and 14.17 cm³ kmol⁻¹ for a single Gd³⁺ ion (⁸S_7/2, S = 7/2, L = 0, g = 2), a Tb³⁺ ion (⁷F₆, S = 3, L = 3, g = 3/2) and a Dy³⁺ ion (⁶H_15/2, S = 5/2, L = 5, g = 4/3), respectively. The overall shape of each variable temperature curve provides further information about the magnetic properties of each compound. The ⁸S_7/2 ground state of Gd^III has zero orbital angular momentum, and the straight line is consistent with its isotropic nature. Conversely, the deviation from linearity observed for dysprosium at low temperatures is to be expected as a result of the large magnetic anisotropy of Dy^III, with its 4f⁹ configuration and a ⁶H_15/2 ground state. Interestingly, at low temperature Tb^III demonstrates a slight increase in the χT product below 10 K. This probably indicates an intermolecular ferromagnetic interaction between neighbouring spin carriers (closest Tb–Tb intermolecular distance 9.9 Å).

Fig. 4 Temperature dependence of the static magnetic susceptibility times the temperature, for all lanthanide-porphyrin complexes described in the text. Spheres indicate dinuclear compounds, open triangles mononuclear ones; Gd compounds are in pink, Dy in green, Tb in blue and the hetero-dinuclear one in black. All measurements were performed on crystalline powder samples under an external static magnetic field of 1 kOe.

The room temperature χT values for each P2_Gd2, P2_Tb2, P2_Dy2 and P2_Dy-Tb dimers are 14.50, 22.69, 28.56 and 25.22 cm³ kmol⁻¹ and are in good agreement with twice the theoretical values of monomers 15.76, 23.64, 28.34 and 25.99 cm³ kmol⁻¹ respectively. For P2_Gd2 the χT product remains basically unchanged from room temperature to 5 K, with only a very slight decrease below 5 K. This agrees well with an S = 7/2 centre with no anisotropy and with very small anti-ferromagnetic interactions. For P2_Tb2 and P2_Dy2 the χT product decreases steadily from 300 K down to 5 K, only to eventually increase again below 5 K. In this case, the marked decrease below 50 K is attributed primarily to the large inherent magnetic anisotropy of Dy^III and Tb^III ions, while the slight low temperature increase may indicate weak ferromagnetic coupling between intramolecular Ln^III ions. Such magnetic interactions can be dipolar or exchange-coupling in nature. The distances between Ln ions are presented in Table 1.

Incoherent dynamic magnetic properties. AC susceptibility studies of all Dy^III and Tb^III complexes revealed no out-of-phase zero field magnetic susceptibility in the χ′′ plot at 1.8 K, presumably due to the presence of ground state quantum tunnelling of the magnetisation (QTM). Under a dc field of 2000 Oe, full frequency and temperature dependent peaks are observed in the χ′′ susceptibility for all Dy^III and Tb^III complexes.

The two monomeric systems, P1_Tb and P1_Dy, show very similar slow magnetic relaxation dynamics, with frequency dependent peaks that shift to higher frequencies on increasing T (Fig. 5a and c). Using the χ′′ peak maxima determined using a Lorentz fit, and the Arrhenius law (τ = τ₀ exp(U_eff/k_BT) where τ = 1/(2πυ) and υ is the frequency corresponding to the maxima of each χ′′ peak), the effective energy barriers were obtained; P1_Tb, U_eff = 7(2) K with a τ₀ = 3.4 × 10⁻⁵ s and P1_Dy, U_eff = 14(4) K with a τ₀ = 2.5 × 10⁻⁶ s. The relaxation barrier and τ₀ for P1_Dy is very close to that reported for the similar 10,15,20-tetraphenylporphyrin dysprosium complex.^9e

Fig. 5 Dynamic magnetic susceptibility of crystalline samples of P1_Tb (a), P2_Tb2 (b), P1_Dy (c), P2_Dy2 (d) and P2_Dy-Tb (e) for different temperatures (colour scale). All measurements were performed under a 2000 Oe constant magnetic field, in a 1–1500 Hz frequency range. (f) A plot of ln(τ) versus 1/T for all compounds, showing the Arrhenius behaviour. Data points are represented with error bars obtained from the fitting of the peak maxima of the frequency dependence in a–e with Lorentzian curves. Solid lines represent linear regression to the data.

Homoleptic dimers P2_Tb2 and P2_Dy2 displayed slightly different magnetic properties to their corresponding monomer complexes: P2_Tb2, U_eff = 15(2) K with τ₀ = 3.6 × 10⁻⁵ s; P2_Dy2, U_eff = 7(2) K with a τ₀ = 4.5 × 10⁻⁵ s. The slight differences in the U_eff between mono- and bimetallic complexes is most likely a result of the slight change in the second coordination sphere where the TIPS pendent is replaced by an acetylene resulting in a less symmetric second coordination sphere. Changes to the second-coordination sphere can often influence the magnetic properties of lanthanide ions.¹⁵P2_Dy-Tb also showed similar behaviour to the above with U_eff = 9(2) K with τ₀ = 2.9 × 10⁻⁵ s confirming slow-magnetic relaxation at low temperatures.

Coherent dynamic magnetic properties. We accessed the quantum coherence properties of the systems by pulsed electron paramagnetic resonance (EPR) techniques. These can provide extremely valuable information on whether the compounds can be used for quantum operation processing via single-electron transport at low temperatures, which for example includes T₁ and T₂ times. By applying different sequences of microwave pulses, we can extract the spin-lattice relaxation time, T₁, and the spin–spin dephasing (or phase-memory) time, T_M. Furthermore, pulsed EPR techniques can provide valuable information on very low spin–spin dipolar and exchange interactions in dimer systems, as is valuable for the molecular spin-valves. We restrict our analysis to the compounds based on Gd^III, P1_Gd and P2_Gd2, because of the extreme broadening and zero-field splitting of Dy^III and Tb^III centres due to spin–orbit interactions with their environment. As the orbital momentum is zero for the ground state, Gd^III systems can be treated as pure spin systems with a total spin of S = 7/2. Interactions and mixing with higher excited multiplets and higher-order spin–orbit interactions will be neglected here. The low symmetry crystal-fields of P1_Gd and P2_Gd2 result in a complicated crystal-field splitting with four Kramers doublet eigenstates, where states with different magnetic quantum numbers m_s = ±1/2, ±3/2, ±5/2, ±7/2 are mixed. Fig. 6 shows the spectra obtained by the derivative of the free-induction-decay (FID) detected absorption at X-band. The broad linewidths, hundreds of Gauss, are typical of Gd^III EPR spectra.^6h Comparison of P1_Gd and P2_Gd2 reveals a very similar spectrum, with the additional Gd^III ion in the P2_Gd2 dimer producing only minimal changes around 2000 and 3300 Oe. The small changes in the spectrum can be due to weak dipolar and exchange interactions between the Gd^III ions, where the latter one is mediated by the alkyne bridge. Due to the broad linewidth, the interactions cannot be resolved by standard CW spectroscopy at X-band frequency. For improved resolution and to determine coherence times, we used pulsed EPR techniques.

Fig. 6 EPR spectra comparing the derivative of the free-induction-decay-detected absorption of P1_Gd (red line) and P2_Gd2 (blue line). The spectra were acquired at X-band frequency (9.4 GHz) at 3 K in a 1 : 1 mixture of dichloromethane and toluene.

The phase-memory times, T_M, were determined using a variation of the Hahn-echo sequence (Fig. 7, ESI for details†). T_M is found to increase steadily on lowering the temperature, with an almost linear trend from 10 to 2 K, reaching 3.6 μs and 3.0 μs for P1_Gd and P2_Gd2, respectively. This indicates only a marginal decoherence time introduced by the dimerisation, and provides times that amply meet the requirements for possible quantum operations. The echo signal decays mono-exponentially versus the free evolution time τ of the Hahn echo (ESI†), indicating no spectral diffusion. Analogously, no detectable dependence of the T₂ times on the magnetic field is observed (ESI†). This means that the dinuclear system will indeed appear composed of two distinct magnetic centres in low temperature transport experiments. Measurements of the single molecule transport and magnetic properties are currently underway.

Fig. 7 Temperature dependence of the quantum coherence time T_M for P1_Gd (red full dots) and P2_Gd2 (blue empty dots), as measured at 344 mT at 9.37 GHz in a 1 : 1 mixture of dichloromethane and toluene.

Conclusions

In conclusion, we have created hetero-dinuclear single molecule magnets based on porphyrin ligands, using a rational stepwise approach. We demonstrated the use of protecting groups of very different polarity for the selective creation of hetero-metallic porphyrin complexes. The protecting group methodology can now be used as a tool for the tuning of functional poly-porphyrin molecules, providing a route to hetero-dinuclear molecular magnets.^7,16 From this perspective, the dimers here created are the first that are sufficiently stable to undergo purification steps such as chromatographic separations.

An extension of these novel synthetic approach would be to rationally assemble poly-heteronuclear lanthanide chains. For example, tetra- and hexa-heteronuclear chains with a pre-defined order of the magnetic centres are planned via stepwise repetition of the oxidative homo-coupling step. Including the optical properties of the porphyrin ligands themselves, future work will include addition of functional groups such as luminescent appendages, light antennae or anchoring groups.

The magnetic properties of the compounds created show the desired slow magnetic relaxation and the coherent properties indicate μs coherence times at low temperatures for Gd^III compounds. These compounds offer a testbed to probe magnetic interactions in molecular systems, and indicate a possible applicative direction for ring-shaped poly-porphyrins¹⁷ in spintronics and in single-molecule EPR investigations.¹⁸ All stringent requirements for the creation of spin-valves are simultaneously met and the new molecules show slow relaxation, magnetic anisotropy, the presence of molecular electronic quantum dots, different spin centres and stability that permits deposition on different surfaces.

Conflicts of interest

There are no conflicts to declare.

Acknowledgements

We thank Dr F. Tuna and Prof. R. Winpenny for SQUID time at Manchester University during a downtime in Oxford. We acknowledge financial support by the European Union (ERC-StG-338258 OptoQMol, ERC-AdG-320969) and Marie Skłodowska-Curie Research Fellowship MC-SpinReMag-707252, the Violette & Samuel Glasstone Research Fellowship, the Royal Society (URF and Grant funds), EPSRC (QuEEN grant), EPSRC (Centre for Advanced ESR spectroscopy (CAESR) at the University of Oxford (EP/L011972/1)), German DV scholarship and Corpus Christi College. We thank the UK National Mass Spectrometry Facility at Swansea University for recording mass spectra.

Notes and references

1. (a) T. Tanaka and A. Osuka, Chem. Soc. Rev., 2015, 44, 943–969; (b) D. Holten, D. F. Bocian and J. S. Lindsey, Acc. Chem. Res. 2002, 35, 1, 57; (c) C. Maeda, T. Taniguchi, K. Ogawa and T. Ema, Angew. Chem., Int. Ed., 2015, 54, 134–138; (d) G. Sedghi, V. M. García-Suárez, L. J. Esdaile, H. L. Anderson, C. J. Lambert, S. Martin, D. Bethell, S. J. Higgins, M. Elliott, N. Bennett, J. E. Macdonald and R. J. Nichols, Nat. Nanotechnol., 2011, 6, 517; (e) M. D. Peeks, C. E. Tait, P. Neuhaus, G. M. Fischer, M. Hoffmann, R. Haver, A. Cnossen, J. R. Harmer, C. R. Timmel and H. L. Anderson, J. Am. Chem. Soc., 2017, 139, 10461.
2. (a) A. Harriman, F. Odobel and J. P. Sauvage, J. Am. Chem. Soc., 1995, 117, 9461; (b) J. Fortage, J. Boixel, E. Blart, L. Hammarstrçm, H. C. Becker and F. Odobel, Chem.–Eur. J., 2008, 14, 3467; (c) J. Cremers, S. Richert, D. Kondratuk, T. D. W. Claridge, C. Timmel and H. L. Anderson, Chem. Sci., 2016, 7, 6961.
3. (a) S. A. Wolf, D. D. Awschalom, R. A. Buhrman, J. M. Daughton, S. von Molnár, M. L. Roukes, A. Y. Chtchelkanova and D. M. Treger, Science, 2001, 294(5546), 1488; (b) M. Slota, A. Keerthi, W. K. Myers, E. Tretyakov, M. Baumgarten, A. Ardavan, H. Sadeghi, C. J. Lambert, A. Narita, K. Müllen and L. Bogani, Nature, 2018, 557, 691.
4. (a) X. Sun, S. Vélez, A. Atxabal, A. Bedoya-Pinto, S. Parui, X. Zhu, R. Llopis, F. Casanova and L. E. Hueso, Science, 2017, 357(6352), 677; (b) S. Sanvito, Chem. Soc. Rev., 2011, 40, 3336; (c) C. Cervetti, A. Rettori, M. G. Pini, A. Cornia, A. Repollés, F. Luis, M. Dressel, S. Rauschenbach, K. Kern, M. Burghard and L. Bogani, Nat. Mater., 2016, 15, 164; (d) E. Ö. Karabulut, F. Malet, A. L. Fetter, G. M. Kavoulakis and S. M. Reimann, New J. Phys., 2016, 18, 015013; (e) D. D. Awschalom, L. C. Bassett, A. S. Dzurak, E. L. Hu and J. R. Petta, Science, 2013, 339, 1174; (f) L. Bogani and W. Wernsdorfer, Nat. Mater., 2008, 7, 179; (g) S. D. Jiang, K. Goß, C. Cervetti and L. Bogani, Sci. China: Chem., 2012, 55(6), 867; (h) L. Bogani, Jpn. J. Appl. Phys., 2011, 109, 07B115; (i) L. Bogani, R. Maurand, L. Marty, C. Sangregorio, C. Altavilla and W. Wernsdorfer, J. Mater. Chem., 2010, 20, 2099–2107; (j) C. Cervetti, E. Heintze and L. Bogani, Dalton Trans., 2014, 43, 4220–4232; (k) E. Heintze, F. El Hallak, C. Clauß, A. Rettori, M. G. Pini, F. Totti, M. Dressel and L. Bogani, Nat. Mater., 2013, 12, 202.
5. (a) S. Thiele, F. Balestro, R. Ballou, S. Klyatskaya, M. Ruben and W. Wernsdorfer, Science, 2014, 344, 1135; (b) F. Troiani, C. Godfrin, S. Thiele, F. Balestro, W. Wernsdorfer, S. Klyatskaya, M. Ruben and M. Affronte, Phys. Rev. Lett., 2017, 118, 257701; (c) A. Candini, S. Klyatskaya, M. Ruben, W. Wernsdorfer and M. Affronte, Nano Lett., 2011, 11, 2634; (d) M. Urdampilleta, S. Klayatskaya, M. Ruben and W. Wernsdorfer, ACS Nano, 2015, 9(4), 4458; (e) M. Ganzhorn, S. Klyatskaya, M. Ruben and W. Wernsdorfer, Nat. Nanotechnol., 2013, 8, 165; (f) M. Misiorny, E. Burzurì, R. Gaudenzi, K. Park, M. Leijnse, M. R. Wegewijs, J. Paaske, A. Cornia and H. S. J. van der Zant, Phys. Rev. B: Condens. Matter Mater. Phys., 2015, 91, 035442; (g) E. Burzurì, Y. Yamamoto, M. Warnock, X. Zhong, K. Park, A. Cornia and H. S. J. van der Zant, Nano Lett., 2014, 14, 3191; (h) R. Frisenda, G. D. Harzmann, J. A. Celis Gil, J. M. Thijssen, M. Mayor and H. S. J. van der Zant, Nano Lett., 2016, 16, 4733; (i) J. Ponce, C. R. Arroyo, S. Tatay, R. Frisenda, P. Gaviña, D. Aravena, E. Ruiz, H. S. J. van der Zant and E. Coronado, J. Am. Chem. Soc., 2014, 136, 8314; (j) D. Schmid-Lorch, T. Häberle, F. Reinhard, A. Zappe, M. Slota, L. Bogani, A. Finkler and J. Wrachtrup, Nano Lett., 2015, 15, 4942–4947.
6. (a) G. Christou, D. Gatteschi, D. N. Hendrickson and R. Sessoli, MRS Bull., 2000, 66; (b) R. Sessoli, D. Gatteschi, A. Caneschi and M. A. Novak, Nature, 1993, 365, 141; (c) J. J. Le Roy, L. Ungur, I. Korobkov, L. F. Chibotaru and M. Murugesu, J. Am. Chem. Soc., 2014, 136, 8003; (d) S. Chen, V. Mereacre, C. E. Anson and A. K. Powell, Dalton Trans., 2016, 45, 98; (e) M. Gysler, F. El Hallak, L. Ungur, R. Marx, M. Hakl, P. Neugebauer, Y. Rechkemmer, Y. Lan, I. Sheikin, M. Orlita, C. E. Anson, A. K. Powell, R. Sessoli, L. F. Chibotaru and J. van Slageren, Chem. Sci., 2016, 7, 4347; (f) T. Lacelle, G. Brunet, R. J. Holmberg, B. Gabidullin and M. Murugesu, Cryst. Growth Des., 2017, 17(10), 5044; (g) J. J. Baldoví, Y. Duan, C. Bustos, S. Cardona-Serra, P. Gouzerh, R. Villanneau, G. Gontard, J. M. Clemente-Juan, A. Gaita-Ariño, C. Giménez-Saiz, A. Proust and E. Coronado, Dalton Trans., 2016, 45, 16653; (h) M. D. Jenkins, Y. Duan, B. Diosdado, J. J. García-Ripoll, A. Gaita-Ariño, C. Giménez-Saiz, P. J. Alonso, E. Coronado and F. Luis, Phys. Rev. B, 2017, 95, 064423; (i) J. D. Hilgar, B. S. Flores and J. D. Rinehart, Chem. Commun., 2017, 53, 7322; (j) J. Xiong, H.-Y. Ding, Y.-S. Meng, C. Gao, X.-J. Zhang, Z.-S. Meng, Y.-Q. Zhang, W. Shi, B.-W. Wang and S. Gao, Chem. Sci., 2017, 8, 1288; (k) C. A. Gould, L. E. Darago, M. I. Gonzalez, S. Demir and J. R. Long, Angew. Chem., Int. Ed., 2017, 56, 10103; (l) M. Ding, G. E. Cutsail, D. Aravena, M. Amoza, M. Rouzières, P. Dechambenoit, Y. Losovyj, M. Pink, E. Ruiz, R. Clérac and J. M. Smith, Chem. Sci., 2016, 7, 6132; (m) E. Kiefl, M. Mannini, K. Bernot, X. Yi, A. Amato, T. Leviant, A. Magnani, T. Prokscha, A. Suter, R. Sessoli and Z. Salman, ACS Nano, 2016, 10, 5663; (n) F. Pointillart, K. Bernot, B. Le Guennic, S. Golhen, O. Cador and L. Ouahab, Angew. Chem., Int. Ed., 2015, 54, 1504; (o) J.-K. Ou-Yang, N. Saleh, G. Fernandez Garcia, L. Norel, F. Pointillart, T. Guizouarn, O. Cador, F. Totti, L. Ouahab, J. Crassous and B. Le Guennic, Chem. Commun., 2016, 52, 14474; (p) L. Bogani, Jpn. J. Appl. Phys., 2011, 109, 07B115.
7. 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.
8. (a) K. Ueta, K. Naoda, S. Ooi, T. Tanaka and A. Osuka, Angew. Chem., Int. Ed., 2017, 56, 7223; (b) C. Stähler, D. Shimizu, K. Yoshida, K. Furukawa, R. Herges and A. Osuka, Chem.–Eur. J., 2017, 23, 7217; (c) R. Wang, A. M. Brugh, J. Rawson, M. J. Therien and M. D. E. Forbes, J. Am. Chem. Soc., 2017, 139(29), 9759; (d) Y. Bai, J.-H. Olivier, H. Yoo, N. F. Polizzi, J. Park, J. Rawson and M. J. Therien, J. Am. Chem. Soc., 2017, 139(46), 16946; (e) B. Basic, O. B. Locos, L. Rintoul, J. Kan, J. Jiang and D. P. Arnold, ChemPlusChem, 2014, 79, 813; (f) W.-K. Wong, L. Zhang, F. Xue and T. C. W. Mak, J. Chem. Soc., Dalton Trans., 1999, 0, 3053; (g) C.-M. Che and S.-Y. Wong, US pat., US20100069345A1, 2005; (h) Y. V. Korovina, N. V. Rusakova, Z. I. Zhilina, Y. V. Ishkovb, S. V. Vodzinsky and V. P. Dotsenko, Mendeleev Commun., 2002, 12(4), 151.
9. (a) J. H. Shi, Z. Li, W. Y. Liu, M. Yu, Z. H. Li and G. F. Liu, Inorg. Nano-Met. Chem., 2012, 43, 316; (b) A. G. Coutsolelos and G. A. Spyroulias, Polyhedron, 1994, 13, 647; (c) W. Kläui, H. Neukomm, H. Werner and G. Huttner, Chem. Ber., 1977, 110, 2283; (d) W. Kläui, Angew. Chem., Int. Ed., 1990, 29, 627; (e) F. Gao, M.-X. Yao, Y.-Y. Li, Y.-Z. Li, Y. Song and J.-L. Zuo, Inorg. Chem., 2013, 52, 6407.
10. (a) Y. Lan, S. Klyatskaya, M. Ruben, O. Fuhr, W. Wernsdorfer, A. Candini, V. Corradini, A. L. Rizzini, U. del Pennino, F. Troiani, L. Joly, D. Klar, H. Wendeh and M. Affronte, J. Mater. Chem. C, 2015, 3, 9794; (b) D. J. Lewis, P. B. Glover, M. C. Solomons and Z. Pikramenou, J. Am. Chem. Soc., 2011, 133, 1033; (c) R. Sato, K. Suzuki, M. Sugawa and N. Mizuno, Chem.–Eur. J., 2013, 19, 12982; (d) J. Gónzalez-Fabra, N. A. G. Bandeira, V. Velasco, L. A. Barrios, D. Aguilà, S. J. Teat, O. Roubeau, C. Bo and G. Aromí, Chem.–Eur. J., 2017, 23, 5117; (e) S. Haas, E. Heintze, S. Zapf, B. Gorshunov, M. Dressel and L. Bogani, Phys. Rev. B: Condens. Matter Mater. Phys., 2014, 89, 174409.
11. P. N. Taylor, J. Huuskonen, G. Rumbles, R. T. Aplin, E. Williams and H. L. Anderson, Chem. Commun., 1998, 909.
12. (a) N. Meundaeng, A. Rujiwatra and T. J. Prior, Transition Met. Chem., 2016, 41(7), 783; (b) R. J. Davey, Cryst. Growth Des., 2002, 2(6), 675; (c) L. Borka and J. K. Haleblain, Acta Pharm. Jugosl., 1990, 40, 71; (d) H. G. Brittain, Polymorphism in Pharmaceutical Solids, CRC Press, 2016.
13. M. D. Peeks, P. Neuhaus and H. L. Anderson, Phys. Chem. Chem. Phys., 2016, 18, 5264.
14. M. Llunell, D. Casanova, J. Girera, P. Alemany and S. Alvarez, SHAPE, Continuous Shape Measures Calculation, version 2.0, Universitat de Barcelona, Barcelona, Spain, 2010.
15. (a) J. J. Le Roy, M. Jeletic, S. I. Gorelsky, I. Korobkov, L. Ungur, L. F. Chibotaru and M. Murugesu, J. Am. Chem. Soc., 2013, 135, 3502; (b) G. Cucinotta, M. Perfetti, J. Luzon, M. Etienne, P.-E. Car, A. Caneschi, G. Calvez, K. Bernot and R. Sessoli, Angew. Chem., Int. Ed., 2012, 51, 1606; (c) R. Liu, C. Zhang, L. Li, D. Liao and J.-P. Sutter, Dalton Trans., 2012, 41, 12139; (d) M. Dörfel, M. Kern, H. Bamberger, P. Neugebauer, K. Bader, R. Marx, A. Cornia, T. Mitra, A. Müller, M. Dressel, L. Bogani and J. van Slageren, Magnetochemistry, 2016, 2, 25; (e) K. Bernot, J. Luzon, L. Bogani, M. Etienne, C. Sangregorio, M. Shanmugam, A. Caneschi, R. Sessoli and D. Gattesch, J. Am. Chem. Soc., 2009, 131, 5573–5579.
16. (a) N. Ishikawa, T. Iino and Y. Kaizu, J. Am. Chem. Soc., 2002, 124, 11440; (b) J. P. Costes, F. Dahan and F. Nicodéme, Inorg. Chem., 2003, 42, 6556; (c) X.-Y. Chen, Y. Bretonniére, J. Pécaut, D. Imbert, J.-C. Bünzli and M. Mazzanti, Inorg. Chem., 2007, 46, 625.
17. S. Richert, I. Kuprov, M. D. Peeks, E. A. Suturina, J. Cremers, H. L. Anderson and C. R. Timmel, Phys. Chem. Chem. Phys., 2017, 19, 16057.
18. S. Baumann, W. Paul, T. Choi, C. P. Lutz, A. Ardavan and A. J. Heinrich, Science, 2015, 350(6259), 417.

---

Footnotes

† Electronic supplementary information (ESI) available: Synthesis and characterization of new compounds. CCDC 1565176–1565188. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/c8sc03762k

‡ These authors contributed equally to this work.

---