Single-molecule magnet behavior in 2,2’-bipyrimidine-bridged dilanthanide complexes

  1. Wen Yu1,
  2. Frank Schramm1,
  3. Eufemio Moreno Pineda1,
  4. Yanhua Lan2,
  5. Olaf Fuhr1,
  6. Jinjie Chen3,
  7. Hironari Isshiki3,
  8. Wolfgang Wernsdorfer2,
  9. Wulf Wulfhekel3 and
  10. Mario Ruben1,4

1Institute of Nanotechnology (INT), Karlsruhe Institute of Technology (KIT), Hermann-von-Helmholtz-Platz 1, D-76344 Eggenstein-Leopoldshafen, Germany
2Université Grenoble Alpes / CNRS, Institut NEEL, 25 rue des Martyrs, F-38000 Grenoble, France
3Physikalisches Institut (PI), Karlsruhe Institute of Technology (KIT), Wolfgang-Gaede-Str. 1, D-76131 Karlsruhe, Germany,

4Institut de Physique et Chimie des Matériaux de Strasbourg (IPCMS), CNRS-Université de Strasbourg, 23 rue du Loess, BP 43, F-67034 Strasbourg Cedex 2, France

  1. Corresponding author email

This article is part of the Thematic Series "Functional nanostructures – structure formation and self-organization, electrochemistry".

Guest Editor: R. J. Behm
Beilstein J. Nanotechnol. 2016, 7, 126–137. doi:10.3762/bjnano.7.15
Received 07 Oct 2015, Accepted 08 Jan 2016, Published 28 Jan 2016


A series of 2,2’-bipyrimidine-bridged dinuclear lanthanide complexes with the general formula [Ln(tmhd)3]2bpm (tmhd = 2,2,6,6-tetramethyl-3,5-heptanedionate, bpm = 2,2’-bipyrimidine, Ln = Gd(III), 1; Tb(III), 2; Dy(III), 3; Ho(III), 4 and Er(III), 5) has been synthesized and characterized. Sublimation of [Tb(tmhd)3]2bpm onto a Au(111) surface leads to the formation of a homogeneous film with hexagonal pattern, which was studied by scanning tunneling microscopy (STM). The bulk magnetic properties of all complexes have been studied comprehensively. The dynamic magnetic behavior of the Dy(III) and Er(III) compounds clearly exhibits single molecule magnet (SMM) characteristics with an energy barrier of 97 and 25 K, respectively. Moreover, micro-SQUID measurements on single crystals confirm their SMM behavior with the presence of hysteresis loops.

Keywords: bipyrimidine; diketone; hysteresis; lanthanide; magnetism


The application of magnetic molecular compounds within molecular electronic devices is combined in the progressive field of spintronics. An anisotropic spin is confined by the ligand field of the metal complex, which is typically of a few nm in dimensions. It is proposed that the ultimate size limits of modern electronics and information processing can be tackled by such device geometry [1-3].

Research on compounds with single-molecule magnet (SMM) characteristics discovered striking advantages that make this class of molecules promising candidates for molecular spintronics [4]: (a) Each molecule behaves as a single spin domain, thus they could act as the base of a high-density information storage or processing unit; (b) They possess an energy barrier to inversion of total spin, leading to slow magnetic relaxation and magnetic hysteresis at low temperatures. Combined with their long coherence times they could open the door to quantum computing [5,6].

After the first SMM was discovered in 1980 [7,8], for the next 15 years the SMM field was dominated by cluster compounds containing mainly high-spin Mn(III). In such compounds, the SMM behavior was due to a combination of the Jahn–Teller distortion of the Mn(III), the ferromagnetic alignment of the Mn(IV), and overall antiferromagnetic exchange between the Mn(III) and Mn(IV) leading to a S = 10 ground state [9]. To date, the largest anisotropy barrier observed in a transition metal SMM is 226 cm−1 in the Fe(II) complex [Fe(C(SiMe3)3)2] [10]. Recently, molecular compounds employing lanthanide ions led to a considerable increase in strength of molecular magnets due to the observation of the SMM character at the single-ion level [11,12].

SMMs based on 4f ions possess larger thermal energy barriers for magnetization reversal caused by their large single-ion magnetic anisotropy, which originates from spin–orbit coupling and crystal-field splitting of the 4f ions. Technological and structural development of lanthanide SMM compounds to access specific surface deposition drive the innovation towards device applications [13-16].

The design of a suitable system that includes quantum bits (qubits) and quantum gates (qugates) is the main challenge to realize quantum computing. There, the electronic spins of magnetically anisotropic lanthanide ions can possibly act as basic units of quantum computing, i.e., as qubits. Universal qugates may be engineered by designing one molecule with two interacting lanthanide ions [17]. We are particularly interested in generating dinuclear lanthanide complexes with a bridging ligand that provides a communication pathway between the lanthanide ions. The 2,2’-bipyrimidine (bpm) ligand has been selected for many transition metal complexes to bridge metal ions and it acts as a connector to influence emission and magnetic properties [18-21]. An adsorption-site-dependent zero-bias (Kondo) resonance was clearly observed in the dinuclear transition metal complexes with the bpm bridging ligand [20]. In the same molecular architecture, the substitution of Ising-like spins, such as lanthanides, can offer a weak exchange coupling to study their possible implementation as qugates. Dinuclear lanthanide complexes with bpm as the bridge were reported before, but only a few studies describe their magnetic properties [22-30]. In this work we report the synthesis, characterization and single-crystal structure determination of five examples of homo-dinuclear complexes of tris-β-diketonate adducts of Gd(III), Tb(III) [30], Dy(III), Ho(III) and Er(III) with a bpm bridging ligand. We employed 2,2’,6,6’-tetramethyl-2,4-heptanedionate (tmhd) ligands as peripheric ligands, providing overall charge-neutral compounds. The magnetic behavior of the compounds was measured using AC, DC and micro-SQUID magnetometry techniques. The homo-dinuclear complexes of Dy(III) and Er(III) show single-molecule magnet behavior featuring hysteresis loops. The [Tb(tmhd)3]2bpm was sublimated on Au(111) surfaces and scanning tunneling microscopy results are presented in this report.



