The Intricate Determination of Magnetic Anisotropy in Quasi-octahedral Vanadium(III): An HF-EPR and Magnetic Study

We report here the synthesis and a preliminary characterization of the tetranuclear complex of formula [Ga3V(LEt)2(dpm)6], Ga3VEt, in which H3LEt = 2-Ethyl-2-(hydroxymethyl)-propane-1,3-diol and Hdpm = dipivaloylmethane, containing a single paramagnetic vanadium(III) center, from a structural, magnetic, and spectroscopic point of view. Structural characterization by X-ray diffraction evidenced that this derivative is isostructural with the star-shaped Single-Molecule Magnet [Fe3V(LEt)2(dpm)6], Fe3VEt, and can, thus, be considered a model to analyze the magnetic anisotropy of the vanadium(III) ion in that system. The observed results confirm the complexity in obtaining a rationalization of the magnetic behavior of this metal ion, with magnetization data and High Field Electron Paramagnetic Resonance (HF-EPR) spectroscopy providing apparently conflicting results. Indeed, the former were rationalized assuming a rhombic distortion of the ligand field and a dominant easy-axis type anisotropy (equivalent to D ≈ −14.1 cm−1, E ≈ 1.2 cm−1), while a simple axial Spin Hamiltonian approach could explain HF-EPR data (|D| ≈ 6.98 cm−1).


Introduction
The discovery in the early 90s of Single-Molecule Magnet (SMM) behavior, i.e., magnetic bistability at the molecular level in polynuclear transition metal complexes [1], sparked intense research in the field of molecular magnetism, aiming at the possible use of these systems as magnetic memory units [2]. The cause of the slow relaxation of the magnetization in these systems was soon attributed to the existence of an anisotropy barrier that amounts to the total energy spanned by the ground spin state, S, of the complexes [3]. For a system with an integer spin state and an easy axis of magnetization, this is equal to |D|S 2 , where D is the parameter describing the zero-field splitting (i.e., magnetic anisotropy) of the S multiplet (D < 0). When thermal energy is much lower than the height of the barrier, a system previously magnetized will retain its magnetization; on the contrary, on increasing temperature, thermal energy becomes large enough to allow fast relaxation of the magnetization, leading to a loss of the information stored in it. This prompted many groups to find ways to improve the magnetic anisotropy barrier rationally [4,5]. This is, however, by no means a simple task: in polynuclear clusters, it requires careful engineering of the anisotropic features of the constituent ions, the relative orientation of their magnetic anisotropy axis, and the exchange coupling among the centers [6]. In this context, one of the most accurately studied families is that of tetranuclear star-shaped complexes generally referred to as Fe 3 M′ (M′ = Fe, Cr, V, see Scheme 1) [7][8][9][10], which also proved very successful in nanostructuration processes [11][12][13].
For this family, some of us showed in the past that changing the central metal M′ allows playing with both factors affecting the magnetic anisotropy barrier height. By changing M′ from iron(III) to chromium(III) and vanadium(III), the global spin state increased from S = 5 to S = 6 to S = 13/2 [8,9]. At the same time, this modification affected the anisotropy parameter D, as a consequence of both the different singleion contribution of the central site and the different ratio between single-ion and Scheme 1 General structure of star-shaped Fe 3 M′ SMMs, depicting the tripodal alkoxide ligands above and below the plane defined by the metal ions. The wireframe colors are gray for carbon and red for oxygen atoms. The alkyl or aryl groups bound to the bridgehead carbon in each chelating alkoxide (Et, Me, Ph, and Bu t in different literature reports) are not shown (colour figure online) 1 3 The Intricate Determination of Magnetic Anisotropy in… global spin states. While the anisotropy in Fe 3 Cr and Fe 3 Fe could be rationalized entirely on the basis of single-ion and dipolar anisotropy as derived by the study of iron(III)-and chromium(III)-doped isostructural Ga 4 diamagnetic clusters [8], the situation remained unclear for Fe 3 V. At that time, Ga 3 V was not available, and we resorted to simple Ligand Field arguments to obtain a rough estimate of the singleion anisotropy of vanadium(III), which is an orbitally degenerate ion in octahedral symmetry ( 3 T 1g ). As deeply discussed in a series of paper by the late Tregenna-Piggott, trigonal distortion may provide either an orbitally non-degenerate ground state, with hard-axis anisotropy, or an orbitally degenerate ground state, whose anisotropic properties cannot be directly mapped onto a Spin Hamiltonian in the absence of further rhombic distortions [14][15][16]. This complicates the analysis of the origin of the anisotropy in exchange coupled systems containing vanadium(III) [17] and, in the case of Fe 3 V, a non-negligible contribution of anisotropic exchange was postulated [9]. The isostructural Ga 3 V derivative has now been obtained, and we report here some preliminary results concerning its structural, High Field Electron Paramagnetic Resonance (HF-EPR), and magnetic properties, which evidence the complexity of the rationalization of this metal ion and suggest that our previous estimate might be wrong.
Elemental analyses were carried out by Medac Laboratories Ltd. (Chobham, Surrey, UK). The equipment employed was a Thermal Scientific Flash EA 1112 Series Elemental Analyser (C and H contents) and a Varian Vista MPX ICP-OES spectrometer for Gallium and Vanadium contents. Infrared spectroscopy data presented as Supplementary Material (400-4000 cm −1 , Figure S1) were recorded from KBr pellets on BIORAD FTS3500GX equipment with a spectral resolution of 2 cm −1 .

