research papers
Quantifying magnetic anisotropy using X-ray and neutron diffraction
aDepartment of Chemistry, Aarhus University, Langelandsgade 140, Aarhus C 8000, Denmark, bLLB, CEA, CE de Saclay, Gif sur Yvette 91191, France, and cDepartment of Chemistry, IIT Bombay, Powai, Mumbai, Maharashtra 400076, India
*Correspondence e-mail: jacobo@chem.au.dk
In this work, the magnetic anisotropy in two iso-structural distorted tetrahedral Co(II) complexes, CoX2tmtu2 [X = Cl(1) and Br(2), tmtu = tetramethylthiourea] is investigated, using a combination of polarized neutron diffraction (PND), very low-temperature high-resolution synchrotron X-ray diffraction and CASSCF/NEVPT2 ab initio calculations. Here, it was found consistently among all methods that the compounds have an easy axis of magnetization pointing nearly along the bisector of the compression angle, with minute deviations between PND and theory. Importantly, this work represents the first derivation of the atomic susceptibility tensor based on powder PND for a single-molecule magnet and the comparison thereof with ab initio calculations and high-resolution X-ray diffraction. Theoretical ab initio theory (AILFT) analysis finds the dxy orbital to be stabilized relative to the dxz and dyz orbitals, thus providing the intuitive explanation for the presence of a negative zero-field splitting parameter, D, from coupling and thus mixing of dxy and . Experimental d-orbital populations support this interpretation, showing in addition that the metal–ligand covalency is larger for Br-ligated 2 than for Cl-ligated 1.
1. Introduction
Molecular magnetism (Kahn, 1993; Gatteschi et al., 2006) is one particular branch of science where researchers have been able to build a level of understanding that allows the prediction (Rinehart & Long, 2011; Chilton et al., 2015) and subsequent realization of a single-molecule magnet with a colossal energy barrier towards magnetic relaxation (Goodwin et al., 2017; Guo et al., 2018, 2017). This was a groundbreaking result for the science of single-molecule magnets (SMMs), as well as a massive boost of faith in traditional scientific behavior.
SMMs are molecules that exhibit a unique magnetic axis, along which an external magnetic field can easily align the molecular magnetization vector either parallel or anti-parallel. We can view this as assuming values of either `0' or `1', in which case the molecule itself becomes a carrier of binary information. Molecules belonging to this class will resist the loss of directionality of the magnetization, and we may consider them as miniature memory bits. Unfortunately, the above-mentioned SMMs are not ideal for building new technology, as they are coordinatively unsaturated lanthanide-ions (and most often dysprosium) which are generally unstable and difficult to handle (Gupta et al., 2016). Instead, focus could be placed on ions that are more abundant and yield more stable compounds, and there is general consensus in the literature that Co(II) is an ideal target. It exhibits a large parameter; has an odd number of electrons meaning it is a Kramers ion, which significantly reduces the detrimental effects of quantum tunneling that shortcuts the relaxation barrier; and most of its complexes (with ≥ 4) are stable under ambient conditions.
At the very heart of SMM function is the phenomenon of magnetic anisotropy. This quantifies the difference in the magnetic response when an external magnetic field is applied in different directions relative to the molecule. One caveat is that, although high-level ab initio calculations are able to calculate both magnetic anisotropy (Neese et al., 2019) and, recently, magnetic relaxation properties (Reta et al., 2021), such calculations are experts tools, and we are still not able to predict the exact magnetic properties of such compounds from a simple set of rules.
The origin of magnetic anisotropy is the presence of orbital angular momentum in the ground state or in a relatively near-excited state. For transition metals, we often quantify this anisotropy using the zero-field splitting parameter (D). We are able to estimate the value of D from high-level theory (Chibotaru, 2015), which provides a magnitude and a direction of the anisotropy and also allows an explanation that is based on the electronic structure (Gomez-Coca et al., 2013). Energy differences between states and the relative magnitude of the eigenvectors of the g-tensor can also be experimentally measured using spectroscopic methods such as EPR and INS (Abragam & Bleaney, 2012; Sigrist et al., 2015), but these approaches provide no geometrical insight into the direction of the anisotropy axes relative to the molecular structure. Thus, the only link that binds theory and experiment together is the comparison of the D values, which in itself is of very limited use as a design criterion. Therefore, to reach a quantitative correlation, we need to measure direction-dependent properties, and this is exactly the focus of this work.
It has been shown on several occasions that four-coordinate Co(II) complexes in distorted tetrahedral geometries exhibit very high negative D-values, and some explanations have been presented by us and others (Zadrozny et al., 2013; Zadrozny & Long, 2011; Tripathi et al., 2019; Vaidya et al., 2018, 2016, 2014; Rechkemmer et al., 2016; Damgaard-Møller et al., 2020c). In the current work, we combine results from neutron and synchrotron diffraction to provide irrefutable experimental evidence for the origin of this behavior. Thus, we have re-examined the previously published Co(II)-complexes, CoX2tmtu2 [tmtu = tetramethylthiourea, X = Cl (1), Br (2)]. Detailed magnetic data of these two complexes were reported by some of us (Vaidya et al., 2018) and the dc magnetic data presented therein is recalled here for comparison. The temperature-dependent powder data [χMT(T)] in the temperature range 2–300 K are shown in Fig. S1. The χMT value of both complexes 1 and 2 decreases gradually from room temperature to 40 K, below this temperature χMT drops precipitously and reaches a final value of 0.8 and 1.2 cm3 K mol−1 at 2.0 K, respectively. The low temperature drop in χMT in both complexes was attributed to the magnetic anisotropy associated with the S = 3/2 ground state of 1 and 2. Consistent with this observation, the of both 1 and 2 does not show any sign of magnetization saturation, even at 70 kOe external magnetic field at 2.0 K in the field-dependent magnetization measurements. Significantly low values for 1 (2.07 NμB) and 2 (2.13 NμB) at this limit (70 kOe and 2.0 K) further strongly corroborate the presence of relatively large magnetic anisotropy in these structurally analogous complexes. To quantify the spin Hamiltonian (SH) parameters associated with these complexes (both 1 and 2), the magnetic data [χM(T) and M(H)] were fitted simultaneously providing the following SH parameters for 1, D = −18.1 cm−1, giso = 2.26; and for 2, D = −16.4 cm−1, giso = 2.33. We emphasize the good agreement between the experimental and simulated magnetic data [χMT(T)] using the computed SH parameters (gx = 2.25, gy = 2.28 and gz = 2.48, D = −18.57 cm−1. ∣E/D∣ = 0.067 for 1; and gx = 2.26, gy = 2.30 and gz = 2.49, D = −17.75 cm−1, ∣E/D∣ = 0.086 for 2). This shows that the incorporation of a transverse component in the simulation of magnetic data does not significantly alter the magnetic anisotropy (D) extracted experimentally in the previous study (Vaidya et al., 2018), where values of −18.1 cm−1 and −16.4 cm−1 for 1 and 2 were obtained using giso and axial anisotropy. This is presumably due to the relatively small ∣E/D∣ value associated with these complexes.
Recently, some of us have shown how polarized neutron diffraction (PND) from single crystals can provide accurate et al., 2018; Ridier et al., 2016; Gukasov & Brown, 2002; Tripathi et al., 2021). Now, a potentially disruptive innovation allows similar information to be obtained from PND on powder (pPND) samples (Kibalin & Gukasov, 2019) and we present the first thorough analysis using this approach. In a parallel research path, we have recently shown how the d-orbital populations obtained from multipole modeling of single-crystal X-ray diffraction (Holladay et al., 1983) provide significant insight into the magnetic properties of transition-metal based SMMs (Craven et al., 2018; Thomsen et al., 2019; Bunting et al., 2018; Damgaard-Møller et al., 2020a,b,c). In the current work, we have, for the first time, combined these advanced diffraction techniques to obtain both the electronic and the magnetic structure, and added theoretical calculations, which allow for an ab initio theory analysis (Atanasov et al., 2012). As a result, we obtain an unrivaled understanding of the molecular magnetic properties of 1 and 2.
tensors, clearly showing the direction and size of the easy axis of magnetization (Klahn2. Experimental methods
2.1. Synthesis
Crystals of 1 and 2 were obtained from synthesis and recrystallization according to the published procedure (Vaidya et al., 2018).
2.2. Synchrotron X-ray diffraction and multipole modeling.
Single crystals of 1 and 2 of suitable size for the synchrotron beam were mounted on a Huber goniometer head via a glass-fiber and using glue. The crystal quality was tested at room temperature and subsequently cooled to low temperature. Data were collected at beamline BL02B1 at SPring-8, which is equipped with a Pilatus3 X CdTe 1M detector and a gaseous He-cooling device. The datasets consist of 180° omega-scans in steps of 0.5° with χ fixed at 0, 20 and 45° and with 2θ = 0 and 20° (scan speeds of 2 and 4° s−1 for 1, and 1 and 2° s−1 for 2), thus six runs for each crystal. The detector frames are converted to Bruker format, followed by a Lorentz and polarization correction, and the raw images are then integrated using SAINT+ (v8.38A). Subsequently, integrated intensities are corrected for absorption and other effects in SADABS (Krause et al., 2015), and resulting unmerged data without application of an error-model are then finally merged in SORTAV (Blessing, 1997). Final hkl files are used to solve and refine structures in SHELXT (Sheldrick, 2015a) and SHELXL (Sheldrick, 2015b), respectively. The resulting structures are imported into XD (Volkov et al., 2006) and a multipole model is incrementally built, until we reached the final model, which for 2 included anharmonic parameters for bromium (Herbst-Irmer et al., 2013). We emphasize here the importance of using radial functions based on Hartree–Fock calculations for the metal center, in order to obtain physically reliable d-orbital populations (Damgaard-Møller et al., 2020c). Further details about the multipole model are deposited in the supporting information; Table 1 contains essential crystallographic results. Fig. 1 shows the molecular structures, while Fig. 2 shows the fractal dimensionality plots (Meindl & Henn, 2008) for both compounds, indicating the high reliability of the final models. For the ensuing discussion of the d-orbital populations, it is important to clarify the choice of local coordinate system. In both complexes, we have defined the z axis to bisect one X—Co—S angle such that it points from Co towards the tmtu ligands. There are several reasons for this choice: (1) the deformation density maps around Co show clear maxima corresponding to these directions; (2) it corresponds to directions obtained from theoretical calculations; (3) the perception of the coordination sphere of Co as a distorted tetrahedron fits best with this direction as the unique axis in approximate D2d point symmetry.
|
2.3. Polarized neutron diffraction
Flipping ratio measurements based on experiments with single-crystal PND, and its intimate connection with the magnetization density were first explained in 1959 by Nathans et al. (1959). The period from 1970 to 1990 saw many experimental studies with particular emphasis on spin-only 3d metal molecular complexes to derive distributions, with substantial contributions from Figgis and Reynolds (Daul et al., 1988; Figgis et al., 1988, 1987, 1983, 1982) and others (Forsyth, 1977). Focus was, at that time, on complexes with little or no orbital angular momentum, for which the induced magnetization density vector aligns perfectly with the external field. The effect of significant orbital angular momentum is, however, that we can no longer assume such collinearity. The direct link between flipping ratios and magnetization density thus becomes invalid and the latter cannot be experimentally recovered. Under such circumstances, Brown and Gukasov showed that flipping ratio measurements can instead be used to retrieve a local atomic susceptibility tensor (Gukasov & Brown, 2002) which models this absence of collinearity.
Recently, their approach has been extended to study the magnetic anisotropy of paramagnetic compounds in polycrystalline samples (Kibalin & Gukasov, 2019; Gukasov & Brown, 2010). Paramagnetic complexes are at the center of SMM research, and this development is therefore of enormous importance. Until recently, the lack of appropriate software made the powder PND technique practically inapplicable to data but a library of dedicated computer code is now available and, having established an initial proof-of-concept, the technique can be considered ready for application to novel systems (Kibalin & Gukasov, 2019). We envisage that this approach will receive a massive boost when the European Spallation Source becomes operational in a few years time.
A detailed description of the theoretical basis for the PND method for both powder and single crystals is provided in the supporting information.
PND measurements for this work for both the powder- and single-crystal studies were performed at the thermal polarized neutron lifting counter diffractometer 6T2 (LLB-Orphée, Saclay). Neutrons were monochromated to a wavelength of 1.4 Å by a vertically focusing graphite crystal and polarized by a supermirror bender. The polarization factor P of the beam was 0.95 for the powder measurements on 1 and 0.78 for the single-crystal measurements on 2. Details of the powder data collection are provided by Kibalin & Gukasov (2019) and details on the data reduction and are supplied in the supporting information. The powder diffraction patterns for 1 were measured at a magnetic field of 1 T and a temperature of 2 K using a position-sensitive detector. Single-crystal measurements for 2 were made in the magnetic field of 1 T at 3 K, and the flipping ratios were collected for five different sample orientations. Flipping ratios for compound 2 were then extracted from the raw images by employing the in-house data reduction suite at the LLB.
The Cryspy (v. 0.5.8) through the Cryspy Editor (v.1.5.6), both available for Python3.X through PyPI (Kibalin & Gukasov, 2019).
of susceptibility tensors both for single-crystal and powder PND data were performed with the newly developed software library2.4. Theoretical calculations
Ab initio calculations were performed using the ORCA software (4.1; Neese, 2012, 2018) with the solid-state geometry of complexes 1 and 2 obtained from synchrotron X-ray diffraction at 20 K. CASSCF(7,5) (Malmqvist & Roos, 1989) and subsequent NEVPT2 (Angeli et al., 2001) correction was performed including all the ten quartet S = 3/2 states and using the Douglas–Kroll–Hess triple-ζ DKH-def2-TZVP basis set (Schäfer et al., 1992, 1994; Weigend & Ahlrichs, 2005). Including the 40 doublet states (S = 1/2 states) did not change the results. The AILFT (Atanasov et al., 2012, 2013) program was used to get an estimated d-orbital splitting of the compounds. The SOC was accounted for on the basis of non-relativistic eigenstates using quasi-degenerate perturbation theory. Relativistic electron densities used in the analysis were obtained from the ground Kramers doublet and first exited Kramers doublet. Projection of the two lowest Kramer doublets onto an S = 3/2 pseudo spin furthermore allowed for the extraction of the SH parameters gx, gy, gz, D and E.
3. Results and discussion
We start with an analysis of the X2S2 coordination geometries in 1 and 2 obtained from the crystal structures at 20 K are summarized in Table 2. The deviation from tetrahedral symmetry is quite complicated to quantify, and the typical shape index descriptor (Pinsky & Avnir, 1998) does not add valuable insight. There is little asymmetry (<0.03 Å) in the bond lengths from cobalt to both X and S, internally in both compounds. Obviously, the bond angles show a strong deviation from tetrahedral geometry. Closer inspection suggests that two of three angles involving X2—Co are significantly larger than 109.47° and the third is quite close to this value, whereas two of three angles involving X1—Co are much smaller. Thus, we may describe the distortion of the initial ideal tetrahedral coordination sphere as one where we view down the X2—Co axis and squeeze S1 and X1, while S2 is more or less left in place. However, we choose instead to view the coordination sphere as a compressed tetrahedron with approximate D2d symmetry, where the unique axis is the one that bisects the X1—Co—S1 angle. This axis is the one we use as the z axis in the definition of the orbital functions in the multipole modeling (vide infra).
The Co
|
3.1. Theoretical results
Before we describe the experimental results (Section 3.2 and onwards), we look at the calculated electronic structures of 1 and 2, obtained using the ORCA program suite from CASSCF(7,5) and subsequent NEVPT2 correction to the energies. Two main results from the calculations will be described in turn in the following: (1) the AILFT orbitals provide an intuitive interpretation and prediction of the magnetic anisotropy; (2) the g, D and E values, which quantify the magnetic anisotropy, and can be compared with the experimental measurement of the tensors. In addition, Section 3.3 contains a comparison of the theoretically calculated electron density (ED) and the experimentally derived ED.
(1) The AILFT analysis of an ab initio calculation provides orbitals independent of the input molecular coordinate system, thus providing the most consistent orbitals. As already discussed, both 1 and 2 are best described as distorted tetrahedral complexes (see above). Therefore, we expect to have the e-type (in Td symmetry) orbitals (i.e. and dz2) stabilized, and the t2-type orbitals (i.e. dxy, dxz and dyz) destabilized. The resulting AILFT orbitals and their energies are depicted in Fig. 3. Several indicators suggest that the chosen assignment is correct. First of all, it matches the results of several other distorted tetrahedral Co(II) complexes, in which the dz2 orbital is stabilized, and has its lobes pointing in the direction of the narrowest L—Co—L bisecting angle (Carl et al., 2015; Rechkemmer et al., 2016; Zadrozny & Long, 2011; Zadrozny et al., 2013). Secondly, this choice stabilizes dxy compared with the other t2-orbitals, and the coupling of and dxy would predict easy-axis type anisotropy in the direction of the z axis, which is also true in the present case. The last indication that these are `good' orbitals is provided by the composition of the CASSCF wavefunction, which in the given coordinate system consist of fairly pure single Slater-determinants (91.6% for the ground state and 84.0% for the first in 1 and 89.7% for the ground state and 76.2% for the first in 2). We note that the orbital ordering reported for both complexes in the earlier report (Vaidya et al., 2018) are slightly different from the orbital ordering reported here. This is because, in the earlier paper, we reported the CASSCF ground state wavefunction which is multideterminant in character, whereas in this manuscript we report AILFT orbital ordering. However, we noticed the same orbital ordering (as we notice in this manuscript) when we compute the AILFT orbitals.
(2) The pseudo-spin 3/2 SH approximation allows for the extraction of the magnetic anisotropy from the ab initio calculations. The eigenvectors for the D- and g-tensor are almost parallel, and the eigenvalues as well as the directions of the vectors with respect to the coordinate system chosen for the AILFT orbitals above are found in Table 3.
|
The D-tensor gives D = −18.57 cm−1, ∣E/D∣ = 0.067 for 1 and D = −17.75 cm−1, ∣E/D∣ = 0.086 for 2. The easy axis of the compounds (i.e. characterized by Dzz or gz) is not parallel with the chosen molecular z axis for the AILFT orbitals, but forms an angle with this z axis of 15° for 1 and 11° for 2 (see Fig. 4). In the case of pure orbitals, the easy axis aligns perfectly along the molecular z axis, and the observed easy axis thus shows the inadequacy of the d-orbital diagram to accurately predict the magnetic anisotropy. This is caused by the multi-configurational nature of the ground state and first-excited state, where only about 90% is a pure Slater determinant, and the remaining 10% of the ground state allows for change in the direction of the easy axis. Nevertheless, the d-orbital scheme gives a useful and intuitive explanation for both sign and magnitude of the magnetic anisotropy, providing valuable guidelines for a synthetic chemist in the search for new transition metal SMMs (Ruamps et al., 2013; Gomez-Coca et al., 2013).
Previous studies have found clear evidence that distortion of the tetrahedral geometry around Co(II) can lead to strong magnetic anisotropy (Rechkemmer et al., 2016; Tripathi et al., 2019; Vaidya et al., 2018, 2017; 2016; Legendre et al., 2021). This has been explained by the appearance of a strong splitting of the d-orbitals of t2-symmetry and, to a smaller extent, of the e-symmetry orbitals, with the great benefit that the energies of the and dxy orbitals approach each other. Perhaps surprisingly, it is often observed that the is destabilized relative to dz2 upon tetrahedral compression. However, regardless of the order of these two latter orbitals, a large and negative D-value is predicted and most often realized. As explained above, the geometrical distortions of the coordination spheres in 1 and 2 are rather complicated and the unique axis of magnetization is extremely difficult to rationalize from structural considerations alone. Similarly, the theoretical result indicates a significant deviation of the anisotropy axis away from the molecular z axis. We have therefore determined the direction and the associated anisotropy of the unique axis of magnetization using PND for both compounds, as described in Section 3.2.
3.2. Experimental results: PND
The experimental susceptibility tensors, given in the normalized reciprocal coordinate system,
have the form
Fig. 4 shows these susceptibility tensors as ellipsoids overlaid on the molecular structures of 1 and 2 and with the structures rotated so as to minimize the root-mean-square distance (RMSD) between the first coordination sphere and Co atoms in the two structures. The final RMSD value obtained from this procedure is 0.0942 Å. It is evident that the experimental easy axis of magnetization in both compounds approximately bisects the X1—Co—S1 angles, with deviations from the Cl1—Co—S1 and Br1—Co—S1 bisectors being 24.6 and 10.8° for 1 and 2, respectively, and an approximate deviation between the two experimental easy axes of 23.3° [Fig. 4(b)]. For comparison, the angle between the two theoretical easy axes of 1 and 2 is calculated to be 12.4°. The extent of axiality, which we define here as the ratio of the largest eigenvalue to the average of the two smallest eigenvalues, is 14.8 for 1 and a staggering 33.1 for 2. An even better metric for the axiality, using the susceptibility tensors for these two structures might however be what we here call an `effective' ∣E/D∣, based on the relationship , where , ; and χ1, χ2 and χ3 are the largest, intermediate and smallest eigenvalues of the tensor. This metric gives 0.15 and 0.018 for 1 and 2, respectively. The theoretically calculated value of ∣E/D∣ for 1 is roughly half of what we find here experimentally. The experimental extent of axiality is quite extreme for 2, and thus the effective ∣E/D∣ is much smaller for this compound than what is found theoretically.
In addition to a comparison based on their eigenvalues, the susceptibility tensors also allow us to compare the easy axis directions obtained from PND with the calculated ones [Fig. 4(b)]. We see a striking agreement between the calculated and measured directions, particularly for 2, with deviations between the theoretical and experimental result of only 12.1 and 1.8° for 1 and 2, respectively.
3.3. Experimental results: electron density
On the basis of this unique result from the PND analysis, we now discuss the experimental and theoretical ED analysis. As mentioned above, the z axis of the local coordinate system that defines the multipolar functions is chosen to match the direction of the AILFT orbitals from ab initio theory (i.e. along the angle bisector of X1—Co—S1). The y axes point between the two halogens, and the x axes approximately bisect the X1—Co—S2 angle. With this convention, we obtained the experimental d-orbital populations (Table 4) from the multipole parameters using the transformation matrix given by Holladay et al. (1983). To discuss the values obtained, it is important to know what they represent. When looking at the simplest description of the system, with pure d-orbitals and electrons placed according to the Aufbau principle (see Fig. 3), we expect the populations to be 1.00 in the t2 antibonding orbitals and 2.00 in the e non-bonding orbitals. This is obviously a too-simplistic description of the true ED around cobalt, as it is affected by metal–ligand and SOC, and perhaps also geometrical distortions and relaxation effects. Next, we discuss the potential impact on the d-orbital populations of the two former effects.
|
Metal–ligand electron transfer/covalent bonding is a well known effect for transition metal compounds. Its influence on d-orbital populations is due to the fact that ligands of the same symmetries form bonding and antibonding MOs with contributions from some of the d-orbitals (Figgis & Hitchman, 2000). In approximate tetrahedral symmetry as in 1 and 2, the t2-orbitals [d(xy), d(yz), d(xz)] have the same symmetry as the symmetry-adapted linear combination of the four σ-donating lone pairs on the ligands. The bonding orbitals of these are fully occupied and primarily constituted by ligand orbitals, and not part of the conventional active space. However, due to the finite ligand contribution to these MOs, they will add modest ED near the metal center, and therefore it is to be expected that the average d-orbital population of the t2-orbitals is larger than 1.000. Furthermore, this effect would add more electrons to the dyz and dxz orbitals than to dxy since the latter is slightly stabilized relative to the two former (see Fig. 3) and thus experience less overlap with the ligand orbitals. This is clearly expressed in the experimental d-orbital populations of 2, which show a significantly larger population of dyz and dxz compared with the dxy orbital, and thus indicates that 2 has a much larger degree of metal–ligand compared with 1. This is backed by a topological analysis of the ED (Table 5) which shows that the total energy density – to some extent a measure of covalency in a chemical bond (Gatti, 2005) – is consequently less positive for 2 than for 1, indicating that the σ-donation from the halide is indeed higher for the Br-complex. However, in the populations for the other models shown in Table 4, the dxy orbital is populated more than the dyz and dxz orbitals. We may ascribe this to the effects of SOC, as explained below.
|
SOC is an effect that couples ground and excited states, and thus in effect leads to withdrawal of electrons from the doubly occupied orbitals, and insertion into some of the singly occupied orbitals. In this specific case, the ground-state Kramers doublet consists of 97% of the spin-orbit free ground state () and 3% of the first-excited spin-orbit free state () in both 1 and 2. Given the spin-orbit free states mentioned earlier, the resulting d-orbital populations of dxy and would be , thus in addition to the electrons from metal–ligand interaction, the population of dxy would be increased due to SOC.
There is thus a competition between the metal–ligand covalency on the one hand, which increases the populations of dxz, dyz more than it does dxy, and SOC on the other hand, which effectively moves electrons from to dxy. We have recently shown how the effect of the SOC can be used to quantify the ZFS from experimental d-orbital populations (Damgaard-Møller et al., 2020c). However, in the present case with much smaller SOC, the quantification is not reliable. We also note that the population of dz2 is less than 2.0, which is very likely due to mixing, which is possible as they have the same symmetry.
To evaluate the expected size of the two opposing effects described above, we obtained d-orbital populations from theory (Table 4), by of theoretical structure factors calculated from the relativistic, SOC ground state of compounds that provide the best possible description of real-world ED (Genoni, 2020; Gao et al., 2020). The relativistic SOC electron densities extracted from the ab initio calculation (i.e. the ED of the lowest two Kramer doublets) were combined into one ED using a Boltzmann averaging at 20 K, and this was modeled using the Hansen–Coppens multipole formalism (Hansen & Coppens, 1978). Table 4 shows that there is a clear, albeit small, difference suggesting that the σ-donation is larger in 2 than in 1, while the identical populations of in 1 and 2 suggest that the effect of SOC is identical, not surprising given the nearly identical energy level differences seen in Fig. 3.
4. Conclusions
We have quantified the magnetic anisotropy in the form of atomic susceptibility tensors for two isostructural four-coordinate Co(II) complexes using PND. The results clearly show highly axial and comparable magnetic anisotropy for both compounds, with unique axes aligned approximately along the direction of the molecular z axis, from Co and bisecting the X1—Co—S1 angle. Exceptionally, this study shows for the first time that PND studies using powder samples can provide atomic susceptibility tensors of comparable quality to more traditional single-crystal studies. Using the z axis direction, an AILFT analysis of the CASSCF results offer an explanation for the magnetic anisotropy showing the increased proximity of dxy and orbitals relative to an ideal tetrahedral coordination. The small and opposite effects of metal–ligand covalency and are identified in the electron density models refined using very high-resolution, low-temperature synchrotron X-ray diffraction data. Combined, the range of experimental results unanimously show that the σ-donation to the Co(II) ion is larger in 2 with Br-ligands than in the Cl-ligated compound 1.
5. Related literature
The following references are cited in the supporting information: Clark et al.(2007); Gukasov et al. (2007); Politzer & Murray (2019).
Supporting information
Combined https://doi.org/10.1107/S2052252521008290/lt5038sup1.cif
for 1 and 2. No structure factors are available as a multipole model was used. DOI:Additional information, figures and tables. DOI: https://doi.org/10.1107/S2052252521008290/lt5038sup2.pdf
C10H24Cl2CoN4S2 | F(000) = 820 |
Mr = 394.28 | Dx = 1.472 Mg m−3 |
Monoclinic, P21/n | Synchrotron radiation, λ = 0.2486 Å |
a = 9.9071 (4) Å | Cell parameters from 9216 reflections |
b = 12.7019 (5) Å | θ = 2.2–18.6° |
c = 14.1556 (6) Å | µ = 0.09 mm−1 |
β = 92.824 (2)° | T = 20 K |
V = 1779.16 (13) Å3 | Block, yellow |
Z = 4 | 0.1 × 0.1 × 0.1 mm |
Synchrotron diffractometer | Rint = 0.075 |
Absorption correction: multi-scan SADABS 2016/2: Krause, L., Herbst-Irmer, R., Sheldrick G.M. & Stalke D., J. Appl. Cryst. 48 (2015) 3-10 | θmax = 24.5°, θmin = 0.8° |
Tmin = 0.903, Tmax = 0.958 | h = −29→33 |
650031 measured reflections | k = −42→39 |
68940 independent reflections | l = −41→47 |
Refinement on F2 | 25757 reflections |
Least-squares matrix: full | 279 parameters |
R[F2 > 2σ(F2)] = 0.014 | 0 restraints |
wR(F2) = 0.031 | w2 = 1/[s2(Fo2)] |
S = 1.02 | (Δ/σ)max < 0.001 |
x | y | z | Uiso*/Ueq | ||
Co(01) | 0.776252 (2) | 0.757076 (2) | 0.513297 (2) | 0.005 | |
Cl(02) | 0.930775 (4) | 0.811110 (3) | 0.626713 (3) | 0.008 | |
Cl(03) | 0.851971 (4) | 0.758476 (3) | 0.365631 (3) | 0.008 | |
S(004) | 0.583861 (4) | 0.859940 (3) | 0.531818 (3) | 0.007 | |
S(005) | 0.717352 (5) | 0.583469 (3) | 0.555147 (3) | 0.007 | |
N(006) | 0.749841 (14) | 0.414117 (10) | 0.447409 (11) | 0.007 | |
N(007) | 0.619969 (14) | 0.546578 (10) | 0.376627 (10) | 0.007 | |
N(008) | 0.584101 (16) | 0.990104 (11) | 0.678718 (11) | 0.008 | |
N(009) | 0.739622 (14) | 1.030992 (11) | 0.567087 (11) | 0.007 | |
C(00A) | 0.694061 (15) | 0.510899 (11) | 0.452117 (11) | 0.006 | |
C(00B) | 0.83262 (2) | 1.086588 (15) | 0.633504 (17) | 0.012 | |
C(00C) | 0.523570 (17) | 0.632889 (13) | 0.383007 (13) | 0.009 | |
C(00D) | 0.641439 (14) | 0.967982 (11) | 0.596661 (11) | 0.006 | |
C(00E) | 0.857743 (19) | 0.379727 (15) | 0.514693 (15) | 0.011 | |
C(00F) | 0.684759 (19) | 0.328028 (13) | 0.392754 (15) | 0.01 | |
C(00G) | 0.56316 (2) | 1.098108 (15) | 0.711391 (15) | 0.011 | |
C(00H) | 0.51175 (2) | 0.909258 (16) | 0.729726 (15) | 0.011 | |
C(00I) | 0.65058 (2) | 0.518460 (15) | 0.279561 (13) | 0.01 | |
C(00J) | 0.783207 (19) | 1.026896 (15) | 0.470084 (15) | 0.011 | |
H(00A) | 0.812005 | 1.1698 | 0.630365 | 0.01809 (4)* | |
H(00B) | 0.935177 | 1.072463 | 0.614941 | 0.01809 (4)* | |
H(00C) | 0.819249 | 1.058432 | 0.704197 | 0.01809 (4)* | |
H(00D) | 0.489833 | 0.637657 | 0.454053 | 0.01304 (3)* | |
H(00E) | 0.438098 | 0.619008 | 0.334631 | 0.01304 (3)* | |
H(00F) | 0.57147 | 0.705823 | 0.365048 | 0.01304 (3)* | |
H(00G) | 0.923554 | 0.445099 | 0.532175 | 0.01623 (3)* | |
H(00H) | 0.914928 | 0.317754 | 0.483524 | 0.01623 (3)* | |
H(00I) | 0.81426 | 0.350943 | 0.577983 | 0.01623 (3)* | |
H(00J) | 0.586331 | 0.352972 | 0.365827 | 0.01533 (3)* | |
H(00K) | 0.674429 | 0.260415 | 0.437738 | 0.01533 (3)* | |
H(00L) | 0.745752 | 0.307593 | 0.33457 | 0.01533 (3)* | |
H(00M) | 0.620014 | 1.110537 | 0.777266 | 0.01694 (3)* | |
H(00N) | 0.457348 | 1.110506 | 0.721658 | 0.01694 (3)* | |
H(00O) | 0.596588 | 1.152778 | 0.659287 | 0.01694 (3)* | |
H(00P) | 0.410907 | 0.900942 | 0.698609 | 0.01690 (3)* | |
H(00Q) | 0.506682 | 0.931939 | 0.802776 | 0.01690 (3)* | |
H(00R) | 0.564354 | 0.835335 | 0.72549 | 0.01690 (3)* | |
H(00S) | 0.744964 | 0.476337 | 0.280102 | 0.01512 (3)* | |
H(00T) | 0.658167 | 0.588969 | 0.237897 | 0.01512 (3)* | |
H(00U) | 0.571103 | 0.469348 | 0.249256 | 0.01512 (3)* | |
H(00V) | 0.824198 | 1.102142 | 0.45128 | 0.01600 (3)* | |
H(00W) | 0.697983 | 1.008512 | 0.422788 | 0.01600 (3)* | |
H(00X) | 0.85946 | 0.967073 | 0.464594 | 0.01600 (3)* |
U11 | U22 | U33 | U12 | U13 | U23 | |
Co(01) | 0.005130 (8) | 0.004028 (8) | 0.004756 (8) | 0.000072 (6) | 0.000123 (6) | −0.000341 (6) |
Cl(02) | 0.008604 (12) | 0.006881 (11) | 0.007834 (12) | −0.000920 (8) | −0.002438 (9) | −0.000228 (9) |
Cl(03) | 0.009311 (12) | 0.008733 (12) | 0.006318 (12) | −0.000732 (9) | 0.002232 (9) | −0.000944 (9) |
S(004) | 0.005286 (11) | 0.005694 (11) | 0.008994 (14) | −0.000098 (8) | 0.000240 (9) | −0.002247 (9) |
S(005) | 0.011675 (14) | 0.004705 (11) | 0.005419 (12) | −0.000954 (10) | −0.000641 (10) | −0.000069 (9) |
N(006) | 0.00783 (3) | 0.00497 (3) | 0.00868 (4) | 0.00077 (3) | −0.00098 (3) | −0.00060 (3) |
N(007) | 0.00708 (3) | 0.00643 (3) | 0.00600 (3) | 0.00112 (3) | −0.00048 (3) | −0.00042 (3) |
N(008) | 0.01021 (4) | 0.00642 (3) | 0.00662 (4) | 0.00055 (3) | 0.00158 (3) | −0.00063 (3) |
N(009) | 0.00661 (3) | 0.00611 (3) | 0.00917 (4) | −0.00124 (3) | 0.00054 (3) | 0.00023 (3) |
C(00A) | 0.00679 (4) | 0.00472 (3) | 0.00608 (4) | 0.00004 (3) | −0.00033 (3) | −0.00028 (3) |
C(00B) | 0.00968 (5) | 0.00916 (5) | 0.01698 (7) | −0.00267 (4) | −0.00303 (5) | −0.00091 (5) |
C(00C) | 0.00807 (4) | 0.00785 (4) | 0.01008 (5) | 0.00191 (3) | −0.00044 (4) | 0.00008 (4) |
C(00D) | 0.00617 (3) | 0.00486 (3) | 0.00645 (4) | −0.00018 (3) | 0.00041 (3) | −0.00038 (3) |
C(00E) | 0.01027 (5) | 0.00865 (5) | 0.01321 (6) | 0.00207 (4) | −0.00269 (4) | 0.00126 (4) |
C(00F) | 0.01179 (5) | 0.00631 (4) | 0.01249 (6) | −0.00050 (4) | −0.00012 (4) | −0.00243 (4) |
C(00G) | 0.01506 (6) | 0.00881 (5) | 0.00996 (5) | 0.00281 (4) | 0.00018 (4) | −0.00305 (4) |
C(00H) | 0.01360 (6) | 0.01130 (5) | 0.00922 (5) | 0.00039 (4) | 0.00365 (4) | 0.00249 (4) |
C(00I) | 0.01169 (5) | 0.01193 (5) | 0.00656 (5) | 0.00171 (4) | −0.00007 (4) | −0.00106 (4) |
C(00J) | 0.01094 (5) | 0.00971 (5) | 0.01172 (6) | −0.00015 (4) | 0.00430 (4) | 0.00209 (4) |
CO(01)—CL(02) | 2.2694 (1) | CO(01)—S(005) | 2.3639 (1) |
CO(01)—CL(03) | 2.2554 (1) | S(004)—C(00D) | 1.7319 (1) |
CO(01)—S(004) | 2.3363 (1) | S(005)—C(00A) | 1.7311 (2) |
CL(02)—CO(01)—CL(03) | 114.1860 (18) | CL(03)—CO(01)—S(004) | 114.2076 (19) |
CL(02)—CO(01)—S(004) | 106.0885 (18) | CL(03)—CO(01)—S(005) | 109.6999 (18) |
CL(02)—CO(01)—S(005) | 105.7281 (18) | S(004)—CO(01)—S(005) | 106.3043 (18) |
C10H24Br2CoN4S2 | V = 1856.9 (16) Å3 |
Mr = 483.18 | Z = 4 |
Monoclinic, P21/n | Synchrotron radiation, λ = 0.24860 Å |
a = 9.798 (4) Å | Cell parameters from 9630 reflections |
b = 12.961 (8) Å | θ = 2.2–17.6° |
c = 14.636 (6) Å | T = 20 K |
β = 92.363 (8)° | Block, yellow |
Synchrotron diffractometer | θmax = 24.5°, θmin = 0.7° |
Absorption correction: multi-scan SADABS 2016/2: Krause, L., Herbst-Irmer, R., Sheldrick G.M. & Stalke D., J. Appl. Cryst. 48 (2015) 3-10 | h = −32→28 |
Tmin = 0.903, Tmax = 0.958 | k = −33→43 |
695086 measured reflections | l = −43→48 |
Rint = 0.068 |
Refinement on F2 | 26535 reflections |
Least-squares matrix: full | 315 parameters |
R[F2 > 2σ(F2)] = 0.014 | 0 restraints |
wR(F2) = 0.029 | w2 = 1/[s2(Fo2)] |
S = 0.99 | (Δ/σ)max = 0.00001 |
x | y | z | Uiso*/Ueq | ||
BR(1) | 0.135757 (5) | 0.252140 (4) | 0.638058 (4) | 0.006 | |
BR(2) | 0.065725 (5) | 0.193389 (5) | 0.369338 (4) | 0.006 | |
CO(1) | 0.227864 (3) | 0.248045 (3) | 0.488470 (2) | 0.005 | |
S(1) | 0.291052 (7) | 0.416212 (5) | 0.447714 (4) | 0.007 | |
S(2) | 0.417238 (6) | 0.143396 (5) | 0.468000 (5) | 0.007 | |
N(1) | 0.25467 (2) | 0.581262 (17) | 0.552539 (14) | 0.007 | |
N(2) | 0.38465 (2) | 0.451494 (17) | 0.621403 (13) | 0.007 | |
N(3) | 0.25896 (2) | −0.026455 (17) | 0.467010 (14) | 0.007 | |
N(4) | 0.37387 (2) | 0.005188 (18) | 0.334528 (14) | 0.009 | |
C(1) | 0.31177 (2) | 0.486789 (17) | 0.548005 (14) | 0.006 | |
C(2) | 0.14593 (3) | 0.61478 (2) | 0.487671 (18) | 0.010 | |
C(3) | 0.31738 (3) | 0.66558 (2) | 0.606841 (19) | 0.010 | |
C(4) | 0.48357 (2) | 0.36775 (2) | 0.615738 (17) | 0.009 | |
C(5) | 0.35058 (3) | 0.47832 (2) | 0.715017 (16) | 0.010 | |
C(6) | 0.34373 (2) | 0.032804 (18) | 0.419994 (14) | 0.006 | |
C(7) | 0.42557 (4) | 0.07976 (3) | 0.26928 (2) | 0.014 | |
C(8) | 0.38020 (3) | −0.10277 (2) | 0.30495 (2) | 0.013 | |
C(9) | 0.25285 (3) | −0.01970 (2) | 0.566503 (18) | 0.010 | |
C(10) | 0.14743 (3) | −0.08541 (2) | 0.42206 (2) | 0.011 | |
H(2A) | 0.185443 | 0.649532 | 0.427366 | 0.01446 (4) | |
H(2B) | 0.083429 | 0.670944 | 0.520498 | 0.01446 (4) | |
H(2C) | 0.082591 | 0.550549 | 0.466308 | 0.01446 (4) | |
H(3A) | 0.326641 | 0.732742 | 0.564212 | 0.01485 (5) | |
H(3B) | 0.420157 | 0.645138 | 0.628704 | 0.01485 (5) | |
H(3C) | 0.255544 | 0.684841 | 0.663556 | 0.01485 (5) | |
H(4A) | 0.517688 | 0.362317 | 0.546888 | 0.01316 (4) | |
H(4B) | 0.438584 | 0.295973 | 0.636167 | 0.01316 (4) | |
H(4C) | 0.567345 | 0.385347 | 0.663133 | 0.01316 (4) | |
H(5A) | 0.425587 | 0.527633 | 0.748571 | 0.01522 (5) | |
H(5B) | 0.346083 | 0.407632 | 0.753457 | 0.01522 (5) | |
H(5C) | 0.253263 | 0.51701 | 0.71454 | 0.01522 (5) | |
H(7A) | 0.535231 | 0.085761 | 0.274863 | 0.02038 (6) | |
H(7B) | 0.383004 | 0.154725 | 0.282164 | 0.02038 (6) | |
H(7C) | 0.394939 | 0.052736 | 0.20193 | 0.02038 (6) | |
H(8A) | 0.298395 | −0.123093 | 0.256903 | 0.01886 (5) | |
H(8B) | 0.378404 | −0.152898 | 0.363578 | 0.01886 (5) | |
H(8C) | 0.475662 | −0.112971 | 0.272135 | 0.01886 (5) | |
H(9A) | 0.350365 | 0.004347 | 0.595796 | 0.01456 (4) | |
H(9B) | 0.226518 | −0.09468 | 0.59216 | 0.01456 (4) | |
H(9C) | 0.17886 | 0.036548 | 0.586479 | 0.01456 (4) | |
H(10A) | 0.128457 | −0.058042 | 0.35323 | 0.01689 (5) | |
H(10B) | 0.055374 | −0.075345 | 0.458827 | 0.01689 (5) | |
H(10C) | 0.170348 | −0.166558 | 0.41867 | 0.01689 (5) |
U11 | U22 | U33 | U12 | U13 | U23 | |
BR(1) | 0.00717 (2) | 0.00656 (2) | 0.00545 (2) | −0.000062 (15) | 0.000980 (15) | −0.000408 (15) |
BR(2) | 0.00690 (2) | 0.00600 (2) | 0.00615 (2) | −0.000521 (15) | −0.001268 (14) | −0.000667 (15) |
CO(1) | 0.005187 (11) | 0.004078 (12) | 0.004913 (11) | 0.000080 (10) | 0.000028 (9) | −0.000383 (10) |
S(1) | 0.010302 (18) | 0.004581 (18) | 0.005170 (16) | −0.000414 (14) | −0.000084 (13) | −0.000283 (13) |
S(2) | 0.005543 (16) | 0.005215 (18) | 0.009608 (19) | −0.000169 (13) | 0.000499 (13) | −0.001355 (15) |
N(1) | 0.00833 (5) | 0.00514 (5) | 0.00743 (5) | 0.00111 (4) | −0.00134 (4) | −0.00089 (4) |
N(2) | 0.00785 (5) | 0.00655 (5) | 0.00588 (5) | 0.00138 (4) | −0.00065 (4) | −0.00023 (4) |
N(3) | 0.00722 (5) | 0.00632 (6) | 0.00807 (5) | −0.00116 (4) | 0.00038 (4) | 0.00030 (4) |
N(4) | 0.01288 (6) | 0.00693 (6) | 0.00615 (5) | 0.00084 (5) | 0.00138 (5) | −0.00047 (4) |
C(1) | 0.00715 (5) | 0.00484 (6) | 0.00566 (5) | 0.00030 (4) | −0.00033 (4) | −0.00031 (4) |
C(2) | 0.01027 (6) | 0.00850 (7) | 0.00995 (7) | 0.00250 (6) | −0.00197 (5) | −0.00007 (6) |
C(3) | 0.01199 (7) | 0.00650 (7) | 0.01108 (7) | −0.00005 (6) | −0.00088 (6) | −0.00252 (6) |
C(4) | 0.00852 (6) | 0.00786 (7) | 0.00985 (6) | 0.00203 (5) | −0.00069 (5) | 0.00028 (6) |
C(5) | 0.01300 (7) | 0.01120 (8) | 0.00619 (6) | 0.00209 (6) | −0.00030 (5) | −0.00078 (6) |
C(6) | 0.00719 (5) | 0.00508 (6) | 0.00633 (5) | −0.00001 (4) | 0.00047 (4) | −0.00030 (4) |
C(7) | 0.01965 (10) | 0.01272 (10) | 0.00867 (7) | 0.00142 (8) | 0.00376 (7) | 0.00307 (7) |
C(8) | 0.01873 (10) | 0.00899 (8) | 0.00998 (7) | 0.00244 (7) | 0.00030 (7) | −0.00326 (6) |
C(9) | 0.01121 (7) | 0.00915 (8) | 0.00900 (6) | 0.00016 (6) | 0.00311 (5) | 0.00089 (6) |
C(10) | 0.00960 (7) | 0.00831 (8) | 0.01569 (9) | −0.00243 (6) | −0.00152 (6) | −0.00049 (7) |
BR(1)—CO(1) | 2.4028 (1) | CO(1)—S(1) | 2.3495 (1) |
BR(2)—CO(1) | 2.4160 (1) | CO(1)—S(2) | 2.3276 (1) |
BR(1)—CO(1)—BR(2) | 113.804 (2) | BR(2)—CO(1)—S(1) | 105.132 (2) |
BR(1)—CO(1)—S(1) | 108.884 (2) | BR(2)—CO(1)—S(2) | 103.895 (2) |
BR(1)—CO(1)—S(2) | 117.482 (2) | S(1)—CO(1)—S(2) | 106.744 (3) |
Acknowledgements
Affiliation with the Center for Integrated Materials Research (iMAT) at Aarhus University is gratefully acknowledged.
Funding information
The following funding is acknowledged: Villum Fonden (grant No. 12391 awarded to JO); Danscatt (grant awarded to JO); Science and Engineering Research Board (grant Nos. CRG/2019/004185, SPR/2019/001145 awarded to MS); Council of Scientific and Industrial Research, India (grant No. 01(2933)/18/EMR-II awarded to MS); Indian Institute of Technology Bombay (grant awarded to MS). The synchrotron experiment was performed on beamline BL02B1 at SPring-8 with the approval of the Japan Synchrotron Radiation Research Institute as a Partner User (proposal No. 2018B0078).
References
Abragam, A. & Bleaney, B. (2012). Electron Paramagnetic Resonance of Transition Ions. Oxford University Press. Google Scholar
Angeli, C., Cimiraglia, R. & Malrieu, J.-P. (2001). Chem. Phys. Lett. 350, 297–305. CrossRef CAS Google Scholar
Atanasov, M., Ganyushin, D., Sivalingam, K. & Neese, F. (2012). Molecular Electronic Structures of Transition Metal Complexes II, edited by D. M. P. Mingos, P. Day & J. P. Dahl, pp. 149–220. Berlin, Heidelberg: Springer. Google Scholar
Atanasov, M., Zadrozny, J. M., Long, J. R. & Neese, F. (2013). Chem. Sci. 4, 139–156. CrossRef CAS Google Scholar
Blessing, R. H. (1997). J. Appl. Cryst. 30, 421–426. CrossRef CAS Web of Science IUCr Journals Google Scholar
Bunting, P. C., Atanasov, M., Damgaard-Moller, E., Perfetti, M., Crassee, I., Orlita, M., Overgaard, J., van Slageren, J., Neese, F. & Long, J. R. (2018). Science, 362, eaat7319. CSD CrossRef PubMed Google Scholar
Carl, E., Demeshko, S., Meyer, F. & Stalke, D. (2015). Chem. Eur. J. 21, 10109–10115. CSD CrossRef CAS PubMed Google Scholar
Chibotaru, L. F. (2015). Molecular Nanomagnets and Related Phenomena, edited by S. Gao, pp. 185–229. Berlin, Heidelberg: Springer. Google Scholar
Chilton, N. F., Goodwin, C. A. P., Mills, D. P. & Winpenny, R. E. P. (2015). Chem. Commun. 51, 101–103. CSD CrossRef CAS Google Scholar
Clark, T., Hennemann, M., Murray, J. S. & Politzer, P. (2007). J. Mol. Model. 13, 291–296. Web of Science CrossRef PubMed CAS Google Scholar
Craven, M., Nygaard, M. H., Zadrozny, J. M., Long, J. R. & Overgaard, J. (2018). Inorg. Chem. 57, 6913–6920. CSD CrossRef CAS PubMed Google Scholar
Damgaard-Møller, E., Krause, L., Lassen, H., Malaspina, L. A., Grabowsky, S., Bamberger, H., McGuire, J., Miras, H. N., Sproules, S. & Overgaard, J. (2020a). Inorg. Chem. 59, 13190–13200. PubMed Google Scholar
Damgaard-Møller, E., Krause, L. & Overgaard, J. (2020b). 21st Century Challenges in Chemical Crystallography II, edited by D. Michael P. Mingos & P. R. Raithby, pp. 145–182. Cham: Springer International Publishing. Google Scholar
Damgaard-Møller, E., Krause, L., Tolborg, K., Macetti, G., Genoni, A. & Overgaard, J. (2020c). Angew. Chem. Int. Ed. 59, 21203–21209. Google Scholar
Daul, C. A., Day, P., Figgis, B. N., Gudel, H. U., Herren, F., Ludi, A. & Reynolds, P. A. (1988). Proc. R. Soc. London Ser. A, 419, 205–219. CAS Google Scholar
Figgis, B. N., Forsyth, J. B. & Reynolds, P. A. (1987). Inorg. Chem. 26, 101–105. CrossRef CAS Google Scholar
Figgis, B. N. & Hitchman, M. A. (2000). Ligand Field Theory and its Applications. New York: Wiley-VCH. Google Scholar
Figgis, B. N., Reynolds, P. A. & Forsyth, J. B. (1988). J. Chem. Soc. Dalton Trans. 640, 117–122. CrossRef Google Scholar
Figgis, B. N., Reynolds, P. A. & Mason, R. (1982). Proc. R. Soc. London Ser. A, 384, 49–55. CAS Google Scholar
Figgis, B. N., Reynolds, P. A. & Mason, R. (1983). J. Am. Chem. Soc. 105, 440–443. CrossRef CAS Google Scholar
Forsyth, J. B. (1977). Phys. Scr. 15, 69–79. CrossRef CAS Google Scholar
Gao, C., Genoni, A., Gao, S., Jiang, S., Soncini, A. & Overgaard, J. (2020). Nat. Chem. 12, 213–219. CSD CrossRef CAS PubMed Google Scholar
Gatteschi, D., Sessoli, R. & Villain, J. (2006). Molecular Nanomagnets. New York: Oxford University Press. Google Scholar
Gatti, C. (2005). Z. Kristallogr. 220, 399–457. Web of Science CrossRef CAS Google Scholar
Genoni, A. (2020). Acta Cryst. A76, 172–179. CrossRef IUCr Journals Google Scholar
Gomez-Coca, S., Cremades, E., Aliaga-Alcalde, N. & Ruiz, E. (2013). J. Am. Chem. Soc. 135, 7010–7018. CAS PubMed Google Scholar
Goodwin, C. A. P., Ortu, F., Reta, D., Chilton, N. F. & Mills, D. P. (2017). Nature, 548, 439–442. CSD CrossRef CAS PubMed Google Scholar
Gukasov, A. & Brown, P. (2010). J. Phys. Condens. Matter, 22, 502201. CrossRef PubMed Google Scholar
Gukasov, A. & Brown, P. J. (2002). J. Phys. Condens. Matter, 14, 8831–8839. Web of Science CrossRef CAS Google Scholar
Gukasov, A., Goujon, A., Meuriot, J.-L., Person, C., Exil, G. & Koskas, G. (2007). J. Phys. Condens. Matter, 397, 131–134. CAS Google Scholar
Guo, F. S., Day, B. M., Chen, Y. C., Tong, M. L., Mansikkamäki, A. & Layfield, R. A. (2017). Angew. Chem. Int. Ed. 56, 11445–11449. CSD CrossRef CAS Google Scholar
Guo, F. S., Day, B. M., Chen, Y. C., Tong, M. L., Mansikkamäki, A. & Layfield, R. A. (2018). Science, 362, 1400–1403. CSD CrossRef CAS PubMed Google Scholar
Gupta, S. K., Rajeshkumar, T., Rajaraman, G. & Murugavel, R. (2016). Chem. Sci. 7, 5181–5191. CSD CrossRef CAS PubMed Google Scholar
Hansen, N. K. & Coppens, P. (1978). Acta Cryst. A34, 909–921. CrossRef CAS IUCr Journals Web of Science Google Scholar
Herbst-Irmer, R., Henn, J., Holstein, J. J., Hübschle, C. B., Dittrich, B., Stern, D., Kratzert, D. & Stalke, D. (2013). J. Phys. Chem. A, 117, 633–641. Web of Science CAS PubMed Google Scholar
Holladay, A., Leung, P. & Coppens, P. (1983). Acta Cryst. A39, 377–387. CrossRef CAS Web of Science IUCr Journals Google Scholar
Kahn, O. (1993). Molecular Magnetism. New York: Wiley-VCH. Google Scholar
Kibalin, I. A. & Gukasov, A. (2019). Phys. Rev. Res. 1, 033100. CrossRef Google Scholar
Klahn, E. A., Gao, C., Gillon, B., Gukasov, A., Fabrèges, X., Piltz, R. O., Jiang, S. D. & Overgaard, J. (2018). Chem. Eur. J. 24, 16576–16581. CSD CrossRef CAS PubMed Google Scholar
Krause, L., Herbst-Irmer, R., Sheldrick, G. M. & Stalke, D. (2015). J. Appl. Cryst. 48, 3–10. Web of Science CSD CrossRef ICSD CAS IUCr Journals Google Scholar
Legendre, C., Damgaard-Moller, E., Overgaard, J. & Stalke, D. (2021). EurJIC, 2021, 3108–3114 CAS Google Scholar
Malmqvist, P.-Å. & Roos, B. O. (1989). Chem. Phys. Lett. 155, 189–194. CrossRef CAS Google Scholar
Meindl, K. & Henn, J. (2008). Acta Cryst. A64, 404–418. Web of Science CrossRef CAS IUCr Journals Google Scholar
Nathans, R., Shull, C. G., Shirane, G. & Andresen, A. (1959). J. Phys. Chem. Solids, 10, 138–146. CrossRef CAS Web of Science Google Scholar
Neese, F. (2012). WIREs Comput. Mol. Sci. 2, 73–78. Web of Science CrossRef CAS Google Scholar
Neese, F. (2018). WIREs Comput. Mol. Sci. 8, e1327. Google Scholar
Neese, F., Atanasov, M., Bistoni, G., Maganas, D. & Ye, S. (2019). J. Am. Chem. Soc. 141, 2814–2824. Web of Science CrossRef CAS PubMed Google Scholar
Pinsky, M. & Avnir, D. (1998). Inorg. Chem. 37, 5575–5582. Web of Science CrossRef PubMed CAS Google Scholar
Politzer, P. & Murray, J. S. (2019). Crystals, 9, 165. Web of Science CrossRef Google Scholar
Rechkemmer, Y., Breitgoff, F. D., van der Meer, M., Atanasov, M., Hakl, M., Orlita, M., Neugebauer, P., Neese, F., Sarkar, B. & van Slageren, J. (2016). Nat. Commun. 7, 10467. Web of Science CSD CrossRef PubMed Google Scholar
Reta, D., Kragskow, J. G. C. & Chilton, N. F. (2021). J. Am. Chem. Soc. 138, 5441–5450. Google Scholar
Ridier, K., Gillon, B., Gukasov, A., Chaboussant, G., Cousson, A., Luneau, D., Borta, A., Jacquot, J.-F., Checa, R., Chiba, Y., Sakiyama, H. & Mikuriya, M. (2016). Chem. Eur. J. 22, 724–735. CSD CrossRef CAS PubMed Google Scholar
Rinehart, J. D. & Long, J. R. (2011). Chem. Sci. 2, 2078–2085. Web of Science CrossRef CAS Google Scholar
Ruamps, R., Batchelor, L. J., Maurice, R., Gogoi, N., Jiménez-Lozano, P., Guihéry, N., de Graaf, C., Barra, A.-L., Sutter, J.-P. & Mallah, T. (2013). Chem. Eur. J. 19, 950–956. CrossRef CAS PubMed Google Scholar
Schäfer, A., Horn, H. & Ahlrichs, R. (1992). J. Chem. Phys. 97, 2571–2577. CrossRef Web of Science Google Scholar
Schäfer, A., Huber, C. & Ahlrichs, R. (1994). J. Chem. Phys. 100, 5829–5835. Google Scholar
Sheldrick, G. M. (2015a). Acta Cryst. A71, 3–8. Web of Science CrossRef IUCr Journals Google Scholar
Sheldrick, G. M. (2015b). Acta Cryst. C71, 3–8. Web of Science CrossRef IUCr Journals Google Scholar
Sigrist, M. L. W., Tregenna-Piggott, P. S., Pedersen, K. A., Sørensen, M., Barra, A.-L., Hauser, J., Liu, S.-X., Decurtins, S., Mutka, H. & Bendix, J. (2015). Eur. J. Inorg. Chem. 2015, 2683–2689. CSD CrossRef CAS Google Scholar
Thomsen, M. K., Nyvang, A., Walsh, J. P. S., Bunting, P. C., Long, J. R., Neese, F., Atanasov, M., Genoni, A. & Overgaard, J. (2019). Inorg. Chem. 58, 3211–3218. Web of Science CSD CrossRef CAS PubMed Google Scholar
Tripathi, S., Vaidya, S., Ahmed, N., Andreasen Klahn, E., Cao, H., Spillecke, L., Koo, C., Spachmann, S., Klingeler, R., Rajaraman, G., Overgaard, J. & Shanmugam, M. (2021). Cell. Rep. Phys. Sci. 2, 100404. CrossRef Google Scholar
Tripathi, S., Vaidya, S., Ansari, K. U., Ahmed, N., Rivière, E., Spillecke, L., Koo, C., Klingeler, R., Mallah, T., Rajaraman, G. & Shanmugam, M. (2019). Inorg. Chem. 58, 9085–9100. CSD CrossRef CAS PubMed Google Scholar
Vaidya, S., Shukla, P., Tripathi, S., Rivière, E., Mallah, T., Rajaraman, G. & Shanmugam, M. (2018). Inorg. Chem. 57, 3371–3386. CSD CrossRef CAS PubMed Google Scholar
Vaidya, S., Singh, S. K., Shukla, P., Ansari, K., Rajaraman, G. & Shanmugam, M. (2017). Chem. Eur. J. 23, 9546–9559. CSD CrossRef CAS PubMed Google Scholar
Vaidya, S., Tewary, S., Singh, S. K., Langley, S. K., Murray, K. S., Lan, Y., Wernsdorfer, W., Rajaraman, G. & Shanmugam, M. (2016). Inorg. Chem. 55, 9564–9578. CSD CrossRef CAS PubMed Google Scholar
Vaidya, S., Upadhyay, A., Singh, S. K., Gupta, T., Tewary, S., Langley, S. K., Walsh, J. P. S., Murray, K. S., Rajaraman, G. & Shanmugam, M. (2014). Chem. Commun. 51, 3739–3742 . CSD CrossRef Google Scholar
Volkov, A., Macchi, P., Farrugia, L. J., Gatti, C., Mallinson, P. R., Richter, T. & Koritsanszky, T. (2006). XD2006. University at Buffalo, State University of New York, NY, USA; University of Milan, Italy; University of Glasgow, UK; CNRISTM, Milan, Italy; Middle Tennessee State University, TN, USA; and Freie Universität, Berlin, Germany. Google Scholar
Weigend, F. & Ahlrichs, R. (2005). Phys. Chem. Chem. Phys. 7, 3297–3305. Web of Science CrossRef PubMed CAS Google Scholar
Zadrozny, J. M. & Long, J. R. (2011). J. Am. Chem. Soc. 133, 20732–20734. CrossRef CAS PubMed Google Scholar
Zadrozny, J. M., Telser, J. & Long, J. R. (2013). Polyhedron, 64, 209–217. CSD CrossRef CAS Google Scholar
This is an open-access article distributed under the terms of the Creative Commons Attribution (CC-BY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.