Solvents and reagents were of commercial grade and used without further purification. 15 were prepared by modification of published procedures [29]. To a mixture of Ln(tmhd)3 [Ln = Gd(III) 1; Tb(III) 2; Dy(III) 3; Ho(III) 4 and Er(III) 5] (1 mmol) and bpm (0.5 mmol) in a 50 mL flask absolute ethanol (20 mL) was added. The suspension was stirred overnight at room temperature, leading to a precipitate that was then collected by filtration. The powder was washed with 10 mL of cold absolute ethanol and dried at 90 °C overnight. The powder was dissolved in a minimum volume of CHCl3 at room temperature, immediately filtered to remove insoluble particles, and then reprecipitated by the addition of cold absolute ethanol. The precipitate was again filtered and dried at 90 °C. X-ray diffraction quality crystals were formed by recrystallization from Et2O for 2 and 3 at −17 °C and by layering EtOH onto a CHCl3 solution of 1, 4 and 5 at room temperature, respectively. Results of elemental analyses and isolated yields are given in Table 1.

Table 1: Elemental analysis and yield (%) for compounds 15.

  Formula Yielda Elemental analysis: found (calculated)
      C H N
1 [Gd(thmd)3]2bpm 33% 56.33 (56.53) 8.05 (7.69) 3.57 (3.56)
2 [Tb(thmd)3]2bpm·Et2O 42% 59.08 (59.06) 8.21 (8.04) 3.61 (3.72)
3 [Dy(thmd)3]2bpm 42% 56.22 (56.15) 8.08 (7.64) 3.57 (3.54)
4 [Ho(thmd)3]2bpm 43% 56.23 (55.98) 8.23 (7.62) 3.55 (3.53)
5 [Er(thmd)3]2bpm 43% 55.71 (55.82) 7.60 (7.60) 3.52 (3.52)

aCalculated based on the lanthanide starting material.

Physical measurements and instrumentation

IR transmission measurements of pressed KBr pellets were recorded at room temperature with a Perkin Elmer Spectrum GX FT-IR system spectrophotometer. Elemental analysis data were collected on an ELEMENTAR Vario Micro Cube. NMR spectra were carried out on a Bruker Ultrashield 500 PLUS spectrometer.

Scanning tunneling microscopy. The STM measurements were realized with a homebuilt scanning tunneling microscope. The sample preparation, other than molecule deposition and Ar sputtering, were carried out in ultrahigh vacuum conditions (≈10−10 mbar). Degassing of the [Tb(thmd)3]2bpm compound was performed carefully by heating to 373 K in a ceramic crucible for hours prior to evaporation. The Au(111) single crystal substrate was cleaned with a standard Ar sputtering and annealing process in a separate preparation chamber. After annealing and cooling down to room temperature, the substrate was transferred to a molecule deposition chamber and was exposed to a molecule flow of about 0.01 monolayer/sec for several ten second steps at a sublimation temperature of about 433 K. The pressure during deposition was ca. 2 × 10−7 mbar. After deposition, the samples were transferred to the STM chamber immediately and cooled down to 5 K. During the measurement, the sample temperature was kept at 5 K.

Magnetic measurements. Magnetic susceptibility measurements were collected using a Quantum Design MPMS®3 and MPMS-XL SQUID magnetometer. DC susceptibility measurements for all compounds were performed at temperatures ranging from 2 to 300 K, using an applied field of 1 kOe. The AC data were collected using an oscillating magnetic field of 3.5 Oe. All data were corrected for diamagnetic contributions from the eicosane and core diamagnetism, estimated using Pascal’s constants [31]. Low temperature (0.03–5 K) magnetization measurements were performed on single crystals using a micro-SQUID apparatus at different sweep rates between 0.280 and 0.002 T s−1 [32]. The applied field was parallel to the experimentally observed easy axis of magnetization.

X-ray data collection and structure solution. Data collection for all complexes was carried out on a STOE StadiVari 25 diffractometer with a Pilatus300 K detector using GeniX 3D HF micro focus X-ray source (Mo Kα, λ = 0.71073 Å). The structures were solved by direct methods and refined employing full-matrix least-squares refinement against F2 using SHELX2014 [33] and OLEX2 [34] packages. H atoms were added at idealized positions on their respective parent atoms. Full crystallographic details can be found in CIF format (see the Cambridge Crystallographic Data Centre database, 1427566, 1419831–1419834). These data can be obtained free of charge via (or from Cambridge Crystallographic Data Centre, 12 Union Road, Cambridge CB21EZ, UK; fax: (+44)1223-336-033; or (See Table 2).

Table 2: Crystallographic information for clusters 15.

  1 2 3 4 5
Formula C75H121Cl3Gd2N4O12 C78H130N4O13Tb2 C78H130Dy2N4O13 C75H121Cl3Ho2N4O12 C74H122Er2N4O12
fw/g mol−1 1691.60 1649.69 1656.85 1706.96 1594.27
Crystal system monoclinic triclinic triclinic monoclinic triclinic
Space group P21 [Graphic 1] [Graphic 2] P21 [Graphic 3]
a 10.8114(3) 10.8604(5) 10.8535(4) 10.8520(5) 10.8332(5)
b 29.6046(10) 13.8827(6) 13.8502(6) 29.5505(10) 12.5416(6)
c 13.8460(4) 14.8544(6) 14.8305(6) 13.7466(6) 16.8646(8)
α/° 90 92.101(3) 92.011(3) 90 106.948(4)
β/° 106.850(2) 93.253(3) 93.342(3) 106.945(4) 101.337(4)
γ/° 90 107.539(3) 107.530(3) 90 105.671(4)
V/Å3 4241.4(2) 2128.66(16) 2118.97(15) 4216.9(3) 2013.84(18)
Z 2 1 1 2 1
ρcalcd/g cm−3 1.325 1.287 1.298 1.344 1.315
T/K 180.15 180.15 180.15 180.15 180.15
μ/mm−1 1.699 1.704 1.806 2.078 2.125
R1(I > 2σ)(I))a 0.0225 0.0175 0.0231 0.0216 0.0429
wR2a 0.0592 0.0432 0.0604 0.0479 0.1166