3
The resulting white suspension was kept stirring at 0 °C for 5 h and was then filtered to give a quantitative yield of the tripodal lithium alkoxide Li 3 L Et . The white solid was dried under vacuum and immediately used as follows. Solid [VCl 3 (thf) 3 ] (280 mg, 0.749 mmol) was slowly extracted for 48 h into a boiling suspension of Li 3 L Et (221 mg, 1.45 mmol) in 1,2-dimethoxyethane (glyme, 120 mL). This produced a light-brown suspension that was added at room temperature to a white suspension of [Ga 2 (OMe) 2 (dpm) 4 ] (1.04 g, 1.11 mmol) in glyme (50 mL). The resulting mixture was heated under reflux for 5 h, cooled down to room temperature, and then filtered, producing a greenish-brown solid mixed with a white contaminant (possibly LiCl and LiOMe). This solid was extracted with thf, and the light yellow extract was submitted to vapor diffusion of glyme over 4 weeks. From the crystallization chamber, brownishgreen crystals (Ga 3 VEt-1) were isolated by filtration and dried under vacuum. Samples of Ga 3 VEt can be handled in air in the solid state for several days without detectable change. Their structural, spectroscopic, and magnetic characterization was performed in air. Crystals are greenish-brown rhombohedra. From a sample under oil, one, ca 0.46 × 0.40 × 0.23 mm, was mounted on a micro-mount and fixed in the cool nitrogen stream on a Bruker D8 Venture diffractometer equipped with a Photon 100 CMOS detector, Mo-Kα radiation, and graphite monochromator. Intensity data were measured at 273(2) K by thin-slice ω-and φ-scans. The total number of reflections recorded, to θ max = 25°, was 146490 of which 2651 were unique [R(int) = 0.056]; 2186 were 'observed′ with I > 2σ I .

Crystal Structure Analysis of Ga 3 VEt
Data were processed using the Bruker APEX2 software [21]. The structure was determined by direct methods in the SHELXS program [22] and refined by full-matrix least-squares methods, on F 2 ′s, in SHELXL [23]. Scattering factors for neutral atoms were taken from the literature [24]. All computer programs used in these analyses were run through WinGX [25,26] at the Universidade Federal do Paraná and the University of East Anglia. Crystal data, data collection, and general structure refinement data are presented in Table S1.

