Multiple Magnetization Reversal Channels Observed in a 3 d-4 f Single Molecule Magnet

The present study discusses the magnetic dynamics of a previously reported cyanide bridged 3d-4f dinuclear DyIIICoIII complex. Following the axial anisotropy suggested by previous Electron Paramagnetic Resonance spectroscopy (EPR) analysis, the complex turned out to show slow relaxation of the magnetization at cryogenic temperature, and this was studied in different temperature and field regimes. The existence of multichannel relaxation pathways that reverse the magnetization was clearly disclosed: a tentative analysis suggested that these channels can be triggered and controlled as a function of applied static magnetic field and temperature. Persistent evidence of a temperature independent process even at higher fields, attributable to quantum tunneling, is discussed, while the temperature dependent dynamics is apparently governed by an Orbach process. The broad distribution of relaxation rates evidenced by the ac susceptibility measurements suggest a relevant role of the intermolecular interactions in this system.


Introduction
Many interdisciplinary efforts are underway to comprehend the magnetic behavior of lanthanide based molecular magnets (MMs); however, this has proven to be a complex and tangled task owing to the large orbital contribution of lanthanides (large magnetic anisotropy), along with weak exchange interactions and significant crystal field effects normally occurring in such systems [1][2][3][4][5].The objective of these dedicated attempts is to individually identify and explain the contribution of the aforementioned processes to the static and dynamic magnetic properties of the lanthanide based molecules.Many of these attempts have been performed in mixed d-f metal ion complexes, as a way to bypass the weakness of the f-f magnetic exchange interactions.Interestingly, these mixed d-f complexes allow the use of the so-called diamagnetic substitution approach, where the transition paramagnetic ion exchange-coupled to the lanthanide ion is substituted with a diamagnetic analogue, at the same time keeping the integrity of the crystal structure.This particular approach provides a qualitative picture of the relative magnitude of the active interactions in the system [6,7] which was particularly precious in the absence of reliable theoretical methods to evaluate them.
On the other hand, more recently, a lot of success has been achieved in understanding the key factors affecting the processes leading to slow relaxation of the magnetization in lanthanide based MMs-much more ubiquitous than initially thought-by use of multiple techniques, including magnetometry, susceptibility, Electron Paramagnetic Resonance spectroscopy (EPR), luminescence and magnetic circular dichroism [8][9][10][11], often coupled to detailed ab initio calculations.This provided breakthroughs in the understanding of the role of the magnetic anisotropy and the crystal field effects in determining high barriers to the relaxation of magnetization, and provided key inputs for the design and synthesis of compounds showing relatively high magnetization blocking temperatures [12][13][14].In this framework, it is interesting to note that most of the systems earlier synthesized for the study of their static magnetic properties have not yet been characterized for the dynamics of their magnetization, despite the few exceptions that have provided clear evidence of their potentialities in this regard [15,16].
The present work is focused on a Dy(H 2 O) 3 (dmf) 4 (µ-CN)Co-(CN) 5 ]¨nH 2 O where 1 ď n ď 1.5 and dmf = N,N 1 -dimethylformamide, complex (hereafter DyCo) which had been previously synthesized by one of us [17] to study the competing effect of exchange and anisotropy effects on the static magnetic behavior in a family of isotructural LnM complexes.We underline that, for this derivative, only the crystal field affects the temperature dependence of the susceptibility, given the diamagnetic nature of the low-spin Co III , only slightly affected by temperature independent paramagnetism [17].The complex has a monoclinic unit cell containing four discrete hetero-dinuclear molecules and crystallizes in the P2 1 /n (i.e., n ˝14) space group [7,17].A detailed view of the coordination sphere of the Dy III center is reported as reference in Figure 1: Dy III coordination number is eight, and analysis of the coordination polyhedron via use of SHAPE [18] provided evidence of a geometry intermediate between square antiprismatic, bicapped trigonal prismatic and trigonal dodecahedron geometry (see Table S1).The analysis of the packing, as well as detailed X-ray investigation on isostructural systems, [19] showed the existence of a 3D network made up of Dy III -Co III dinuclear entities where the five terminal CN groups of the {Co(CN) 6 } fragment are connected through hydrogen bonds to the Ln-coordinated water molecules via crystallization water molecules.The static magnetic properties revealed that this network of hydrogen bonding interaction might result in non-negligible inter-dimer coupling interactions [17].
and magnetic circular dichroism [8][9][10][11], often coupled to detailed ab initio calculations.This provided breakthroughs in the understanding of the role of the magnetic anisotropy and the crystal field effects in determining high barriers to the relaxation of magnetization, and provided key inputs for the design and synthesis of compounds showing relatively high magnetization blocking temperatures [12][13][14].In this framework, it is interesting to note that most of the systems earlier synthesized for the study of their static magnetic properties have not yet been characterized for the dynamics of their magnetization, despite the few exceptions that have provided clear evidence of their potentialities in this regard [15,16].
The present work is focused on a Dy(H2O)3(dmf)4(μ-CN)Co-(CN)5]•nH2O where 1 ≤ n ≤ 1.5 and dmf = N,N′-dimethylformamide, complex (hereafter DyCo) which had been previously synthesized by one of us [17] to study the competing effect of exchange and anisotropy effects on the static magnetic behavior in a family of isotructural LnM complexes.We underline that, for this derivative, only the crystal field affects the temperature dependence of the susceptibility, given the diamagnetic nature of the low-spin Co III , only slightly affected by temperature independent paramagnetism [17].The complex has a monoclinic unit cell containing four discrete hetero-dinuclear molecules and crystallizes in the P21/n (i.e., n° 14) space group [7,17].A detailed view of the coordination sphere of the Dy III center is reported as reference in Figure 1: Dy III coordination number is eight, and analysis of the coordination polyhedron via use of SHAPE [18] provided evidence of a geometry intermediate between square antiprismatic, bicapped trigonal prismatic and trigonal dodecahedron geometry (see Table S1).The analysis of the packing, as well as detailed X-ray investigation on isostructural systems, [19] showed the existence of a 3D network made up of Dy III -Co III dinuclear entities where the five terminal CN groups of the {Co(CN) 6 } fragment are connected through hydrogen bonds to the Lncoordinated water molecules via crystallization water molecules.The static magnetic properties revealed that this network of hydrogen bonding interaction might result in non-negligible interdimer coupling interactions [17].EPR spectroscopy on a powder sample, performed at 4.2 K was interpreted under the assumption that the observed properties were only due to the anisotropic ground doublet of Dy III [20], and suggested that the DyCo boasts an almost axial anisotropy with an effective g tensor (g eff x, g eff y, g eff z) = (0.1(1), 0.1(1), 18(1)) [21].Since axial anisotropy is often accompanied, in Dy-containing systems [22][23][24], by slow relaxation of the magnetization, this prompted us to study the magnetic dynamics in this complex, which is herein reported as a function of frequency, static magnetic field and temperature.Interestingly, multiple relaxation pathways triggered by wide range of temperature and fields were observed.The different contributions to the slow relaxation of magnetization were then investigated by modeling the experimental data assuming various relaxation mechanisms.Interestingly it was observed that a temperature independent contribution, usually attributed to quantum tunneling of magnetization and thus often ineffective in applied field, survived up to quite a high field.EPR spectroscopy on a powder sample, performed at 4.2 K was interpreted under the assumption that the observed properties were only due to the anisotropic ground doublet of Dy III [20], and suggested that the DyCo boasts an almost axial anisotropy with an effective g tensor (g eff x , g eff y , g eff z ) = (0.1(1), 0.1(1), 18(1)) [21].Since axial anisotropy is often accompanied, in Dy-containing systems [22][23][24], by slow relaxation of the magnetization, this prompted us to study the magnetic dynamics in this complex, which is herein reported as a function of frequency, static magnetic field and temperature.Interestingly, multiple relaxation pathways triggered by wide range of temperature and fields were observed.The different contributions to the slow relaxation of magnetization were then investigated by modeling the experimental data assuming various relaxation mechanisms.Interestingly it was observed that a temperature independent contribution, usually attributed to quantum tunneling of magnetization and thus often ineffective in applied field, survived up to quite a high field.

