K.
Randall McClain‡
a,
Colin A.
Gould‡
b,
Khetpakorn
Chakarawet
b,
Simon J.
Teat
c,
Thomas J.
Groshens
a,
Jeffrey R.
Long
*bde and
Benjamin G.
Harvey
*a
aUS Navy, Naval Air Warfare Center, Weapons Division, Research Department, Chemistry Division, China Lake, California 93555, USA. E-mail: benjamin.g.harvey@navy.mil
bDepartment of Chemistry, University of California, Berkeley, Berkeley, California 94720, USA. E-mail: jrlong@berkeley.edu
cAdvanced Light Source, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
dMaterials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
eDepartment of Chemical and Biomolecular Engineering, University of California, Berkeley, Berkeley, California 94720, USA
First published on 18th October 2018
A series of dysprosium(III) metallocenium salts, [Dy(CpiPr4R)2][B(C6F5)4] (R = H (1), Me (2), Et (3), iPr (4)), was synthesized by reaction of DyI3 with the corresponding known NaCpiPr4R (R = H, iPr) and novel NaCpiPr4R (R = Me, Et) salts at high temperature, followed by iodide abstraction with [H(SiEt3)2][B(C6F5)4]. Variation of the substituents in this series results in substantial changes in molecular structure, with more sterically-encumbering cyclopentadienyl ligands promoting longer Dy–C distances and larger Cp–Dy–Cp angles. Dc and ac magnetic susceptibility data reveal that these structural changes have a considerable impact on the magnetic relaxation behavior and operating temperature of each compound. In particular, the magnetic relaxation barrier increases as the Dy–C distance decreases and the Cp–Dy–Cp angle increases. An overall 45 K increase in the magnetic blocking temperature is observed across the series, with compounds 2–4 exhibiting the highest 100 s blocking temperatures yet reported for a single-molecule magnet. Compound 2 possesses the highest operating temperature of the series with a 100 s blocking temperature of 62 K. Concomitant increases in the effective relaxation barrier and the maximum magnetic hysteresis temperature are observed, with 2 displaying a barrier of 1468 cm−1 and open magnetic hysteresis as high as 72 K at a sweep rate of 3.1 mT s−1. Magneto-structural correlations are discussed with the goal of guiding the synthesis of future high operating temperature DyIII metallocenium single-molecule magnets.
Recently, the DyIII metallocenium salt [Dy(Cpttt)2][B(C6F5)4] (Cpttt = 1,2,4-tri(tert-butyl)cyclopentadienyl) was reported to show magnetic hysteresis at temperatures up to 60 K, as observed using a sweep rate of 2.2 mT s−1,5 a substantial increase over the previous record of 28 K.4c Significantly, this molecule represents an outlier in its class: while the majority of mononuclear single-molecule magnets show only waist-constricted magnetic hysteresis and little or no remanent magnetization due to quantum tunneling of the magnetization,6 [Dy(Cpttt)2]+ retains polarization at zero field and high temperatures.7 Magnetic and computational studies suggest that the impressive magnetic properties of [Dy(Cpttt)2]+ arise from its unique coordination environment—in particular, the rigid η5-cyclopentadienyl ligands constrain metal–ligand vibrational modes, thereby reducing pathways for through-barrier relaxation.5b This represents a new approach to limiting through-barrier magnetic relaxation and enabling higher operating temperatures. It is thus essential to develop a deeper understanding of how structure influences the magnetic properties of [Dy(CpX)2]+ complexes in order to determine how coordination environment can be tuned to further increase operating temperatures.8
The design of mononuclear 4f-element single-molecule magnets can be guided by a simple electrostatic model that compares the shape of lanthanide MJ electron distributions to the crystal field.9 A number of studies of mononuclear, D5h-symmetric complexes with strongly-donating axial ligands have demonstrated the validity of this qualitative approach, with molecules in this class possessing thermal barriers to magnetic relaxation as large as 1260 cm−1.3b–f Importantly, the electrostatic model provides a rationale for the high anisotropy engendered by the axial cyclopentadienyl ligands in [Dy(Cpttt)2]+ and outlines a clear route for improvement. Increasing the Cp–Dy–Cp angle of this bent metallocene from 152.56(7)° toward 180° and shortening the Dy–Cp bond distances should increase the axiality of the crystal field, leading to larger thermal barriers to magnetic relaxation and higher potential operating temperatures.
Equally important in the design of new [Dy(CpX)2]+ molecules are the metal–ligand vibrational modes. Phonons—vibrational modes of the lattice that couple to the electronic states of a molecule—provide the necessary energy for magnetic relaxation10 and determine viable relaxation processes at a given temperature. Computational studies indicate that high-temperature magnetic relaxation in [Dy(Cpttt)2]+ is induced by Cp–H bending modes,5b leading to the suggestion that substituting the hydrogen positions in this molecule with rigid functional groups could result in more constrained vibrations, further slowing magnetic relaxation and resulting in magnetic hysteresis at higher temperatures.
Herein, we report a series of DyIII metallocenium cation salts, [Dy(CpiPr4R)2][B(C6F5)4] (R = H (1), Me (2), Et (3), iPr (4)), in which subtle variation of the substituents on the cyclopentadienyl rings produces molecules with a large range of Dy–C distances and Cp–Dy–Cp angles. This series includes DyIII metallocenium salts of two new cyclopentadienyl ligands, CpiPr4R (R = Me, Et). By varying the substituents, it is also possible to tune the metal–ligand vibrational modes in these complexes. Modifications result in a record high thermal barrier to magnetic relaxation (Ueff) of 1468 cm−1 and a record high 100 s magnetic blocking temperature of 62 K for the methyl-substituted complex, 2. A careful comparison of the structural and magnetic properties in this series elucidates important design principles for the synthesis of further new DyIII metallocenium single-molecule magnets.
In order to complete a stepwise series of metallocenium salts in which the R group in CpiPr4R varies by only one methyl or methylene group in adjacent members of the series, it was necessary to synthesize the heretofore unknown tetraisopropyl substituted Cp ligand salts NaCpiPr4R (R = Me, Et). This was achieved by exploitation of the fulvene route used to prepare NaCpiPr5. Here, tetraisopropylfulvene and methyl-substituted tetraisopropylfulvene12 reacted with excess LiAlH4 and were then hydrolyzed to produce the corresponding cyclopentadiene isomer mixtures. These were then converted to the sodium salts NaCpiPr4Me and NaCpiPr4Et with sonication in the presence of excess NaNH2 (Scheme 1). In all cases, the sodium salts were preferred for their moderate solubility and ready accessibility as ether-free compounds, allowing facile workup in subsequent steps. For example, the lithium salt of CpiPr5 is obtained as the diethyl etherate and is soluble in pentane,12 whereas the sodium salt can be freed of coordinated THF by heating under vacuum and is soluble in arenes but insoluble in alkanes.
With the Cp transfer reagents in hand, we turned to the synthesis of the neutral dysprosium compounds, Dy(CpiPr4R)2X (R = H, Me, Et, iPr) (X = Cl, Br, I). The most common strategy for transfer of Cp ligands to LnIII centers involves the reaction of MCpX (commonly M = Li, Na, K) with a lanthanide halide. For Sm, Eu, and Yb, which have relatively stable divalent oxidation states, the initial reaction can occur between MCpX and LnI2, followed by oxidation to obtain the LnIII complex; an example of this approach is the well-known synthesis of (Cp∗ = pentamethylcyclopentadienyl).14
In many cases, the lanthanide iodide is preferred over the other halides, especially in the case of very bulky ligand environments.15 For instance, the reaction of YCp″2X (Cp″ = 1,3-bis(trimethylsilyl)cyclopentadienyl, X = Cl, I) with KCp″ in refluxing toluene to produce YCp″3 was most successful when X = I.16 Our own exploratory reactions with YCl3 and excess NaCpiPr5 (2.2 equiv.) in toluene at 130 °C indicated YCpiPr5Cl2 as the dominant product, and thus we opted for DyI3 or DyI3(THF)3 as the most promising DyIII sources. In some cases, the THF adduct allows for substitutions that are not possible with the unsolvated precursor, such as in the reaction of YbCp″2I(THF) with KCp″ to form YbCp″3.16 Conversely, the presence of THF can result in unproductive side reactions.17 In our hands, reaction of DyI3(THF)3 with NaCpiPr5 in toluene at 130 °C resulted in isolation of THF ring-opened products, as has been observed previously during attempts to synthesize from the di-substituted SmIII salt.14 Therefore, DyI3 was used in all subsequent reactions.
In general, DyI3 and a Cp transfer reagent NaCpiPr4R (R = H, Me, Et, iPr) (2.5 equiv.) were heated together at 160 °C in toluene for three days to generate the putative intermediate complexes Dy(CpiPr4R)2I (R = H, Me, Et, iPr), which were isolated from excess ligand and NaI by extraction with pentane (Scheme 2). A pressure flask was used to obtain the high temperature necessary for complete di-substitution.
Crystallization of Dy(CpiPr4R)2I from pentane at low temperature (−35 °C) was largely fruitless owing to the high solubility of the compounds. Instead, the crude products were used directly in the subsequent halide abstraction step. The addition of a sub-stoichiometric quantity (0.8 equiv.) of [H(SiEt3)2][B(C6F5)4]—freshly prepared before use18—to pale yellow benzene solutions of Dy(CpiPr4R)2I (R = H, Me, Et, iPr) resulted in a rapid color change to a more vibrant color ranging from yellow to orange. In most cases, yellow-orange powdery to oily precipitate was also observed, which changed over the course of several days, typically becoming more solid and increasing in volume. It is known that reactions involving [SiEt3]+ in benzene are complex, forming biphasic mixtures.18 This may be the source of the oily precipitates observed initially in some cases. It was assumed that the initial challenging reaction of DyI3 with NaCpiPr4R was not quantitative, therefore a sub-stoichiometric (versus DyI3) quantity of [H(SiEt3)2][B(C6F5)4] was employed in the reactions. Here, we wanted to avoid the presence of excess [SiEt3]+ in the reaction mixtures, which can render purification more difficult. After stirring for several days at room temperature, the desired DyIII metallocene salts [Dy(CpiPr4R)2][B(C6F5)4] (R = H (1), Me (2), Et (3), iPr (4)) were obtained as crude powders (Scheme 2) and subsequently recrystallized 2–3 times by layering dichloromethane solutions of the salts with pentane. In this way, analytically pure material and crystals suitable for X-ray diffraction analysis were obtained. After drying under reduced pressure, compounds 1–4 were obtained as yellow (1), orange-yellow (2–3) or orange crystals (4) and found to be indefinitely stable at room temperature under argon; however at ambient temperatures in air they decay rapidly with an accompanying color change to green/blue.
To obtain solid-state dilutions of 1-4, it was necessary to synthesize the diamagnetic YIII analogues, [Y(CpiPr4R)2][B(C6F5)4] (R = H (Y1), Me (Y2), Et (Y3), iPr (Y4)) (see ESI† for experimental procedures). Utilizing an analogous route to that employed for 1–4, YI3 was converted to the putative Y(CpiPr4R)2I (R = H, Me, Et, iPr) complex by reaction with the appropriate NaCpiPr4R (R = H, Me, Et, iPr) (2.5 equiv.) at high temperature, followed by iodide abstraction with [H(SiEt3)2][B(C6F5)4] (Scheme 2). The YIII metallocenium complexes were obtained as colorless (Y1) to yellow crystals (Y2–Y4) after multiple recrystallizations followed by prolonged drying under vacuum. Like the DyIII metallocenes salts, Y1–Y4 are stable under argon at ambient temperature, but decay rapidly in air as indicated by the appearance of a green/blue color.
NMR spectra were obtained for the diamagnetic YIII analogues Y1–Y4 to assess their structural dynamics in solution (Fig. S15–S26†). Y1 exhibited sharp, fully resolved signals at 25 °C in the 1H-NMR spectra. Two methine (–CH(CH3)2) septets and four methyl (–CH(CH3)2) doublets were observed, consistent with restricted Cp–iPr bond rotation giving rise to separate methyl group (inner and outer CH3) environments.19 At 25 °C, Y4 displayed somewhat broadened peaks, which sharpened at −10 °C to three methyl (–CH(CH3)2) doublets (δ = 1.07, 1.66, 1.69) and one methine (–CH(CH3)2) septet. The two closely spaced downfield doublets probably correspond to the inner methyl group of the rac and meso isomers, while the upfield doublet likely corresponds to the outer methyl group of both isomers. Similar features have been associated with rac and meso isomers for [(CpiPr5)2M] (M = Ca, Sr, Ba).20 The 13C spectrum of Y4 also indicates rac and meso isomers, with four methyl (–CH(CH3)2) and two methine (–CH(CH3)2) peaks. The 1H spectra of Y2 and Y3 vary substantially with temperature. For Y2, one methyl singlet (Cp–CH3), two overlapping methine (–CH(CH3)2) septets and broadened peaks representing the four methyl (–CH(CH3)2) environments are observed at 25 °C, however, when heated to 75 °C, the four methyl (–CH(CH3)2) doublets are sharp and fully resolved. In the case of Y3 at 25 °C, sharp peaks are observed for the ethyl (Cp–CH2CH3) group and one of the methine (–CH(CH3)2) signals, while one methine (–CH(CH3)2) group is broadened and the methyl (–CH(CH3)2) peaks are observed as a single broad resonance. Heating to 85 °C results in sharp peaks for all groups, although three of the four methyl (–CH(CH3)2) doublets are overlapping.
The solid-state structures of the metallocene salts 1–4 were elucidated by single-crystal X-ray diffraction analysis (Fig. 1). Compounds 1 and 2 crystallize with two [Dy(CpiPr4R)2][B(C6F5)4] molecules in the asymmetric unit, while compounds 3 and 4 crystallize with only one molecule in the asymmetric unit. The structural analyses were complicated by pseudomerohedral twinning for 1 and 2 and positional disorder in 2 and 4. Two of the cyclopentadienyl ligands in the structure of 2 are disordered, each over two positions, and the DyIII center in the structure of 4 displays positional disorder over four sites (see ESI† for refinement details). Twinning and disorder were found to be consistent over several tested crystal samples.
The Dy–Cp separations, as assessed by the distances from the Dy atom to the centroid of each cyclopentadienyl ligand, lie in the ranges 2.27(1)–2.30(1) Å for 1, 2.273(3)–2.382(3) Å for 2, 2.297(4)–2.306(4) Å for 3, and 2.245(15)–2.392(16) Å for 4. Because of the disorder present in each structure and the similarity in the resulting bond distances (within statistical error), it was not possible to derive a meaningful structural trend from these distances alone. Nonetheless, a trend can be drawn from the major component in each of the crystal structures: the average Dy–Cp(centroid) distances for the major component are 2.29(1), 2.298(5), 2.302(6), and 2.340(7) Å for 1–4, respectively (Table 1). Although the average Dy–Cp(centroid) distance does not vary by more than the statistical error in some adjacent members of the series (e.g., 1 and 2), the change is statistically significant across the series as a whole. Similarly, the average Dy–C distances for the major component are 2.587(5), 2.600(3), 2.596(2), and 2.625(3) Å for 1–4, respectively. The increase of both the average Dy–Cp(centroid) and Dy–C distances suggests elongation of the interaction between the Dy atom and the cyclopentadienyl rings as the series progresses from 1 to 4. In addition, the Cp–Dy–Cp angle increases across the series, lying in the ranges 146.5(8)–148.0(8), 155.9(3)–157.3(4), 161.1(2), and 157.8(2)–167.9(15)° in 1–4, respectively (Table 1). The angle between the planes formed by the cyclopentadienyl rings follows a similar trend with values in the range of 148.1–150.5, 162.3–162.8, 165.0, and 172.5° for 1–4, respectively. While the cyclopentadienyl rings in 4 are nearly parallel, the Dy atom moves off-center, leading to a reduced Cp–Dy–Cp angle. These structural trends are consistent with the expected steric environment imposed by the cyclopentadienyl ligands; that is, more sterically-encumbered substituents prompt the two cyclopentadienyl ligands in [Dy(CpiPr4R)2]+ to be positioned further away from each other, increasing linearity and lengthening the average Cp–Dy bond distance.
1 | 2 | 3 | 4 | [(Cpttt)2Dy]+e | |
---|---|---|---|---|---|
a Average value for all positions in the crystal structure. b Average value for the highest occupancy component in the crystal structure. c T b is defined as the 100 s magnetic blocking temperature. d T max is the maximum hysteresis temperature. e Values are an average of those reported in ref. 5. | |||||
Cp–Dy–Cpa (°) | 147.2(8) | 156.6(3) | 161.1(2) | 162.1(7) | 152.70(7) |
Dy–Cpb (Å) | 2.29(1) | 2.298(5) | 2.302(6) | 2.340(7) | 2.316(3) |
U eff (cm−1) | 1285 | 1468 | 1380 | 1334 | 1223 or 1277 |
T b (K) | 17 | 62 | 59 | 56 | 53 |
τ tunnel (s) | 439 | 2452 | 447 | 1187 | n/a |
T max (K) | 32 | 72 | 66 | 66 | 60 |
2 K Mr (µB) | 2.6 | 4.4 | 3.9 | 4.2 | 4.2 |
2 K Hc (T) | 1.4 | 2.4 | 1.4 | 1.8 | 2.4 |
The solid-state structures of the yttrium metallocene salts Y1–Y4 were also investigated by single-crystal X-ray diffraction (Fig. S31–S34†). As in 1–4, positional disorder complicated the structural analyses. Compounds Y1 and Y2 crystalize with two molecules in the asymmetric unit; both the YIII center and one cyclopentadienyl ligand display positional disorder over two sites.
Similar structural trends to those seen in 1–4 are observed for the series of YIII congeners. Namely the average Cp–Y bond distance and Cp–Y–Cp angle generally increase across the series from Y1 to Y4. The average Y–Cp(centroid) distances for the major component are 2.273(10), 2.305(6), 2.294(3), and 2.331(2) Å for Y1–Y4, respectively, and the average Cp–Y–Cp angles are 145.7(5), 157.1(4), 161.3(1) and 180°. While the Dy atom in 4 moves off-center of the cyclopentadienyl rings, the Y atom in Y4 sits on an inversion center in the asymmetric unit, leading to a 180° Cp–Y–Cp angle. While both rac and meso isomers of Y4 are observed in the 1H- and 13C-NMR spectra, only the meso isomer is observed in the crystal structure.
In order to investigate the effect of the structural changes on the magnetic properties of each complex, dc magnetic susceptibility data were collected for 1–4 from 2 to 300 K under an applied magnetic field of 1 kOe (Fig. S35, S37, S39 and S41†). The room temperature χMT values for 2 and 3 are 14.17 and 14.28 emu K mol−1, respectively, agreeing well with the expected value of 14.17 emu K mol−1 for a free DyIII ion. The room temperature χMT values for 1 and 4 are slightly higher at 14.62 and 14.57 emu K mol−1, respectively, likely due to the presence of weak temperature-independent paramagnetism (TIP). In the presence of a stronger applied field, TIP should be diminished, and accordingly χMT decreases to 14.18 and 14.25 emu K mol−1 for 1 and 4, respectively, under a 10 kOe dc field (Fig. S36 and S42†). For all compounds, χMT gradually decreases as the temperature is lowered followed by a precipitous drop at low temperatures. This sudden decrease in χMT is indicative of magnetic blocking, wherein the magnetic moments of individual crystallites in the sample become pinned along a preferred axis and no longer respond to an external magnetic field.
Field-cooled and zero-field cooled measurements collected under an applied magnetic field of 1 kOe were compared for 1–4 (Fig. 2), as the temperature at which these measurements diverge (Tirrev) can provide an estimate of the magnetic blocking temperature.2 Divergence in these curves occurs at 28, 65, 60, and 60 K for 1 through 4, respectively. Such a large range in Tirrev provided an initial indication that the structural changes across this series induce major changes in the magnetic behavior. This result also suggested that compound 2 might exhibit a higher operating temperature than [Dy(Cpttt)2]+, which shows a divergence in field-cooled and zero-field cooled data at 61 K.5 However, Tirrev varies with both the applied magnetic field used to conduct the experiment and the sweep rate,2,21 and thus we turned to measurements of the magnetic relaxation time to more accurately determine the operating temperature.
High-temperature magnetic relaxation in 1–4 was probed by ac magnetic susceptibility measurements from 70 to 114 K (Fig. 3, S43, S46 and S48†). Magnetic relaxation times, τ, were extracted from a simultaneous fit of in-phase (χ′) and out-of-phase (χ″) components of the magnetic susceptibility to a generalized Debye model (Fig. S44, S45, S47 and S49†). Despite the two molecules observed in the crystal structures of 1 and 2, and the positional disorder in 2 and 4, only one peak was observed in the ac susceptibility data for each compound. In this high temperature regime, Arrhenius plots of inverse temperature versus the natural log of τ are linear, indicative of a thermally-activated Orbach relaxation process (Fig. 4, red circles).22 All data could be fit to the equation τ−1 = τ0−1e−Ueff/kT, yielding thermal barriers to magnetic relaxation (Ueff) of 1285, 1468, 1380, and 1334 cm−1 for 1 through 4, respectively. These Ueff values surpass those reported for [Dy(Cpttt)2]+, 1277 or 1223 cm−1,5a,5b and the previous record of Ueff = 1260 cm−1 for [Dy(OtBu)2(C5H5N)5]+.3f
The values of Ueff for 1–4 can be rationalized through structural trends. The magnitude of Ueff should depend on both the metal–ligand bond distances and the geometry of the ligand field (i.e., axial or equatorial). Such correlations have been quantified via ab initio calculations for a hypothetical two-coordinate DyIII complex with alkyl ligands.23 In that study, it was found that U increased linearly with the L–M–L angle at low angles (90 to 150°), tapering its growth at high angles, and flattening above 170°. In contrast, U increased linearly as the M–L bond distance decreased across all bond lengths calculated.
In the series 1–4, the most bent metallocene is 1, with a Cp–Dy–Cp angle of 147.2(8)°. Complex 1 also shows the smallest thermal barrier to magnetic relaxation. Linearity increases by nearly 10° in 2 (156.6(3)°) and even further in 3 (161.1(2)°) and 4 (162.1(7)°), resulting in larger Ueff values for the high-angle set 2–4. Within these three compounds, Ueff should vary less substantially with bond angle, and instead bond distance should have the largest impact on the relaxation barrier. Indeed, Ueff increases upon going from 4 to 2, as the average Dy–C distance for the major disordered component in the crystal structure decreases.
This magneto-structural correlation is significant for the design of future DyIII metallocene cations. While steric bulk in the cyclopentadienyl ligand increases the linearity of the complex in 1–4, it also tends to promote longer Dy–C distances. These factors have opposing effects on the magnitude of Ueff, with the angle having a larger impact at low angles and the bond distance having a larger impact at high angles. Clearly, cyclopentadienyl ligands that can achieve a balance between promoting a large Cp–Dy–Cp angle and short Dy–C bond distances are worth pursuing.
Increasing the Cp–Dy–Cp angle and decreasing the Dy–C bond distances increases the axiality of the ligand field surrounding the DyIII ion in these metallocenium salts, resulting in larger thermal barriers to magnetic relaxation. This is consistent with electrostatic models that have been used to guide the synthesis of lanthanide-based single-molecule magnets.9 Studies on molecules of the type (E = equatorial ligand) have also found that increasing the axiality of the ligand field at DyIII results in larger thermal barriers, with weaker equatorial ligands leading to larger Ueff.24
To probe magnetic relaxation behavior at lower temperatures, dc magnetic relaxation measurements were conducted for 1–4 from 2 to 64 K (Fig. 4, blue circles). Relaxation times were extracted from magnetization versus time plots by fitting the data with a stretched exponential function (for full details see the ESI, Fig. S50–S68 and Tables S3–S6†). The temperature dependence of τ at intermediate temperatures (20 to 64 K) is indicative of Raman relaxation, while at low temperatures the relaxation time tends toward a constant value, indicative of quantum tunneling of the magnetization. The full range of ac magnetic susceptibility and dc magnetic relaxation data could be fit to the equation τ−1 = τ0−1e−Ueff/kT + CTn + τtunnel−1 (Fig. 4 and Table 1).
Hundred-second magnetic blocking temperatures for 1–4 were extracted from the dc magnetic relaxation data, yielding Tb = 17, 62, 59, and 56 K for 1 through 4, respectively. The values for complexes 2–4 all surpass the 100 s blocking temperature of 53 K measured for [Dy(Cpttt)2]+.5 Indeed, the value for 2 represents the highest magnetic blocking temperature yet measured for a single-molecule magnet. While τ = 100 s has been chosen as the benchmark relaxation time for comparing blocking temperatures, 2 shows longer magnetic relaxation times than [Dy(Cpttt)2]+ across the entire temperature range measured (2–114 K).
Ab initio spin dynamics calculations performed on [Dy(Cpttt)2]+ suggest that magnetic relaxation is predominantly moderated by local molecular vibrations.5b In particular, four C–H bending modes of the Cpttt ligand were identified as giving rise to the initial MJ = 15/2 to 13/2 transition that leads to Orbach relaxation. These calculations led to the proposal that changing ring hydrogens in Cpttt to a different substituent could lead to changes in magnetic relaxation time, a hypothesis that has yet to be experimentally verified.5b Substituting the ring hydrogen in 1 with alkyl substituents in 2–4 results in up to a 45 K increase in the blocking temperature. This result clearly illustrates that substitution of the cyclopentadienyl ligands in dysprosium metallocene cations is a powerful tool for tuning magnetic relaxation and Tb.
While changes in ligand substitution modify molecular vibrations, molecular structure and anisotropy are also affected. It is therefore likely that a variety of factors influence changes in Tb. Clearly, spin dynamics calculations are worth pursuing to determine if the contributions of molecular vibrations to magnetic relaxation in this series can be generalized to other molecules in this class. It is worth noting that the 100 s blocking temperature for 1–4 occurs in the regime in which Raman relaxation dominates. While spin dynamics calculations can provide insight into the contribution of molecular vibrations to high-temperature, Orbach relaxation, models have not yet been developed for through-barrier relaxation processes such as Raman relaxation or quantum tunneling of the magnetization.5b As such, developing a more quantitative magneto-structural correlation for the operating temperature in this series will require further development of theories on the impact of molecular vibrations on through-barrier magnetic relaxation. A well-characterized series like 1–4 that varies only slightly in substitution, but substantially in molecular structure and magnetic relaxation could lay the experimental groundwork for such theoretical studies. Insights from these calculations could be significant to the design of new cyclopentadienyl ligands to maximize operating temperatures in DyIII metallocene cations.
The rates of quantum tunneling of the magnetization, τtunnel, extracted from dc relaxation measurements are among the slowest reported for mononuclear single-molecule magnets. This result is consistent with the drastically reduced quantum tunneling of the magnetization for [Dy(Cpttt)2]+.5 While the reduction of quantum tunneling in [Dy(Cpttt)2]+ has been ascribed to the high axiality of the ligand field (which gives rise to nearly pure MJ states in the lowest eight Kramers doublets5a), this axiality cannot fully account for the observed behavior. Other highly anisotropic mononuclear single-molecule magnets, such as [Dy(OtBu)2(C5H5N)5]+, possess nearly pure MJ ground doublets and show fast quantum tunneling of the magnetization.3f
Another possible explanation for the reduced quantum tunneling in [Dy(Cpttt)2]+ is the constrained metal–ligand vibrational modes.5b,25 A lack of low-energy, first-coordination sphere metal–ligand vibrations in [Dy(Cpttt)2]+—such as those created by the weakly bound, equatorial pyridine ligands in [Dy(OtBu)2(C5H5N)5]+—may limit the phonons available for magnetic relaxation at low temperatures. Limiting low-energy phonons has been shown to slow magnetic relaxation at low fields for Ho atoms on an MgO surface,26 as well as for TbPc2 (Pc2− = phthalocyanine dianion) molecules on carbon nanotubes.27 The hypothesis that local molecular coordination is responsible for limiting quantum tunneling in [Dy(Cpttt)2]+ is supported by an increase in the zero-field step in magnetic hysteresis measurements performed on solution samples, which differ from bulk samples in accessible metal–ligand vibrational modes.5b
The trend in τtunnel for 1–4 deviates from that observed for Tb; namely, compounds 1 and 3 show the fastest quantum tunneling of the magnetization, while 2 and 4 show much larger τtunnel values (Table 1). As the rate of quantum tunneling can be influenced by intermolecular dipole interactions, largely dictated by crystal packing,28 magnetically dilute samples of 1–4 were synthesized in order to enable a more accurate comparison of τtunnel.
Dilute samples 1–4@Y were prepared by co-crystallization of 1–4 in a 1:9 ratio with the analogous YIII compound Y1–Y4. Dc magnetic susceptibility measurements of the dilute samples matched those of 1–4 at temperatures above 50 K. Relaxation times extracted from dc magnetic relaxation measurements of 1–4@Y also matched those determined for 1–4 at the highest temperatures measured for each compound. This consistency suggests that the local molecular structure of 1–4—responsible for the high-temperature magnetic relaxation dynamics—is unchanged upon dilution, despite differences in the structures of Y1–Y4. The value of τtunnel, extracted from dc magnetic relaxation measurements of 1–4@Y at 2 K, differs from the undiluted samples, increasing to 1783, 11200, 1125, and 1956 s, respectively. While the values of τtunnel for 1–4@Y are substantially larger than those determined for the undiluted samples, the same overall trend is observed for the series, with the exception that 1@Y shows slower tunneling than 3@Y.
The trend in τtunnel breaks from that observed for the 100 s blocking temperature. While 2@Y shows the slowest quantum tunneling, 1@Y shows a value of τtunnel close to that of 4@Y and larger than that observed for 3@Y. Substitution of the ring hydrogen in 1 with alkyl substituents in 2–4 does not have a consistent influence on the rate of quantum tunneling, as it did on the rate of Raman relaxation, which dictates Tb. This underscores the need to develop a deeper understanding of the role of molecular vibrations in through-barrier magnetic relaxation. Though molecular symmetry has been shown to influence magnetic relaxation in single-molecule magnets,29 there is no apparent magneto-structural correlation between molecular symmetry and quantum tunneling of the magnetization or blocking temperature in the series 1–4, which all possess C1 symmetry.
Magnetic hysteresis data were collected for 1–4 at several temperatures in the range 2–72 K by sweeping the field between 7 and −7 T. For fields below 2 T, a magnetic field sweep rate of 3.1(4) mT s−1 was employed, similar to that used in hysteresis measurements of [Dy(Cpttt)2]+.5 A step is observed at zero field at 2 K, resulting in remanent magnetization values (Mr) of 2.6, 4.4, 3.9, and 4.2 µB for complexes 1 through 4. The magnitude of the remanent magnetization and the coercive field (Hc) for each complex mirrors the trend in τtunnel (Fig. 5a, Table 1).
The hysteresis loops remain open at temperatures as high as 32, 72, 66 and 66 K for 1 through 4, respectively (Fig. 5b). The maximum temperatures in the case of complexes 2–4 are higher than the 60 K maximum determined for [Dy(Cpttt)2]+,5 and this observation is consistent with the trend observed in Tb. Notably, the maximum hysteresis temperature for compound 2 is the highest that has yet been reported for a single-molecule magnet, approaching the temperature of liquid nitrogen (77 K).
Compounds 1–4 provide experimental support for the hypothesis that modifying cyclopentadienyl ring substituents can tune metal–ligand vibrational modes, impacting magnetic relaxation.5b In the series investigated here, substituent modification notably results in a 45 K increase in the magnetic blocking temperature across the series, with complex 2 exhibiting the highest 100 s blocking temperature and hysteresis temperature yet measured for a single molecule magnet. Expanding the magneto-structural correlations presented in this study with ab initio calculations is a clear next step towards developing generalizable design principles for the synthesis of [Dy(CpX)2]+ single molecule magnets that function well above 77 K.
A report of a dysprosium(III) metallocenium salt, [(CpiPr5)(Cp*)Dy][B(C6F5)4], showing a 100 s magnetic blocking temperature of 65 K and open magnetic hysteresis loops as high as 80 K was published on the same date as this work. The magneto-structural correlations outlined here can help to explain the magnetic properties observed in [(CpiPr5)(Cp*)Dy][B(C6F5)4]: the sterically encumbered CpiPr5 ligand promotes a high Cp–Dy–Cp angle of 162.5°, while the less-hindered Cp* allows for short Dy–Cp(centroid) bond distances of 2.296(1) and 2.284(1) Å for the CpiPr5 and Cp* ligand, respectively. The combination of a large Cp–Dy–Cp angle and short Dy–Cp bond distances enables a large thermal barrier to magnetic relaxation of Ueff = 1541 cm−1 and a 100 s magnetic blocking temperature 3 K higher than that observed for 2.30
Footnotes |
† Electronic supplementary information (ESI) available: IR and NMR spectroscopy, synthesis of Y1–Y4, and magnetic characterization. Structures of compounds 1–4 are available under CCDC 1865386–1865389 and compounds Y1–Y4 under CCDC 1873921–1873924. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/c8sc03907k |
‡ These authors contributed equally to this work. |
This journal is © The Royal Society of Chemistry 2018 |