3
The Intricate Determination of Magnetic Anisotropy in… Figure S2 presents the first coordination sphere of the central (vanadium) ion, while Figures S3 and S4 show different views of the tetranuclear structure, with the atom numbering scheme. The asymmetric unit of [V{(OCH 2 ) 3 CEt} 2 {Ga(dpm) 2 } 3 ] (Figure S5) contains two different metal centers in special positions: a Ga III ion that lies on a twofold axis and a V III ion lying on several symmetry elements: one threefold axis (parallel to c), three twofold axes (one parallel to b, the second one parallel to a, and the third bisecting the angle between the first two), and three glide planes. The asymmetric unit ( Figure S5) also contains a complete β-diketonate (dpm − ) ligand coordinated to Ga III , and one-third of a tripodal alkoxide ligand (including one of the disordered orientations of its ethyl group) bridging the vanadium and gallium ions through the O3 donor atom.

Notes on the Structure
A high degree of disorder affects the tert-butyl groups of the β-diketonate ligands. In the asymmetric unit, the disorder involving the carbon atoms attached to C4 (C1-C3) was modeled with two orientations and occupancy factors of 0.56 and 0.44. In the minor component, the C100, C200, and C300 atoms were refined isotropically; these were the only non-hydrogen atoms refined in this way. Disorder was also evident for C15 in the ethyl group of the tripodal alkoxide; this atom was refined with an occupancy factor of 0.33 as a consequence of the threefold symmetry axis passing through C14 and C13. The distance restraint DFIX was applied to all disordered atoms mentioned above: C1, C2, C3, C15, C100, C200, and C300. All B alerts in the CheckCIF file refer to these disordered portions of the asymmetric unit.
Hydrogen atoms were included in idealized positions and were refined isotropically; their U iso values were set to ride on the U eq values of the parent carbon atoms.
No hydrogen bond or other intra/intermolecular interaction of considerable strength was found in the structure.
Replacement of the central vanadium atom by gallium was investigated, but no improvement was seen; the best refinement was reached with a Ga 3 V system.

Magnetic Characterization
The magnetization of the sample was measured using an MPMS SQUID magnetometer (Quantum Design, Inc.). Microcrystalline powders were pressed in a pellet to avoid preferential orientation of the crystallites due to magnetic torque. Raw data were corrected for the sample holder contribution, measured in the same field and temperature range, and for the intrinsic diamagnetism of the sample, estimated using Pascal′s constants [27]. Fits and simulations of magnetic data were obtained using PHI [28].

High Field Electron Paramagnetic Resonance (HF-EPR)
HF-EPR measurements were performed on a home-built multi-frequency spectrometer operating in a double-pass configuration. A 110 GHz frequency source (Virginia Diodes Inc.) alone or with multipliers up to the third harmonic), as well as 95 and 115 GHz Gunn oscillators (Radiometer Physics GmbH) together with multipliers up to the fourth harmonic were used. The main magnetic field was supplied by a 16 T superconducting magnet associated with a Variable Temperature Insert (Cryogenic). A hot electron InSb bolometer (QMC Instruments) was used to measure the light transmitted through the sample. The spectrometer is designed to operate in the perpendicular mode, which means that the oscillating microwave field B mw is perpendicular to the static B dc and the selection rule for the transition is thus ΔM S = ± 1. The measurements were done on the microcrystalline powder sample pressed into a pellet in order to limit torqueing effects (actually, the same pellet of ca. 12 mg which was used for SQUID measurements), in the 190-460 GHz frequency range at 5 and 15 K. Spectral simulation and analysis were performed using EasySpin [29].

