# Oxovanadium(IV) Thiocarboxylate Paddlewheels Containing Ancillary Group 10 Metals: A Comparative Study on Pd and Pt Derivatives

Olga Mironova, Giacomo Bellini, Alessio Nicolini, Manuel Imperato, Antonio

Ranieri, Marco Borsari, Matteo Briganti, Rodolphe Clérac, Mathieu

Rouzières, Enrico Salvadori, et al.

► To cite this version:

Olga Mironova, Giacomo Bellini, Alessio Nicolini, Manuel Imperato, Antonio Ranieri, et al. Oxovanadium(IV) Thiocarboxylate Paddlewheels Containing Ancillary Group 10 Metals: A Comparative Study on Pd and Pt Derivatives. Inorganic Chemistry, 2025, 65 (1), pp.638-651. 10.1021/acs.inorgchem.5c04835. [hal-05466821](https://hal.science/hal-05466821v1)

HAL Id: hal-05466821

https://hal.science/hal-05466821v1

Submitted on 20 Jan 2026

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_2.png" /></figure>

<figure><img src="image_3.png" /></figure>

# Oxovanadium(IV) Thiocarboxylate Paddlewheels Containing Ancillary Group 10 Metals: A Comparative Study on Pd and Pt Derivatives

Published as part of *Inorganic Chemistry* special issue "Quantum Materials from an Inorganic Chemistry Perspective".

Olga Mironova,* Giacomo Bellini, Alessio Nicolini, Manuel Imperato, Antonio Ranieri, Marco Borsari, Matteo Briganti, Rodolphe Clérac, Mathieu Rouzières, Enrico Salvadori, Maria Chiara Pagliero, Mario Chiesa, and Andrea Cornia*

<figure><img src="image_4.png" /></figure>

Cite This: *Inorg. Chem.* 2026, 65, 638–651

<figure><img src="image_5.png" /></figure>

Read Online

## ACCESS

<figure><img src="image_6.png" /></figure>

Metrics & More

<figure><img src="image_7.png" /></figure>

Article Recommendations

<figure><img src="image_8.png" /></figure>

Supporting Information

**ABSTRACT:** Vanadyl-containing paddlewheel structures have recently joined the pool of molecular spin systems showing respectable coherence times ($T_m$). We extended the investigation of [PtVO(SOCR)$_4$] (R = Me, Ph) by synthesizing the corresponding Pd derivatives. Crystals of [PdVO(SOCMe)$_4$] and [PdVO(SOCPh)$_4$]·DCM are isomorphous to their Pt analogues and comprise staggered and square dimers, respectively. For both Group 10 metals, M···M' contacts in staggered dimers promote a stronger antiferromagnetic coupling than M···S' contacts in square dimers, with Pd complexes showing a slightly more effective interaction consistent with DFT predictions. In solution, the paddlewheels exist in monomeric form and undergo a quasi-reversible one-electron reduction, more favorable for M = Pd than Pt and for R = Ph than Me. X-band EPR spectra in frozen CD$_2$Cl$_2$/toluene-$d_8$ show the eight-line hyperfine pattern characteristic of $^{51}$V ($I = 7/2$), with identical spin-Hamiltonian parameters across the series. Unlike the interaction with $^{195}$Pt ($I = 1/2$, 34%), the superhyperfine coupling with $^{105}$Pd ($I = 5/2$, 22%) is unresolved. $T_m$ values are essentially the same for Pt and Pd derivatives and follow the same substituent dependence (Ph > Me). Therefore, the ancillary Group 10 metal affects redox potentials and, to a lesser extent, solid-state magnetic properties, but leaves quantum coherence unchanged.

<figure><img src="image_9.png" /></figure>

## 1. INTRODUCTION

Lantern structures, also known as paddlewheels (PWs), are robust and versatile molecular complexes widely studied in metal–organic chemistry, e.g., in catalysis and the design of metal–organic frameworks.¹⁻³ Among the multitude of homo- or heterobimetallic PWs so far synthesized, examples of vanadyl (VO$^{2+}$)-based lanterns are surprisingly scarce. All of them belong to the class of heterobimetallic PWs containing, in addition to a 3d metal (Tr), a Group 10 metal (M), specifically Pt or Pd.⁴⁻⁶ In these structures, the two metals are bridged by bidentate ligands carrying different donor atoms, like (O,S), (N,S), or (N,O). Their binding mode aligns with the concept of soft–hard Lewis acids and bases (HSAB), where hard and soft donor atoms preferentially coordinate to Tr and M, respectively.⁷

For almost a decade, Prof. Doerrer's group has been advancing this approach and has applied it to the synthesis of heterobimetallic PWs using thioacetate and thiobenzoate (O,S)

ligands,⁶,⁸⁻¹³ with other groups recently joining the work.¹⁴⁻¹⁶ In the course of this time, many complexes with the general formula [PtTr(SOCR)$_4$(L)$_n$] (Tr = Cr, Mn, Fe, Ni, Co, and Zn; R = Me, Ph; L = axial ligand) have been comprehensively studied, and major regularities between their structure and magnetic properties have been established. Significantly, the Doerrer's group was the first to isolate two VO$^{2+}$ complexes [PtVO(SOCR)$_4$] with R = Me (1Pt) or R = Ph (2Pt), which were designed as metalloligands to enhance the performance of lanthanide complexes.¹⁷

Received: October 15, 2025

Revised: November 24, 2025

Accepted: December 5, 2025

Published: December 19, 2025

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>
We have recently highlighted the potential of these VO<sup>2+</sup>-based lanterns in the field of quantum information science.<sup>18</sup> Vanadyl complexes are, in general, well-established molecular spin qubits with consistently long coherence times (T<sub>m</sub>). Their S = 1/2 electronic spin is coupled with the nuclear spin (I = 7/2) of <sup>51</sup>V, a naturally abundant isotope (99.75%), yielding an electronuclear spin qudit.<sup>19–21</sup> Among strategies to enhance their performance is reducing spin–lattice relaxation effects using rigid ligand scaffolds.<sup>22–24</sup> Macrocyclic vanadyl complexes with phthalocyanines<sup>25,26</sup> and porphyrins,<sup>20,21,27</sup> for instance, have T<sub>m</sub> values surpassing those of bis-ligand systems like acetylacetonates.<sup>28–30</sup> Further sources of decoherence arising from neighboring electronic or nuclear spin carriers<sup>26,31–33</sup> can be mitigated by diluting the qubits in a diamagnetic host (e.g., a frozen solution or a crystalline matrix<sup>23,25</sup>) and utilizing O- or S-donor ligands (I = 0) as well as deuterated or nuclear-spin-free solvents.<sup>24,26,30,34</sup> Combining these approaches, we have recently isolated a quasi-macrocyclic bis(β-diketonato) vanadyl complex with T<sub>m</sub> = 13 µs in a frozen deuterated toluene-dichloromethane glass at 10 K.<sup>35</sup> 1Pt and 2Pt also contain a V<sup>4+</sup> ion coordinated solely by O donors and, in the same experimental conditions, exhibit quite respectable coherence times of 6 and 11 µs, respectively.<sup>18</sup> From a design perspective, the square planar coordination geometry of the d<sup>8</sup> Group 10 metal ion serves as an important structural element, imposing an idealized 4-fold symmetry to the PWs. Furthermore, the rigid structure of the M(SOCR)<sub>4</sub><sup>2−</sup> metalloligands offers an alternative to organic macrocyclic systems and enables multiple routes for qubit functionalization and processing without the direct involvement of the paramagnetic center, e.g., by chemi- or physisorption on surfaces.

In this study,<sup>36</sup> we continue our exploration of vanadyl PWs by focusing on Pd-based derivatives. While the influence of R substituents on the structure and properties of [PtTr-(SOCR)<sub>4</sub>(L)<sub>n</sub>] complexes has been extensively examined in various contexts, the impact of replacing Pt with Pd remains relatively unexplored.<sup>15</sup> Our initial idea was that a lighter Group 10 metal would reduce spin-orbit coupling (SOC) effects, conferring improved quantum coherence properties on vanadyl PWs. Transitioning from Pt to Pd presents a tricky synthetic challenge due to the dual coordination nature of Pd, which can interact with both soft and hard donor atoms.<sup>15,37–39</sup> By developing specialized synthetic protocols, we were able to isolate both [PdVO(SOCMe)<sub>4</sub>] (1Pd) and [PdVO(SOCPh)<sub>4</sub>] (2Pd), and to study their solid-state and solution behavior in direct comparison to the Pt predecessors. Importantly, these Pd PWs are isostructural to the corresponding Pt analogues and form isomorphic crystals containing dimeric assemblies supported by M···M' or M···S' contacts. Thus, they also provide a unique opportunity for directly investigating how the strength of magnetic coupling varies depending on the diamagnetic heavy metal ion involved. In solution, Pd derivatives exhibit distinct behavior characterized by an increased propensity for reduction. Moreover, their X-band EPR spectra in a frozen CD<sub>2</sub>Cl<sub>2</sub>/toluene-d<sub>8</sub> matrix lack the superhyperfine line splitting introduced by the <sup>195</sup>Pt isotope in 1Pt and 2Pt; however, their quantum coherence closely mirrors the trends observed for Pt counterparts, remaining largely unaffected by replacing Pt with the lighter Pd.

## 2. RESULTS AND DISCUSSION

**2.1. Synthesis.** In general, for the selective synthesis of [MTr(SOCR)<sub>4</sub>(L)<sub>n</sub>] complexes, the 3d metal salt is added to

the corresponding metalloligand [M(SOCR)<sub>4</sub>]<sup>2−</sup>, obtained in situ in water, which causes almost immediate precipitation of the desired product. However, the palladium analogues of 1Pt and 2Pt cannot be obtained by strictly following the procedures developed earlier.<sup>15</sup> The reason is that the higher lability of Pd over Pt in terms of the HSAB principle favors the formation of undesired products and also speeds up the reactions. We found that the target thioacetate derivative 1Pd cannot be obtained in water, presumably due to the prevailing formation of binary palladium(II) sulfides.<sup>15</sup> Instead, we could access 1Pd by mixing the reagents (K<sub>2</sub>PdCl<sub>4</sub>, 4 equiv of KSOCMe, and VOSO<sub>4</sub>·3.6H<sub>2</sub>O) in anhydrous methanol under an atmosphere of purified nitrogen. The reaction mixture turned into a brown suspension after 18 h of stirring, then the solvent was evaporated, and an emerald-green product was extracted with dichloromethane (DCM). The layering of this concentrated solution with *n*-hexane afforded dark green plates of 1Pd in 67% yield. As proof of the necessity of air-free and nonaqueous conditions, the synthesis failed when performed in degassed water under a N<sub>2</sub> atmosphere or in dry methanol in the air.

In contrast to the difficulties encountered in the synthesis of 1Pd, thiobenzoate 2Pd was accessed by only slightly modifying the “classical” procedure used for the Pt analogue. In fact, the greater “hardness” of Pd vs Pt causes a higher propensity to form homometallic lanterns with ligands coordinated through both soft and hard donors, as found in homoleptic [Pd<sub>2</sub>(pyt)<sub>4</sub>]<sup>37</sup> (Hpyt = pyridine-2-thiol) and heteroleptic [Pd<sub>2</sub>((C<sub>6</sub>H<sub>4</sub>)-PPh<sub>2</sub>)<sub>2</sub>(SOCPh)<sub>2</sub>]<sup>38</sup>. To reduce the formation of these homobimetallic species, a diluted solution of K<sub>2</sub>PdCl<sub>4</sub> in water was added dropwise over 10 min to a stirred aqueous solution of NaSOCPh, resulting in a clear orangish-red solution. The addition of VOSO<sub>4</sub>·3.6H<sub>2</sub>O dissolved in water caused the formation of an other precipitate, which was isolated by filtration. When this precipitate was treated with DCM, a red solution flowed down, leaving a grassy-green precipitate on the filter. The green product turned out to be poorly soluble in DCM, which is usually a good solvent for this type of PWs, but could be dissolved in hot THF. Upon cooling, green crystals formed and were structurally characterized as 2Pd·THF. In the red fraction, some green plates also formed upon standing, which were analyzed as a DCM solvate (2Pd·DCM). Even upon gentle drying, this solvate loses the lattice DCM completely. Larger crystals of both solvates can be grown by subjecting solutions of 2Pd to cycles of gentle heating and slow cooling.

Solutions of 2Pd in THF remain stable as long as they are kept under an inert atmosphere, but display signs of sensitivity to moisture upon standing in the air. For instance, a THF solution over crystals gradually turns brown in a few days, and signals attributed to diamagnetic species appear in the <sup>1</sup>H NMR spectrum (see below). Presumably, this is caused by an increase in the polarity of the solvent, which facilitates structural rearrangement. However, the slow degradation rate and the nonwetting<sup>40</sup> of the crude product allow for carrying out the synthesis in the air up to the recrystallization step. In contrast, solutions of 1Pd remain visibly unchanged after weeks of storage in the air.

The IR spectra of 1Pd and 2Pd closely mirror those of the corresponding Pt analogues, with differences in the peak positions not exceeding 10 cm<sup>-1</sup> (Figure S5). In particular, the band attributed to V=O stretching shifts from 984 in 1Pt to 995 cm<sup>-1</sup> in 1Pd, and from 989 in 2Pt to 1002 cm<sup>-1</sup> in 2Pd.

**2.2. X-ray Crystallography.** In a square planar coordination environment, Pd(II) and Pt(II) have very similar ionic radii
<table><caption>Table 1. Selected Bond Distances and Angles in 1Pd, 1Pt, 2Pd·DCM, 2Pd·THF, and 2Pt·DCM<sup>a</sup></caption><thead><tr><th rowspan="2">Compound</th><th colspan="3">Distances (Å)</th><th colspan="2">Angles (deg)</th></tr><tr><th>V=O</th><th>V-M</th><th>M···M'/M···S'</th><th>V-M···M'</th><th>τ<sub>intra</sub><sup>b</sup></th></tr></thead><tbody><tr><td>1Pd</td><td>1.5802(19)</td><td>2.9112(5)</td><td>3.1868(4)/3.97 (avg)</td><td>176.048(13)</td><td>14.9</td></tr><tr><td>1Pt<sup>c</sup></td><td>1.592(2)</td><td>2.8635(6)</td><td>3.1747(4)/3.94 (avg)</td><td>177.168(11)</td><td>15.2</td></tr><tr><td>2Pd·DCM</td><td>1.5827(13)</td><td>2.8355(3)</td><td>3.9280(3)/3.0737(5)</td><td>139.972(9)</td><td>21.3</td></tr><tr><td>2Pd·THF</td><td>1.5777(15)</td><td>2.8377(4)</td><td>4.2113(4)/3.0936(6)</td><td>134.557(10)</td><td>21.0</td></tr><tr><td>2Pt·DCM<sup>c</sup></td><td>1.581(4)</td><td>2.7823(10)</td><td>3.8408(5)/3.1266(14)</td><td>143.62(2)</td><td>20.6</td></tr></tbody></table>

<sup>a</sup>Dots indicate intermolecular contacts; primed symbols are used for atoms of the neighboring PW in a dimer. <sup>b</sup>Intramolecular twisting angle, evaluated as the average O–V–Pt–S dihedral angle (with O and S belonging to the same thiocarboxylate ligand). <sup>c</sup>Data taken from ref 17 (structures collected at 100 K).

<figure><img src="image_12.png" /><figcaption>**Figure 1.** Molecular structures of dimers {1Pd}₂ (staggered) and {2Pd}₂ (square, in 2Pd·DCM). Metals and heteroatoms are displayed as 50% probability ellipsoids, while carbon atoms are drawn in ball and stick style; hydrogens are omitted for clarity. Primed atoms are related to unprimed ones by 2-fold rotation in (a) and inversion in (b).</figcaption></figure>

(0.64 and 0.60 Å, respectively).<sup>41</sup> By consequence, compounds 1Pd and 2Pd·DCM are isomorphous with their Pt analogues (Figure S1). 1Pd crystallizes in the monoclinic C2/c space group and the unit cell contains eight [PdVO(SOCMe)₄] molecules in general position. The selected geometrical parameters gathered in Table 1 show that the V–Pd distance (2.91 Å) is slightly longer than the V–Pt one in 1Pt (2.86 Å). However, the two derivatives have virtually identical intramolecular twisting angles ($τ_{intra}$ = 14.9 vs 15.2°).<sup>17</sup> Molecules form approximately coaxial dimers (Figure 1a) with short M···M' contacts and virtually straight V–M···M' angles. Differences in intermolecular M···M' and (average) M···S' distances lie well below 0.05 Å in the two derivatives, which both have an intermolecular twisting angle of 32.6° (evaluated by averaging the smallest set of S–M···M'–S' dihedrals). This indicates that the conformation is virtually identical in the two compounds and closer to staggered (45°) than eclipsed (0°).<sup>6</sup>

Both 2Pd·DCM and 2Pd·THF crystallize in the triclinic P\bar{1} space group and are isomorphous with each other and with 2Pt·DCM.<sup>17</sup> The unit cell contains two [PdVO(SOCPh)₄] molecules and two solvent molecules (DCM or THF) in general positions. The V–Pd distance (2.84 Å) is again longer than the V–Pt distance in 2Pt·DCM (2.78 Å), although the intramolecular twisting angle is almost identical in 2Pd·DCM/THF and 2Pt·DCM (21.0–21.3 vs 20.6°). As in 2Pt·DCM, the PWs assemble into noncoaxial (square) dimers (Figure 1b),

with a displacement of one molecule relative to the other and formation of two reciprocating intermolecular M···S' contacts (~3.1 Å). By consequence, the M···M' distances in these square dimers (3.84–4.21 Å) are considerably larger than in the staggered dimers (3.17–3.19 Å), and the displacement angles (V–M···M') are correspondingly much more acute (134.6–143.6°). A significant difference between the three thiobenzoate derivatives is the displacement angle, which increases in the order 2Pd·THF < 2Pd·DCM < 2Pt·DCM with concomitant shortening of the M···M' separation (Table 1).

Statistics collected on the structurally characterized [MTr-(SOCR)₄(L)ₙ] complexes (Figure S4) show that the Tr–M distance varies within narrow limits and does not depend on the nature of the transition metal or the type of complex. On the other hand, the contact distance M···M' is expectedly much shorter in staggered dimers than in square dimers. Predictably, a bulky axial ligand (such as quinuclidine<sup>42</sup>) or a large cocrystalized compound (acridine<sup>43,44</sup>) prevents the formation of dimers, but the key factors ruling the way PWs assemble into dimers have not been fully ascertained.<sup>45</sup> Thioacetate PWs show a greater propensity to form staggered dimers than thiobenzoate ones, which may be caused by a higher steric hindrance in the latter, but other factors must play a role. For instance, the three known staggered dimers with thiobenzoate ligands are found in crystals of [PtCo(SOCPh)₄(H₂O)], [PtNi(SOCPh)₄(H₂O)],<sup>8</sup> and [PtZn(SOCPh)₄(H₂O)]<sup>12</sup> obtained from THF. However,
the cobalt(II) complex gives square dimers when the compound is recrystallized from DCM. By contrast, 2Pd assembles into noncoaxial dimers upon recrystallization from both solvents.

**2.3. Magnetic Measurements.** Magnetic characterization was performed on crystalline samples of 1Pd and 2Pd·THF (Figure 2). For 1Pd, the χ<sub>M</sub>T value (χ<sub>M</sub> is the magnetic

<figure><img src="image_13.png" /><figcaption>**Figure 2.** Temperature (T) dependence of the molar magnetic susceptibility (χ<sub>M</sub>) per vanadyl unit for 1Pd (a) and 2Pd·THF (b) measured in applied fields (H) of 0.1 and 1 T for temperatures below and above 15 K, respectively. Solid red lines are the best fits of the experimental data to a Heisenberg S = 1/2 spin dimer model. Insets: Field dependence of the magnetization (M) per vanadyl unit plotted as M vs HT<sup>-1</sup> for 1Pd (a) and 2Pd·THF (b) at the indicated temperatures. Solid red lines are the best fits of the experimental data to a Heisenberg S = 1/2 spin dimer model for 1Pd (a) and a simple S = 1/2 Brillouin function for 2Pd·THF (b).</figcaption></figure>

susceptibility per mole of vanadyl units and T is temperature) is 0.37 cm³ K mol⁻¹ at 300 K and remains approximately constant down to about 50 K. Upon further cooling, it progressively drops, reaching 0.03 cm³ K mol⁻¹ at 1.9 K (Figure 2a). The χ<sub>M</sub>T vs T data can be perfectly fitted to an S = 1/2 spin dimer model in the low-field approximation with g = 1.99(5) and J = 5.23(4) cm⁻¹ (Table 2), which indicates antiferromagnetic interactions between vanadium(IV) ions in neighboring complexes (J values are based on the JS₁·S₂ convention). Isothermal magnetization (M) vs applied field (H) data from 1.9 to 10 K exhibit a clear inflection point. They are consistent with a nonmagnetic S = 0 state and field-induced population of the excited S = 1 state, and can be fitted to give g = 1.97(5) and J = 5.16(4) cm⁻¹, in perfect agreement with the analysis of χ<sub>M</sub>T vs T data (Table 2). Considering that the closest V···V contact in the unit cell of 1Pd is shorter than the intradimer V···V' separation (Figure S2),

<table><caption>**Table 2.** Superexchange Coupling Constants *J* within Dimers of 1M and 2M from DC Magnetic Measurements and DFT Calculations<sup>*a*</sup></caption><thead><tr><th>Compound</th><th>Expt.</th><th>DFT</th><th>ref.</th></tr></thead><tbody><tr><td>1Pd</td><td>5.23(4)<sup>b</sup>, 5.16(4)<sup>c</sup></td><td>7.44</td><td>This work</td></tr><tr><td>2Pd</td><td>0.38(1)<sup>b</sup></td><td>0.79</td><td>This work</td></tr><tr><td>1Pt</td><td>4.70<sup>b</sup></td><td>6.55</td><td>17, 18</td></tr><tr><td>2Pt</td><td>0.644(2)<sup>b,c</sup></td><td>0.47</td><td>18</td></tr></tbody></table>

<sup>*a*</sup>*J* uses JS₁·S₂ convention and is expressed in cm⁻¹. <sup>*b*</sup>From the analysis of χ<sub>M</sub>T vs T data. <sup>*c*</sup>From the analysis of isothermal M vs H data.

interdimer couplings may, in principle, contribute to the observed magnetic behavior. The validity of our model is, however, supported by the quality of the fit and the intradimer *J* values estimated by DFT calculations (see below and ref 18).

The χ<sub>M</sub>T value for 2Pd·THF is 0.37 cm³ K mol⁻¹ at 300 K and undergoes only a slight decrease upon cooling (Figure 2b). χ<sub>M</sub>T vs T data indeed provide g = 1.99(5) and J = 0.38(1) cm⁻¹, a value indicative of a very weak antiferromagnetic interaction (Table 2). Consistent with this, the isothermal M vs H plots from 1.9 to 8 K conform to simple Brillouin behavior and can be fitted with g = 2.01(5) (Figure 2b).

The described dc magnetic analysis indicates that the magnetic coupling is much larger in the staggered thioacetate dimers than in the square thiobenzoate ones, in agreement with previous findings on the Pt derivatives.<sup>17,18</sup> In the staggered dimers (1M), the magnetic coupling is mediated by the metallophilic M···M' interaction<sup>16</sup> and occurs between vanadium(IV) ions that are 8.9–9.0 Å apart from each other. The coupling is, however, weaker than found in staggered dimers containing Co²⁺ (22–25 cm⁻¹) and Ni²⁺ ions (25–120 cm⁻¹).<sup>8,10,16</sup> The difference is clearly rooted in the fact that the unpaired electron in vanadyl complexes is mainly localized in the 3d<sub>xy</sub> metal orbital.<sup>30</sup> The vanadium 3d<sub>xy</sub> and palladium 4d<sub>xy</sub> orbitals are in a cofacial configuration and are engaged in a direct δ-type interaction, which is less effective than σ- and π-type interactions. At the same time, the 3d<sub>xy</sub> orbital of vanadium has a σ-nonbonding character and can only delocalize on the thiocarboxylate ligands through weak π-type interactions with the O atoms, making the superexchange paths also inefficient, as previously reported for the Pt analogues.<sup>18</sup> Overall, in both Pt and Pd derivatives, there will be a limited spin delocalization on the Group 10 metal via both direct and superexchange pathways (see below). As a consequence, the intradimer magnetic interaction is minimal in magnitude with respect to other PW complexes where the unpaired electrons are in σ-antibonding 3d orbitals with components along the Tr–M direction (e.g., 3d<sub>z²</sub>).<sup>8,16</sup> In the square thiobenzoate dimers (2M), the intradimer interaction is very weak but can still be clearly detected as a deviation from the straightforward magnetic response of two independent S = 1/2 spins. Fe²⁺, Co²⁺, and Ni²⁺ derivatives with a similar dimeric structure<sup>8,9</sup> also exhibit a much-weakened coupling with respect to the corresponding staggered dimers. However, in that case, the magnetic analysis was complicated by single-ion zero-field splitting effects and orbital contributions to the susceptibility, and an intradimer *J* value could not be reliably determined.

Finally, while the small *J* values in 2M may be biased by interdimer coupling (Figure S3),<sup>18</sup> the results for 1M clearly indicate a significantly stronger coupling in the Pd derivative, which is fully consistent with DFT predictions (see below and Table 2). It is noted that in the series of four eclipsed dimers
<figure><img src="image_14.png" /><figcaption>**Figure 3.** CV curves of **1M** (a) and **2M** (b). Glassy carbon (GC) working electrode, Ag wire quasi-reference electrode, -15 °C; for **1Pd** and **2Pd**: complex concentration 0.5 mM, 100 mM TBACl in DCM, scan rate 1 V s⁻¹; for **1Pt** and **2Pt**: complex concentration 1.0 mM, 50 mM TBACl in DCM, scan rate 0.05 V s⁻¹.¹⁸ The experimental conditions were adjusted for each sample to reach the best signal resolution.</figcaption></figure>

<table><caption>**Table 3.** *E*°' Values and Corresponding Thermodynamic Contributions in **1M** and **2M**ª</caption><thead><tr><th>Compound</th><th><i>E</i><sup>O'</sup>(V)<sup>d</sup></th><th><i>E</i><sub>LUMO</sub> (eV)</th><th>Δ<i>S</i><sup>O'</sup><sub>rc</sub> (J K<sup>-1</sup> mol<sup>-1</sup>)</th><th>Δ<i>H</i><sup>O'</sup><sub>rc</sub> (kJ mol<sup>-1</sup>)</th><th><i>T</i>Δ<i>S</i><sup>O'</sup><sub>rc</sub>/F (V)<sup>d</sup></th><th>-Δ<i>H</i><sup>O'</sup><sub>rc</sub>/F (V)</th></tr></thead><tbody><tr><td>1Pd<sup>b</sup></td><td>-1.632</td><td>-1.964</td><td>-98 ± 4</td><td>132.2 ± 6.9</td><td>-0.262</td><td>-1.370</td></tr><tr><td>1Pt<sup>c</sup></td><td>-1.752</td><td>-1.813</td><td>-167 ± 5</td><td>125.9 ± 7.3</td><td>-0.446</td><td>-1.305</td></tr><tr><td>2Pd<sup>b</sup></td><td>-1.418</td><td>-2.123</td><td>-139 ± 4</td><td>101.0 ± 5.3</td><td>-0.372</td><td>-1.046</td></tr><tr><td>2Pt<sup>c</sup></td><td>-1.517</td><td>-2.116</td><td>-140 ± 3</td><td>110.3 ± 3.7</td><td>-0.374</td><td>-1.143</td></tr></tbody></table>

⁴Potential values and thermodynamic parameters are referred to the Fc⁺/Fc redox couple. ⁵0.5 mM solution in DCM with 100 mM TBACl as base electrolyte, scan rate 1 V s⁻¹. ⁶1 mM solution in DCM with 50 mM TBACl as base electrolyte, scan rate 0.05 V s⁻¹ (see ref 18). ⁷At -15 °C (T = 258 K); the error associated with the *E*°' values is ±0.002 V.

[MTr(SOCMe)₄L] (Tr = Co, Ni; M = Pd, Pt; L = 3-nitropyridine) reported by Larsen et al.,¹⁵ replacement of Pt with Pd yields a remarkable decrease of 3d–3d coupling from 10–25 to 1–2 cm⁻¹, a feature that so far lacks a sound theoretical explanation.

## 2.4. UV–Vis and NMR Investigations in Solution.

Solutions of **1Pd** and **2Pd** were investigated by UV–vis spectroscopy and ¹H NMR. The UV–vis absorption spectra in DCM solution are very similar for Pt and Pd compounds, with somewhat larger extinction coefficients in the Pd derivatives and the same weak double-peaked band in the 550–900 nm interval assigned to the vanadyl group absorption (Figure S6).

Like **1Pt**, **1Pd** is ¹H NMR silent as a consequence of the strong paramagnetic effect of vanadium(IV) on the nearby methyl protons.¹⁸ The ¹H NMR spectra of **2Pd** in both CD₂Cl₂ (Figure S7) and THF-d₈ (Figure S8) show a broad but clearly resolved peak at 8.70 ppm (Δν₁/₂ = 0.17 kHz), which is accompanied by an exceedingly broad signal at ca. 4 ppm (Δν₁/₂ > 2.5 kHz). Peaks at similar chemical shift values (8.50 and 4.4 ppm) are present in the ¹H NMR spectrum of **2Pt** in CD₂Cl₂;¹⁸ their significantly smaller line width (89 Hz and 0.9 kHz) allowed integration and assignment to the *p*- and *m*-Ar protons, respectively. When THF-d₈ was not dried with molecular sieves prior to use, the solution of **2Pd** gradually turned brownish on standing, and new, much narrower peaks emerged in the aromatic region of the ¹H NMR spectrum. One set of signals can be tentatively assigned to the complex “[Pd(SOCPh)₂]”, while a second set is similar to the NMR pattern of thiobenzoic acid. At the same time, the water peak becomes broader and shifts from 1.52 to 2.65 ppm, indicating the solvation of vanadium(IV) by water molecules (Figure S8). We contend that the presence of

water triggers a structural rearrangement of the complex, which decomposes into separate V- and Pd-containing species. The latter then form insoluble brown precipitates of [Pd(SOCPh)₂]ₙ on the walls of the NMR tube.

The ¹H DOSY NMR spectrum of **2Pd** in CD₂Cl₂ (Figure S9) allowed a direct estimation of the molecular weight (MW) in solution based on external calibration curves (ECCs).⁴⁶,⁴⁷ Following the same approach employed for **2Pt**,¹⁸ we obtained MW = 763 ± 236 and 872 ± 294 g mol⁻¹, using the “dissipated spheres and ellipsoids” and merged calibration curves, respectively (details can be found in the *Experimental Section*).⁴⁷ A third ECC reported by Byers et al.,⁴⁸ which accounts for the presence of 3d transition elements, afforded MW = 812 ± 368 g mol⁻¹. These results strongly suggest that **2Pd**, like **2Pt**, is monomeric in solution (expected MW = 722.08 g mol⁻¹).

## 2.5. Electrochemistry.
Cyclic voltammetry (CV) curves recorded at -15 °C on **1Pd** and **2Pd** and on their Pt analogues are shown in Figure 3.¹⁸ Both Pd PWs display one efficient quasi-reversible one-electron reduction/oxidation process, leading to a well-defined cathodic peak (A and D) with the corresponding anodic counterpart (B and E). The latter is, however, not very intense at low scan rates (ν < 0.200 V s⁻¹) and becomes evident and well-formed only by increasing ν (ν > 0.250 V s⁻¹). This suggests that the reduced forms of **1Pd** and **2Pd** are somehow unstable and therefore only partially reoxidizable in the reverse scan. By contrast, the signals of the corresponding Pt complexes are stable and reversible at all scan rates. Table 3 presents the formal potentials (*E*°' vs Fc⁺/Fc) in the series of four compounds, as well as the calculated thermodynamic
<figure><img src="image_15.png" /></figure>

<figure><img src="image_16.png" /></figure>

<figure><img src="image_17.png" /></figure>

<figure><img src="image_18.png" /><figcaption>**Figure 4.** X-band CW- (a and b) and ESE-detected EPR spectra (c and d) of **1Pd** (a and c) and **2Pd** (b and d). CW spectra (a and b) were recorded at *T* = 77 K and *ν* = 9.4 GHz. ESE-detected spectra (c and d) were recorded at *T* = 10 K and *ν* = 9.7 GHz. The experimental/simulated spectra are drawn in yellow/black (a and c) or in blue/black (b and d) for **1Pd** and **2Pd**, respectively.</figcaption></figure>

<table><caption>**Table 4.** Spin-Hamiltonian Parameters for **1Pd** and **2Pd** from X-Band CW-EPR Spectra and DFT Calculations<sup>*a*</sup></caption><thead><tr><th rowspan="2">Parameters</th><th colspan="2">1Pd</th><th colspan="2">2Pd</th></tr><tr><th>EPR<sup>b</sup></th><th>DFT</th><th>EPR<sup>b</sup></th><th>DFT</th></tr></thead><tbody><tr><td>g<sub>x,y</sub></td><td>1.982(1)</td><td>1.9751</td><td>1.982(1)</td><td>1.9744</td></tr><tr><td>g<sub>z</sub></td><td>1.935(2)</td><td>1.9123</td><td>1.935(2)</td><td>1.9147</td></tr><tr><td>V<sub>A<sub>x,y</sub></sub></td><td>196(2)</td><td>-272.802, -272.740</td><td>195(2), 200(2)</td><td>-267.458, -267.485</td></tr><tr><td>V<sub>A<sub>z</sub></sub></td><td>517(3)</td><td>-609.372</td><td>517(3)</td><td>-602.460</td></tr><tr><td>Pd<sub>A<sub>x,y</sub></sub></td><td>-</td><td>2.601</td><td>-</td><td>2.593</td></tr><tr><td>Pd<sub>A<sub>z</sub></sub></td><td>-</td><td>5.338</td><td>-</td><td>5.084</td></tr></tbody></table>

<sup>*a*</sup>V<sub>A<sub>x,y,z</sub></sub> and Pd<sub>A<sub>x,y,z</sub></sub> describe the hyperfine and superhyperfine coupling with <sup>51</sup>V and <sup>105</sup>Pd nuclei, respectively, and are given in MHz; the experimental V<sub>A<sub>x,y,z</sub></sub> parameters are absolute values. <sup>*b*</sup>Frozen CD<sub>2</sub>Cl<sub>2</sub>/toluene-d<sub>8</sub> (1:1 v/v), 1 mM, 77 K.

contributions to E°, extracted from the temperature-dependent measurements (see Supplementary Note 1).

The E° values span the range from -1.75 to -1.42 V and increase in the order **1Pt** < **1Pd** < **2Pt** < **2Pd**. In particular, each Pd-based PW is easier to reduce than the corresponding Pt derivative by ca. 100-120 mV. This trend correlates with the higher second ionization energy of Pd vs Pt (19.43 vs 18.56 eV).<sup>49</sup> Furthermore, for the same Group 10 metal, the thiobenzoate derivative (**2M**) is easier to reduce than the thioacetate one (**1M**) by ca. 210-240 mV. This suggests that the R substituent on the thiocarboxylate ligand and the Group 10 metal act in a largely independent fashion in modulating the reduction potential of the PWs. According to the DFT calculations (see below), the LUMO is localized on the heavy metal in all the derivatives; therefore, the quasi-reversible peak

can be attributed to the reduction M<sup>2+</sup>/M<sup>+</sup>. At more negative potentials, the thioacetate derivative **1Pd** displays a further cathodic peak (-2.190 V; C) which is not observed for **1Pt**, while **2Pd** shows two one-electron reduction peaks at -1.969 (F) and -2.249 V (G). These peaks have no evident anodic counterpart and are therefore due to irreversible reduction processes.

**2.6. EPR Measurements.** Continuous-wave (CW) and electron spin echo (ESE) detected EPR spectra of **1Pd** and **2Pd** were recorded at X-band frequency on 1 mM frozen solutions in CD<sub>2</sub>Cl<sub>2</sub>/toluene-d<sub>8</sub> (1:1 v/v). The experimental EPR spectra (Figure 4) are typical of isolated vanadyl moieties and confirm that both complexes are present in solution as monomeric species. The spectra are dominated by the anisotropic hyperfine
<figure><img src="image_19.png" /><figcaption>Figure 5. (a) Temperature dependence of $T_1$ (squares) and $T_m$ (crosses) for 1Pd (yellow) and 2Pd (blue). Fit of the dependence of $T_1$ vs temperature for 1Pd (b) and 2Pd (c) using eq 2.</figcaption></figure>

coupling of the electron spin ($S = 1/2$) with the nuclear spin of $^{51}$V ($I = 7/2$) that yields a characteristic eight-line pattern.

The spectra of 1Pd and 2Pd can be satisfactorily simulated (Figure 4a,b, black lines) assuming collinear $\overline{\mathbf{g}}$ and $\overline{\mathbf{v}}\overline{\mathbf{A}}$ matrices in the spin Hamiltonian (SH) of eq 1:

$$ \hat{H} = \mu_B \mathbf{B} \cdot \overline{\mathbf{g}} \cdot \hat{\mathbf{S}} + \hat{\mathbf{S}} \cdot \overline{\mathbf{v}}\overline{\mathbf{A}} \cdot \hat{\mathbf{i}}_V \quad (1) $$

where $S = 1/2$ and $I_V = 7/2$. The SH parameters are in line with those reported for 1Pt and 2Pt$^{18}$ and are listed in Table 4. The same set of simulation parameters also successfully reproduces the ESE-detected spectra (Figure 4c,d, black lines). The SH parameters are consistent with the unpaired electron mainly localized in the 3d$_{xy}$ orbital of vanadium(IV). However, while the spectra of 1Pt and 2Pt show additional satellite lines arising from the hyperfine coupling with the naturally occurring fraction of $^{195}$Pt nuclei ($I = 1/2$, 34% abundance),$^{18}$ no such hyperfine splitting is detected for $^{105}$Pd ($I = 5/2$, 22% abundance) due to the low nuclear g factor of Pd compared to that of Pt ($g_{Pd} = -0.257$, $g_{Pt} = +1.219$). A simulation analysis (see Supplementary Note 2) indicates a maximum coupling on the order of 5 MHz and suggests a lower spin density transfer to the Pd atom compared to the Pt analogue,$^{18}$ consistent with DFT estimates (see below).

The spin relaxation properties of 1Pd and 2Pd in solution were studied using the same standard pulse EPR protocols employed for the Pt analogues$^{18}$ in order to ensure a meaningful comparison. All data were measured at a static magnetic field corresponding to the maximum intensity of the ESE-detected EPR spectrum (ca. 330 mT), which corresponds to the so-called powder position, where all possible molecular orientations are excited, providing the average contribution of all magnetic interactions (Zeeman and hyperfine).$^{50}$ The temperature dependence of $T_1$ and $T_m$ was measured in the temperature interval 10–100 K, the upper temperature limit being dictated by the softening of the matrix. At all temperatures, the echo decay displays Electron Spin Echo Envelope Modulation (ESEEM) oscillations due to the deuterated solvent (see Figure S11).

Figure 5a reports the $T_1$ and $T_m$ values obtained by fitting the experimental traces with a biexponential function (see Figures S11 and S12, and Table S2). The $T_m$ values were extracted by fitting the center of the echo decay traces. To test the robustness of the analysis, the time traces for both 1Pd and its Pt analogue$^{18}$ were also analyzed using mono- and stretched exponential functions (Table S3). At low temperatures ($T < 50$ K), the weighted average of the biexponential components agrees well with the monoexponential fit (e.g., 4.0 µs at 10 K vs a

biexponential average of 3.5 µs for 1Pd). Above 50 K, the decay is monoexponential. As a matter of fact, the biexponential fit converges to identical values for the slow and fast components.

Furthermore, given the pronounced ESEEM modulations in the low-temperature decay traces, the fits were also performed considering the maxima of the modulation (Figures S13 and S14). This does not alter the overall trend, but leads to a systematic increase in $T_m$ by ~10% when comparing the weighted average value obtained from the biexponential fit with the corresponding values obtained from mono- and stretched-exponential fits. Taking into account other sources of uncertainty,$^{50}$ we estimate the overall uncertainty in both $T_1$ and $T_m$ of the order of 20%. In the temperature interval 10–90 K, the slow component of a biexponential fit affords $T_m$ values of 4.5–0.7 µs in 1Pd and 12–1.3 µs in 2Pd. As previously observed for the 1Pt and 2Pt analogues,$^{18}$ 2Pd generally displays slightly longer relaxation times than 1Pd. It has been shown in multiple systems$^{50,51}$ that the rotation of the methyl groups has a detrimental effect on $T_m$. The bigger and sterically hindered phenyl substituents, on the other hand, have a smaller effect. The effect of the Group 10 metal is negligible as shown by the value of $T_m$ measured at 10 K, which goes from 6.1 to 4.5 µs and from 11 to 12 µs upon replacing Pt with Pd in 1M and 2M, respectively.

Figure 5a also shows a comparison of the dependence of $T_1$ on temperature for 1Pd and 2Pd. The experimental points correspond to the slow component of a biexponential fit, which is usually taken as a measure of the actual $T_1$.26 In line with $T_m$, 2Pd displays longer $T_1$ in the whole temperature range. For instance, at 10 K, the maximum values are of ca. 4.5 and 14 ms in 1Pd and 2Pd, respectively. The temperature dependence of the relaxation times proves that the values for Pd are very similar to those of Pt (4.0 and 11 ms for 1Pt and 2Pt, respectively, at 10 K), indicating that the replacement of Pt with Pd does not substantially affect the spin–lattice relaxation, even though the SOC constant of Pd is approximately three times smaller than that of Pt.$^{52,53}$ A plausible explanation is related to the minute spin density delocalization (ca. 1%) on the Group 10 metal.

To assess the effect of the lighter Group 10 metal on molecular vibrations and, in turn, on the spin–lattice relaxation, the temperature dependence of $1/T_1$ was analyzed using the same model employed for the Pt counterparts.$^{18}$ This is important to ensure a meaningful comparison between the Pd and Pt complexes. The model (eq 2) assumes that low-energy optical vibrations are responsible for Raman-like processes and also includes a direct mechanism contribution active at low temperature:$^{54,55}$
$$ \frac{1}{T_1} = A_{\text{dir}} \frac{e^{\hbar \omega_{\text{mw}} / k_B T}}{e^{\hbar \omega_{\text{mw}} / k_B T} - 1} + A_{\text{loc}} \frac{e^{\hbar \omega_{\text{loc}} / k_B T}}{e^{\hbar \omega_{\text{loc}} / k_B T} - 1}^2 \quad (2) $$

where $\omega_{mw}$ is the working microwave frequency ($\omega_{mw}/2\pi = 9.74$ GHz), and $\omega_{loc}$ is the frequency of an average effective vibration responsible for the so-called Raman process. The experimental data were fitted, leaving as free parameters only $\omega_{loc}$ and the coefficients $A_{\text{dir}}$ and $A_{\text{loc}}$ which quantify the effectiveness of each relaxation process. The fitted curves are shown in Figure Sb,c, and the parameters are reported in Table S4. For the two complexes, the model yields an average effective vibration with a frequency of 110–120 cm$^{-1}$. This result agrees with DFT simulations. The computed low-energy (<200 cm$^{-1}$) vibrational frequencies for 1Pd and 2Pd show a strong similarity with those calculated for the Pt analogues (see Tables S6, S7, and ref 18); even in this case, the thiobenzoate derivative presents a larger number of vibrations below 50 cm$^{-1}$, which are localized on the phenyl rings. Furthermore, in the examined frequency range, the vibrational modes of the thioacetate complex 1Pd have slightly lower energy than those of 1Pt, and the same trend is visible when comparing the vibrational frequencies above 80 cm$^{-1}$ for the two thiobenzoate derivatives.

The observation of an electron spin echo, per se, demonstrates the presence of electron spin coherence. However, in literature, this is further substantiated through Rabi nulation experiments, which were performed at 70 K, demonstrating the possibility of placing the spins in any arbitrary superposition of states (Figure S15).

## 2.7. DFT Calculations.
DFT calculations were first performed to evaluate the sign and magnitude of intradimer magnetic interactions in the solid state (Table 2). Within the Broken Symmetry (BS) approach, the simulations for 1Pd showed an antiferromagnetic intradimer coupling with $J = 7.44$ cm$^{-1}$, hence slightly larger than calculated for 1Pt at the same level of theory (6.55 cm$^{-1}$).¹⁸ A much smaller antiferromagnetic interaction is expected in 2Pd (0.79 cm$^{-1}$) and 2Pt (0.47 cm$^{-1}$),¹⁸ again with a small increase of antiferromagnetic coupling when replacing Pt with Pd. Although DFT calculations slightly overestimate most $J$ values, the experimental data are largely consistent with these theoretical predictions (Table 2). Since intradimer geometrical parameters ($M \cdots M'$ distances, twisting angles, etc.) are almost the same in the Pd and Pt derivatives (Table 1), the reason behind the observed trend lies in the electronic structure of Pd. The 4d orbitals are less diffuse than the 5d orbitals and lead to a larger overlap and a larger antiferromagnetic contribution. Worth noting is the opposite trend that would be suggested by the extent of spin delocalization on the Group 10 metal. According to Löwdin spin population analysis, the spin density on Pd is 0.012 e$^{-}$ in 1Pd and 0.010 e$^{-}$ in 2Pd. These values are lower than those computed for the platinum derivatives (0.014 and 0.011 e$^{-}$ on Pt in 1Pt and 2Pt, respectively¹⁸), in agreement with the EPR findings and the longer V–Pd distances, which lead to a less effective spin delocalization through V–Pd δ interaction. Even if full validation of such small variations would require a proper benchmark on a series of different compounds, the overall trend nicely correlates with the observed exchange constants.

DFT calculations were also used to relate the reduction potentials for the two compounds to the energies of their LUMO orbitals, which are –1.964 eV (1Pd) and –2.123 eV (2Pd). The reversed sign energy of the LUMO represents the electron affinity within Koopmans’ theorem. 2Pd should then be easier to reduce than 1Pd, as found for the two Pt analogues¹⁸

and confirmed experimentally. In addition, the greater propensity of 1Pd and 2Pd to undergo reduction with respect to their Pt counterparts agrees with their computed larger electron affinity (Table 3). It is also worth stressing that the LUMO of 1Pd and 2Pd, as that of their Pt counterparts,¹⁸ is localized on the Group 10 metal and shows a large $4d_{x^2-y^2}$ component (37.4 and 16.9% in 1Pd and 2Pd, respectively, according to Löwdin population analysis).

Finally, to account for the absence of a detectable hyperfine coupling with ¹⁰⁵Pd, we computed the ¹⁰⁵Pd hyperfine tensor ($Pd\AA$) in each Pd derivative (Table 4) and its composition in terms of Fermi contact, dipolar, and SOC contributions (Table S5). The overall hyperfine coupling amounts to ca. 5 MHz and is thus much smaller than in the Pt analogues.¹⁸ This estimate explains the absence of a detectable hyperfine coupling with ¹⁰⁵Pd and is consistent with the simulation analysis presented in the Supplementary Note 2. As a direct consequence of the lower SOC in Pd, the computed $Pd\AA$ tensor is also much less anisotropic than $Pd\AA$. For the same reason, the fraction of isotropic coupling arising from SOC is very different: 4.9–5.2% in 1Pd and 2Pd, and 15.3–16.7% in 1Pt and 2Pt (the remaining fraction originates from Fermi contact interaction).

## 3. CONCLUSIONS

Thiocarboxylate PW structures [PdVO(SOCR)₄] with R = Me (1Pd) and Ph (2Pd) contain an O-bonded vanadyl ion and an S-bonded Pd²⁺ ion, and have proven synthetically accessible despite the dual coordination nature of Pd²⁺, which can interact with both soft and hard donor atoms. Molecular and crystal structures are closely similar to those of the corresponding Pt derivatives 1Pt and 2Pt,¹⁷ with formation of staggered dimers supported by short Pd···Pd' contacts (1Pd) or square dimers featuring two reciprocating Pd···S' interactions (2Pd). The Group 10 metal has a small but significant effect on the strength of magnetic coupling within the staggered dimers via the metallophilic contact. The interaction is invariably antiferromagnetic, with $J = 4.7$ and $5.2$ cm$^{-1}$ in 1Pt and 1Pd, respectively (JS₁·S₂ convention). The slightly larger coupling promoted by the Pd···Pd' bridge is correctly reproduced by DFT calculations within the BS approach.

The dimers dissociate into monomers in organic solution, as evidenced by ¹H DOSY NMR and CW-EPR spectra. CV demonstrated that all four complexes undergo a quasi-reversible one-electron reduction, which is mostly centered on the Group 10 metal according to DFT studies. This reduction is experimentally found (and theoretically predicted to be) easier in Pd than in Pt derivatives, and in thiobenzoates than in thioacetates.

The X-band EPR spectra recorded in a frozen CD₂Cl₂/ toluene-d₈ matrix display a well-resolved hyperfine pattern due to the interaction of the $S = 1/2$ spin with the ⁵¹V nucleus ($I = 7/2$, ~100% abundance). The minority fraction of magnetically active ¹⁰⁵Pd isotope ($I = 5/2$, 22% abundance) gives no detectable hyperfine splitting of the X-band EPR lines, unlike the ¹⁹⁵Pt isotope in the Pt derivatives ($I = 1/2$, 34% abundance).¹⁸ Simulation analysis shows that this is due to the combined effect of a decreased spin delocalization on Pd (confirmed by DFT) and the approximately 5-fold smaller nuclear g-factor of ¹⁰⁵Pd vs ¹⁹⁵Pt. Coherence times measured by X-band pulsed EPR spectroscopy are comparable to those found in the Pt analogues in the same conditions, with thiobenzoate derivatives displaying systematically longer values (11–12 vs 4–6 µs at 10 K).¹⁸ Apparently, the largely decreased SOC constant
of the outermost d electrons of Pd has no effect on quantum coherence, presumably because of the minute (ca. 1%) spin density delocalized on the proximal Group 10 metal.

Our results demonstrate that thiocarboxylate-based PWs containing vanadyl ion and a diamagnetic Group 10 metal can be chemically tailored by modifying both the R substituent on the thiocarboxylate ligand and the ancillary heavier metal. This enables a fine-tuning of properties like molecular size, hydrophobicity, solubility, redox potentials, and intradimer magnetic coupling, while preserving quantum coherence. Further possibilities can be envisaged, like exploiting the different radial extension of 4d vs 5d orbitals to modulate surface-molecule interaction energy and electronic coupling.

## 4. EXPERIMENTAL SECTION

**4.1. General Methods.** Anhydrous DCM, n-hexane, acetone, D₂O (99.9 atom % D), CD₂Cl₂ (99.8 atom % D), THF-d₈ (99.5 atom % D), and toluene-d₈ (99.5 atom % D) were used as received. MeOH was predried with activated 3 Å molecular sieves (MS), then degassed by purging with dry nitrogen flux, distilled from Mg(OMe)₂, and stored over activated 3 Å MS in a nitrogen-filled Schlenk vessel.⁵⁶ THF was dried over Na, then distilled over Na with the addition of benzophenone under a dry nitrogen atmosphere, and stored over activated 4 Å MS prior to use. **1Pd** was synthesized under nitrogen using standard Schlenk techniques. Compound **2Pd** is air-sensitive in solution and was recrystallized and stored inside a nitrogen-filled MBraun glovebox. Both deuterated and regular solvents for use in the glovebox were degassed by three freeze-pump-thaw cycles and stored over activated 4 Å MS. The water content in commercially available VOSO₄·nH₂O was determined by elemental analysis (EA) as n = 3.6. KSOCMe (98%) and K₂PdCl₄ were used in commercially available form. NaSOCPh was prepared by reaction of HSOCPh (93%, contaminated with benzoyl disulfide according to ¹H NMR) with NaHCO₃ in water, and its purity was checked by ¹H NMR in D₂O, EA, and IR spectroscopy.

EA (CHN) was performed using a Thermo Fisher Scientific Flash 2000 analyzer. IR spectra were measured on a Jasco 4700 FT-IR spectrometer in ATR mode with a resolution of 1 cm⁻¹. ¹H NMR spectra of **1Pd** and **2Pd** were collected in CD₂Cl₂ or THF-d₈ (C = 0.01 M) at 298 K on an Avance 400 spectrometer from Bruker Biospin. Chemical shifts (δ, ppm) are expressed downfield from TMS and referenced to the residual proton resonances of the solvent (CD₂Cl₂: δ = 5.32; THF-d₈: δ = 1.72 [CHD(3,4)])⁵⁷ or to TMS as an internal standard. ¹H DOSY NMR measurements on **2Pd** were carried out with a ledbpgp2s sequence (Bruker library) using bipolar gradient pulses.⁵⁸ The diffusion time (or big delta) and the gradient length (or small delta) were tuned to 0.03 s and 800 µs, respectively, using the monodimensional sequence ledbpgp2s1d (Bruker library). The signal decay was fitted with a single exponential function using Bruker Dynamic Center software (version 2.8.3). Diffusion coefficients were normalized using TMS as an internal reference following the procedure described by Stalke et al.⁴⁷ Alternatively, the residual proton signal of CD₂Cl₂ was used as a reference for the ECC developed by Byers et al.⁴⁸ The absolute errors on estimated *MWs* were calculated as reported by Stalke et al., who showed that these errors are largely determined by the uncertainties on ECC parameters.⁵⁹ Electronic absorption spectra were collected on ca. 3 × 10⁻³ M solutions in DCM using a Jasco V-770 UV-vis-NIR spectrometer operating in double-beam mode (optical path length *l* = 0.1 cm).

**4.2. Synthesis of [PdVO(SOCMe)₄] (1Pd).** Solid VOSO₄·3.6H₂O (14.8 mg, 0.0649 mmol), K₂PdCl₄ (21.4 mg, 0.0656 mmol), and KSOCMe (29.8 mg, 0.261 mmol) were introduced into a Schlenk vessel, which was then evacuated and refilled with nitrogen. MeOH (5 mL) was added via a cannula under a nitrogen atmosphere, and the mixture was stirred for 18 h. Upon the course of the reaction, it turned into a brown suspension. The solvent was evaporated on a vacuum line, then the Schlenk vessel was opened to the air, and the emerald-green soluble product was extracted with 4 mL of DCM, while other products

were filtered off. Evaporation gave a light green powder (20.7 mg, 67% yield based on Pd). EA (C₈H₁₂O₅PdS₄V, 473.80) calcd. C 20.28, H 2.55%; found C 21.08, H 2.69%. To obtain X-ray quality crystals, the compound was redissolved in DCM (0.5 mL), and n-hexane (6 mL) was layered above the solution. Upon standing, large emerald-green crystals grew up. EA (crystals, C₈H₁₂O₅PdS₄V, 473.80) calcd. C 20.28, H 2.55%; found C 20.28, H 2.37%. IR (Figure S5, $\tilde{\nu}_{max}$ cm⁻¹): 1540 s, 1516 m sh, 1502 m, 1473 m, 1414 m br, 1346 m, 1238 w, 1144 s, 1100 w sh, 998 s, 930 w sh, 723 s, 672 w sh, 539 w, 509 s, 447 s. UV-vis (Figure S6, DCM, 3.5 × 10⁻³ M): $\lambda_{max}$ nm (ε, M⁻¹cm⁻¹) ~608 (26), 714 (26).

**4.3. Synthesis of [PdVO(SOCPh)₄] (2Pd).** NaSOCPh (102.2 mg, 0.6380 mmol) and K₂PdCl₄ (51.8 mg, 0.159 mmol) were separately dissolved in 3 and 5 mL of water, respectively. Upon vigorous stirring, the solution of K₂PdCl₄ was dropwise added to the solution of NaSOCPh over 10 min in such a manner that no precipitate formed. The solution turned orangish-red. VOSO₄·3.6H₂O (36.4 mg, 0.160 mmol) was dissolved in water (2 mL) and added immediately after the addition of K₂PdCl₄ was over, resulting in the formation of an other precipitate. The reaction was allowed to proceed for 1 h, and the precipitate was collected on a glass filter and extensively washed with water. The precipitate on the filter was then washed with small portions of DCM that caused the dissolution of a red component, leaving only a grassy-green solid on the filter (acetone can be used as an alternative to DCM). The filter was dried in a vacuum until a constant mass. Crude yield 63.0 mg (56.5% based on Pd). The green precipitate was dissolved in hot THF, the solution was passed through the same filter, and evaporated to dryness. The subsequent recrystallization steps were carried out inside the glovebox. THF (0.5 mL) was added to the precipitate, and the suspension was allowed to stand for a few days. The compound recrystallized in the form of grassy-green plates. A slightly brownish mother liquid was decanted, and a fresh portion of THF was added. The crystals were identified by X-ray diffraction as THF solvate [PdVO(SOCPh)₄]·THF (2Pd·THF). After a few days, a green mother liquid was decanted, the crystals were briefly washed with 1 mL of n-hexane and dried in a flow of nitrogen, giving 41.3 mg (36%) of the pure product. EA better agrees with formula [PdVO(SOCPh)₄]·0.2THF (C₂₈.₈H₂₁.₆O₅.₂PdS₄V, 736.50): calcd. C 46.97, H 2.95%; found C 46.76, H 3.23%. A portion of the product was recrystallized from DCM in the same fashion, by leaving the powder under a small amount of solvent. The crystals obtained from this batch were identified by X-ray diffraction as DCM solvate [PdVO(SOCPh)₄]·DCM (2Pd·DCM). The mother liquid was decanted, and the crystals were washed with n-hexane and dried in a flow of nitrogen. EA results agree with a solvent-free compound (C₂₈H₂₀O₃PdS₄V, 722.08): calcd. C 46.57, H 2.79%; found 46.68, H 2.82%. IR (Figure S5, $\tilde{\nu}_{max}$ cm⁻¹): 1591 w, 1504 m, 1462 m, 1435 m, 1338 w, 1307 w, 1214 s, 1172 s, 1096 w, 1072 w, 1006 s, 955 s, 928 m sh, 894 w sh, 841 w, 770 s, 724 s, 682 s, 642 s, 614 w, 578 s, 512 s, 503 m, 484 w, 438 s. UV-vis (Figure S6, DCM, 3.1·10⁻³ M): $\lambda_{max}$ nm (ε, M⁻¹cm⁻¹) ~340 (8.8 × 10⁴), 620 (51), 711 (50).

¹H NMR (Figure S7, 400.13 MHz, CD₂Cl₂) δ 8.70 (br s, 4H; *p*-C₆H₅), ~4 (vbr s, 8H; *m*-C₆H₅).

¹H NMR (Figure S8, 400.13 MHz, THF-d₈) δ 8.69 (br s, 4H; *p*-C₆H₅), ~4 (vbr s, 8H; *m*-C₆H₅).

**4.4. X-ray Crystallography.** Single-crystal X-ray diffraction data on **1Pd**, **2Pd**·THF, and **2Pd**·DCM were collected on a Bruker-Nonius XAPEX diffractometer, equipped with a Mo Kα generator, an area detector, and a Kryoflex cryostat for data collection at 200 K. One of the components of an epoxy glue was used for crystal handling, selection, and mounting without significant loss of any lattice solvent. Matrix frames and data collection were done with APEX2 v1.0-22 software, while data reduction was performed with SAINT v7.06A program and was followed by scaling and multiscan absorption correction using SADABS v2.10.⁶⁰ The structures were solved by direct methods (SIR92)⁶¹ and refined by full matrix least-squares methods on *F*₀² using SHELXL-2018/3⁶² program and WINGX⁶³ v2020.2 suite. Unless otherwise noted, all non-hydrogen atoms were assigned anisotropic displacement parameters while H atoms were treated isotropically, with *U* = 1.5*U*<sub>eq</sub> (C) for methyl groups and *U* = 1.2*U*<sub>eq</sub> (C) for the other H atoms. In **1Pd**, *CH*₃ groups were subject to rotating group refinement (AFIX 37). In **2Pd**·THF and **2Pd**·DCM, all aromatic H atoms were
fully refined, while methylene hydrogens were placed in geometrically idealized positions and allowed to ride on the parent C atom. Atom C31 of the lattice THF molecule in 2Pd·THF required split-atom refinement with a common isotropic displacement parameter for the two components (ca. 60:40).

Deposition Numbers 2443816–2443818 contain the supplementary crystallographic data for this paper. These data are provided free of charge by the joint Cambridge Crystallographic Data Centre and Fachinformationszentrum Karlsruhe Access Structures service.

**4.5. Magnetic Measurements.** The magnetic measurements were carried out with a Quantum Design MPMS-XL SQUID magnetometer and a PPMS-9 susceptometer housed at the Center de Recherche Paul Pascal (Pessac, France). The MPMS-XL instrument works between 1.85 and 400 K with applied static fields (H) ranging from -7 to 7 T. In order to avoid loss of lattice solvent, crystals of 2Pd·THF stored under the mother liquor were collected by filtration and sealed in a polypropylene bag (typical size ≈ 2 × 1 × 0.02 cm) under an argon atmosphere just before measurement. The magnetic susceptibility was evaluated as $χ = M/H$, where M is the magnetization. Data reduction for 1Pd and 2Pd·THF (12.76 and 25.40 mg, respectively) was carried out assuming $MW = 473.8$ and $794.2$ g mol⁻¹, corresponding to the formulas [PdVO(SOCMe)₄] and [PdVO(SOCPh)₄]·THF, respectively, and a diamagnetic correction estimated as $-0.5$ MW·10⁻⁶ cm³ mol⁻¹. Note that the use of an MW corresponding to the unsolvated 2Pd yields unacceptably low g values in the fit. The ac susceptibility measurements down to 1.85 K were performed with the MPMS-XL and PPMS-9 systems using an oscillating field from 1 to 6 Oe for frequencies between 1 Hz and 10 kHz in a zero-dc field and up to 1 T. In the available temperature and frequency ranges, the samples do not display any observable slow relaxation of the magnetization.

**4.6. Electrochemistry.** A potentiostat/galvanostat PARSTAT mod. 2273A (EG&G PAR, Oak Ridge, USA) was used for CV measurements set up in a nitrogen-filled glovebox equipped with a cooling bath of dry n-heptane and a coldfinger cryostat for low-temperature operations. Experiments were carried out in the temperature range from -30 to -5 °C in steps of 5 °C and at different scan rates ($v = 0.02–5$ V s⁻¹) using a cell for small-volume samples (2 mL). A 1 mm diameter GC disk (PAR), a Pt wire, and an Ag wire were used as working, counter, and quasi-reference electrodes, respectively. The GC electrode was cleaned following a previously reported procedure.⁶⁴,⁶⁵ The measurements were carried out in DCM on 0.5 mM (1 mM) solutions of 1Pd and 2Pd (1Pt and 2Pt) using 100 mM (50 mM) TBACl as a supporting electrolyte. A careful feedback correction minimized the ohmic drop between the working and the reference electrodes. The ferrocenium/ferrocene (Fc⁺/Fc) redox couple (in DCM, $E^{°}$ = 0.460 V vs SCE and 0.946 ± 0.005 V vs the Ag wire quasi-reference electrode used in our experiments) was used to calibrate the potential of the quasi-reference electrode.⁶⁶ For the reversible CV signals, the formal potential value ($E^{°}$) of the electron transfer process was calculated as the semisum of the cathodic and anodic peak potentials, $E^{°}$ = $E_{1/2}$ = ($E_{pc}$ + $E_{pa}$)/2. The peak currents of all the signals were found to be proportional to the square root of v (not shown), indicating diffusion-controlled electron transfer processes. All the reported $E^{°}$ values are referred to the Fc⁺/Fc redox couple. Experiments with the addition of ferrocene standard were repeated at each temperature. Variable temperature experiments were carried out in an isothermal cell configuration in which the temperature of the reference and working electrodes was always the same. For this experimental configuration, the reduction entropy referenced to the Fc⁺/Fc redox couple ($ΔS^{°}_{rc}$) is given by eq 3:

$$ \Delta S^{°}_{rc} = S^{°}_{red} - S^{°}_{ox} = nF(\partial E^{°}/\partial T)_{P,ni} \quad (3) $$

Thus, $ΔS^{°}_{rc}$ can be calculated from the slope of the $E^{°}$ vs T plot, when it is linear. Under the same conditions, the enthalpy change (also referenced to the Fc⁺/Fc redox couple, $ΔH^{°}_{rc}$) was obtained from the Gibbs–Helmholtz equation, namely from the $E^{°}$/T vs 1/T plot. Each experiment was repeated at least five times. The CV curves remained stable over time in the investigated temperature range, and $E^{°}$ values were found to be reproducible within ± 0.002 V.

**4.7. EPR Measurements.** All EPR experiments were performed on frozen solutions using a CD₂Cl₂/toluene-d₈ (1:1 v/v) mixture as a solvent. The samples were prepared in an argon-filled glovebox with oxygen and water levels below 0.5 ppm. The final concentration of the 1Pd and 2Pd samples was 1 mM. The EPR tubes were sealed with PTFE tape before removal from the glovebox.

CW-EPR spectra were recorded at X-band frequency ($ν$ ≈ 9.40 GHz) on a Bruker EMX spectrometer equipped with an SHQ cavity. Low-temperature measurements were obtained using a finger Dewar working at 77 K. The experimental parameters were: microwave power 1 mW, modulation amplitude 0.4 mT, conversion time 20 ms.

Pulsed EPR spectra were also recorded at X-band frequency ($ν$ ≈ 9.74 GHz) on a Bruker Eleksys E580 spectrometer equipped with a dielectric ring resonator (ER4118X-MD5) housed in a cryogen-free variable temperature cryostat (Cryogenic Ltd.). During the measurements, the resonator was overcoupled to minimize ringdown effects due to the application of the intense microwave pulses. The ESE-detected EPR spectra were measured at T = 10 K, integrating the whole spin echo generated by the Hahn echo sequence ($π/2–τ–π–τ–echo$, with $τ$ = 200 ns and $π/2$ = 16 ns) as a function of the applied magnetic field. Spin coherence time traces ($T_m$) were measured by integrating the top of the spin echo (40 ns) generated by the Hahn echo sequence as a function of $τ$. Attention was paid to set the shot repetition time (SRT) at least five times larger than $T_1$ in order to avoid artifacts related to partial saturation of the echo signal. Spin–lattice relaxation time ($T_1$) traces were recorded by integrating the whole spin echo generated by the inversion recovery sequence ($π–t_w–π/2–τ–π–τ–echo$, with $τ$ = 200 ns and $π/2$ = 16 ns) as a function of the waiting time $t_w$. Rabi oscillations were measured by integrating the whole spin echo generated by the sequence $t_{nut}–t_w–π/2–τ–π–τ–echo$ as a function of the length of the nutation pulse $t_{nut}$ (from 0 to 2046 ns in steps of 2 ns), while keeping $t_w$ = 4 µs and $τ$ = 200 ns fixed. The linear dependence with respect to the applied microwave field was demonstrated by recording Rabi traces as a function of the microwave power attenuation at 1, 7, and 14 dB for 1Pd and 0, 6, and 12 dB for 2Pd. The primary $T_1$ and $T_m$ time traces were fitted according to a biexponential model, described by eqs 4 and 5.

$$ I(t_w) = y_0 + A_{\text{fast}} e^{-\frac{t_w}{T_{1,\text{fast}}}} + A_{\text{slow}} e^{-\frac{t_w}{T_{1,\text{slow}}}} \quad (4) $$

$$ I(2\tau) = y_0 + A_{\text{fast}} e^{-\frac{2\tau}{T_{m,\text{fast}}}} + A_{\text{slow}} e^{-\frac{2\tau}{T_{m,\text{slow}}}} \quad (5) $$

The $T_m$ values were extracted by fitting the center of the echo decay traces. To test the robustness of the analysis, the time traces for both 1Pd and its Pt analogue were analyzed using also mono- (eq 6) and stretched exponential (eq 7) functions.

$$ I(2\tau) = y_0 + A e^{-\frac{2\tau}{T_m}} \quad (6) $$

$$ I(2\tau) = y_0 + A e^{-(\frac{2\tau}{T_m})^{\beta}} \quad (7) $$

All simulations were performed using the EasySpin 6.0.6 package working in Matlab.⁶⁷

**4.8. DFT Calculations.** All the calculations were performed with the ORCA 6.0.1 quantum chemistry package.⁶⁸ The X-ray structures of 1Pd and 2Pd were optimized in the gas phase by DFT. Solvent (toluene) effects were included within the CPCM model.⁶⁹ We used the PBE0 hybrid exchange-correlation functional⁷⁰ along with D3 atom pairwise dispersion corrections with BJ damping.⁷¹ Def2-TZVPP basis set⁷² was employed for all the atoms. For palladium, we used an effective core potential⁷³ to describe the inner core electrons. For the optimization runs, “TightOpt” convergence criteria were set. Afterward, we ran analytical frequency calculations on the final optimized structures to verify that the geometry is a true local minimum of the potential energy surface.

The calculation of single-ion SH parameters, i.e., g-tensor and hyperfine coupling, was performed at the same level of theory on the final optimized structure with Douglas-Kroll-Hess (DKH) relativistic corrections and basis sets (DKH-def2-TZVPP). SARC basis sets were employed for palladium.⁷⁴
The intradimer exchange interaction was computed in dimer models. The models were built using two molecular neighboring units extracted from the crystal structure, as in our previous work on the Pt derivatives.¹⁸ The B3LYP functional⁷⁵ was chosen with def2-TZVPP basis for all atoms.

The isotropic magnetic-coupling constants were extracted within the BS approach,⁷⁶ as developed by Noodleman. Such a method enables the calculation of magnetic interactions using DFT without the need to employ computationally expensive multiconfigurational approaches. Therefore, single-point calculations on both the triplet and the BS state were performed. Afterward, we mapped the spin states on the Heisenberg–Dirac–Van Vleck Hamiltonian $\hat{H}_{\text{HDVV}} = J\hat{S}_1\hat{S}_2$, where $\hat{S}_1 = \hat{S}_2 = 1/2$, and computed the value of the superexchange constant $J$ with the following formula:

$$J = (E_{\text{HS}} - E_{\text{BS}})/(2\hat{S}_1\hat{S}_2) \qquad (8)$$

In eq 8, $E_{\text{HS}}$ is the energy of the high-spin (HS) triplet state (total spin $S=1$) and $E_{\text{BS}}$ is the energy of the BS state (nontotal symmetric $S=0$ singlet).

## ASSOCIATED CONTENT

### Supporting Information

The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.inorgchem.5c04835.

Additional crystallographic information and structural diagrams, IR, UV-vis, and NMR spectra, additional electrochemical data, details on EPR spectroscopy and magnetic measurements (PDF)

Computational data including Cartesian coordinates of DFT optimized structures (XYZ)

### Accession Codes

Deposition Numbers 2443816–2443818 contain the supplementary crystallographic data for this paper. These data can be obtained free of charge via the joint Cambridge Crystallographic Data Centre (CCDC) and Fachinformationszentrum Karlsruhe Access Structures service.

## AUTHOR INFORMATION

### Corresponding Authors

Andrea Cornia – Dipartimento di Scienze Chimiche e Geologiche e UdR INSTM, Università degli Studi di Modena e Reggio Emilia, Modena 41125, Italy; ● orcid.org/0000-0001-9765-3128; Email: acornia@unimore.it

Olga Mironova – Dipartimento di Scienze Chimiche e Geologiche e UdR INSTM, Università degli Studi di Modena e Reggio Emilia, Modena 41125, Italy; ● orcid.org/0000-0002-8650-7771; Email: mironovaoa.nsk@gmail.com

### Authors

Giacomo Bellini – Dipartimento di Scienze Chimiche e Geologiche e UdR INSTM, Università degli Studi di Modena e Reggio Emilia, Modena 41125, Italy

Alessio Nicolini – Dipartimento di Scienze Chimiche e Geologiche e UdR INSTM, Università degli Studi di Modena e Reggio Emilia, Modena 41125, Italy; ● orcid.org/0000-0002-4742-5458

Manuel Imperato – Dipartimento di Scienze Chimiche e Geologiche e UdR INSTM, Università degli Studi di Modena e Reggio Emilia, Modena 41125, Italy; Dipartimento di Scienze Fisiche, Informatiche e Matematiche, Università degli Studi di Modena e Reggio Emilia, Modena 41125, Italy; ● orcid.org/0000-0002-4963-0025

Antonio Ranieri – Dipartimento di Scienze della Vita, Università degli Studi di Modena e Reggio Emilia, Modena 41125, Italy

Marco Borsari – Dipartimento di Scienze Chimiche e Geologiche e UdR INSTM, Università degli Studi di Modena e Reggio Emilia, Modena 41125, Italy; ● orcid.org/0000-0002-3612-4764

Matteo Briganti – Dipartimento di Chimica “Ugo Schiff” e UdR INSTM, Università degli Studi di Firenze, Sesto Fiorentino, FI 50019, Italy; ● orcid.org/0000-0001-8576-3792

Rodolphe Clérac – Univ. Bordeaux, CNRS, Centre de Recherche Paul Pascal, CRPP, UMR 5031, Pessac 33600, France; ● orcid.org/0000-0001-5429-7418

Mathieu Rouzières – Univ. Bordeaux, CNRS, Centre de Recherche Paul Pascal, CRPP, UMR 5031, Pessac 33600, France

Enrico Salvadori – Dipartimento di Chimica e NIS Centre, Università degli Studi di Torino, Torino 10125, Italy; ● orcid.org/0000-0003-4394-9438

Maria Chiara Pagliero – Dipartimento di Chimica e NIS Centre, Università degli Studi di Torino, Torino 10125, Italy

Mario Chiesa – Dipartimento di Chimica e NIS Centre, Università degli Studi di Torino, Torino 10125, Italy; ● orcid.org/0000-0001-8128-8031

Complete contact information is available at: https://pubs.acs.org/10.1021/acs.inorgchem.5c04835

### Author Contributions

O.M.: synthesis and investigation, methodology, visualization, writing—original draft; G.B.: synthesis and investigation; A.N.: DOSY NMR characterization, formal analysis, methodology; M.I.: methodology, validation; A.R.: electrochemical measurements, formal analysis; M.Bo.: electrochemical measurements, formal analysis, validation; M.Br.: DFT calculations; R.C.: magnetic measurements, supervision, validation; M.R.: magnetic measurements, formal analysis; E.S.: EPR measurements, supervision; M.C.P.: EPR measurements, formal analysis; M.C.: EPR measurements, supervision, validation; A.C.: conceptualization, funding acquisition, project administration, crystallography, methodology, writing—original draft. All authors contributed to writing—review and editing, and writing relevant sections of the article.

### Notes

The authors declare no competing financial interest.

## ACKNOWLEDGMENTS

This work was funded by the Italian Ministry of University and Research under the program Dipartimenti di Eccellenza 2023-2027 through project CH4.0 (Department of Chemistry, University of Turin), as well as by the European Research Council through the ERC-2022-SYG project CASTLe (no. 101071533). M.C. and E.S. acknowledge the NQSTI project PRECESS, D13C24001190001. R.C. and M.R. thank the University of Bordeaux, the Region Nouvelle Aquitaine, Quantum Matter Bordeaux, the Association Française de Magnétisme Moléculaire (https://asso-am2.fr), Idex Bordeaux (Research Program GPR Light), and the Centre National de la Recherche Scientifique (CNRS) for their continuous support.

## REFERENCES

(1) Köberl, M.; Cokoja, M.; Herrmann, W. A.; Kühn, F. E. From Molecules to Materials: Molecular Paddle-Wheel Synthons of Macro-
molecules, Cage Compounds and Metal-Organic Frameworks. *Dalton Trans.* **2011**, *40* (26), 6834–6859.

(2) Gosselin, A. J.; Rowland, C. A.; Bloch, E. D. Permanently Microporous Metal-Organic Polyhedra. *Chem. Rev.* **2020**, *120* (16), 8987–9014.

(3) Senkovska, I.; Bon, V.; Abylgazina, L.; Mendt, M.; Berger, J.; Kieslich, G.; Petkov, P.; Luiz Fiorio, J.; Joswig, J. O.; Heine, T.; Schaper, L.; Bachetzky, C.; Schmid, R.; Fischer, R. A.; Pöppl, A.; Brunner, E.; Kaskel, S. Understanding MOF Flexibility: An Analysis Focused on Pillared Layer MOFs as a Model System. *Angew. Chem., Int. Ed.* **2023**, *62* (33), No. e202218076.

(4) Kitano, K. N. I.; Tanaka, R.; Kimura, T.; Tsuda, T.; Shimizu, S.; Takagi, H.; Nishioka, T.; Shiomi, D.; Ichimura, A.; Kinoshita, I.; Isobe, K.; Ooi, S. Lantern-Type Dinuclear $Cr^{III}Pt^{II}$ and $V^{IV}Pt^{II}$ Complexes Bridged by Pyridine-2-Thiolate. Synthesis and Characterization. *J. Chem. Soc., Dalton Trans.* **2000**, No. 6, 995–1000.

(5) Funahashi, Y.; Kamimura, T.; Mizutani, M.; Wasada-Tsutsui, Y.; Sakai, M.; Yashiro, H.; Ozawa, T.; Jitsukawa, K.; Masuda, H. Micro-Columnar Crystals with $(-O=V-Pd-)_n$ Infinite Chain Structures Assembled by Hetero-Dinuclear Metal Complexes. *Adv. Mater.* **2006**, *11–12*, 273–276.

(6) Beach, S. A.; Doerrer, L. H. Heterobimetallic Lantern Complexes and Their Novel Structural and Magnetic Properties. *Acc. Chem. Res.* **2018**, *51* (5), 1063–1072.

(7) Skipper, H. E.; May, C. V.; Rheingold, A. L.; Doerrer, L. H.; Kamenetska, M. Hard-Soft Chemistry Design Principles for Predictive Assembly of Single Molecule-Metal Junctions. *J. Am. Chem. Soc.* **2021**, *143* (40), 16439–16447.

(8) Dahl, E. W. W.; Baddour, F. G. G.; Fiedler, S. R. R.; Hoffert, W. A. A.; Shores, M. P. P.; Yee, G. T. T.; Djukic, J. P. P.; Bacon, J. W. W.; Rheingold, A. L. L.; Doerrer, L. H. H. Antiferromagnetic Coupling across a Tetrametallic Unit through Noncovalent Interactions. *Chem. Sci.* **2012**, *3* (2), 602–609.

(9) Baddour, F. G.; Fiedler, S. R.; Shores, M. P.; Bacon, J. W.; Golen, J. A.; Rheingold, A. L.; Doerrer, L. H. Pt···Pt vs Pt···S Contacts between Pt-Containing Heterobimetallic Lantern Complexes. *Inorg. Chem.* **2013**, *52* (23), 13562–13575.

(10) Baddour, F. G.; Fiedler, S. R.; Shores, M. P.; Golen, J. A.; Rheingold, A. L.; Doerrer, L. H. Heterobimetallic Lantern Complexes That Couple Antiferromagnetically through Noncovalent Pt···Pt Interactions. *Inorg. Chem.* **2013**, *52* (9), 4926–4933.

(11) Guillet, J. L.; Bhowmick, I.; Shores, M. P.; Daley, C. J. A.; Gembicky, M.; Golen, J. A.; Rheingold, A. L.; Doerrer, L. H. Thiocyanate-Ligated Heterobimetallic {PtM} Lantern Complexes Including a Ferromagnetically Coupled 1D Coordination Polymer. *Inorg. Chem.* **2016**, *55* (16), 8099–8109.

(12) Baddour, F. G.; Hyre, A. S.; Guillet, J. L.; Pascual, D.; Lopez-De-Luzuriaga, J. M.; Alam, T. M.; Bacon, J. W.; Doerrer, L. H. Pt-Mg, Pt-Ca, and Pt-Zn Lantern Complexes and Metal-Only Donor-Acceptor Interactions. *Inorg. Chem.* **2017**, *56* (1), 452–469.

(13) Beach, S. A.; Zuckerman, L. A.; Portillo, R. I.; Shores, M. P.; Rheingold, A. L.; Doerrer, L. H. Heterobimetallic {PtMn} and {PtFe} Lantern Complexes with Exceptionally Long Metallophilic Contacts. *Inorg. Chim. Acta* **2019**, *493*, 81–90.

(14) Gennari, M.; Givaja, G.; Castillo, O.; Hermosilla, L.; Gómez-García, C. J.; Duboc, C.; Lledós, A.; Mas-Ballesté, R.; Zamora, F. On the Road to MM'X Polymers: Redox Properties of Heterometallic Ni···Pt Paddlewheel Complexes. *Inorg. Chem.* **2014**, *53* (19), 10553–10562.

(15) Larsen, E. M. H.; Bonde, N. A.; Weihe, H.; Ollivier, J.; Vosch, T.; Lohmiller, T.; Holldack, K.; Schnegg, A.; Perfetti, M.; Bendix, J. Experimental Assignment of Long-Range Magnetic Communication through Pd & Pt Metallophilic Contacts. *Chem. Sci.* **2023**, *14* (2), 266–276.

(16) Larsen, E. M. H.; Vosch, T.; Bendix, J. Ductile Metallophilic Interactions and Their Mediation of Magnetic Interactions in $[Ni(Pt(SAc)₄)(PyNO₂)]_2$ Complexes. *Inorg. Chem.* **2025**, *64* (11), 5366–5374.

(17) Beach, S. A.; Guillet, J. L.; Lagueux, S. P.; Perfetti, M.; Livesay, B. N.; Shores, M. P.; Bacon, J. W.; Rheingold, A. L.; Arnold, P. L.; Doerrer, L. H. Hetero-trimetallic {LnOVPt} Complexes with Antiferromagnetic Ln-V Coupling and Magnetic Memory. *Chem. Commun.* **2020**, *56* (75), 11062–11065.

(18) Imperato, M.; Nicolini, A.; Borsari, M.; Briganti, M.; Chiesa, M.; Liao, Y. K.; Ranieri, A.; Raza, A.; Salvadori, E.; Sorace, L.; Cornia, A. Quantum Spin Coherence and Electron Spin Distribution Channels in Vanadyl-Containing Lantern Complexes. *Inorg. Chem. Front.* **2024**, *11* (1), 186–195.

(19) Atzori, M.; Garlatti, E.; Allodi, G.; Chicco, S.; Chiesa, A.; Albino, A.; De Renzi, R.; Salvadori, E.; Chiesa, M.; Carretta, S.; Sorace, L. Radiofrequency to Microwave Coherent Manipulation of an Organometallic Electronic Spin Qubit Coupled to a Nuclear Qudit. *Inorg. Chem.* **2021**, *60* (15), 11273–11286.

(20) Chicco, S.; Chiesa, A.; Allodi, G.; Garlatti, E.; Atzori, M.; Sorace, L.; De Renzi, R.; Sessoli, R.; Carretta, S. Controlled Coherent Dynamics of [VO(TPP)], a Prototype Molecular Nuclear Qudit with an Electronic Ancilla. *Chem. Sci.* **2021**, *12* (36), 12046–12055.

(21) Gimeno, I.; Urtizberea, A.; Román-Roche, J.; Zueco, D.; Camón, A.; Alonso, P. J.; Roubeau, O.; Luis, F. Broad-Band Spectroscopy of a Vanadyl Porphyrin: A Model Electronuclear Spin Qudit. *Chem. Sci.* **2021**, *12* (15), 5621–5630.

(22) Bader, K.; Dengler, D.; Lenz, S.; Endeward, B.; Jiang, S. D. D.; Neugebauer, P.; Van Slageren, J. Room Temperature Quantum Coherence in a Potential Molecular Qubit. *Nat. Commun.* **2014**, *5* (1), 5304.

(23) Atzori, M.; Morra, E.; Tesi, L.; Albino, A.; Chiesa, M.; Sorace, L.; Sessoli, R. Quantum Coherence Times Enhancement in Vanadium(IV)-Based Potential Molecular Qubits: The Key Role of the Vanadyl Moiety. *J. Am. Chem. Soc.* **2016**, *138* (35), 11234–11244.

(24) Yu, C.-J.; Graham, M. J.; Zadrozny, J. M.; Niklas, J.; Krzyaniak, M. D.; Wasielewski, M. R.; Poluektov, O. G.; Freedman, D. E. Long Coherence Times in Nuclear Spin-Free Vanadyl Qubits. *J. Am. Chem. Soc.* **2016**, *138* (44), 14678–14685.

(25) Atzori, M.; Tesi, L.; Morra, E.; Chiesa, M.; Sorace, L.; Sessoli, R. Room-Temperature Quantum Coherence and Rabi Oscillations in Vanadyl Phthalocyanine: Toward Multifunctional Molecular Spin Qubits. *J. Am. Chem. Soc.* **2016**, *138* (7), 2154–2157.

(26) Bader, K.; Winkler, M.; Van Slageren, J. Tuning of Molecular Qubits: Very Long Coherence and Spin-Lattice Relaxation Times. *Chem. Commun.* **2016**, *52* (18), 3623–3626.

(27) Yamabayashi, T.; Atzori, M.; Tesi, L.; Cosquer, G.; Santanni, F.; Boulon, M. E. E.; Morra, E.; Benci, S.; Torre, R.; Chiesa, M.; Sorace, L.; Sessoli, R.; Yamashita, M. Scaling Up Electronic Spin Qubits into a Three-Dimensional Metal-Organic Framework. *J. Am. Chem. Soc.* **2018**, *140* (38), 12090–12101.

(28) Eaton, G. R.; Eaton, S. S. Solvent and Temperature Dependence of Spin Echo Dephasing for Chromium(V) and Vanadyl Complexes in Glassy Solution. *J. Magn. Reson.* **1999**, *136* (1), 63–68.

(29) Tesi, L.; Lunghi, A.; Atzori, M.; Lucaccini, E.; Sorace, L.; Totti, F.; Sessoli, R. Giant Spin-Phonon Bottleneck Effects in Evaporable Vanadyl-Based Molecules with Long Spin Coherence. *Dalton Trans.* **2016**, *45* (42), 16635–16643.

(30) Tesi, L.; Lucaccini, E.; Cimatti, I.; Perfetti, M.; Mannini, M.; Atzori, M.; Morra, E.; Chiesa, M.; Caneschi, A.; Sorace, L.; Sessoli, R. Quantum Coherence in a Processable Vanadyl Complex: New Tools for the Search of Molecular Spin Qubits. *Chem. Sci.* **2016**, *7* (3), 2074–2083.

(31) Winpenny, R. E. P. Quantum Information Processing Using Molecular Nanomagnets as Qubits. *Angew. Chem., Int. Ed.* **2008**, *47* (42), 7992–7994.

(32) Graham, M. J.; Yu, C. J.; Krzyaniak, M. D.; Wasielewski, M. R.; Freedman, D. E. Synthetic Approach To Determine the Effect of Nuclear Spin Distance on Electronic Spin Decoherence. *J. Am. Chem. Soc.* **2017**, *139* (8), 3196–3201.

(33) Atzori, M.; Chiesa, A.; Morra, E.; Chiesa, M.; Sorace, L.; Carretta, S.; Sessoli, R. A Two-Qubit Molecular Architecture for Electron-Mediated Nuclear Quantum Simulation. *Chem. Sci.* **2018**, *9* (29), 6183–6192.
(34) Zadrozny, J. M.; Niklas, J.; Poluektov, O. G.; Freedman, D. E. Millisecond Coherence Time in a Tunable Molecular Electronic Spin Qubit. ACS Cent. Sci. **2015**, *1* (9), 488–492.

(35) Imperato, M.; Nicolini, A.; Boniburini, M.; Sartini, D.; Benassi, E.; Chiesa, M.; Gigli, L.; Liao, Y. K.; Raza, A.; Salvadori, E.; Sorace, L.; Cornia, A. Dual Structure of a Vanadyl-Based Molecular Qubit Containing a Bis(β-Diketonato) Ligand. *Inorg. Chem.* **2024**, *63* (17), 7912–7925.

(36) Mironova, O.; Bellini, G.; Nicolini, A.; Imperato, M.; Ranieri, A.; Borsari, M.; Briganti, M.; Clérac, R.; Rouzières, M.; Salvadori, E.; Pagliero, M. C.; Chiesa, M.; Cornia, A. Oxovanadium(IV) thiocarboxylate paddlewheels containing ancillary Group 10 metals: a comparative study on Pd and Pt derivatives. *ChemRxiv*, 2025.

(37) Umakoshi, K.; Kinoshita, I.; Shun'ichiro, O. O. I. Dinuclear Palladium(II) Complex of Pyridine-2-Thiol. *Inorg. Chim. Acta* **1987**, *127* (2), L41–L42.

(38) Estevan, F.; Ibáñez, S.; Ofori, A.; Hirva, P.; Sanaú, M.; Úbeda, M. A. Benzoato and Thiobenzoato Ligands in the Synthesis of Dinuclear Palladium(III) and -(II) Compounds: Stability and Catalytic Applications. *Eur. J. Inorg. Chem.* **2015**, *2015* (17), 2822–2832.

(39) Halder, P.; Santalucia, D. J.; Park, S. V.; Berry, J. F. From Pincer to Paddlewheel: C-H and C-S Bond Activation at Bis(2-Pyridylthio)-Methane by Palladium(II). *Inorg. Chem.* **2019**, *58* (4), 2270–2274.

(40) Mironova, O. A.; Ryadun, A. A.; Sukhikh, T. S.; Pushkarevsky, N. A.; Konchenko, S. N. Synthesis and Photophysical Properties of Rare Earth (La Nd. Gd Ho) Complexes with Silanediamido Ligands Bearing a Chelating Phenylbenzothiazole Chromophore. *New J. Chem.* **2023**, *47* (7), 3406–3416.

(41) Shannon, R. D. Revised Effective Ionic Radii and Systematic Studies of Interatomic Distances in Halides and Chalcogenides. *Acta Crystallogr.* **1976**, A32 (5), 751–767.

(42) Guillet, J. L.; Ozumerzifon, T. J.; Shores, M. P.; Moore, C. E.; Rheingold, A. L.; Royappa, A. T.; Doerrer, L. H. Metalphilic Interactions and Magnetism in Heterobimetallic {Pt,M} Lantern Complexes. *Polyhedron* **2024**, *250*, 116788.

(43) Rheingold, A. L.; Skipper, H. E. CCDC 1977722: *Experimental Crystal Structure Determination*; 2020.

(44) Rheingold, A. L.; Doerrer, L. H.; Skipper, H. CCDC 2089337: *Experimental Crystal Structure Determination*; 2021.

(45) Mironova, O.; Bellini, G.; Nicolini, A.; Cornia, A. Multiple Dimerization Modes in Thiocarboxylate Paddlewheel Complexes: A Comprehensive View of Energy Landscapes from DFT Calculations and Statistics. *Chem. - Eur. J.* **2025**, No. e03004.

(46) Neufeld, R.; Stalke, D. Accurate Molecular Weight Determination of Small Molecules via DOSY-NMR by Using External Calibration Curves with Normalized Diffusion Coefficients. *Chem. Sci.* **2015**, *6* (6), 3354–3364.

(47) Bachmann, S.; Neufeld, R.; Dzemski, M.; Stalke, D. New External Calibration Curves (ECCs) for the Estimation of Molecular Weights in Various Common NMR Solvents. *Chem. - Eur. J.* **2016**, *22* (25), 8462–8465.

(48) Crockett, M. P.; Zhang, H.; Thomas, C. M.; Byers, J. A. Adding Diffusion Ordered NMR Spectroscopy (DOSY) to the Arsenal for Characterizing Paramagnetic Complexes. *Chem. Commun.* **2019**, *55* (96), 14426–14429.

(49) Kramida, A.; Ralchenko, Y.; Reader, J. and NIST ASD Team. *NIST Atomic Spectra Database* (version 5.12); National Institute of Standards and Technology: Gaithersburg, MD, 2024.

(50) Jackson, C. E.; Ngendahimana, T.; Lin, C. Y.; Eaton, G. R.; Eaton, S. S.; Zadrozny, J. M. Impact of Counter Ion Methyl Groups on Spin Relaxation in $[V(C_6H_4O_2)_3]^{2-}$. *J. Phys. Chem. C* **2022**, *126* (16), 7169–7176.

(51) Du, J. L.; Eaton, G. R.; Eaton, S. S. Effect of Molecular Motion on Electron Spin Phase Memory Times for Copper(II) Complexes in Doped Solids. *Appl. Magn. Reson.* **1994**, *6*, 373–378.

(52) Griffith, J. S. *The Theory of Transition-Metal Ions*; Cambridge University Press: London, 1961.

(53) Koseki, S.; Matsunaga, N.; Asada, T.; Schmidt, M. W.; Gordon, M. S. Spin-Orbit Coupling Constants in Atoms and Ions of Transition

Elements: Comparison of Effective Core Potentials, Model Core Potentials, and All-Electron Methods. *J. Phys. Chem. A* **2019**, *123* (12), 2325–2339.

(54) Eaton, S. S.; Harbridge, J.; Rinard, G. A.; Eaton, G. R.; Weber, R. T. Frequency Dependence of Electron Spin Relaxation for Three S = 1/2 Species Doped into Diamagnetic Solid Hosts. *Appl. Magn. Reson.* **2001**, *20* (1), 151–157.

(55) Lunghi, A.; Sanvito, S. The Limit of Spin Lifetime in Solid-State Electronic Spins. *J. Phys. Chem. Lett.* **2020**, *11* (15), 6273–6278.

(56) Armarego, W. L. Purification of Laboratory Chemicals, 8th ed.; Butterworth-Heinemann: Oxford, UK, 2017.

(57) Fulmer, G. R.; Miller, A. J. M.; Sherden, N. H.; Gottlieb, H. E.; Nudelman, A.; Stoltz, B. M.; Bercaw, J. E.; Goldberg, K. I. NMR Chemical Shifts of Trace Impurities: Common Laboratory Solvents, Organics, and Gases in Deuterated Solvents Relevant to the Organometallic Chemist. *Organometallics* **2010**, *29* (9), 2176–2179.

(58) Wu, D. H.; Chen, A.; Johnson, C. S. An Improved Diffusion-Ordered Spectroscopy Experiment Incorporating Bipolar-Gradient Pulses. *J. Magn. Reson., Ser. A* **1995**, *115* (2), 260–264.

(59) Bachmann, S.; Gernert, B.; Stalke, D. Solution Structures of Alkali Metal Cyclopentadienides in THF Estimated by ECC-DOSY NMR-Spectroscopy (Incl. Software). *Chem. Commun.* **2016**, *52* (87), 12861–12864.

(60) Bruker AXS Inc, APEX2, SAINT, And SADABS, Bruker Advanced X-Ray Solutions; Bruker AXS Inc: Madison, Wisconsin, USA, 2012.

(61) Altomare, A.; Cascarano, G.; Giacovazzo, C.; Guagliardi, A.; Burla, M. C.; Polidori, G.; Camalli, M. SIR92 – a program for automatic solution of crystal structures by direct methods. *J. Appl. Crystallogr.* **1994**, *27*, 435–435.

(62) Sheldrick, G. M. Crystal Structure Refinement with *SHELXL*. *Acta Crystallogr., Sect. C: Struct. Chem.* **2015**, C71 (1), 3–8.

(63) Farrugia, L. J. WinGX and ORTEP for Windows: An Update. *J. Appl. Crystallogr.* **2012**, *45* (4), 849–854.

(64) Meneses, A. B.; Antonello, S.; Arévalo, M. C.; Maran, F. Double-Layer Correction for Electron-Transfer Kinetics at Glassy Carbon and Mercury Electrodes in N,N-Dimethylformamide. *Electroanalysis* **2006**, *18* (4), 363–370.

(65) Nicolini, A.; Galavotti, R.; Barra, A. L.; Borsari, M.; Caleffi, M.; Luo, G.; Novitchi, G.; Park, K.; Ranieri, A.; Rigamonti, L.; Roncaglia, F.; Train, C.; Cornia, A. Filling the Gap in Extended Metal Atom Chains: Ferromagnetic Interactions in a Tetrairon(II) String Supported by Oligo-α-Pyridylamido Ligands. *Inorg. Chem.* **2018**, *57* (9), 5438–5448.

(66) Antonello, S.; Arrigoni, G.; Dainese, T.; De Nardi, M.; Parisio, G.; Perotti, L.; René, A.; Venzo, A.; Maran, F. Electron Transfer through 3D Monolayers on $Au_{25}$ Clusters. *ACS Nano* **2014**, *8* (3), 2788–2795.

(67) Stoll, S.; Schweiger, A. EasySpin, a Comprehensive Software Package for Spectral Simulation and Analysis in EPR. *J. Magn. Reson.* **2006**, *178* (1), 42–55.

(68) Neese, F. Software Update: The ORCA Program System—Version 5.0. *Wiley Interdiscip. Rev.: Comput. Mol. Sci.* **2022**, *12* (5), No. e1606.

(69) Garcia-Ratés, M.; Neese, F. Effect of the Solute Cavity on the Solvation Energy and Its Derivatives within the Framework of the Gaussian Charge Scheme. *J. Comput. Chem.* **2020**, *41* (9), 922–939.

(70) Adamo, C.; Barone, V. Toward Reliable Density Functional Methods without Adjustable Parameters: The PBE0 Model. *J. Chem. Phys.* **1999**, *110* (13), 6158–6170.

(71) Grimme, S.; Ehrlich, S.; Goerigk, L. Effect of the Damping Function in Dispersion Corrected Density Functional Theory. *J. Comput. Chem.* **2011**, *32* (7), 1456–1465.

(72) Weigend, F.; Ahlrichs, R. Balanced Basis Sets of Split Valence, Triple Zeta Valence and Quadruple Zeta Valence Quality for H to Rn: Design and Assessment of Accuracy. *Phys. Chem. Chem. Phys.* **2005**, *7* (18), 3297–3305.

(73) Andrae, D.; Häußermann, U.; Dolg, M.; Stoll, H.; Preuß, H. Energy-Adjusted Ab Initio Pseudopotentials for the Second and Third Row Transition Elements. *Theor. Chim. Acta* **1990**, *77* (2), 123–141.
(74) Rolfes, J. D.; Neese, F.; Pantazis, D. A. All-Electron Scalar Relativistic Basis Sets for the Elements Rb-Xe. *J. Comput. Chem.* **2020**, *41* (20), 1842–1849.

(75) Becke, A. D. Density-Functional Thermochemistry. III. The Role of Exact Exchange. *J. Chem. Phys.* **1993**, *98* (7), 5648–5652.

(76) Bencini, A.; Totti, F. A Few Comments on the Application of Density Functional Theory to the Calculation of the Magnetic Structure of Oligo-Nuclear Transition Metal Clusters. *J. Chem. Theory Comput.* **2009**, *5* (1), 144–154.

<figure><img src="image_20.png" /></figure>
Supporting Information

# Oxovanadium(IV) thiocarboxylate paddlewheels containing ancillary Group 10 metals: a comparative study on Pd and Pt derivatives

Olga Mironova,ª* Giacomo Bellini,ª Alessio Nicolini,ª Manuel Imperato,ª,ᵇ Antonio Ranieri,ᶜ

Marco Borsari,ª Matteo Briganti,ᵈ Rodolphe Clérac,ᵉ Mathieu Rouzières,ᵉ Enrico Salvadori,ᵉ

Maria Chiara Pagliero,ᵉ Mario Chiesa,ᵉ and Andrea Corniaª*

## Table of contents

<table><tr><td>1. X-Ray diffraction</td><td>S2</td></tr><tr><td>2. IR, UV-Vis, and NMR spectroscopy</td><td>S5</td></tr><tr><td>3. Cyclic voltammetry</td><td>S8</td></tr><tr><td>Supplementary Note 1</td><td>S8</td></tr><tr><td>4. EPR spectroscopy</td><td>S9</td></tr><tr><td>Supplementary Note 2</td><td>S12</td></tr><tr><td>5. DFT calculations</td><td>S14</td></tr><tr><td>6. References</td><td>S16</td></tr></table>

<footnote>ªDipartimento di Scienze Chimiche e Geologiche e UdR INSTM, Università degli Studi di Modena e Reggio Emilia, via G. Campi 103, 41125 Modena, Italy. *E-mail: acornia@unimore.it, mironovaoa.nsk@gmail.com</footnote>

<footnote>ᵇDipartimento di Scienze Fisiche, Informatiche e Matematiche, Università degli Studi di Modena e Reggio Emilia, via G. Campi 213/A, 41125 Modena, Italy</footnote>

<footnote>ᶜDipartimento di Scienze della Vita, Università degli Studi di Modena e Reggio Emilia, via G. Campi 103, 41125 Modena, Italy</footnote>

<footnote>ᵈDipartimento di Chimica “Ugo Schiff” e UdR INSTM, Università degli Studi di Firenze, via della Lastruccia 3, 50019 Sesto Fiorentino, FI, Italy</footnote>

<footnote>ᵉUniv. Bordeaux, CNRS, Centre de Recherche Paul Pascal, CRPP, UMR 5031; 33600 Pessac, France</footnote>

<footnote>ᵉDipartimento di Chimica e NIS Centre, Università degli Studi di Torino, via P. Giuria 7, 10125 Torino, Italy</footnote>
# 1. X-Ray diffraction

<table><caption>**Table S1.** Crystal data and refinement parameters for **1Pd**, **2Pd·DCM**, and **2Pd·THF**.</caption><thead><tr><th></th><th>1Pd</th><th>2Pd·DCM</th><th>2Pd·THF</th></tr></thead><tbody><tr><td>Radiation</td><td>MoKα (λ = 0.71073 Å)</td><td>MoKα (λ = 0.71073 Å)</td><td>MoKα (λ = 0.71073 Å)</td></tr><tr><td>Chemical formula</td><td>C<sub>8</sub>H<sub>12</sub>O<sub>5</sub>PdS<sub>4</sub>V</td><td>C<sub>29</sub>H<sub>22</sub>Cl<sub>2</sub>O<sub>5</sub>PdS<sub>4</sub>V</td><td>C<sub>32</sub>H<sub>28</sub>O<sub>6</sub>PdS<sub>4</sub>V</td></tr><tr><td>Formula weight</td><td>473.76</td><td>806.94</td><td>794.12</td></tr><tr><td>T (K)</td><td>200(2)</td><td>200(2)</td><td>200(2)</td></tr><tr><td>Crystal size (mm<sup>3</sup>)</td><td>0.35 × 0.07 × 0.02</td><td>0.41 × 0.23 × 0.18</td><td>0.20 × 0.14 × 0.10</td></tr><tr><td>Crystal system</td><td>Monoclinic</td><td>Triclinic</td><td>Triclinic</td></tr><tr><td>Space group</td><td>C2/c</td><td>P<span style="text-decoration: overline">1</span></td><td>P<span style="text-decoration: overline">1</span></td></tr><tr><td>α (Å)</td><td>27.3395(11)</td><td>11.7722(4)</td><td>11.8395(9)</td></tr><tr><td>b (Å)</td><td>8.4421(4)</td><td>11.7823(4)</td><td>11.8559(8)</td></tr><tr><td>c (Å)</td><td>15.7903(7)</td><td>12.3804(4)</td><td>12.7320(9)</td></tr><tr><td>α (deg)</td><td>90</td><td>70.6953(13)</td><td>74.128(3)</td></tr><tr><td>β (deg)</td><td>121.7797(14)</td><td>75.3529(14)</td><td>68.833(3)</td></tr><tr><td>γ (deg)</td><td>90</td><td>87.6764(14)</td><td>86.828(3)</td></tr><tr><td>V (Å<sup>3</sup>)</td><td>3098.1(2)</td><td>1566.24(9)</td><td>1601.2(2)</td></tr><tr><td>Z</td><td>8</td><td>2</td><td>2</td></tr><tr><td>ρ<sub>calcd</sub> (g cm<sup>-3</sup>)</td><td>2.031</td><td>1.711</td><td>1.647</td></tr><tr><td>2θ<sub>min</sub>/2θ<sub>max</sub> (deg)</td><td>6.418/56.094</td><td>6.922/58.038</td><td>7.072/57.324</td></tr><tr><td>Reflections collected/independent</td><td>16023/3729</td><td>34447/8213</td><td>30096/8018</td></tr><tr><td>No of parameters/ restraints</td><td>176/0</td><td>459/0</td><td>476/0</td></tr><tr><td>R1/wR2 (all data)</td><td>0.0339/0.0756</td><td>0.0318/0.0786</td><td>0.0447/0.0868</td></tr><tr><td>R1/wR2 (I > 2σ(I))</td><td>0.0268/0.0708</td><td>0.0276/0.0743</td><td>0.0308/0.0807</td></tr><tr><td>GOF</td><td>1.075</td><td>1.027</td><td>1.093</td></tr><tr><td>Largest diff. peak/hole (eÅ<sup>-3</sup>)</td><td>1.21/-0.83</td><td>0.93/-1.21</td><td>0.93/-0.82</td></tr></tbody></table>
<figure><img src="image_21.png" /><figcaption>(b)</figcaption></figure>

<figure><img src="image_22.png" /></figure>

**Figure S1.** Overlayed structures of staggered dimers in **1Pd** and **1Pt** (a), and of square dimers in **2Pd**·DCM and **2Pt**·DCM (b). Pd- and Pt-containing paddlewheels are drawn in orange and blue, respectively. Capped sticks representation, hydrogen atoms omitted.

<figure><img src="image_23.png" /></figure>

**Figure S2.** Shortest V···V contacts in the unit cell of **1Pd**.
<figure><img src="image_24.png" /><figcaption>**Figure S3.** Shortest V···V contacts in the expanded unit cell of 2Pd·THF.</figcaption></figure>

<figure><img src="image_25.png" /><figcaption>**Figure S4.** Statistics collected on the first-row transition metal (Tr) paddlewheels [MTr(SCOR)₄(L)ₙ] with R = Me (a) and Ph (b) using CCDC 2025.2.0.¹ M···M contacts displayed as circles, Tr···M as triangles. Restriction of M···M contact was set as 0–5 Å. The dimers are classified according to Doerrer's group classification.²</figcaption></figure>
## 2. IR, UV-Vis, and NMR spectroscopy

<figure><img src="image_26.png" /><figcaption>**Figure S5.** IR spectra of 1Pd (a), 2Pd·THF (b), and the corresponding alkali metals salts and Pt analogues.</figcaption></figure>
<figure><img src="image_27.png" /></figure>

<figure><img src="image_28.png" /><figcaption>**Figure S6.** Electronic absorption spectra of **1Pd** (a) and **2Pd** (b) compared with the spectra of the corresponding Pt analogues in the region of vanadyl group absorption. DCM, 0.1 cm optical path length.</figcaption></figure>

<figure><img src="image_29.png" /><figcaption>**Figure S7.** ¹H NMR spectrum of **2Pd·THF** in CD₂Cl₂. TMS was added as an internal standard for DOSY experiments. THF and *n*-hexane are residuals from crystallization and crystal washing, respectively.</figcaption></figure>
<figure><img src="image_30.png" /><figcaption>**Figure S8.** ¹H NMR spectrum of 2Pd in wet THF-d₈, showing degradation in the presence of water. The red lower spectrum was recorded on a freshly prepared solution, and the green upper spectrum was collected 3 days later.</figcaption></figure>

<figure><img src="image_31.png" /><figcaption>**Figure S9.** ¹H DOSY NMR spectrum of 2Pd·THF in CD₂Cl₂ (0.01 M, 298 K, 400.13 MHz).</figcaption></figure>
### 3. Cyclic voltammetry

**Supplementary Note 1: thermodynamic parameters.** Temperature-dependent electrochemical measurements showed that the $E^{\circ}$ values of both **1Pd** and **2Pd** undergo a monotonic linear decrease with increasing temperature from -30 to -5 °C (Figure S10). The calculated values of $\Delta S^{\circ'}_{rc}$ and $\Delta H^{\circ'}_{rc}$ are reported in Table 3 together with those of the two Pt derivatives. The $E^{\circ}$ values consist of two additive contributions: a predominant enthalpic term ($-\Delta H^{\circ'}_{rc}/F$) and a smaller entropic term ($T\Delta S^{\circ'}_{rc}/F$), both negative in sign. The enthalpic and entropic contributions increase in the order **1Pd** < **1Pt** < **2Pt** < **2Pd** and **1Pt** < **2Pt** ~ **2Pd** < **1Pd**, respectively. It is apparent that the ordering of $E^{\circ}$ values is the same as that of the enthalpic contribution, except for **1Pt**. This derivative has a moderately negative value of $-\Delta H^{\circ'}_{rc}/F$ but is the hardest to reduce because of an unusually large and unfavourable entropic contribution. Overall, the ligand dependence of $E^{\circ}$ primarily mirrors the different enthalpic contributions to the formal potential of thioacetate and thiobenzoate derivatives, while the M dependence of $E^{\circ}$ reflects a more complex interplay between enthalpic and entropic terms.

<figure><img src="image_32.png" /></figure>

<figure><img src="image_33.png" /><figcaption>**Figure S10.** Plot of $E^{\circ}$ vs. T for **1Pd** (a) and **2Pd** (b) (0.5 mM solutions in DCM with 100 mM TBACl as base electrolyte).</figcaption></figure>
4. EPR spectroscopy

<figure><img src="image_34.png" /><figcaption>**Figure S11.** Echo decay traces and corresponding fits obtained using Equation 5 for 1Pd (a) and 2Pd (b).</figcaption></figure>

<figure><img src="image_35.png" /><figcaption>**Figure S12.** Inversion recovery traces and corresponding fits obtained using Equation 4 for 1Pd (a) and 2Pd (b).</figcaption></figure>
<figure><img src="image_36.png" /><figcaption>**Figure S13.** Different fit methods for the **1Pd** sample. Top, fit of the center of the ESEEM modulation using from left to right **Equations 6** (mono-exponential decay), **5** (bi-exponential decay), and **7** (stretched-exponential decay); bottom fit of the maxima of the ESEEM modulations using the same models. The corresponding *T*<sub>m</sub> values extracted from the different methods are listed in **Table S3**.</figcaption></figure>

<figure><img src="image_37.png" /><figcaption>**Figure S14.** Different fit methods for the **1Pt** sample. The data are taken from Ref. ³. Top, fit of the center of the ESEEM modulation using from left to right **Equations 6** (mono-exponential decay), **5** (bi-exponential decay), and **7** (stretched-exponential decay); bottom fit of the maxima of the ESEEM modulations using the same models. The corresponding *T*<sub>m</sub> values extracted from the different methods are listed in **Table S3**.</figcaption></figure>
<table><caption>**Table S2.** *T*<sub>1</sub> and *T*<sub>m</sub> values extracted from the biexponential fit of the experimental inversion recovery and echo decay traces (Equations 4 and 5). The fits are presented in **Figures S11** and **S12**.</caption><thead><tr><th rowspan="2">Compound</th><th rowspan="2">T (K)</th><th colspan="4">T<sub>1</sub> (μs)</th><th colspan="4">T<sub>m</sub> (μs)</th></tr><tr><th>A<sub>slow</sub></th><th>T<sub>1</sub> slow</th><th>A<sub>fast</sub></th><th>T<sub>1</sub> fast</th><th>A<sub>slow</sub></th><th>T<sub>m</sub> slow</th><th>A<sub>fast</sub></th><th>T<sub>m</sub> fast</th></tr></thead><tbody><tr><td rowspan="6">1Pd</td><td>10</td><td>52</td><td>4.5×10<sup>3</sup></td><td>48</td><td>8.0×10<sup>2</sup></td><td>75</td><td>4.5</td><td>25</td><td>0.38</td></tr><tr><td>20</td><td>51</td><td>9.3×10<sup>2</sup></td><td>49</td><td>219</td><td>65</td><td>3.6</td><td>35</td><td>0.44</td></tr><tr><td>30</td><td>55</td><td>2.1×10<sup>2</sup></td><td>45</td><td>62</td><td>59</td><td>3.0</td><td>41</td><td>0.45</td></tr><tr><td>50</td><td>52</td><td>26</td><td>48</td><td>11</td><td>-</td><td>2.3</td><td>-</td><td>2.3</td></tr><tr><td>70</td><td>50</td><td>6.7</td><td>50</td><td>2.7</td><td>-</td><td>1.7</td><td>-</td><td>1.7</td></tr><tr><td>90</td><td>21</td><td>3.6</td><td>79</td><td>1.6</td><td>-</td><td>0.74</td><td>-</td><td>0.74</td></tr><tr><td>100</td><td>84</td><td>2.0</td><td>16</td><td>1.0</td><td>-</td><td>0.70</td><td>-</td><td>0.70</td></tr><tr><td rowspan="6">2Pd</td><td>10</td><td>51</td><td>1.4×10<sup>4</sup></td><td>49</td><td>2.5×10<sup>3</sup></td><td>61</td><td>12</td><td>39</td><td>3.0</td></tr><tr><td>20</td><td>65</td><td>1.6×10<sup>3</sup></td><td>35</td><td>3.7×10<sup>2</sup></td><td>85</td><td>9.0</td><td>15</td><td>1.5</td></tr><tr><td>30</td><td>68</td><td>3.0×10<sup>2</sup></td><td>32</td><td>89</td><td>-</td><td>5.2</td><td>-</td><td>5.2</td></tr><tr><td>50</td><td>36</td><td>60</td><td>64</td><td>25</td><td>-</td><td>2.5</td><td>-</td><td>2.5</td></tr><tr><td>70</td><td>47</td><td>14</td><td>53</td><td>5.9</td><td>-</td><td>2.3</td><td>-</td><td>2.3</td></tr><tr><td>90</td><td>65</td><td>4.8</td><td>35</td><td>1.8</td><td>-</td><td>1.3</td><td>-</td><td>1.3</td></tr><tr><td>100</td><td>-</td><td>2.5</td><td>-</td><td>2.5</td><td>-</td><td>0.565</td><td>-</td><td>0.565</td></tr></tbody></table>

<table><caption>**Table S3.** Best fit parameters for **1Pd** and **1Pt** using different fitting functions (Equations 5-7). The fitting curves are shown in **Figures S13** and **S14**.</caption><thead><tr><th rowspan="3">Compound</th><th rowspan="3">T (K)</th><th colspan="8">T<sub>m</sub> (μs)</th></tr><tr><th colspan="2">monoexponential</th><th colspan="4">biexponential</th><th colspan="2">stretched exponential</th></tr><tr><th>A<sub>slow</sub></th><th>T<sub>m,slow</sub></th><th>A<sub>fast</sub></th><th>T<sub>m,fast</sub></th><th>T<sub>m,av</sub></th><th>T<sub>m,str</sub></th><th>β</th></tr></thead><tbody><tr><td colspan="9">Fit considering the "center" of the ESEEM oscillations</td></tr><tr><td rowspan="2">1Pd</td><td>10</td><td>4.0</td><td>74</td><td>4.5</td><td>26</td><td>0.38</td><td>3.4</td><td>3.4</td><td>0.8</td></tr><tr><td>70</td><td>1.7</td><td>-</td><td>1.7</td><td>-</td><td>1.7</td><td>1.7</td><td>2.0</td><td>1.7</td></tr><tr><td rowspan="2">1Pt</td><td>10</td><td>5.7</td><td>72</td><td>6.1</td><td>28</td><td>0.169</td><td>4.4</td><td>3.1</td><td>0.7</td></tr><tr><td>70</td><td>2.0</td><td>62</td><td>2.3</td><td>38</td><td>0.182</td><td>1.5</td><td>1.4</td><td>0.7</td></tr><tr><td colspan="9">Fit considering the maxima of the ESEEM oscillations</td></tr><tr><td rowspan="2">1Pd</td><td>10</td><td>3.6</td><td>20</td><td>7.0</td><td>80</td><td>3.0</td><td>3.8</td><td>3.5</td><td>0.9</td></tr><tr><td>70</td><td>1.4</td><td>-</td><td>1.4</td><td>-</td><td>1.4</td><td>1.4</td><td>1.6</td><td>1.2</td></tr><tr><td rowspan="2">1Pt</td><td>10</td><td>5.0</td><td>11</td><td>12.6</td><td>89</td><td>4.3</td><td>5.2</td><td>4.8</td><td>1.0</td></tr><tr><td>70</td><td>2.0</td><td>-</td><td>2.0</td><td>-</td><td>2.0</td><td>2.0</td><td>2.0</td><td>1.0</td></tr></tbody></table>
<table><caption>**Table S4.** Best fit parameters in Equation 1 for 1Pd and 2Pd.</caption><thead><tr><th>Compound</th><th>A<sub>dir</sub> [Hz]</th><th>A<sub>loc</sub> [MHz]</th><th>ħω<sub>loc</sub> [cm<sup>-1</sup>]</th></tr></thead><tbody><tr><td>1Pd</td><td>12</td><td>1.5</td><td>122</td></tr><tr><td>2Pd</td><td>3.62</td><td>0.62</td><td>107</td></tr></tbody></table>

<figure><img src="image_38.png" /><figcaption>**Figure S15.** Rabi frequency as a function of the relative intensity of the oscillating field $B_1$ for 1Pd (yellow) and 2Pd (blue).</figcaption></figure>

**Supplementary Note 2: simulation analysis of ¹⁰⁵Pd hyperfine coupling.** In our previous study on 1Pt and 2Pt,³ clear hyperfine splittings of the ⁵¹V EPR signals were observed, arising from the interaction of the unpaired electron with the ¹⁹⁵Pt nucleus (I = 1/2, 33.83% natural abundance, $g_n^{Pt}$ = 1.2190). From the measured hyperfine tensors — [in 1Pt: $^{Pt}A_z$ = 80.0(5) MHz, $^{Pt}A_{x,y}$ = 47(2) MHz; in 2Pt: $^{Pt}A_z$ = 74.0(5) MHz, $^{Pt}A_{x,y}$ = 42(2) MHz] — a spin density transfer to the Pt atoms of approximately $\rho_{Pt}$ = 2.7% was determined for both complexes.

In contrast, the CW-EPR spectra of the Pd analogues 1Pd and 2Pd show no detectable hyperfine structure attributable to the magnetically active ¹⁰⁵Pd isotope (I = 5/2, 22.33% natural
abundance, $g_n^{\text{Pd}} = -0.2570$). Assuming the same spin density transfer as in the Pt complexes, the expected $^{105}\text{Pd}$ hyperfine coupling constants can be estimated by scaling the elements of the $^{\text{Pt}\bar{\text{A}}}$ tensor according to the ratio of the nuclear g-factors. Applying the described approach to the thioacetate derivatives **1Pt** and **1Pd**, one obtains:

$$ [\text{Pt}\bar{\text{A}}_x, \text{Pt}\bar{\text{A}}_y, \text{Pt}\bar{\text{A}}_z] = [\text{Pt}\bar{\text{A}}_x, \text{Pt}\bar{\text{A}}_y, \text{Pt}\bar{\text{A}}_z] \times \left|\frac{g_n^{\text{Pd}}}{g_n^{\text{Pt}}}\right| = [47, 47, 80] \times \left|\frac{-0.2570}{1.2190}\right| = [9.9, 9.9, 16.9] \text{ MHz} $$

The corresponding spectral simulation, performed under this assumption, is shown in **Figure S16**. The simulated spectrum (red trace) displays clear $^{105}\text{Pd}$ hyperfine satellites, which are not observed in the experimental data (Figure S16b).

<figure><img src="image_39.png" /></figure>

<figure><img src="image_40.png" /><figcaption>**Figure S16.** Computer simulation of the X-band CW-EPR spectrum of **1Pd** based on the assumption that the spin density transfer on Pd is the same as that observed in the Pt analogue (see text). The figure shows both the full spectrum (a) and a magnification (b) of the low field transitions ($^{51}\text{V}: m_I = -7/2$ and $m_I = -5/2$, red square in **Figure S16a**). The following parameters were used in the simulation: $v = 9.457776$ GHz, $\bar{\text{g}} = [1.982, 1.982, 1.935]$, $\bar{\text{A}} = [197, 197, 518]$ MHz, $^{\text{Pd}\bar{\text{A}}} = [9.9, 9.9, 16.9]$ MHz. A Gaussian line shape with full width at half maximum (FWHM) of 0.5 mT was employed. This ensures that broadening effects are not related to the low abundance $^{105}\text{Pd}$ nuclei. In **Figure S16b**, the stick diagram highlights the $^{105}\text{Pd}$ hyperfine transitions</figcaption></figure>

To further investigate the extent of $^{105}\text{Pd}$ hyperfine coupling, we performed a series of spectral simulations starting from the DFT-computed hyperfine tensor for **1Pd**, namely $^{\text{Pd}\bar{\text{A}}} = [2.6, 2.6, 5.3]$ MHz. The results are presented in **Figure S17**. When no hyperfine interaction with $^{105}\text{Pd}$ is included (Figure S17a), the simulated spectral linewidth does not match the experimental lineshape. The best agreement is achieved using the DFT-calculated tensor (Figure S17b), whereas larger coupling values lead to excessive broadening in the spectral wings. Based on this analysis, a maximum $^{105}\text{Pd}$ hyperfine coupling of approximately 5 MHz can be inferred, which implies a reduced spin density transfer with respect to the Pt complex.
<figure><img src="image_41.png" /></figure>

<figure><img src="image_42.png" /></figure>

<figure><img src="image_43.png" /></figure>

<figure><img src="image_44.png" /><figcaption>**Figure S17.** Computer simulation of the X-band CW-EPR spectrum of **1Pd** assuming different $^{105}$Pd hyperfine tensors. Only the low field transitions ($^{51}$V: $m_I = -7/2$ and $m_I = -5/2$) are shown.

(a) $^{Pd}\overline{\text{A}} = 0$ MHz, (b) $^{Pd}\overline{\text{A}} = [2.6, 2.6, 5.3]$ MHz, (c) $^{Pd}\overline{\text{A}} = [2.6, 2.6, 5.3] \times 1.5 = [3.9, 3.9, 7.95]$ MHz, (d) $^{Pd}\overline{\text{A}} = [2.6, 2.6, 5.3] \times 2 = [5.2, 5.2, 10.6]$ MHz. All other parameters are the same as those reported in the caption of **Figure S16**.</figcaption></figure>

## 5. DFT calculations

<table><caption>**Table S5.** Contributions to the $^{105}$Pd hyperfine coupling tensor for **1Pd** and **2Pd** as computed by DFT. All values are in MHz.</caption><thead><tr><th>1Pd</th><th>PdA<sub>x</sub></th><th>PdA<sub>y</sub></th><th>PdA<sub>z</sub></th><th>Iso</th></tr></thead><tbody><tr><td>Fermi contact</td><td>3.3302</td><td>3.3302</td><td>3.3302</td><td>3.3302</td></tr><tr><td>spin-dipole</td><td>-0.4582</td><td>-0.4595</td><td>0.9177</td><td>-</td></tr><tr><td>spin-orbit</td><td>-0.2709</td><td>-0.2701</td><td>1.0900</td><td>0.1830</td></tr><tr><td>Total</td><td>2.6011</td><td>2.6005</td><td>5.3378</td><td>3.5131</td></tr></tbody></table>
<table><caption>**Table S6.** Calculated vibrational frequencies below 200 cm⁻¹ for **1Pd** and **1Pt**.</caption><thead><tr><th>2Pd</th><th>PdA<sub>x</sub></th><th>PdA<sub>y</sub></th><th>PdA<sub>z</sub></th><th>Iso</th></tr></thead><tbody><tr><td>Fermi contact</td><td>3.2554</td><td>3.2554</td><td>3.2554</td><td>3.2554</td></tr><tr><td>spin-dipole</td><td>-0.4135</td><td>-0.4149</td><td>0.8284</td><td>-</td></tr><tr><td>spin-orbit</td><td>-0.2487</td><td>-0.2481</td><td>1.0002</td><td>0.1678</td></tr><tr><td>Total</td><td>2.5932</td><td>2.5924</td><td>5.0839</td><td>3.4231</td></tr></tbody></table>

<table><caption>**Table S7.** Calculated vibrational frequencies below 200 cm⁻¹ for **2Pd** and **2Pt**.</caption><thead><tr><th rowspan="2">Normal Mode</th><th colspan="2">Frequency (cm<sup>-1</sup>)</th></tr><tr><th>1Pd</th><th>1Pt</th></tr></thead><tbody><tr><td>6</td><td>53.88</td><td>60.82</td></tr><tr><td>7</td><td>58.17</td><td>61.28</td></tr><tr><td>8</td><td>58.90</td><td>62.04</td></tr><tr><td>9</td><td>59.27</td><td>69.08</td></tr><tr><td>10</td><td>93.48</td><td>103.65</td></tr><tr><td>11</td><td>94.71</td><td>107.87</td></tr><tr><td>12</td><td>95.51</td><td>108.15</td></tr><tr><td>13</td><td>104.07</td><td>112.45</td></tr><tr><td>14</td><td>106.17</td><td>116.42</td></tr><tr><td>15</td><td>107.29</td><td>119.92</td></tr><tr><td>16</td><td>116.98</td><td>124.84</td></tr><tr><td>17</td><td>125.25</td><td>128.02</td></tr><tr><td>18</td><td>146.30</td><td>140.60</td></tr><tr><td>19</td><td>162.39</td><td>164.95</td></tr><tr><td>20</td><td>164.07</td><td>168.43</td></tr><tr><td>21</td><td>164.67</td><td>170.69</td></tr><tr><td>22</td><td>187.44</td><td>184.71</td></tr><tr><td>23</td><td>188.34</td><td>188.33</td></tr><tr><td>24</td><td>191.58</td><td>198.18</td></tr><tr><td>25</td><td>199.24</td><td>-</td></tr></tbody></table>

<table><thead><tr><th rowspan="2">Normal Mode</th><th colspan="2">Frequency (cm<sup>-1</sup>)</th></tr><tr><th>2Pd</th><th>2Pt</th></tr></thead><tbody><tr><td>6</td><td>17.48</td><td>7.64</td></tr><tr><td>7</td><td>19.10</td><td>8.03</td></tr><tr><td>8</td><td>19.45</td><td>11.97</td></tr><tr><td>9</td><td>27.90</td><td>14.80</td></tr><tr><td>10</td><td>29.11</td><td>15.34</td></tr></tbody></table>
<table><tr><td>11</td><td>29.74</td><td>30.15</td></tr><tr><td>12</td><td>33.27</td><td>31.93</td></tr><tr><td>13</td><td>40.49</td><td>35.72</td></tr><tr><td>14</td><td>46.92</td><td>43.75</td></tr><tr><td>15</td><td>67.44</td><td>68.83</td></tr><tr><td>16</td><td>77.31</td><td>77.41</td></tr><tr><td>17</td><td>77.81</td><td>80.97</td></tr><tr><td>18</td><td>83.28</td><td>96.08</td></tr><tr><td>19</td><td>109.74</td><td>106.46</td></tr><tr><td>20</td><td>109.91</td><td>117.70</td></tr><tr><td>21</td><td>112.78</td><td>124.44</td></tr><tr><td>22</td><td>124.04</td><td>125.44</td></tr><tr><td>23</td><td>124.31</td><td>128.22</td></tr><tr><td>24</td><td>127.32</td><td>131.45</td></tr><tr><td>25</td><td>137.24</td><td>135.07</td></tr><tr><td>26</td><td>150.45</td><td>155.89</td></tr><tr><td>27</td><td>176.28</td><td>168.43</td></tr><tr><td>28</td><td>176.50</td><td>170.61</td></tr><tr><td>29</td><td>176.86</td><td>171.92</td></tr><tr><td>30</td><td>185.40</td><td>175.07</td></tr><tr><td>31</td><td>187.31</td><td>194.46</td></tr><tr><td>32</td><td>-</td><td>196.90</td></tr></table>

## 6. References

(1) Groom, C. R.; Bruno, I. J.; Lightfoot, M. P.; Ward, S. C. The Cambridge Structural Database. Acta Crystallogr. B. Struct. Sci. Cryst. Eng. Mater. **2016**, B72 (2), 171–179. https://doi.org/10.1107/S2052520616003954.

(2) Beach, S. A.; Doerrer, L. H. Heterobimetallic Lantern Complexes and Their Novel Structural and Magnetic Properties. Acc. Chem. Res. **2018**, 51 (5), 1063–1072. https://doi.org/10.1021/acs.accounts.7b00585.

(3) Imperato, M.; Nicolini, A.; Borsari, M.; Briganti, M.; Chiesa, M.; Liao, Y. K.; Ranieri, A.; Raza, A.; Salvadori, E.; Sorace, L.; Cornia, A. Quantum Spin Coherence and Electron Spin Distribution Channels in Vanadyl-Containing Lantern Complexes. Inorg. Chem. Front. **2024**, 11, 186–195. https://doi.org/10.1039/d3qi01806g.