aR1 = ∑||F0| − |Fc|| / ∑|F0|, wR2 = [∑w(F02Fc2)2 / ∑w(F0)2]1/2.

Results and Discussion

Synthesis and crystal structure

2,2'-bipyrimidine was chosen as the bridging ligand for its general chemical stability, for its strong σ-donor and π-acceptor characteristics and its symmetric shape. In combination these properties stabilize the symmetric coordination on both 1,4-diimine coordination sites. For homo-dinuclear complexes, the electronic or magnetic interaction between two metal centers remains unperturbed by asymmetry. Additionally, the coordination of the peripheral ligands tris(tetramethylheptanedionato) lanthanide(III) for f-series metal ions results in charge-neutral homo-dinuclear compounds; an aspect potentially important for later deposition of the compounds in device environments by sublimation techniques.

Compounds 15 were synthesized by adding Ln(tmhd)3 (Ln = Gd(III) 1; Tb(III) 2; Dy(III) 3; Ho(III) 4 and Er(III) 5) to 2,2'-bipyrimidine and stirring the mixture in EtOH overnight. The molecular structure of compounds 15 is isostructural and the complexes crystallize in two different unit cells depending on the lanthanide source. Compounds 2, 3, and 5 crystallize in the triclinic [Graphic 4] space group, while 4 crystallizes in the monoclinic P21 space group; a behavior that can be attributed to the lanthanide contraction along the series of f-elements.

The molecular structure of 3 is given here in detail as a representative case for all compounds 15 (see Figure 1). The bpm moiety of the dilanthanide complexes resides on a crystallographic inversion center, whilst the individual ions are localized between the two aromatic rings of the bpm bridging ligand. The two metal centers are equivalent by symmetry. Each metal ion is eight-coordinated by six oxygen atoms from the tmhd and two nitrogen atoms of the bpm, resulting in an O6N2 donor set. The coordination polyhedron around the lanthanide ion can be most closely described as a square antiprism (see Supporting Information File 1, Table S1) (D4d, CShM of 0.605) [35,36]. The coordination polyhedron, as obtained from the single crystal analysis, is depicted in Supporting Information File 1, Figure S1. The square planes of the polyhedron are defined by the O1, O2, N1, N2 and the O3, O4, O5, O6 atoms, respectively. The lanthanide ion sits closer to the plane that is defined by four oxygen atoms. The normal vector distances between the Ln atoms and the O3, O4, O5, O6 planes are about 1.13 Å in all five compounds 15, while the O1, O2, N1, N2 planes vary between 1.42 Å and 1.45 Å (Supporting Information File 1, Table S1). The Ln(III) atoms are bound to the bpm ligand slightly above and below the plane defined by the aromatic pyrimidine rings of the bpm ligand, creating a chair-like conformation (Figure 1). The point group of the molecules varies between C1 and Ci depending on the degree of distortion within the molecule. The two coordination spheres within each dinuclear compound are more or less inversion-symmetry-related to each other, resulting in the formation of an achiral meso-form of the polyhedra (see Supporting Information File 1, Figure S1).


Figure 1: Molecular structure of complex 3 obtained from single-crystal diffraction. Dy(III) ions are marked in purple, oxygen in red, nitrogen in blue and carbon in dark grey. Hydrogen and methyl groups are omitted for clarity. Thermal ellipsoid of C, N and O (50% probability), whilst Dy is represented employing 99% probability.

The bond lengths between the lanthanide ions and the coordinating oxygen and nitrogen atoms of the ligands are summarized in Table 3. They decrease along the lanthanide series, with a corresponding decrease in ionic radii referred to as the lanthanide contraction. The distances between the two lanthanide ions have the same trend and are also showed in Table 3. These values are often observed in other known bpm-bridged Ln(III) dimers [28,29].

Table 3: Range of bond lengths between Ln and coordinating atoms, Ln···Ln distance in complexes 15.

  Ln–O [Å] Ln–N [Å] Ln···Ln [Å]
1 2.277(3)–2.363(3) 2.613(4)–2.657(4) 6.8752(2)
2 2.265(2)–2.343(2) 2.599(2)–2.616(2) 6.8124(3)
3 2.260(2)–2.330(2) 2.581(2)–2.604(2) 6.7841(3)
4 2.243(6)–2.327(5) 2.583(8)–2.596(8) 6.7933(2)
5 2.241(4)–2.296(4) 2.578(4)–2.587(4) 6.7545(3)

Magnetic properties

Static susceptibility measurements. The static magnetic behavior for all complexes was investigated on polycrystalline samples from 2 to 300 K under an applied DC magnetic field (H) of 1 kOe, whilst magnetization as a function of applied field was investigated in the field and temperature range of 0–7 T and 2–5 K, respectively.

We firstly explored the magnetic conduct of the isotropic gadolinium molecule, 1, where some insight whether any interaction (i.e. ferromagnetic or antiferromagnetic) operates within the complex. The room temperature χMT value (where χM is the molar magnetic susceptibility) for 1 is in good agreement with the expected values for two noninteracting Gd(III) ions 15.6 cm3 K mol−1 (where 15.8 cm3 K mol−1 is the expected value for two noninteracting Gd(III) ions with gGd = 2.00; S = 7/2). The χMT stays practically constant along the whole temperature range, where a very small decrease at the lowest temperature occurs (leading to a value of 15.1 cm3 K mol−1) (Figure 2). This small downturn could correspond to antiferromagnetic exchange or dipolar interactions between the Gd(III)···Gd(III) ions (see below). Similarly, the anisotropic χMT(T) behavior of the anisotropic analogues was also explored, yielding χMT values close to the expected for the sum of noninteracting Ln(III) ions. The χMT values of 23.4, 28.4, 28.2 and 24.6 cm3 K mol−1 were found for 2 to 5, respectively, whilst expected values were 23.6, 28.3, 28.1 and 24.9 cm3 K mol−1 (for two: Tb(III), gJ = 3/2, J = 6; Dy(III), gJ = 4/3, J = 15/2; Ho(III), gJ = 5/4, J = 8 or Er(III), gJ = 6/5, J = 15/2). For 3 the χMT(T) remains almost constant down to about 100 K, where it decreases reaching a value of 20.5 cm3 K mol−1 at 2 K (see Figure 3). A similar profile is observed for compounds 2, 4 and 5, where χMT(T) of 19.11, 17.8 and 15.4 cm3 K mol−1 at the lowest temperature (i.e. 2 K) are correspondingly obtained. The decrease in χMT(T) at low temperatures observed for all these compounds is due to a gradual depopulation of the Stark levels, i.e. depopulation of the crystal field split mJ sublevels of the ground J multiplet, although magnetic exchange between the Ln(III) sites in these compounds could also cause such behavior (see Figure 3).