Synthesis and Structure
Previous works described several propeller-like complexes with the general formula [M 4 (L R ) 2 (dk) 6 ], in which Hdk is a β-diketone and H 3 L R is a tripodal alcohol, containing different alkyl or aryl (R) groups in the tripodal ligand [8,10,19,30]. Synthetically speaking, the nature of the R group influences the solubility of the intermediate species in the three-step procedure (see below), and of the final complex, modifying the reaction yield. The choice of ethyl as the R group in the tripodal alkoxide was made in this case to allow comparison with the heterometallic M 3 M′ species available in the literature, and in particular with Fe 3 V.
The complex [Ga 3 V(L Et ) 2 (dpm) 6 ] (Ga 3 VEt), in which H 3 L Et = 2-ethyl-2-(hydroxymethyl)-propane-1,3-diol and Hdpm = dipivaloylmethane, was prepared by the three-step route developed earlier for the isostructural Fe 3 M′Et compounds in which M′ = Cr III and V III [9,19]: In agreement with our previous findings [9], step 2 was confirmed to be crucial to the success of the synthesis and isolation of a chemically pure heterometallic M 3 M′ species. This step leads to the formation of the central portion ("core", formulated as V(L Et ) 3− 2 , Fig. 1, left) of the star-shaped compound, around which the propeller "blades", Ga(dpm) 2 + , assemble themselves in the last reaction step. We found that, for each tetra(hetero)nuclear compound, different reaction conditions had to be adjusted to maximize the global reaction yield, depending on the chemical nature of the central ion. As far as the preparation of Ga 3 VEt is concerned, longer and slower extraction of the vanadium(III) starting material into the boiling lithium alkoxide suspension in step 2 gives the best yields of the product, as described in the Experimental section. The addition of the gallium(III) dimer [Ga 2 (OMe) 2 (dpm) 4 ] to the suspension of the "core" in dimethoxyethane (step 3), followed by heating under reflux, gives a greenish-(light)brown solution from which microcrystals of Ga 3 VEt and a white powder, probably a mixture of LiCl and LiOMe, are isolated by filtration. This impure product is then recrystallized by vapor diffusion of dimethoxyethane into a thf solution, producing beautiful greenish-(light)brown rhombohedra.
Despite the well-formed appearance of the crystals, they produce weak X-ray diffraction patterns, even at low temperature, because of static and dynamic disorder in the ethyl group of the tripodal alkoxides and in the tert-butyl groups of the β-diketonate (dpm − ) ligands. This disorder is common to several members of the star-shaped M 4 family reported to date [8,30,32], and, in the case of Ga 3 VEt, it hampered the obtainment of an accurate X-ray structural description of the molecule. After many attempts involving distinct synthesis and recrystallization conditions, we have been able to confirm that the tetranuclear molecules found in the crystals of Ga 3 VEt (Figs. 1 right, S3 and S4) are isostructural, as expected, with the 'star' structures reported for Fe 3 CrEt [19] and Fe 3 VEt [9] in the same trigonal space group R-3c:h, and that the bond distances and angles in all three complexes are similar (Tables 1, S1-S3). The lower quality of the X-ray diffraction for Ga 3 VEt compared to Fe 3 CrEt [19] and Fe 3 VEt [9] provides a higher incertitude in the geometrical parameters. However, it allows us to establish the shape and symmetry of the molecule.
Further evidence in favor of the composition of the crystalline product came from elemental analysis (see Experimental) and spectroscopic data, as described below  Figures S3 and S4 and in the Supplementary Material. Metal analysis results indicate ca. 95% chemically pure samples as far as the Ga 3 V formulation is concerned. This is particularly relevant considering that the modular synthesis described in steps 1-3 competes with the self-assembly of homometallic species (in this case, Ga 4 Et), as discussed earlier [8,9,19]. This success in the desired synthesis of the heterometallic Ga 3 V entity probably relates both to the chelate effect stabilizing the vanadium "core" assembled in step 2, and the similar labilities of vanadium(III) and gallium(III) towards ligand substitution reactions [33][34][35][36][37]; as a consequence, none of the homometallic Ga 4 or V 4 structures is particularly favored when compared to the heterometallic species. The loss of vanadium(III) from the "core" to give ca. 5% Ga 4 (see elemental analysis results) could come from the lower affinity of the softer V III ion (as compared to Ga III ) towards the hard oxygen-donor atoms of the tripodal alkoxide [38,39]. Figures 1 right, S3, and S4 confirm the star-shaped, propeller-like molecular structure of a new component of the [M 3 M′(L Et ) 2 (dpm) 6 ] series. In addition, data presented in Tables 1, S2 and S3 give evidence in favor of the desired M = Ga III and M′ = V III identities for the peripheral and central metal ions respectively, and reassert that the synthetic route summarized in steps 1-3 is an efficient pathway to pure heterometallic M 3 M′ species.
In Table 1 (last column), all M central and M peripheral bond distances to the O-donor atoms of the tripodal alkoxide and β-diketonate ligands support a larger