Low Field and Temperature Dynamics
Ac susceptibility as a function of frequency was first performed at the lowest experimentally attainable temperature of 2 K, in the presence of a variable static magnetic field.No out-of-phase component of the susceptibility (χ 11 ) is observed in zero applied magnetic field (Figure 2a), even if it is quite evident (Figure S1) that a non-zero signal might be present for frequencies above 10 KHz (i.e., outside available experimental range).On application of a static field of 100 Oe a non-zero χ 11 signal with a clear frequency-dependence is observed, with a maximum around 50 Hz, indicating slow relaxation of magnetization (Figure S1).This behavior is quite common for lanthanide based systems, as a consequence of suppression of quantum tunneling of magnetization (QTM) by an applied magnetic field [25].Figure 2a shows the behavior of χ 11 (ν) as a function of applied static magnetic field up to 5 KOe; it is evident that on increasing field, the magnetization is relaxing slower and at the same time the relaxation times distribution narrows.A more quantitative evaluation of the field dependence of the relaxation rate and its distribution was achieved by fitting the data using the generalized Debye equation [11]: where τ is the magnetic relaxation time, χ T and χ S are the isothermal and adiabatic susceptibility, respectively, while α is a parameter describing the distribution of relaxation times.It is interesting to note that, before the relaxation becomes almost field independent at high fields, a decrease followed by a rise in the relaxation rate is observed from 1.25 to 2.25 KOe, Figure 2b.The field dependence of the α parameter obtained by the Debye fit reveals, as expected, a broad distribution of relaxation times (0.27-0.39) which is essentially field independent above 1 KOe (Figure S2).The observed field dependence of the relaxation rate up to 2 KOe is consistent with the rate being dominated by quantum tunneling processes in zero field, which are then quenched by the application of a field: this in turn increases the efficiency of direct process, resulting in a decrease of τ at higher fields.Following this interpretation, a phenomenological fit of the data was attempted in this regime using the expression [26]: which provided best fit parameters B 1 = 1275 ˘83 Hz, B 2 = (8.8˘0.7) ˆ10 ´5 Oe ´2, A 1 = (5 ˘0.43) 10 ´13 Hz(K) ´1 Oe ´4 and the field independent contribution (either Orbach or Raman), C = 12.5 Hz.