Figure 2: (a) Experimental χMT(T) and fitting (red trace) for compound 1 (parameters for fitting are described in the text); (b) Experimental M(H,T) (open circles) and Brillouin functions at given temperatures employing g = 2.00. The M(H,T) can similarly be reproduced employing the small J used to reproduced χMT(T).


Figure 3: Molar magnetic susceptibility (χMT) vs T plot for 25 under 0.1 T DC field and molar magnetization (Mβ) as a function of applied magnetic field (H) at 2 and 5 K (inset) for (a) 2, (b) 3, (c) 4 and (d) 5.

Molar magnetization (Mβ) as function of applied magnetic field at 2 K in the field range of 0–7 T was additionally investigated for all systems. The Mβ(H) for compound 1 shows a saturation value of 13.9 μB, in very good agreement with the expected behavior obtained from the Brillouin function for noninteracting Gd(III) ions (Figure 2b). Mβ(H) measurements of complexes 25 at the maximum field (7 T) and the lowest temperature (2 K) yield values of 10.4, 11.0, 10.8 and 12.6 μB, respectively (Figure 3 inset).

The very small downturn in the χMT and the excellent agreement of the Mβ(H) data rules out any strong antiferromagnetic interaction operating within the gadolinium complex. The isotropic nature of compound 1 allows us to simulate the χMT using the simple Heisenberg Hamiltonian, taking into account a single exchange interaction [Graphic 5]. A small interaction is found sufficient to reproduce the small downturn in χMT(T) (Figure 2a), giving JGd = −0.003(2) cm−1 between neighboring Gd(III) ions with a fixed g-value of g = 2.0. The calculation of the dipolar interaction [37] between the adjacent Gd(III) ions (i.e., Gd(1)···Gd(2)) at a distance of 6.8752(2) Å, leads to Ddip = −0.002 cm−1 (for a −2J Hamiltonian), demonstrating that the decrease in χMT(T) is purely dipolar. A small zero field splitting of each individual Gd(III) could also be responsible for this behavior; however, due to the isotropic nature of Gd(III), this would be very small.