3
The Intricate Determination of Magnetic Anisotropy in… ion (V III ; effective atomic radius 0.640 Å) at the center of the molecular "core", and a smaller cation (Ga III , effective atomic radius 0.620 Å) in the peripheral positions. This size ordering is the inverse of that observed for the Cr III /Fe III pairs in Fe 3 CrEt [19], which is also in agreement with the ionic radii magnitudes. During structure refinement, attempts to replace the central vanadium(III) ion in Ga 3 VEt with gallium(III) to check for the possible formation of Ga 4 Et led to higher R-factors and significant change (increase or decrease) in the anisotropic thermal parameters, indicating the worsening of the structural model. This result supports the presence of vanadium in the molecular "core". Taking the small difference in temperature of data collection into consideration, data in Tables 1 and S2, together with unit cell volumes [13730.5(8) Å 3 for Fe 3 CrEt; 13740.4(5) Å 3 for Fe 3 VEt [9,19] and 13514.9(9) Å 3 for Ga 3 VEt], indicate no large differences in the molecular shape or global size on changing the nature of the central/peripheral metal ions across the fourth period of the periodic classification. These changes are well known to be crucial, however, for the easyor hard-axis nature of the magnetic anisotropy in the entire molecular assembly, as well as to the height of the magnetization reversal barrier [8,30].
An interesting observation can be made at a more localized level: all bond angles involving the peripheral ions are significantly closer to the octahedral values in Ga 3 VEt than in Fe 3 CrEt or Fe 3 VEt (Table S2), leading to a more regular geometry for the peripheral units. This distinguishes Ga 3 VEt from the other two systems in which iron(III) ions are bound to the diketonate ligands. Additionally, the observation of the bond angles about the central atom in Table S2 shows an interesting similarity between the two 'V-centered stars′, lending further support to the  ,  Tables S2 and S3], resulting in the trigonally distorted coordination sphere, as illustrated in Figure S2. The helical pitch angle, defined as the dihedral angle between the M 4 and the M central (O) 2 M peripheral planes, is larger in Ga 3 VEt than in Fe 3 CrEt; this is consistent with "propeller blades" that are more tilted relative to the C 3 axis, and a greater extent of structural distortion involving the central ion [30]. The pitch angles shown in Table 1 fit the trend to increasing values on going from Cr-to V-and then to Fe-centered stars, [M 3 M′(L R ) 2 (dpm) 6 ], a trend that has been correlated with increasing axial zero-field splitting parameters |D| in this family of tetranuclear complexes [19,30,42].
In the unit cell, the triangular molecules are arranged in pairs whose components are related by an inversion center ( Figure S6). As a consequence of the high crystal symmetry, although each tetranuclear molecule is chiral (Λ or Δ enantiomeric forms), the crystals themselves are non-chiral [43,44]. Additionally, each central V III ion of the pair of molecules sits on four crystallographic proper rotation axes, one threefold (main) and three twofold (perpendicular to and related by the threefold symmetry axis) that, in turn, bisect the propeller blades. Symmetry requirements determine the threefold positional disorder of the ethyl groups in the tripodal ligands, and the resulting (disordered) arrangement has higher symmetry than the isolated molecules themselves.