Low Field and Temperature Dynamics
Ac susceptibility as a function of frequency was first performed at the lowest experimentally attainable temperature of 2 K, in the presence of a variable static magnetic field.No out-of-phase component of the susceptibility (χ′′) is observed in zero applied magnetic field (Figure 2a), even if it is quite evident (Figure S1) that a non-zero signal might be present for frequencies above 10 KHz (i.e., outside available experimental range).On application of a static field of 100 Oe a non-zero χ′′ signal with a clear frequency-dependence is observed, with a maximum around 50 Hz, indicating slow relaxation of magnetization (Figure S1).This behavior is quite common for lanthanide based systems, as a consequence of suppression of quantum tunneling of magnetization (QTM) by an applied magnetic field [25].Figure 2a shows the behavior of χ′′(ν) as a function of applied static magnetic field up to 5 KOe; it is evident that on increasing field, the magnetization is relaxing slower and at the same time the relaxation times distribution narrows.A more quantitative evaluation of the field dependence of the relaxation rate and its distribution was achieved by fitting the data using the generalized Debye equation [11]: where τ is the magnetic relaxation time, χT and χS are the isothermal and adiabatic susceptibility, respectively, while α is a parameter describing the distribution of relaxation times.It is interesting to note that, before the relaxation becomes almost field independent at high fields, a decrease followed by a rise in the relaxation rate is observed from 1.25 to 2.25 KOe, Figure 2b.The field dependence of the α parameter obtained by the Debye fit reveals, as expected, a broad distribution of relaxation times (0.27-0.39) which is essentially field independent above 1 KOe (Figure S2).The observed field dependence of the relaxation rate up to 2 KOe is consistent with the rate being dominated by quantum tunneling processes in zero field, which are then quenched by the application of a field: this in turn increases the efficiency of direct process, resulting in a decrease of τ at higher fields.Following this interpretation, a phenomenological fit of the data was attempted in this regime using the expression [26]:  Existence of multiple processes contributing to relaxation of magnetization in rare-earth based complexes is consistent with the findings in the literature [23,27,28].It is however to be noted that, in Existence of multiple processes contributing to relaxation of magnetization in rare-earth based complexes is consistent with the findings in the literature [23,27,28].It is however to be noted that, in the present case, the data clearly indicates that at fields higher than 2 KOe the direct process stops increasing with the field; alternatively, one may assume an exact cancellation between processes the rate of which increases with the field (direct process) and those which should slow down on increasing field (quantum tunneling).Since a small static field was sufficient to trigger slow relaxation of magnetization in this complex, a detailed temperature study was performed in the presence of a static field of 300 Oe to further probe the contribution of different processes to the dynamics.The temperature dependence of the frequency at which the maximum in the imaginary susceptibility component occurs (Figure S3) suggests a single molecular magnetic behavior, as expected for complexes with an axial anisotropic center [29] as is the case here on the basis of previously reported EPR spectroscopy [21].
It is however evident, by a more accurate analysis of the temperature dependence of χ(ν), that at this field two different relaxation phenomena are actually active (Figure S4).At 2 K indication of one fast (>1.6 ˆ10 ´5 s) and one slow (8 ˆ10 ´3 s) process can be clearly seen: on increasing temperature the fast process rapidly moves beyond the upper limit of the instrumental frequency range, and simultaneously the slow process become faster and eventually shifts to frequencies out of experimental range at 4 K.The relaxation time of the slower process, derived by the fit to the Debye Equation ( 1), obeys the Arrhenius law at high temperature.Deviation from linear behavior is evident below 2.2 K, where the relaxation tends toward a temperature independent regime (Figure 3), suggesting the existence of a quantum regime.The respective relaxation time as function of temperature was then modeled using a combination of Arrhenius and quantum tunneling: This provided as best fit parameters, τ 0 = 2.38 ˆ10 ´9 s and an effective barrier U eff = 34 K (23.87 ˘0.9 cm ´1), with a temperature independent contribution D of 72.36 ˘10 Hz which is of the correct order of magnitude compared to what predicted by the fit of field dependent data.Interestingly, a tentative fit performed assuming a combination of Raman (CT n ) and a temperature independent contribution did not provide a meaningful result, the obtained n exponent being far larger than expected (11.8).
Magnetochemistry 2016, 2, 27 4 of 10 the present case, the data clearly indicates that at fields higher than 2 KOe the direct process stops increasing with the field; alternatively, one may assume an exact cancellation between processes the rate of which increases with the field (direct process) and those which should slow down on increasing field (quantum tunneling).Since a small static field was sufficient to trigger slow relaxation of magnetization in this complex, a detailed temperature study was performed in the presence of a static field of 300 Oe to further probe the contribution of different processes to the dynamics.The temperature dependence of the frequency at which the maximum in the imaginary susceptibility component occurs (Figure S3) suggests a single molecular magnetic behavior, as expected for complexes with an axial anisotropic center [29] as is the case here on the basis of previously reported EPR spectroscopy [21].
It is however evident, by a more accurate analysis of the temperature dependence of χ(ν), that at this field two different relaxation phenomena are actually active (Figure S4).At 2 K indication of one fast (>1.6 × 10 −5 s) and one slow (8 × 10 −3 s) process can be clearly seen: on increasing temperature the fast process rapidly moves beyond the upper limit of the instrumental frequency range, and simultaneously the slow process become faster and eventually shifts to frequencies out of experimental range at 4 K.The relaxation time of the slower process, derived by the fit to the Debye Equation ( 1), obeys the Arrhenius law at high temperature.Deviation from linear behavior is evident below 2.2 K, where the relaxation tends toward a temperature independent regime (Figure 3), suggesting the existence of a quantum regime.The respective relaxation time as function of temperature was then modeled using a combination of Arrhenius and quantum tunneling: (3) This provided as best fit parameters, τ0 = 2.38 × 10 −9 s and an effective barrier Ueff = 34 K (23.87 ± 0.9 cm −1 ), with a temperature independent contribution D of 72.36 ± 10 Hz which is of the correct order of magnitude compared to what predicted by the fit of field dependent data.Interestingly, a tentative fit performed assuming a combination of Raman (CT n ) and a temperature independent contribution did not provide a meaningful result, the obtained n exponent being far larger than expected (11.8).