Dynamic magnetic behavior. To probe the presence of the slow relaxation of the magnetization for all anisotropic lanthanide analogues, the dynamic behavior for these was studied employing AC susceptibility measurements in the temperature range of 2–20 K and the frequency range of 0.1–1500 Hz with and without applied DC field. No SMM behavior was observed in compound 4 with or without an applied DC field. On the contrary, compound 3 shows a clear temperature dependent in-phase (χ'M(T)) and out of phase (χ"M(T)) component under zero field (Figure 4 and Supporting Information File 1, Figure S2). Likewise, compounds 2 and 5 show a frequency dependent behavior; however, no maximum is observed in either χ"M(T) or χ"M(ν), probably due to fast quantum tunneling (QT). The application of a small DC field is well known to suppress fast QT, which removes the zero-field degeneracy of the Kramers doublets at both sides of the barrier, allowing the observation of the dynamic magnetic behavior in the thermally activated regime [38,39]. We have therefore studied the optimal magnetic field (0–2 kOe), allowing the observation of the maximum in the frequency and temperature range for compounds 5 under a DC field of 800 Oe (Supporting Information File 1, Figure S3). Application of applied fields on compound 2 reveals the out-of-phase component; nevertheless, no maximum is achieved with either field in the given temperature and frequency range.


Figure 4: Experimental dynamic magnetic behavior for 3. (a) χ"M(T) measured under a zero field (3.5 Oe AC field); (b) Experimental data and fits (solid lines) for χ"M(ν); (c) Cole–Cole plots (χ"M vs χ'M) showing a single relaxation process above 4 K, where solid lines are fits of the AC susceptibility data to a modified Debye function; (d) Arrhenius treatment of χ"M data for the high temperature process.

The out-of-phase (χ"M(T)) contribution of compound 3 reveals a single maximum at about 18 K for the highest frequency (ν = 1488 Hz), which shifts to lower temperatures with lower frequencies (Figure 4a). The peak becomes broader for lower frequencies where nonthermally activated relaxation pathways dominate. In the χ"M(ν), very small shifts of the maximum are present between 2 and 4 K, characteristic of the nonthermally activated pathways taking place. Above 4 K, the maximum shifts to higher frequencies with increasing temperature, indicating a purely, thermally activated regime (Figure 4b and Supporting Information File 1, Figure S2). Fitting the χ'M(ν) and χ"M(ν) to a single relaxation process in the temperature range using a Debye [40] model leads to magnetization relaxation times (τ) that were treated using the Arrhenius law, τ = τ0exp(Ueff/kT), yielding an energy barrier (Ueff) of 97 ± 3 K, where τ0 = (6.73 ± 0.4) × 10−7 s. The nonthermally activated regime is clearly marked in the Arrhenius plot at low temperatures range between 4 and 2 K, where the relaxation time is temperature independent with a tunneling frequency of 1.1 Hz, corresponding to a tunneling time τQTM of 0.145 ms. Cole–Cole fittings give a small distribution of relaxation times for an α-parameter of 0.02 < α < 0.2 (from 5 to 18 K), where a value of 0 would indicate no distribution (Figure 4c).

Compound 5 shows no SMM behavior at zero field; however, application of a small DC field lifts the degeneracy of the barrier, leading to the observation of the out-of-phase component. As evidenced by the clear frequency dependence of the signal (Figure 5a,b), compound 5 shows SMM features. A maximum at about 3.3 K (ν = 1488 Hz) is present in the χ"M(T), shifting to lower temperatures with frequency. The peak becomes broader at lower frequencies where nonthermally activated relaxation pathways dominate. A maximum is also shown in the χ"M(ν), shifting to higher frequencies with increasing temperature. We have similarly fitted both frequency dependent profiles, χ'M(ν) and χ"M(ν), to a single relaxation process in the temperature range employing a Debye model. The treatment of the relaxation times with the Arrhenius law gives an energy barrier (Ueff) of 25 ± 2 K, τ0 = (5.11 ± 7) × 10−8 s, whilst Cole–Cole fittings render a very narrow distribution of α parameters of 0.006 < α < 0.01 (from 2 to 3.3 K) (Figure 5).


Figure 5: Experimental dynamic magnetic behavior for 5. (a) χ"M(T) measured under an applied DC field of 800 Oe (3.5 Oe AC field); (b) Experimental data and fittings (solid lines) for χ"M(ν); (c) Cole–Cole plots (χ"M vs χ'M) showing a single relaxation process, where solid lines are fits of the AC susceptibility data to a modified Debye function; (d) and Arrhenius treatment of χ"M data for the high temperature process.

Magnetic studies at low temperature. We have studied single crystals for complexes 2, 3 and 5 at mK temperatures employing a micro-SQUID apparatus. No hysteresis loop was obtained for compounds 2 (Supporting Information File 1, Figure S5) and 5 (Supporting Information File 1, Figure S6) even down to 0.03 K. The magnetization of Er2 complex 5 has a sharp increase at μ0H = 0, indicative of a very fast tunneling rate of magnetization in this complex. This fact is in agreement with the AC magnetic susceptibility measurements discussed in the earlier section, where the fast tunneling can be slowed down as soon as a small field is applied. In contrast, well-resolved two-step hysteresis loops were obtained for complex 3 (Figure 6). The width of the hysteresis loop increases with decreasing temperature and increasing sweep rate, which is characteristic of SMM behavior. The loops are very typical for two antiferromagnetically coupled Ising-like spins. Around zero field, the loops have a S-shape with two sharp tunnel steps at positive and negative field, suggesting the presence of antiferromagnetic interactions between the two Dy(III) ions. At higher fields, around ±0.5 T, the loops have a broad step, which is strongly field-sweep-rate dependent and is due to a direct relaxation process between the ferromagnetic and antiferromagnetic spin states. The loops also show a small hysteresis at μ0H = 0, which comes from the fact that some of the molecules did not tunnel to the antiferromagnetic ground state, but remain in the ferromagnetic state. A similar observation has also been reported in an antiferromagnetically coupled Dy2 complex [41,42]. The mean exchange field (Hex) can be determined from the inflexion points at about ±0.046 T (inset of Figure 6), yielding the interaction constant J = −0.0019 cm−1 using Hex = 2·J·mS/gJ·μB, where mS = 15/2 and gJ = 4/3. The strength of interactions between the two ions obtained from single-crystal measurement is comparable to that which resulted from the analysis of bulk magnetic data on the Gd compound 1.


Figure 6: Single-crystal measurements of M(H)/MS vs µ0H measured on a micro-SQUID array for 3 (a) at 0.03 K at field sweep rates from 0.002–0.28 T s−1; inset: zoom view of the loops around zero field; and (b) at a sweep rate of 0.14 T s−1 at temperatures from 0.03–5 K.

A good agreement with the dynamic behavior hysteresis loops were observed for compound 3, confirming a better performing SMM compared to compound 5 (Supporting Information File 1, Figure S5). In the light of the observation of the hysteresis loops up to 5 K (Figure 6b), we also tested if hysteresis was also observable, employing a conventional SQUID magnetometer in the temperature range of 2–5 K. Hysteresis loops were indeed observed from 2–3.3 K in the range of ±0.3 T (Supporting Information File 1, Figure S4), confirming that 3 performs rather well as an SMM even at a very slow sweeping rate, i.e. 0.0003 T s−1.

Magnetic axes. The dysprosium dimer compound, i.e. 3 exhibits a well-defined SMM behavior as evidenced in the χ"M(T), χ"M(ν) and hysteresis loops in the Mβ(H) results. This behavior has been associated with a well-defined mJ = ±15/2 ground state doublet and highly axial g-tensors (ideally gx = gy = 0; gz = 20) of the 6H15/2 manifold. To gain some insight into the orientation of the magnetic axes we employ an electrostatic method [43], which employs an electrostatic minimization of the [Graphic 6] Sievers electron density [44] and a minimal valence bond model. The method gives the directionality of the magnetic axes for compound 3, being collinear between the two Dy(III) sites, as expected from symmetric considerations (Figure 7) and almost perpendicularly aligned to the bipyrimidine plane. In such a case, the dipolar field between the two Dy(III) ions can easily be calculated [4] to be 0.03 T, which is less than 0.046 T observed in the micro-SQUID measurements. This means that the main contribution of the interactions between the ions is dipolar but some exchange coupling contributes as well. Undoubtedly, this co-parallel alignment is responsible for the SMM profile observed. Quantum tunneling is commonly accelerated by the noncollinear magnetic arrangement [45] and small exchange interactions [45-48], introducing further relaxation pathways. In compound 3 this is not the case, therefore the magnetic behavior could be associated with small exchange interactions, colinearity of the magnetic axis and a well-defined ground state of the individual Dy(III) ions.


Figure 7: Magnetic axes obtained through the electrostatic method: (a) side and (b) top view for compound 3. Color code: same as in Figure 1.

Scanning tunneling microscopy of [Tb(tmhd)3]2bpm

Scanning tunneling microscopy (STM) studies were performed at 5 K to investigate the structure and behavior of the compound at the single-molecule level. Several repetitions of the experiments suggested that the original molecular structure of 2 is substantially altered when deposited on the Au(111) surface. The experimental results of the deposition are summarized in Figure 8. There, the molecular phase forms a hexagonal surface pattern in a monolayer arrangement on the lower terrace of the Au(111) crystal plane (Figure 8a). Apart from the hexagonal molecular phase, small numbers of bright protrusions form a cluster-type arrangement on the upper terrace of the gold crystal (Figure 8b). The shape of these clusters as well as of the molecular phase differs substantially from the projected molecular structure, which can be derived from the single crystal measurement of the bulk compound (see Figure 1).


Figure 8: STM topography of 2 after deposition on Au(111). Image sizes are (a) 60 × 60 nm2 and (c) 20 × 20 nm2. White scale bars in the lower right corner in the images (a) and (c) indicate the length of 5 nm. (b) Cross section along the thin white line in (a). Set point: (a) 2.5 V, 10 pA, (c) 2 V, 40 pA. Length of the respective lines: black, 3.4 nm; red, 3.0 nm; green, 1.4 nm; blue, 2.0 nm; violet, 1.0 nm.

Obviously, some of the molecules are decomposed on the Au substrate, but the majority of the molecules seem to self-organize as trimolecular nodes. The cross section of the topography is shown in Figure 8b. The heights of the small clusters range from 350 to 450 pm. Heights of 400 pm and 280 pm have been reported on mononuclear ruthenium tris(β-diketonato) complexes on Ag(111) for two different conformations [49]. Thus, the height of the film (450 pm) is reasonable for the height of 2, indicating a slightly distorted surface coordination in comparison to the literature. The small clusters in Figure 8a (left side) are likely decomposition products of [Tb(tmhd)3]2bpm (2), such as Tb(thmd)3, which may form at the gold surface. A magnified view of the molecular film is shown in Figure 8c. The pronounced hexagonal lattice of molecules is very regular and consists of triangular repetition. There, three bright protrusions form a triangle, which is highlighted in violet. The centers of the three dots lay in a distance range of about 1 nm, which corresponds well with the distance of the tert-butyl groups of the tmhd ligands between two adjacent molecules of the single crystal data. The average size of the nearly round bright protrusions is about 0.7 nm, which could describe the entire inner diameter of the Tb(tmhd)3 subunit within one dinuclear molecule. The inter-triangular distance from center to center (blue line) is 2 nm.

The underlying cause of the formation of the hexagonal structure remains unclear. However, size and shape analyses suggest that three molecules of 2 interlock with each other and form a piano stool subunit within the surface pattern. Further methodical investigation is currently underway. The STM study implies that the [Tb(tmhd)3]2bpm (2) molecule seems to be transferable onto the metal substrate with a certain amount of decomposition.


The synthesis, molecular structure and magnetic properties of a series of five isostructural dilanthanide complexes 1–5 with a bipyrimidine bridging ligand were reported. The magnetic characteristics of all complexes were investigated by DC and AC SQUID-measurements, leading to the observation of SMM behavior in the Dy2- and Er2-containing systems, 3 and 5. Clear out-of-phase components were observed for the Dy2-analogue 3 under zero field, whilst application of a small bias DC field to the Er2-analogue 5 slows down the quantum tunneling rate revealing the SMM properties. Furthermore, micro-SQUID studies show hysteresis loops for these complexes demonstrating that they retain their magnetization below a certain temperature, i.e. up to 5 K for the Dy2 complex 3. To gain some insight into the molecular orientation of the magnetic axes of the two lanthanide ions, we have employed an electrostatic method, which gives a parallel alignment of the axes of the Dy(III) ions. These results demonstrate that SMM behavior can be achieved by linking two lanthanide metal ions (which exhibit single-ion magnetic anisotropy) with a 2,2’-bipyrimidine bridging ligand. The observed SMM character, with hysteresis loops observed as high as 5 K, make this class of bipyrimidine-bridged dilanthanide complexes promising systems to be sublimated onto surfaces. In this way, it is possible to study their magnetic behavior as single molecules or in thin film compositions. Due to its charge-neutral character, it was possible to sublime [Tb(tmhd)3]2bpm (2) onto a Au(111) surface. Preliminary results from scanning tunneling microscopy at 5 K suggest the formation of a homogeneous molecular film. Based on this result, it can be envisioned that the presented class of molecules, eventually equipped with linker-substituents, will function as active molecular entities that could be combined into spintronic hybrid device environments [14,50].

Supporting Information

Supporting Information File 1: Additional experimental information.
Format: PDF Size: 844.3 KB Download


This work has been supported by “Simple Quantum Operations using Surface-deposited Molecular Magnets” released in the call “Functional Nanostructures” from the German Baden-Würtemberg Foundation. W.W. and M.R. acknowledge support from the EU with the FP7 FET-Proactive project MoQuaS N. 610449 and with the AgenceNationale de la Recherche project MolQuSpin, N. ANR-13-BS10.


  1. Bogani, L.; Wernsdorfer, W. Nat. Mater. 2008, 7, 179–186. doi:10.1038/nmat2133
    Return to citation in text: [1]
  2. Sanvito, S. Chem. Soc. Rev. 2011, 40, 3336–3355. doi:10.1039/c1cs15047b
    Return to citation in text: [1]
  3. Žutić, I.; Fabian, J.; Das Sarma, S. Rev. Mod. Phys. 2004, 76, 323–410. doi:10.1103/RevModPhys.76.323
    Return to citation in text: [1]
  4. Gatteschi, D.; Sessoli, R.; Villain, J. Molecular Nanomagnets; Oxford University Press: Oxford, United Kingdom, 2006.
    Return to citation in text: [1] [2]
  5. Ladd, T. D.; Jelezko, F.; Laflamme, R.; Nakamura, Y.; Monroe, C.; O'Brien, J. L. Nature 2010, 464, 45–53. doi:10.1038/nature08812
    Return to citation in text: [1]
  6. Leuenberger, M. N.; Loss, D. Nature 2001, 410, 789–793. doi:10.1038/35071024
    Return to citation in text: [1]
  7. Sessoli, R.; Gatteschi, D.; Caneschi, A.; Novak, M. A. Nature 1993, 365, 141–143. doi:10.1038/365141a0
    Return to citation in text: [1]
  8. Sessoli, R.; Tsai, H. L.; Schake, A. R.; Wang, S.; Vincent, J. B.; Folting, K.; Gatteschi, D.; Christou, G.; Hendrickson, D. N. J. Am. Chem. Soc. 1993, 115, 1804–1816. doi:10.1021/ja00058a027
    Return to citation in text: [1]
  9. Aromí, G.; Brechin, E. K. Synthesis of 3d Metallic Single-Molecule Magnets. In Single-Molecule Magnets and Related Phenomena; Winpenny, R., Ed.; Structure and Bonding, Vol. 122; Springer: Berlin, Germany, 2006; pp 1–67. doi:10.1007/430_022
    Return to citation in text: [1]
  10. Zadrozny, J. M.; Xiao, D. J.; Atanasov, M.; Long, G. J.; Grandjean, F.; Neese, F.; Long, J. R. Nat. Chem. 2013, 5, 577–581. doi:10.1038/nchem.1630
    Return to citation in text: [1]
  11. Rinehart, J. D.; Long, J. R. Chem. Sci. 2011, 2, 2078–2085. doi:10.1039/c1sc00513h
    Return to citation in text: [1]
  12. Woodruff, D. N.; Winpenny, R. E. P.; Layfield, R. A. Chem. Rev. 2013, 113, 5110–5148. doi:10.1021/cr400018q
    Return to citation in text: [1]
  13. Kyatskaya, S.; Mascarós, J. R. G.; Bogani, L.; Hennrich, F.; Kappes, M.; Wernsdorfer, W.; Ruben, M. J. Am. Chem. Soc. 2009, 131, 15143–15151. doi:10.1021/ja906165e
    Return to citation in text: [1]
  14. Stepanow, S.; Honolka, J.; Gambardella, P.; Vitali, L.; Abdurakhmanova, N.; Tseng, T.-C.; Rauschenbach, S.; Tait, S. L.; Sessi, V.; Klyatskaya, S.; Ruben, M.; Kern, K. J. Am. Chem. Soc. 2010, 132, 11900–11901. doi:10.1021/ja105124r
    Return to citation in text: [1] [2]
  15. Thiele, S.; Balestro, F.; Ballou, R.; Klyatskaya, S.; Ruben, M.; Wernsdorfer, W. Science 2014, 344, 1135–1138.
    Return to citation in text: [1]
  16. Cornia, A.; Mannini, M.; Sainctavit, P.; Sessoli, R. Chem. Soc. Rev. 2011, 40, 3076–3091. doi:10.1039/c0cs00187b
    Return to citation in text: [1]
  17. Aguilà, D.; Barrios, L. A.; Velasco, V.; Roubeau, O.; Repollés, A.; Alonso, P. J.; Sesé, J.; Teat, S. J.; Luis, F.; Aromí, G. J. Am. Chem. Soc. 2014, 136, 14215–14222. doi:10.1021/ja507809w
    Return to citation in text: [1]
  18. Barquín, M.; Garmendia, M. J. G.; Bellido, V. Transition Met. Chem. 1999, 24, 584–588. doi:10.1023/A:1006940930679
    Return to citation in text: [1]
  19. Wagner, S.; Kisslinger, F.; Ballmann, S.; Schramm, F.; Chandrasekar, R.; Bodenstein, T.; Fuhr, O.; Secker, D.; Fink, K.; Ruben, M.; Weber, H. B. Nat. Nanotechnol. 2013, 8, 575–579. doi:10.1038/nnano.2013.133
    Return to citation in text: [1]
  20. Zhang, L.; Bagrets, A.; Xenioti, D.; Korytár, R.; Schackert, M.; Miyamachi, T.; Schramm, F.; Fuhr, O.; Chandrasekar, R.; Alouani, M.; Ruben, M.; Wulfhekel, W.; Evers, F. Phys. Rev. B 2015, 91, 195424. doi:10.1103/PhysRevB.91.195424
    Return to citation in text: [1] [2]
  21. Herrera, J.-M.; Pope, S. J. A.; Adams, H.; Faulkner, S.; Ward, M. D. Inorg. Chem. 2006, 45, 3895–3904. doi:10.1021/ic0521574
    Return to citation in text: [1]
  22. Jaradat, Q.; Akasheh, T. Polyhedron 1987, 6, 175–179. doi:10.1016/S0277-5387(00)80786-X
    Return to citation in text: [1]
  23. Demir, S.; Zadrozny, J. M.; Nippe, M.; Long, J. R. J. Am. Chem. Soc. 2012, 134, 18546–18549. doi:10.1021/ja308945d
    Return to citation in text: [1]
  24. Irfanullah, M.; Iftikhar, K. Inorg. Chim. Acta 2013, 394, 373–384. doi:10.1016/j.ica.2012.08.015
    Return to citation in text: [1]
  25. Irfanullah, M.; Iftikhar, K. J. Photochem. Photobiol., A: Chem. 2011, 224, 91–101. doi:10.1016/j.jphotochem.2011.09.011
    Return to citation in text: [1]
  26. Irfanullah, M.; Iftikhar, K. J. Lumin. 2010, 130, 1983–1993. doi:10.1016/j.jlumin.2010.05.015
    Return to citation in text: [1]
  27. Zucchi, G.; Jeon, T.; Tondelier, D.; Aldakov, D.; Thuéry, P.; Ephritikhine, M.; Geffroy, B. J. Mater. Chem. 2010, 20, 2114. doi:10.1039/b921740a
    Return to citation in text: [1]
  28. Zucchi, G.; Le Goff, X. F. Polyhedron 2013, 52, 1262–1267. doi:10.1016/j.poly.2012.06.021
    Return to citation in text: [1] [2]
  29. Baker, M. H.; Dorweiler, J. D.; Ley, A. N.; Pike, R. D.; Berry, S. M. Polyhedron 2009, 28, 188–194. doi:10.1016/j.poly.2008.10.019
    Return to citation in text: [1] [2] [3]
  30. Richards, G.; Osterwyk, J.; Flikkema, J.; Cobb, K.; Sullivan, M.; Swavey, S. Inorg. Chem. Commun. 2008, 11, 1385–1387. doi:10.1016/j.inoche.2008.09.011
    Return to citation in text: [1] [2]
  31. Bain, G. A.; Berry, J. F. J. Chem. Educ. 2008, 85, 532. doi:10.1021/ed085p532
    Return to citation in text: [1]
  32. Wernsdorfer, W. Classical and Quantum Magnetization Reversal Studied in Nanometer-Sized Particles and Clusters. In Handbook of Advanced Magnetic Materials; Liu, Y.; Sellmyer, D. J.; Shindo, D., Eds.; Springer: Berlin, Germany, 2007; pp 77–127. doi:10.1007/1-4020-7984-2_3
    Return to citation in text: [1]
  33. Sheldrick, G. M. Acta Crystallogr., Sect. A: Found. Adv. 2008, A64, 112–122. doi:10.1107/S0108767307043930
    Return to citation in text: [1]
  34. Dolomanov, O. V.; Bourhis, L. J.; Gildea, R. J.; Howard, J. A. K.; Puschmann, H. J. Appl. Crystallogr. 2009, 42, 339–341. doi:10.1107/S0021889808042726
    Return to citation in text: [1]
  35. Alvarez, S.; Alemany, P.; Casanova, D.; Cirera, J.; Llunell, M.; Avnir, D. Coord. Chem. Rev. 2005, 249, 1693–1708. doi:10.1016/j.ccr.2005.03.031
    Return to citation in text: [1]
  36. Casanova, D.; Alemany, P.; Bofill, J. M.; Alvarez, S. Chem. – Eur. J. 2003, 9, 1281–1295. doi:10.1002/chem.200390145
    Return to citation in text: [1]
  37. Abragam, A.; Bleaney, B. Electron Paramagnetic Resonance of Transition Ions; Oxford University Press: Oxford, United Kingdom, 2012.
    Return to citation in text: [1]
  38. Ako, A. M.; Mereacre, V.; Hewitt, I. J.; Clérac, R.; Lecren, L.; Anson, C. E.; Powell, A. K. J. Mater. Chem. 2006, 16, 2579–2586. doi:10.1039/b604611h
    Return to citation in text: [1]
  39. Palacios, M. A.; Titos-Padilla, S.; Ruiz, J.; Herrera, J. M.; Pope, S. J. A.; Brechin, E. K.; Colacio, E. Inorg. Chem. 2014, 53, 1465–1474. doi:10.1021/ic402597s
    Return to citation in text: [1]
  40. Aubin, S. M. J.; Sun, Z. M.; Pardi, L.; Krzystek, J.; Folting, K.; Brunel, L.-C.; Rheingold, A. L.; Christou, G.; Hendrickson, D. N. Inorg. Chem. 1999, 38, 5329–5340. doi:10.1021/ic990613g
    Return to citation in text: [1]
  41. Chow, C. Y.; Bolvin, H.; Campbell, V. E.; Guillot, R.; Kampf, J. W.; Wernsdorfer, W.; Gendron, F.; Autschbach, J.; Pecoraro, V. L.; Mallah, T. Chem. Sci. 2015, 6, 4148–4159. doi:10.1039/C5SC01029B
    Return to citation in text: [1]
  42. Sulway, S. A.; Layfield, R. A.; Tuna, F.; Wernsdorfer, W.; Winpenny, R. E. P. Chem. Commun. 2012, 48, 1508–1510. doi:10.1039/C1CC14643B
    Return to citation in text: [1]
  43. Chilton, N. F.; Collison, D.; McInnes, E. J. L.; Winpenny, R. E. P.; Soncini, A. Nat. Commun. 2013, 4, No. 2551. doi:10.1038/ncomms3551
    Return to citation in text: [1]
  44. Sievers, J. Z. Phys. B: Condens. Matter 1982, 45, 289–296. doi:10.1007/BF01321865
    Return to citation in text: [1]
  45. Pineda, E. M.; Chilton, N. F.; Marx, R.; Dörfel, M.; Sells, D. O.; Neugebauer, P.; Jiang, S.-D.; Collison, D.; van Slageren, J.; McInnes, E. J. L.; Winpenny, R. E. P. Nat. Commun. 2014, 5, No. 5243. doi:10.1038/ncomms6243
    Return to citation in text: [1] [2]
  46. Pineda, E. M.; Chilton, N. F.; Tuna, F.; Winpenny, R. E. P.; McInnes, E. J. L. Inorg. Chem. 2015, 54, 5930–5941. doi:10.1021/acs.inorgchem.5b00746
    Return to citation in text: [1]
  47. Sessoli, R.; Powell, A. K. Coord. Chem. Rev. 2009, 253, 2328–2341. doi:10.1016/j.ccr.2008.12.014
    Return to citation in text: [1]
  48. Sorace, L.; Benelli, C.; Gatteschi, D. Chem. Soc. Rev. 2011, 40, 3092–3104. doi:10.1039/c0cs00185f
    Return to citation in text: [1]
  49. Munery, S.; Ratel-Ramond, N.; Benjalal, Y.; Vernisse, L.; Guillermet, O.; Bouju, X.; Coratger, R.; Bonvoisin, J. Eur. J. Inorg. Chem. 2011, 2698–2705. doi:10.1002/ejic.201100116
    Return to citation in text: [1]
  50. van der Molen, S. J.; Calame, M.; Reed, M.; Naaman, R.; Scheer, E.; Neaton, J. B.; Nitzan, A.; Natelson, D.; Tao, N. J.; van der Zant, H.; Mayor, M.; Ruben, M. Nat. Nanotechnol. 2013, 8, 385–389.
    Return to citation in text: [1]

© 2016 Yu et al; licensee Beilstein-Institut.
This is an Open Access article under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
The license is subject to the Beilstein Journal of Nanotechnology terms and conditions: (

Back to Article List