Magnetic Properties
The temperature dependence of the magnetic susceptibility of Ga 3 VEt is reported in the form of χT vs. T plot in Fig. 2.
At room temperature, the observed χT value is 1.04 cm 3 K mol −1 , in agreement with the expected value for S = 1 with g close to 2. However, on cooling, the χT vs. T product steadily decreases to reach 0.8 cm 3 K mol −1 at 40 K. Below this temperature, the χT value drops to 0.35 cm 3 K mol −1 at 2 K. It is evident that such behavior cannot be rationalized on the basis of a simple S = 1 Spin Hamiltonian: Indeed, this model is unable to reproduce the linear increase of χT in the hightemperature range. This behavior is indicative of a partially quenched orbital contribution to the magnetic moment of Ga 3 VEt which, together with spin-orbit coupling (SOC) results in close-lying excited states. The latter are effective in inducing a sizeable temperature independent paramagnetism (TIP) contribution and, thus, a linear behavior of χT for T > 50 K [45]. Indeed, in pure octahedral symmetry, vanadium(III) is a 3 T 1g state (having two electrons in the t 2g d orbital set); this can be split by the combined action of SOC and geometrical distortion of the ligand field from purely octahedral symmetry. This can be described using the following Hamiltonian: where σ is the combined orbital reduction factor and λ is the SOC constant, while B 0 2 and B 2 2 are the axial and rhombic distortion parameters of the ligand field. As largely discussed in the literature [15,46] and anticipated in the Introduction, when B 0 2 is positive, the T 1g orbital state is split in such a way that a singlet orbital ground state is lying lowest, and the ensuing SOC provides an essentially isolated zero-field-split spin triplet state (see Fig. 3). The situation is different when B 0 2 is negative, since this provides an E orbital ground state; while this cannot be mapped on a Spin Hamiltonian, it can be further split by The Intricate Determination of Magnetic Anisotropy in… rhombic distortion, providing two orbital singlet states, close in energy, each of them with rhombic anisotropy. It is evident, however, that the rhombic distortion is only possible in cases where the point symmetry of the vanadium(III) ion is lower than axial. Magnetic data, including isothermal magnetization curves, were fitted using (2), as implemented in PHI [28]. The best-fit parameters, which provide a perfect simultaneous reproduction of both M vs. H and χT vs. T curves, were: B 0 2 = −439 ± 1 cm −1 , B 2 2 =148.4 ± 0.3 cm −1 , and = −0.734 ± 0.003 with λ fixed to 75 cm −1 . With these values, there is a ground singlet with two close-lying states at 1.67 and 14.97 cm −1 ; further excited states are at 173.9, 188.2, and 195.8 cm −1 , and at 812, 813, and 817 cm −1 . It is clear that the levels are now grouped in sets of pseudotriplets: the ground one has essentially easy-axis type, with an effective zero-field splitting of ca. −14.1 cm −1 . Accordingly, Hamiltonian (1) could only be employed to model the low-temperature magnetic properties, i.e., the isothermal magnetization curves, providing as best-fit values D = −14.1 ± 0.3 cm −1 , E = 1.21 ± 0.02 cm −1 , and g = 1.72 ± 0.01. While the latter value may appear quite low, it is very close to the value derived from HF-EPR of Fe 3 VEt [9]. However, this approach is clearly not reproducing the thermal behavior of the magnetic susceptibility above 50 K (see dotted line in Fig. 3), which is much influenced by the presence of low lying excited states.