Higher Field Dynamics
Although relaxation of magnetization was observed to be slowest at about 1.25 KOe, a detailed temperature study to observe the evolution of magnetic behavior was performed at higher static applied magnetic field of 2 KOe, i.e., in a regime where competing contributions of various relaxation phenomena were observed at low temperature (Figure 2b).The motivation for this choice lies in

Higher Field Dynamics
Although relaxation of magnetization was observed to be slowest at about 1.25 KOe, a detailed temperature study to observe the evolution of magnetic behavior was performed at higher static applied magnetic field of 2 KOe, i.e., in a regime where competing contributions of various relaxation phenomena were observed at low temperature (Figure 2b).The motivation for this choice lies in exploring the possibility of more easily identifying the multiple relaxation pathways in this compound, already identified in the analysis at 300 Oe.As expected, the data at 2 KOe reveal strong frequency dependence over the whole temperature range (Figure 4a).For temperatures below 2.7 K a clear maximum (identified as 1 in figure) exists at low frequency; starting from 2.7 K a broadening of χ 11 (ν) is observed, suggesting the presence of a secondary relaxation (identified as 2 in figure).The second observation is confirmed on increasing temperature by the clearer definition and gradual shift of maximum 2. At intermediate temperatures 3.4-4 K signature of both the peaks is quite clear (Figure 4a).On further increasing temperature up to 8 K both the maxima completely disappeared and the magnitude of out-of-phase susceptibility decreased to almost zero at 8 K, indicating fast magnetic dynamics in these experimental conditions (Figure S5).exploring the possibility of more easily identifying the multiple relaxation pathways in this compound, already identified in the analysis at 300 Oe.As expected, the data at 2 KOe reveal strong frequency dependence over the whole temperature range (Figure 4a).For temperatures below 2.7 K a clear maximum (identified as 1 in figure) exists at low frequency; starting from 2.7 K a broadening of χ′′(ν) is observed, suggesting the presence of a secondary relaxation (identified as 2 in figure).The second observation is confirmed on increasing temperature by the clearer definition and gradual shift of maximum 2. At intermediate temperatures 3.4-4 K signature of both the peaks is quite clear (Figure 4a).On further increasing temperature up to 8 K both the maxima completely disappeared and the magnitude of out-of-phase susceptibility decreased to almost zero at 8 K, indicating fast magnetic dynamics in these experimental conditions (Figure S5).The temperature dependences of the relaxation rates obtained by the fit of χ(ν) are reported as Arrhenius plot in Figure 4b: it is evident that for T > 2.5 K the behavior is similar for the two phenomena, being essentially linear.However, the slope is different, indicating that the anisotropy barrier is different for the two processes.Finally, at low temperature, where only process 1 is visible, a deviation from the linear behavior is observed, indicating a quantum regime.Keeping the above observations in mind, the data was tentatively analyzed using Equation (3), providing the following best fit parameters: τ0(1) = (6.2± 1.5) × 10 −10 s, Ueff(1) = (39.9± 0.8) K, τ0(2) = (1.4 ± 1) × 10 −10 s, Ueff(2) = 70.6K and a temperature independent term for 1 of (27 ± 2) Hz.The assumption of two-phonon Raman process contributing to the relaxation in place of Orbach one did not provide reasonable fit for any of the two phenomena.
The slower relaxation phenomenon was further investigated as a function of external applied field at 4 K.In the lower field range an increase in χ′′ is observed above 10 KHz suggesting the possible presence of a fast relaxation process (as confirmed in Figure S6).An increase in static field magnitude gave rise to a relatively slow process at a frequency of 20 Hz, with further slowing down of magnetization reversal with increasing field as shown in Figure 5a.Upon comparing this fielddependent phenomenon to the ones observed at 2 KOe, it is quite clear that this is actually the same relaxation phenomenon observed in the 2 KOe data.Interestingly, upon increasing the magnetic field up to 30 KOe, the fraction of molecules relaxing via this process decreases till it fades away above 20KOe (Figure S7).It is significant that while this relaxation channel is quenched, an additional peak appears for fields above 8 KOe, indicating the opening of another field triggered relaxation channel.This is a very slow process and it continues to gain intensity and to slow down while digressing The temperature dependences of the relaxation rates obtained by the fit of χ(ν) are reported as Arrhenius plot in Figure 4b: it is evident that for T > 2.5 K the behavior is similar for the two phenomena, being essentially linear.However, the slope is different, indicating that the anisotropy barrier is different for the two processes.Finally, at low temperature, where only process 1 is visible, a deviation from the linear behavior is observed, indicating a quantum regime.Keeping the above observations in mind, the data was tentatively analyzed using Equation (3), providing the following best fit parameters: τ 0 (1) = (6.2˘1.5) ˆ10 ´10 s, U eff (1) = (39.9˘0.8) K, τ 0 (2) = (1.4˘1) ˆ10 ´10 s, U eff (2) = 70.6K and a temperature independent term for 1 of (27 ˘2) Hz.The assumption of two-phonon Raman process contributing to the relaxation in place of Orbach one did not provide reasonable fit for any of the two phenomena.
The slower relaxation phenomenon was further investigated as a function of external applied field at 4 K.In the lower field range an increase in χ 11 is observed above 10 KHz suggesting the possible presence of a fast relaxation process (as confirmed in Figure S6).An increase in static field magnitude gave rise to a relatively slow process at a frequency of 20 Hz, with further slowing down of magnetization reversal with increasing field as shown in Figure 5a.Upon comparing this field-dependent phenomenon to the ones observed at 2 KOe, it is quite clear that this is actually the same relaxation phenomenon observed in the 2 KOe data.Interestingly, upon increasing the magnetic field up to 30 KOe, the fraction of molecules relaxing via this process decreases till it fades away above 20 KOe (Figure S7).It is significant that while this relaxation channel is quenched, an additional peak appears for fields above 8 KOe, indicating the opening of another field triggered relaxation channel.This is a very slow process and it continues to gain intensity and to slow down while digressing outside the experimental range (less than 0.1 Hz) at higher fields (Figure S7).Over the whole field range, the Debye modeling of the data yield a range of the α parameter distribution 0.25-0.45,which indicates a broad distribution of relaxation times.The field dependence of the relaxation rate obtained by these data (Figure 5b) clearly points to a non-linear behavior: for field less than 2.5 KOe a large dependence is observed, followed by a smoother increase above this field.Such behavior suggests that quantum tunneling processes dominates the relaxation, with superimposed contributions from field independent processes like the two-phonon Orbach process (Raman one having been discarded on the basis of the temperature dependence of the relaxation rate).In contrast, the direct term does not appear to play a major role here, since the relaxation time continues to increase up to the highest measured field.The reproduction of relaxation times at 4 K using the semi-quantitative technique employed for the 2 K field dependent data verified the aforementioned arguments.The fit procedure performed on the basis of Equation ( 2), while pointing to a well definite value for the field-independent term (92.2 ± 5 Hz), indicated that the absolute value of the parameters phenomenologically describing the QT could not be determined, but only their ratio: best fit was obtained with a ratio B2/B1= 5.71 × 10 −8 , with no appreciable distinction among curves for B1 > 5 × 10 5 Hz, which could then be considered as a lower limit.Despite the phenomenological nature of this set of parameters, one can conclude, in agreement with experimental observation, that a significant contribution of quantum tunneling survived even at the highest measured fields of 0.5 T.

