Yuan-Zhu
Zhang‡
ab,
Silvia
Gómez-Coca‡
a,
Andrew J.
Brown
a,
Mohamed R.
Saber
ac,
Xuan
Zhang
a and
Kim R.
Dunbar
*a
aDepartment of Chemistry, Texas A & M University, College Station, TX 77842-3012, USA. Fax: +1-979-845-7177; E-mail: dunbar@chem.tamu.edu
bDepartment of Chemistry, South University of Science and Technology of China, 518055, Shenzhen, P. R. China
cDepartment of Chemistry, Fayoum University, Fayoum 63514, Egypt
First published on 17th June 2016
The air-stable mononuclear Co(II) compounds [CoII(Tpm)2][ClO4]2 (1, Tpm = tris(pyrazol-1-yl)methane), [CoII(Tpm)2][BPh4]2·2MeCN (2) with trigonal antiprismatic geometry (trigonally elongated octahedral geometry) are reported. Magnetic and theoretical studies reveal that the complexes exhibit single-molecule magnet behavior with uniaxial anisotropy and a huge energy difference between ground and first excited Karmers' doublets (∼200 cm−1). Under applied DC fields, compounds 1 and 2 exhibit frequency and temperature dependence of the imaginary susceptibility. The fit of the data to an Orbach relaxation process yields effective energy barriers of 30.6(1) and 44.7(6) cm−1 for 1 and 2, respectively, but there is no real state at that energy. The inclusion of tunneling, direct and Raman relaxation processes leads to the conclusion that the inclusion of an Orbach process is not required to provide a good fit to the data. More interestingly, a detailed study of the dependence of the relaxation time with field shows that for these Kramers' ions, tunneling is the predominant process at low temperature and that differences in the counteranion allow for a tuning of the Raman process at higher temperatures. These findings underscore the fact that large uniaxial anisotropy can be achieved in hexacoordinate Co(II) trigonal antiprismatic complexes which is an unexplored geometry in mononuclear single molecule magnets.
In 2010, Long and co-workers reported SMM behavior for a mononuclear Fe(II) complex in a trigonal pyramidal geometry,7 and later reported that a linear two-coordinate Fe(I) complex exhibits the remarkably high energy barrier of −226 cm−1.8 Subsequently, a number of other mononuclear transition metal SMMs have been isolated, among which are those containing Co(II) ions. For a series of mononuclear tetra-coordinate (tetrahedral or trigonal pyramidal) Co(II) complexes, the anisotropy trends reveal a remarkable range of zero-field splitting parameters from D = −161 to +16 cm−1 depending on the degree of distortion and ligand field contributions.2,9–19 In contrast, Co(II) SMMs with other coordination environments, e.g., tri-,20 penta-,21–23 hexa-,24–33 hepta-34,35 and octa-coordinated,36 are less common.
Given that most Co(II) based single chain magnets37,38 and SMMs39 are composed of octahedral Co(II) building units, which is the most common geometry for this metal ion, the development of additional mononuclear Co(II) SMMs requires a deeper understanding of magneto–structural correlations. Most examples in this vein exhibit large positive D values. For example, Pardo and co-workers first reported the highly distorted octahedral Co(II) compound CoII(dmphen)2(NCS)2 with an overall rhombic (C2V) distortion.24 The compound exhibits strong positive axial magnetic anisotropy with D = +98 cm−1. Later, Trávníček and co-workers demonstrated that a Co(II) center in a compressed octahedral geometry also gives rise to a large positive D value of +48(2) cm−1.30 For these cases, the barriers are assumed to be governed by the transverse anisotropy parameter, E (xy or easy-plane), which is often of a much smaller magnitude than the axial parameter, Dz. Also Luis, Ruiz and co-workers reported a slightly elongated octahedral compound CoII(acac)2(H2O)2, (acac = acetylacetonate), with a D value of +57 and +63.3 cm−1 obtained from the fit of the magnetization curve and CASPT2-RASSI calculations respectively. In this case, the slow relaxation of the spin under an applied field was explained by direct (induced by the hyperfine interaction, I = 7/2) and Raman spin-phonon processes.28
In terms of trigonal molecules, anions such as [(TpR)FeIII(CN)3]− (TpR = tris/tetra(pyrazoyl)borate) and [(triphos)ReII(CN)3]− exhibit significant uniaxial anisotropy40 and have been widely used for the elaboration of SMMs41–50 and SCMs.51–53 Recently, Gao and co-workers reported a family of CoIICoIII3 chiral star-like compounds where the sole Co(II) ion is located in a distorted trigonal prismatic geometry. These authors demonstrated that the magnetic anisotropy and energy barrier can be finely tuned in this family of compounds by the selection of the counter cation and the peripheral ligand substituents, achieving a huge anisotropy (D = −115 cm−1), and an impressive barrier of 109 K under a zero dc field.29,31 Ruiz and co-workers characterized another trigonal prismatic Co(II) complex with D values of −72 and −141 cm−1 determined by a fit of magnetization curve and from CASSCF-RASSI calculations respectively.2 Another interesting case is the work of Novikov and co-workers who reported a trigonal prismatic Co(II) cage complex with relaxation of the spin being observed at zero dc field and D values of −82, −109 and −110 cm−1 determined by a fit of magnetic susceptibility, NMR spectroscopy and from CASSCF/NEVPT2 calculations respectively.32 An ongoing interest in this coordination geometry is apparent from recent literature.32,54
In light of these aforementioned results, we have focused on Co(II) compounds with trigonal antiprismatic geometry that are less sterically rigid but which exhibit a similar splitting of the d-orbitals (Fig. 1a) as the trigonal prismatic geometry which was previously calculated at the DFT level.55 The presence of three electrons in two nearly degenerate eg orbitals is expected to lead to a large anisotropy, in this case it produces a uniaxial magnetic anisotropy as confirmed by EPR analysis 50 years ago.56 Thus far, however, no SMMs have been reported with this geometry. We postulated that fac-tripodal ligands could enforce a trigonal antiprismatic geometry, thus we undertook a systematic study of Co(II) complexes equipped with tris(pyrazol-1-yl)methane (Fig. 1b), a ligand that has been used in spin crossover compounds.57–65 Herein, we report two mononuclear SMMs, viz., [CoII(Tpm)2][ClO4]2 (1) and [CoII(Tpm)2][BPh4]2·2MeCN (2), which have the same magnetic unit, [CoII(Tpm)2]2+, but different chemical surroundings. Magnetic studies indicate that both compounds possess a large energy difference between the ground and excited states (large energy barrier) and field-induced relaxation of the magnetization.
Compound | 1 | 2 |
---|---|---|
a I ≥ 2σ(I): R1 = ∑ ∥Fo| − |Fc∥ / ∑ |Fo|, wR2 = {∑ [w(Fo2 − Fc2)2] / ∑[w(Fo2)2]}1/2. | ||
Crystal system | Monoclinic | Monoclinic |
Space group | P21/n | P21/n |
a, Å | 10.500(8) | 9.7482(8) |
b, Å | 7.818(6) | 17.1636(14) |
c, Å | 16.997(13) | 18.9136(16) |
β, deg | 103.529() | 98.5041() |
V, Å3 | 1356.5(18) | 3129.7(4) |
D c, g cm−3 | 1.680 | 1.282 |
Z | 2 | 2 |
μ, mm−1 | 0.900 | 0.330 |
F(000) | 698 | 1266 |
GooF | 1.078 | 1.039 |
R 1 | 0.0602 | 0.0323 |
wR2a | 0.1635 | 0.0762 |
Compound | 1 | 2 |
---|---|---|
Symmetry transformations used to generate equivalent atoms: A: −x + 2,−y + 2,−z. | ||
Co1–N2 | 2.112(3) | 2.110(1) |
Co1–N4 | 2.107(3) | 2.104(1) |
Co1–N6 | 2.107(3) | 2.101(1) |
N2–N4 | 2.815(9) | 2.861(1) |
N2–N6 | 2.845(5) | 2.842(1) |
N4–N6 | 2.848(2) | 2.814(1) |
N2–N4A | 3.144(4) | 3.094(1) |
N2–N6A | 3.114(3) | 3.107(1) |
N4–N6A | 3.110(10) | 3.125(1) |
N2–Co1–N4 | 83.68(13) | 85.53(4) |
N2–Co1–N6 | 84.89(14) | 84.90(4) |
N4–Co1–N6 | 84.91(13) | 83.99(4) |
Bite angle avg. | 84.49(13) | 84.81(4) |
N2–Co1–N2A | 180.0 | 180.0 |
N2–Co1–N4A | 96.32(13) | 94.47(4) |
N2–Co1–N6A | 95.11(14) | 95.10(4) |
N4–Co1–N6A | 95.09(13) | 96.10(4) |
Co1⋯C | 3.117(8) | 3.113(1) |
Co⋯Co | 7.818(6) | 9.748(1) |
Fig. 2 The molecular structures of 1 and 2 (thermal ellipsodes are at 30% level and all counter anions and hydrogen atoms have been omitted for clarity). |
The nearest Co⋯Co distance in 2 (9.748(1) Å) is significantly longer than that in 1 (7.818(6) Å), due to the larger [BPh4]− anions (Fig. S1 and S2, ESI†). Two different orientations of the molecules in the crystal packing (the pseudo-C3 axes along the Co1–C10 vectors) are present with an angle of 40.8° for 1 and 69.8° for 2, respectively.
Fig. 3 Variable-temperature dc magnetic susceptibility data in an applied field of 1 kOe for 1 (top) and 2 (bottom). Insets: M vs. H plots at 2, 4, 6 K. Solid lines are the best simulations obtained by the program PHI based on eqn (1). |
To analyze the experimental data, the PHI program was employed using the Hamiltonian in eqn (1).73 This Hamiltonian is divided into three terms, the first two being related to the crystal field Hamiltonian (using the Stevens notation)74 with the third one representing the Zeeman Hamiltonian with an anisotropic g tensor.
(1) |
To further verify the obtained values, ab initio calculations using the experimental geometry and the CASSCF/NEVPT2 approach with the Orca code were performed.67 The calculated energy of the first spin-free states is shown in Table S15.† The energy difference between the ground and first excited states is very small (73.9 and 53.1 cm−1, for 1 and 2 respectively, at NEVPT2 level) while the second excited states is almost 3000 cm−1 above the ground state. These findings are in accord with the first two states being nearly degenerate indicating that the Jahn–Teller effect is small and the first order SOC should be important.
The SOC mixing of the ground and excited states gives rise to four Kramers' doublets (KD) relatively close in energy (Table 3). For a Fe(II) system with 5E ground state in which the distortions can change the energy and lift the orbital degeneracy, a more complicated Spin Hamiltonian (SH) with a larger number of parameters has been proposed.75 This point notwithstanding, the effective SH employed in the experimental part gives us a qualitative estimation of the large energy difference between states (|2D| equal to 184 and 186 cm−1 for 1 and 2 respectively), which is corroborated by the theoretical calculations (232.8 cm−1 for 1 and 2). Although this result does not provide information about the proximity in energy of the third and fourth KD, from the theoretical calculations one can conclude that the population of these two KD at 300 K is smaller than 5% (Table S11†) and that the effective Hamiltonian employed in the experimental part provides a reasonable estimation. Using the effective Hamiltonian implemented in Orca, the D and E values match very well with the experimental ones (Table S12†). The analysis of the principal contribution to the D value revealed that the origin is from the first excited state while the major contribution to E emanates from the second excited state (Table S13†).
KD | 1 | 2 | ||||||
---|---|---|---|---|---|---|---|---|
E | g x | g y | g z | E | g x | g y | g z | |
1 | 0 | 0.71 | 0.73 | 8.93 | 0.0 | 0.72 | 0.74 | 8.92 |
2 | 232.8 | 0.68 | 1.50 | 4.78 | 236.3 | 0.38 | 1.21 | 4.81 |
3 | 509.1 | 0.86 | 1.02 | 1.21 | 504.2 | 0.74 | 0.85 | 0.95 |
4 | 809.8 | 0.06 | 0.08 | 3.05 | 817.2 | 0.06 | 0.08 | 3.06 |
5 | 3290 | 0.13 | 0.14 | 6.38 | 2592 | 0.07 | 0.07 | 6.37 |
6 | 3377 | 2.12 | 3.58 | 3.85 | 2675 | 2.12 | 3.64 | 3.77 |
In addition, the effective gi values (Table 3) confirm the axial character of the ground state and agree nicely with the gi values reported 50 years ago for the analogous compound [Co(Tpm)2](NO3)2.56
For compound 1, variable-frequency ac data were measured at 1.8 K under applied dc fields in the range of 0 to 5000 Oe. No χ′′ signal was observed at zero dc field, while a peak appeared at ∼22 Hz when a 250 Oe was applied. When a larger dc field is applied, this peak increases until 750 Oe and then decreases at higher fields. At 750 Oe a new signal appears at lower frequencies, an indication of several relaxation processes being operative at lower temperatures (Fig. S4†). Due to the superposition of the peaks and the number of points for one of the relaxation processes, not all of the data could be fit. Also, for fields larger than 1500 Oe, the α value is ∼0.3 indicating a large distribution of relaxation times (Fig. S5 and Table S1†), which may be related to the closer distance between molecules for 1 (7.818 Å). These values were discarded for the subsequent analysis.
Under a 500 Oe dc field, variable-temperature (2.0–15.0 K) ac susceptibilities exhibit highly frequency-dependent peaks in both the in-phase and out-of phase signals (Fig. S6†). At 100 Hz, there is an indication of two peaks, the second of which has a maximum at ∼5 K. To gain more insight into the field dependence, variable-frequency ac data were measured under applied dc fields in the range of 0 to 10000 Oe at 5.0 K (Fig. S7†). From 125 to 625 Oe, the maximum in χ′′ is nearly constant with a slight shift to higher frequencies. From 625 to 3000 Oe, when the field is increasing, the maximum in χ′′ moves to lower frequencies. For larger fields, the maximum in χ′′ moves to higher frequencies. The relaxation times obtained from the fittings of the Cole–Cole plots with a modified Debye model were plotted as a function of the dc fields (Fig. S8, S9 and Table S2†). The value of τ−1 first increases up to 625 Oe, then decreases from 625 to 3000 Oe and finally increases for fields larger than 3000 Oe. The decrease up to 3000 Oe is ascribed to a tunneling process and the increase at higher fields is attributed to a direct process due to its dependence with field.
The variable-frequency ac data at different temperatures were collected under 500 and 3000 Oe (Fig. S8† and 4, respectively). At 500 Oe the maximum in χ′′ moves to higher frequencies with increasing temperature. The increase is very slight until 4.5 K and then becomes more pronounced at higher temperatures. Two processes are observed at 3000 Oe up to 3.5 K, but at higher temperatures only one highly frequency dependent peak is observed. The Cole–Cole plots were fitted to a modified Debye model with α values less than 0.1 at 500 Oe and at 3000 Oe above 4.0 K, indicative of a relatively narrow distribution of relaxation pathways (Fig. S11, S12, Tables S3 and S4†).
For compound 2, similar dynamic susceptibility behavior was observed. Variable-frequency ac data were measured at 1.8 K under applied dc fields in the range of 0 to 1000 Oe (Fig. S13†). Only a low frequency signal without an obvious maximum was observed. Under a 500 Oe dc field, variable-temperature (2.0–15.0 K) ac susceptibility data showed highly frequency-dependent peaks in both the in-phase and out-of phase signals (Fig. S14†). At 100 Hz, a maximum is observed at 6 K. To gain more insight into the field dependence, the variable-frequency ac data were measured under applied dc fields in the range of 0 to 10000 Oe at 6.0 K (Fig. S15†). The relaxation times obtained from the fittings of the Cole–Cole plots with a modified Debye model were plotted as a function of the dc fields (Fig. S16, S17 and Table S5†). It was found that the relaxation times increase up to a dc field of 1500 Oe and then decrease at higher dc fields. The first decrease of τ−1 can be attributed to tunneling whereas the subsequent increase is in accord with a direct process.
The variable-frequency ac data at different temperatures were collected under fields of 300, 500 and 1500 Oe (Fig. S18, S20† and 3, respectively). In this case, at the optimum field of 1500 Oe, only one highly frequency dependent peak is observed over the range of temperatures measured. The Cole–Cole plots were fit to a modified Debye model which gave small α values of less than 0.1 at above 4.0 K, indicative of a relatively narrow distribution of relaxation pathways (Fig. S19, S21, S22 and Tables S6–S8†).
For both complexes, the dependence of τ−1 with the temperature at different fields was analyzed. A fit with the Arrhenius law (τ = τ0 exp(Ueff/kBT)), corresponding to the thermally activated relaxation via an Orbach process (Fig. S23†) was performed leading to effective energy barriers (Ueff) and pre-exponential factors (τ0) of 30.6(7)–33.6(5) cm−1/2.0(2)–3.3(7) × 10−7 s for 1; 42.5(6)–44.7(6) cm−1/1.0(1)–1.5(2) × 10−7 s for 2, respectively. These values, however, are far from the energy difference between ground and excited states (∼200 cm−1). This discrepancy and the absence of a real state necessary for an Orbach process at 30–40 cm−1 are clear indications of the presence of other relaxation processes.
Previously the spin relaxation of other mononuclear SMMs has been fitted considering direct, Raman, tunneling and/or Orbach processes. For linear Fe(II) complexes, the dependence of τ with field was fitted to obtain direct and tunneling parameters. The dependence of τ with temperature was fitted by considering direct, Raman, tunneling and Orbach processes but by fixing the previously obtained parameters.79 A similar procedure was employed for a trigonal prismatic Co(II) cage complex.32 For a previously reported octahedral Co(II) mononuclear SMM with a positive D value, the spin relaxation was explained by including the hyperfine coupling and the nuclear spin-lattice interaction, which allows for the one phonon (or direct) relaxation process.28 In that case, the variable temperature dynamic susceptibility of a magnetically diluted sample was fitted with a combination of direct and Raman (two phonon) processes. In a similar way, direct and Raman processes have been considered in order to fit the spin relaxation of other 3d compounds.80,81 For a recently reported four-coordinate Co(II) compound the spin relaxation at zero dc field was fitted considering Raman and Orbach; in this case they fixed the thermal energy barrier to the value obtained from the far infrared transmission spectra to be 230 cm−1 and then fit the rest with the Raman contribution.82 In other cases the data have been fit only with a Raman term when the energy barrier did not match the energy difference between states.26,34,83–86 In another system, the dependence of τ was fit with Orbach, Raman and direct terms although in one of these cases the Raman term was found to be zero.87 while in the other the direct was zero.35
In the current study, the relaxation times were re-evaluated considering direct, Raman, tunneling and Orbach processes. For both compounds, all attempts to fit the dependence of τ−1 with field according to eqn (S1)† (which considers direct, tunneling and a constant to include processes that are independent with field) were not successful. Clearly there is a more complex dependence of τ−1 with field. A fit of the dependence of τ−1 with field including all the possible relaxation pathways (eqn (S2), Fig. S24 and Table S9†) would result in an over-parameterization of the curves, especially without fixing the direct and tunneling parameters. We, therefore, adopted another approach.
The more common relaxation processes for mononuclear compounds are direct, tunneling (usually dominant at low temperatures), Raman and Orbach (usually predominant at higher temperatures).88 First we analyzed the high temperature region. The large energy difference between states should lead to a very slow relaxation time (slower than the measurable relaxation times with the ac measurements that are possible with our SQUID). If the energy difference between the states is larger than the Debye temperature, the Raman term will be the principal one at high temperature with a negligible Orbach process.88 In fact, an attempt to fit the high temperature regime with both processes with an energy barrier of 200 cm−1 as a starting parameter, led to the discovery that the value does not vary and the data can be fit exclusively with the Raman term. For other mononuclear SMMs, it has been observed that both processes coexist even for energy barriers larger than 200 cm−1. In our case, however, if we examine the contributions of Orbach and Raman terms we see that the Raman is the one which describes the shape of the curve (Fig. S25†).
Secondly, we analyzed the lower temperature region. For both compounds, τ−1 is larger at lower fields; the values at 500 Oe vs. 3000 Oe for 1, and 300, 500 and 1500 Oe for 2 were compared. These results indicate that the predominant process is tunneling because the opposite trend is expected for a direct process (τ−1 is proportional to H4 for Kramers' ions).88 Taking this aspect into consideration and to avoid the over-parameterization of the curves, we modelled the dependence of τ−1 with temperature using eqn (2) and S3,† which considers just tunneling and Raman terms.
τ−1 = B + CTn | (2) |
Fig. 5 shows the best fit using eqn (2) with constraining the B parameter to the minimum value of τ−1. The obtained parameters are B = 184 and 0 s−1 at 500 and 3000 Oe respectively, C = 0.0117 s−1 K−n and n = 6.65 for 1 and B = 17.3, 9.9 and 0.9 s−1 at 300, 500 and 1500 Oe respectively, C = 0.0011 s−1 K−n and n is 7.18 for 2. From these parameters, one can deduce that, for an optimal field, the tunneling is essentially negligible and that the exponent for the Raman term in both cases is very close to 7. For Kramers ions, this exponent should be 9 (or 5 in the presence of low-lying states),74 but different values have been reported and justified by the presence of both acoustic and optical phonons.26,30,32,82 The Raman coefficient for 1 is one order of magnitude larger than for 2, a reflection of how a change in counterion modifies the spin-lattice Raman process.
Fig. 5 Dependence of τ−1 with T for 1 and 2 under different dc fields. Solid lines are the best simulation of the curves using eqn (2). |
As revealed by the experimental magnetic data and the theoretical calculations, the anisotropy for both compounds is very similar so the differences in the dynamic susceptibility behavior between 1 and 2 are attributed to the different counter anion. When we move from a smaller anion in 1 (ClO4−) to a larger one in 2 (BPh4−) there is an (i) increase in the distance between Co(II) centers which reduces the dipolar interactions, and (ii) a change in the crystal packing and chemical environment of the SMM unit.
In the case of 1 there is, in general, a larger the distribution of relaxation times (α) and two different processes were observed at lower temperatures and different dc fields, most likely due to the dipolar interactions because of the closest contact being 7.818 Å. In both of the current complexes we have seen from the dependence with field that tunneling is the predominant process at low temperature (whereas for other Co(II) mononuclear SMMs the low temperature behavior is attributed to direct). In addition, the change in the crystal packing and chemical environment due to the larger anion in 2 leads to a decrease in the Raman parameter (C). A deeper understanding of this process will allow us to tune this effect with the appropriate election of the chemical environment. It is not unreasonable to expect that further modifications of this geometry may result in a compound with an extremely high barrier as observed for other geometries.
Footnotes |
† Electronic supplementary information (ESI) available: Packing diagrams, detailed static and dynamic magnetic measurements, analysis of dynamic magnetic measurements and computational details. CCDC 1422285 and 1422286. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/c6sc02035f |
‡ Authors have contributed equally. |
This journal is © The Royal Society of Chemistry 2016 |