HF-EPR Spectroscopy
HF-EPR spectra, measured at frequencies between 190 and 460 GHz and at 5 K, are reported in Fig. 4. They clearly show a main transition that moves toward zero field on decreasing frequency. No spectrum is observed below 210 GHz, indicating a quite large zero-field splitting. A second weaker transition is observed at a somewhat lower field and disappears for frequencies lower than 285 GHz. No other identifiable transitions are observed in the investigated field range. Measuring at 5 K in this frequency range allows getting a high polarization of the spin system with mostly only the lowest energy level populated (the energy of a frequency of 210 GHz corresponds to the one of a ~ 11 K temperature). Furthermore, the intensity of both transitions decreases when increasing the temperature to 15 K, indicating that they are issued from the lowest energy levels.
By plotting the resonance field of the main transition as a function of frequency, a behavior that is linear at high fields and parabolic at low fields is observed (Fig. 5). This, together with the lineshape, clearly identifies the transition as a perpendicular one, while the resonant field points to a formally forbidden one. If we assume a Spin Hamiltonian description of our system, the observed transition is the one between the levels labeled as m S = |0 > and m S = | + 1 > in zero field (m S = | + 1 > and m S = | -1 > in applied field, thus formally forbidden). We can fit this dependence using the following relation [47]: This provides a nice agreement with the experimental data using the following best-fit parameters: D = 209.35 GHz (6.98 cm −1 ), g = 1.9196. It is well known [48,49] that if D/hν < ¾ the forbidden transition at half field can provide two signals in a powder spectrum, one corresponding to a canonical orientation of the field, the other one to an off-axis resonant field.
We then plot in Fig. 6 the angular dependence of the resonant fields for two of the experimentally used frequencies, calculated on the basis of Spin Hamiltonian (1) and the D value obtained by the fit to Eq. (3). This clearly points to the presence of one canonical maximum and one minimum which are observed at 224.4 GHz, while one minimum (non-canonical, θ = 55.5°) and two (canonical) maxima are observed at 460 GHz. It is, however, worth noting that the θ = 90° transition has zero probability for perpendicular polarization of the microwave radiation (i.e., the experimental setup which we used, selecting ΔM S = ± 1 transitions), whereas x,y 2 H 2 + D 2 . Fig. 6 Angular dependence of the HF-EPR resonance fields of Ga 3 VEt at two frequencies in the zx plane, using Spin Hamiltonian (1) and parameters obtained by best fit of (3). Size of the symbols is proportional to the transition probability its probability is highest for parallel polarization (i.e., with oscillating B mw parallel to B dc , selecting ΔM S = 0 transitions). On the contrary, the θ = 0° transition has zero probability in parallel polarization and maximum probability for perpendicular polarization (see Fig. 6). By considering the lineshape of the observed transitions and their resonance fields, the observed spectrum should then be assigned to the parallel polarization, which also results in zero probability for the "allowed" transitions. With this approach, reasonable simulations could be obtained for the experimental spectra at different frequencies (see Figure S7). Unfortunately, due to the low sensitivity of the observed "forbidden" transition to the sign of D, no conclusion can be reached on this point. Finally, it is worth noting that assuming the same model used to rationalize magnetic data did not lead to reasonable reproduction of the experimental spectra.