Discussion
Manifestation of multiple relaxation channels is not unusual in lanthanide based complexes [25,30].However identification and contribution of the individual processes has proven to be quite tricky, due to the complex energetic structure of these ions.Furthermore, analysis of multiple magnetization reversal channels normally requires that each channel is due to a unique or a combination of known relaxation mechanisms like temperature dependent (Orbach, Raman or direct process), temperature independent/field-dependent QTM processes and so on.The case of DyCo reported here however presents a different scenario, where the two peaks identified above behave almost identically on increasing temperature, as seen in Figure 4b.Both the channels apparently undergo reversal of magnetization via the Orbach mechanism, giving an energy barrier of approximately 40 and 70 K, with channel 1 showcasing a small quantum tunneling contribution to The field dependence of the relaxation rate obtained by these data (Figure 5b) clearly points to a non-linear behavior: for field less than 2.5 KOe a large dependence is observed, followed by a smoother increase above this field.Such behavior suggests that quantum tunneling processes dominates the relaxation, with superimposed contributions from field independent processes like the two-phonon Orbach process (Raman one having been discarded on the basis of the temperature dependence of the relaxation rate).In contrast, the direct term does not appear to play a major role here, since the relaxation time continues to increase up to the highest measured field.The reproduction of relaxation times at 4 K using the semi-quantitative technique employed for the 2 K field dependent data verified the aforementioned arguments.The fit procedure performed on the basis of Equation ( 2), while pointing to a well definite value for the field-independent term (92.2 ˘5 Hz), indicated that the absolute value of the parameters phenomenologically describing the QT could not be determined, but only their ratio: best fit was obtained with a ratio B 2 /B 1 = 5.71 ˆ10 ´8, with no appreciable distinction among curves for B 1 > 5 ˆ10 5 Hz, which could then be considered as a lower limit.Despite the phenomenological nature of this set of parameters, one can conclude, in agreement with experimental observation, that a significant contribution of quantum tunneling survived even at the highest measured fields of 0.5 T.

