We resolved to find a set of global parameters to describe the relaxation processes in 1-Dy by simultaneous fitting of the three datasets to Eqn 3. As the high-field process is only observable in the 12 K dataset, we do not have sufficient information to determine the temperature depdencence of this process (which would be expected to have a linear or power-law T-dependence for the Direct or Raman-II relaxation processes, respectively ), and so we simply assume a linear temperature dependence. We used values from the individual fits as starting parameters, and performed the first optimisation with Q and p fixed before allowing all parameters to freely refine. This gives global best-fit parameters: Ueff = 924(17) K, τ0 = should correlate with Ueff. We hypothesized that in solution the N-Dy-N angle in 1-Dy may increase, as has been observed for the Yb congener. ) that accounts for skew as well as a distribution in relaxation times. and Table ). Furthermore, magnetic hysteresis loops are more open in frozen solution than in the solid state (but remain waist-restricted, Figures ). The low-temperature frozen solution relaxation data show a comparable QTM rate to the solid-state sample, whereas the high-temperature fits encompass a major component (65%) that has rates in line with the solid-state data, and a minor component (35%) that has considerably slower dynamics (Figure ). Whilst some of these relaxation data could be explained by precipitation of 1-Dy, the minor component and the hysteresis data are both in accord with a sample that relaxes measurably slower than the solid-state material. These observations are consistent with the frozen solution sample containing a broader and more asymmetric distribution of molecular geometries, with some molecules in the distribution having larger N-Dy-N angles and hence larger anisotropy and slower magnetic relaxation rates.