Discussion
The results reported in the previous sections deserve some comments. It is first to be noted that the analysis of magnetic and HF-EPR data provides conflicting indications: for the first set of techniques, the results can only be rationalized by assuming that vanadium(III) in Ga 3 VEt has easy-axis type anisotropy, and its global behavior cannot be traced back to a simple spin triplet, but rather to a 3 E ligand field state partially split by rhombic distortion. The presence of a rhombic term in the Hamiltonian is not at odds with the observed trigonal crystallographic symmetry. Indeed, the latter is only a consequence of the statistical disorder of the ethyl pendant on the tris-alkoxo ligand, while the molecular symmetry of Ga 3 VEt is not axial. The results obtained by the magnetic characterization provides a magnitude of single-ion anisotropy for V III in the star-shaped compound which would allow rationalizing the global anisotropy of Fe 3 VEt without invoking any anisotropic exchange, at variance with what we reported in ref. [9]. Indeed, in the strong exchange approximation, the Zero-Field Splitting tensor of Fe 3 VEt is related to the microscopic anisotropic parameters (single-ion anisotropy tensor and spin-spin interactions) through: Here, d i ′s and d ij ′s are the projection coefficients on the global anisotropy of spin state S of the single-ion anisotropies (D Fe(i) , D V ) and of the spin-spin interaction (D Fe(i)-Fe(i) , D Fe(i)-V ), respectively, which can be calculated recursively [50] (see Table 2). A similar relation hold for the g tensors: It is now firmly established that peripheral Fe III ions in star-shaped complexes have hard-axis anisotropy with substantial rhombicity (D Fe ≈ + 0.7 cm −1 , E Fe ≈ 0.08 cm −1 ). Furthermore, different HF-EPR and magnetometry studies pointed out that the easy-axis direction of the anisotropy tensor of Fe III is approximately perpendicular to the plane of the star and that the local intermediate axis is oriented along the Fe-M direction (which is often a symmetry imposed direction [8,[51][52][53]). If one substitutes the values and orientations expected for Fe III , the dipolar contribution calculated in the point dipole approximation and D V = −14.1 cm −1 , which models the splitting in the low-energy region for vanadium(III), the estimated anisotropy for the ground S = 13/2 state of Fe 3 VEt from Eq. (4) is −0.290 cm −1 , in good agreement with the experimentally determined value in ref [9] (D(Fe 3 V Et ) = −0.314 cm −1 ). It is further to be noted that the experimentally determined [9] g value for Fe 3 VEt requires g V = 1.7 which is exactly the value one obtains by modeling the low-temperature magnetic data of Ga 3 VEt with a ground pseudotriplet (see above).
On the other hand, HF-EPR data can be rationalized assuming a standard spin triplet Spin Hamiltonian, resulting from a positive value of the trigonal distortion of the ligand field, and, thus, in an orbital singlet ground state. This would confirm our previous hypothesis, based essentially on the literature data [15,16,46], of positive zero-field splitting for vanadium(III) in the Fe 3 VEt cluster. It is, however, to be noted that appropriate simulation of the HF-EPR spectra required to assume, as discussed, parallel polarization of the microwave radiation. This is not the experimental setup used, thus raising further issues about these data. Similar problems with the radiation polarization were recently reported for an experimental study of a Hobased polyoxometalate (POM) single-molecule magnet [54].

Conclusions
We have reported a preliminary characterization of a Ga 3 VEt cluster both from a structural, magnetic, and spectroscopic point of view. The observed results do not allow to obtain a firm conclusion as to the anisotropy of central vanadium(III) in Ga 3 VEt and, thus, in the isostructural SMM Fe 3 VEt. Indeed, while the magnetic measurements point to an easy-axis type anisotropy with a relevant rhombicity, HF-EPR spectra can only be interpreted on the assumption of an easy-plane type anisotropy, of smaller magnitude but with some rhombicity. Also, given the small amount of sample available for HF-EPR measurements, some weak signals, associated with the allowed transitions and giving more stringent constraints on the anisotropy parameters, may have been missed. This confirms the difficulties in getting a reliable characterization of vanadium(III)-containing complexes. Despite that we have succeeded in stabilizing it with tripodal ligand and in exploiting its anisotropy in Fe 3 VEt with a triplication of magnetic blocking temperature compared to the archetypal Fe 4 SMM, a rational design of optimized SMMs based on this metal ion is still challenging. In perspective, it might be useful to extend the HF-EPR study to higher fields and frequencies using a larger amount of sample to obtain more sound conclusions. This would allow to obtain clearer indication as to the contribution of the anisotropic exchange to the anisotropy of this system. In this perspective single crystal techniques, such as cantilever torque magnetometry, which proved very powerful in the investigation of the lanthanide-based variant of tetranuclear star-shaped SMM [55], might also prove useful.