Discussion
Manifestation of multiple relaxation channels is not unusual in lanthanide based complexes [25,30].However identification and contribution of the individual processes has proven to be quite tricky, due to the complex energetic structure of these ions.Furthermore, analysis of multiple magnetization reversal channels normally requires that each channel is due to a unique or a combination of known relaxation mechanisms like temperature dependent (Orbach, Raman or direct process), temperature QTM processes and so on.The case of DyCo reported here however presents a different scenario, where the two peaks identified above behave almost identically on increasing temperature, as seen in Figure 4b.Both the channels apparently undergo reversal of magnetization via the Orbach mechanism, giving an energy barrier of approximately 40 and 70 K, with channel 1 showcasing a small quantum tunneling contribution to relaxation.
Observation of multiple barriers to relaxation in a lanthanide based molecular magnet in a non-varying field is usually attributed to the fact that a certain percentage of the molecules in the cluster behave independently resulting in their own unique dynamics.This could be due to lattice disorder leading to unique species within the complex or the fact that a bunch of molecules resulted in different anisotropy, as first reported for polynuclear SMMs such as Mn 12 and its derivatives [31].In the present case, as mentioned above, the neutral dinuclear units of DyCo are interacting via hydrogen bonds to form a 3D network, involving water molecules from crystallization [7,17].It is then highly probable for the complex to undergo structure distortion upon partial loss of crystallization water over a length of time, leading to formation of two unique entities within the complex and thus in the observation of two different channels for magnetization reversal [32].In this sense, since the two processes were observed to contribute almost equally to the total relaxation of magnetization in the complex, as indicated by the relative integrals of the maxima related to the two channels, one can speculate that quantitatively almost equal number of molecules of each species contribute to their relaxation pathways.Furthermore, the relatively broad distribution of relaxation times pictured by the large α values suggests a large disorder involving the two species.
A more astute rationalization of the field and temperature dependence of the two peaks is shown in Figure S8.The data portrayed depicts behavior of the two peaks at 2 and 4 K, both at 0.3 and 2 KOe.Comparison of the behavior of the temperature dependence of the relaxation rate confirms that the slow process identified in the low temperature and field dynamics section is peak 1 in Figure 4a, whereas peak 2 represents a phenomenon induced by high field and high temperature (see Figure 4b).Indeed, similar slope (i.e., similar effective barriers) are observed for the measures at the two fields, with only a slight increase in the effective barrier at higher field and a corresponding decrease in the temperature independent contribution [32].In particular, in the present case, an increase in static field magnitude from 300 Oe to 2 KOe leads to an apparent increase in the effective barrier from 34.4 to 40 K: however, since we included explicitly temperature independent contributions and direct terms in the fits, the two values should ideally be identical [33].The observed difference should then not be considered as physically significant but possibly portrays the different weights of other unresolved contributions to the relaxation, varying with the field.Finally, we note that the observation of the slow process on increasing the field at 4 K may be attributed, in this framework, to a collective process resulting from intermolecular dipolar interactions: a definite proof of this would require, however, the analysis of isostructural YCo sample doped with Dy III , which is beyond the scope of this paper.

Synthesis and Structural Characterization
DyCo was synthesized following the procedure reported in [17].The purity of the microcrystalline powder used for magnetic dynamic characterization was checked by X-ray powder diffraction carried on a Bruker D8 advance powder diffractometer equipped with a Cu source (K α , λ = 1.54 Å).X-ray powder spectra of the DyCo complex was observed to superimpose with the spectrum obtained from the analysis of the molecular structure obtained by single crystal X-ray diffraction (Figure S9).

Magnetic Characterization
The AC magnetic susceptibility measurements were carried out on a Quantum Design MPMS SQUID magnetometer and Quantum Design PPMS in AC mode at both zero and applied external DC field in the presence of 5 Oe oscillating magnetic field.The later was used for high frequencies of 10 Hz to 10 KHz whereas the was employed for more sensitive measurements 0.1 Hz to 1 KHz.All the measurements were executed on polycrystalline sample pressed in a pellet (diameter 5 mm) to avoid individual field alignment of the crystallites.Raw data were reduced to paramagnetic molar susceptibilities using molecular formula obtained by X-ray data.This also allowed for evaluation of the intrinsic diamagnetism of the sample by use of corresponding Pascal's constants.As a further purity check, dc magnetic data were measured on the same pellet (Figure S10), and turned out to be in agreement with previously reported characterization [17].

Conclusions
We reported here a new example of slowly relaxing Dy III based complex, which has been selected on the basis of previously reported static magnetic characterization.The analysis of its dynamic magnetic properties clearly demonstrates the need for detailed field and temperature dependent studies to disentangle the different contributions to the relaxation.In particular, three different processes have been identified, of which one is attributed to a collective process resulting from the network of intermolecular interactions characterizing these systems, and two single molecular ones, which may arise from the partial loss of crystallization solvent.Despite this problematic issue, we feel that this family of cyano-bridged 3d-4f complexes should be further investigated, in view of the interesting photo-switchable properties of some of its class, up to relatively high temperature [34,35] which may open perspective applications in the fields of photo-controlled switches and information storage.

Figure 2 .
Figure 2. (a) Frequency dependence of the out-of-phase χ′′ component of the molar AC susceptibility of DyCo measured at 2 K at variable external static magnetic field (red 0 Oe to blue 5 KOe).The solid lines are Debye fits at respective fields; (b) Field dependent behavior of the relaxation time observed at 2 K; the solid line is a semi-quantitative fit to the data obtained with parameters reported in the text.

Figure 2 .
Figure 2. (a) Frequency dependence of the out-of-phase χ 11 component of the molar AC susceptibility of DyCo measured at 2 K at variable external static magnetic field (red 0 Oe to blue 5 KOe).The solid lines are Debye fits at respective fields; (b) Field dependent behavior of the relaxation time observed at 2 K; the solid line is a semi-quantitative fit to the data obtained with parameters reported in the text.

Figure 3 .
Figure 3. Arrhenius plot of the temperature dependence of the magnetic relaxation time τ of DyCo, in the presence of a dc magnetic field 300 Oe.The solid line is a fit corresponding to Equation (3) with the best fit parameters reported in the text.

Figure 3 .
Figure 3. Arrhenius plot of the temperature dependence of the magnetic relaxation time τ of DyCo, in the presence of a dc magnetic field 300 Oe.The solid line is a fit corresponding to Equation (3) with the best fit parameters reported in the text.

Figure 4 .
Figure 4. (a) Observation of two identifiable temperature dependent relaxation phenomena in DyCo, in the presence of a static magnetic field of 2 KOe; (b) Temperature behavior of the relaxation rates extracted for the two relaxation phenomena from the Debye fit of the ac data at 2 KOe.The solid lines are the best fit to Equation (3) with the parameters explained in the text.Black squares are the relaxation rate observed at 300 Oe, reported here for comparison's sake.

Figure 4 .
Figure 4. (a) Observation of two identifiable temperature dependent relaxation phenomena in DyCo, in the presence of a static magnetic field of 2 KOe; (b) Temperature behavior of the relaxation rates extracted for the two relaxation phenomena from the Debye fit of the ac data at 2 KOe.The solid lines are the best fit to Equation (3) with the parameters explained in the text.Black squares are the relaxation rate observed at 300 Oe, reported here for comparison's sake.
Debye modeling of the data yield a range of the α parameter distribution 0.25-0.45,which indicates a broad distribution of relaxation times.

Figure 5 .
Figure 5. (a) Frequency dependence of the out-of-phase χ′′ component of the susceptibility measured at 4 K at variable external static magnetic field.The solid lines are guide for eyes; (b) Field dependent behavior of the relaxation time observed at 4 K, extracted from the Debye fit of the AC data, the solid lines represent the best attainable fit.

Figure 5 .
Figure 5. (a) Frequency dependence of the out-of-phase χ 11 component of the susceptibility measured at 4 K at variable external static magnetic field.The solid lines are guide for eyes; (b) Field dependent behavior of the relaxation time observed at 4 K, extracted from the Debye fit of the AC data, the solid lines represent the best attainable fit.