Open Access Article
This Open Access Article is licensed under a
Creative Commons Attribution 3.0 Unported Licence

Modelling local structural and electronic consequences of proton and hydrogen-atom uptake in VO2 with polyoxovanadate clusters

Sourav Chakraborty , Eric Schreiber , Karla R. Sanchez-Lievanos , Mehrin Tariq , William W. Brennessel , Kathryn E. Knowles and Ellen M. Matson *
Department of Chemistry, University of Rochester, Rochester, NY 14627, USA. E-mail:

Received 23rd May 2021 , Accepted 24th August 2021

First published on 25th August 2021


We report the synthesis and characterisation of a series of siloxide-functionalised polyoxovanadate–alkoxide (POV–alkoxide) clusters, [V6O6(OSiMe3)(OMe)12]n (n = 1−, 2−), that serve as molecular models for proton and hydrogen-atom uptake in vanadium dioxide, respectively. Installation of a siloxide moiety on the surface of the Lindqvist core was accomplished via addition of trimethylsilyl trifluoromethylsulfonate to the fully-oxygenated cluster [V6O7(OMe)12]2−. Characterisation of [V6O6(OSiMe3)(OMe)12]1− by X-ray photoelectron spectroscopy reveals that the incorporation of the siloxide group does not result in charge separation within the hexavanadate assembly, an observation that contrasts directly with the behavior of clusters bearing substitutional dopants. The reduced assembly, [V6O6(OSiMe3)(OMe)12]2−, provides an isoelectronic model for H-doped VO2, with a vanadium(III) ion embedded within the cluster core. Notably, structural analysis of [V6O6(OSiMe3)(OMe)12]2− reveals bond perturbations at the siloxide-functionalised vanadium centre that resemble those invoked upon H-atom uptake in VO2 through ab initio calculations. Our results offer atomically precise insight into the local structural and electronic consequences of the installation of hydrogen-atom-like dopants in VO2, and challenge current perspectives of the operative mechanism of electron–proton co-doping in these materials.


Vanadium dioxide (VO2) is a versatile material that can access a variety of crystal phases, each with its own chemical, electronic, and physical properties.1 These properties can be further tuned through the incorporation of impurities/dopants, improving the performance of VO2 in energy-saving smart window technologies and post lithium-ion battery cathode materials.1–6 In particular, the incorporation of hydrogen atoms (H-atoms) as dopants at exterior/surface sites of VO2 has been shown to tune the material's optoelectronic properties by stabilizing the conductive rutile phase of the material (VO2(R); Fig. 1).7–11 H-atom equivalents can be introduced to monoclinic VO2 (VO2(M)) post-synthetically via hydrothermal,12 hydrogen spillover,13,14 and most recently, electron–proton co-doping techniques.
image file: d1sc02809j-f1.tif
Fig. 1 Modelling proton and hydrogen-atom uptake in VO2 (top) using POV–alkoxide clusters (bottom). Silylium cations are used as surrogates for protons to improve cluster stability and facilitate crystallographic resolution.

Electron–proton co-doping of VO2(M) takes advantage of acidic protons and an applied voltage as the source of H-atoms, and is proposed to occur through the initial injection of an electron, followed by protonation of the lattice.8,15 Mechanistic analysis is complicated by poor resolution of structural perturbations that occur upon hydrogen uptake;8,13,15 the small size of the hydrogen nucleus precludes definitive resolution of changes to the proximal lattice structure of the material. As such, insight into the mechanism of H-atom uptake is largely based on computational modelling, which does not explicitly take into account the role that the coupled uptake of protons and electrons might play in lowering the kinetic barrier of the reduction of this insulating material at room temperature. Indeed, recent work from Mayer and coworkers has demonstrated that electrochemical H-doping of metal oxides proceeds via proton-coupled electron transfer, whereby the presence of acidic protons in the electrolyte facilitate reduction.16 These results call the original proposal for the post synthetic formation of H-doped VO2 into question.

One popular approach that offers atomic-level insight into the mechanism and local electronic consequences of hydrogen uptake involves the use of molecular metal oxide cluster complexes as models for bulk solids.17–20 Unlike the majority of POVs reported that are isolated with high-valent vanadium(V) centres,21–23 Lindqvist-type polyoxovanadate–alkoxide (POV–alkoxide) clusters can be isolated in a reduced state, where all six vanadium ions adopt tetravalent oxidation states, allowing our group to study these compounds as models for VO2.24–26 Additionally, the vanadium centres within these clusters adopt edge-sharing pseudo-octahedral coordination geometries that are analogous to the immediate chemical environment of vanadium(IV) ions in all known phases of VO2.27–30 While these clusters possess delocalised electronic structures that can mimic local electronic communication between adjacent vanadium centres, they lack the extended band structure and Mott physics of the bulk material, limiting the effectiveness of our model to understanding the local consequences of surface reactivity and atomistic changes within VO2.

Previously, we reported the reactivity of a reduced POV–alkoxide cluster with protons in an attempt to access models of proton uptake in VO2.31 Upon addition of an organic acid to the reduced cluster, [nBu4N]2[VIV6O7(OEt)12] (Et = C2H5), formation of an oxygen-deficient assembly was observed. The proposed mechanism invokes the formation of a transient hydroxide-functionalised species, [V6O6(OH)(OEt)12]1−, which undergoes rapid intermolecular disproportionation. The hydroxide functionalised assembly is intriguing, as its isolation would provide a handle to interrogate the structural and electronic consequences, as well as the mechanism of H-doping in VO2 in acidic medium.

For decades, trialkylsilylium cations have been used as surrogates for protons in organic synthesis. Researchers have noted that silylium ions are capable of mediating analogous chemical transformations to protons, specifically in examples involving the reduction of carbon–oxygen multiple bonds.32–34 More recently, this analogy has been extended to the use of organosilanes as hydrogen-like reagents in inorganic syntheses.35,36 Drawing inspiration from these works, we hypothesised that the use of SiR3+ as a bulky surrogate for a proton would improve the stability of the siloxide-equivalent of the unobserved hydroxide intermediate, thus allowing for its isolation and analysis.

Here, we present the synthesis and characterisation of a siloxide-functionalised POV–alkoxide cluster, [VIV6O6(OSiMe3)(OMe)12]1− (Me = CH3), generated via addition of trimethylsilyl trifluoromethylsulfonate (TMSOTf) to [VIV6O7(OMe)12]2− (1-V6IVO72−; Fig. 1). Notably, binding of the silylium ion to a surface vanadyl site enables reduction of the cluster by a single electron. The added reducing equivalent introduces a site-differentiated VIII ion into the Lindqvist core adjacent to the silylium dopant. This result suggests that cation coordination to the surface of the vanadium(IV) oxide assembly is required for reduction, challenging the established sequential addition of electrons then protons believed to be necessary to access H-doped VO2. Collectively, this work provides atomically precise insights into the local structural and electronic impacts of H-atom dopants in VO2, as well as mechanistic inferences that suggest proton-coupled electron transfer may be a valid alternative pathway for H-atom uptake in VO2.

Results and discussion

We first investigated the reactivity of TMSOTf with 1-V6IVO72−. We opted to use the methoxide-bridged variant of the POV–alkoxide cluster to minimise steric clashes between surface ligands and the trimethylsilylium ion. Addition of one equivalent of TMSOTf to a bright blue solution of complex 1-V6IVO72− at low temperature resulted in an immediate colour change to green (Scheme 1). Following workup, analysis of the reaction mixture by 1H NMR spectroscopy revealed three new paramagnetically shifted and broadened resonances (δ = 27.50, 24.12 and −9.79 ppm, Fig. S1). This pattern is consistent with the expected reduction of symmetry of the parent cluster (pseudo-Oh → pseudo-C4v) that would accompany surface functionalisation.37–39 Additional characterisation of the product by electrospray ionization mass spectrometry (ESI-MS, (−) mode) revealed a peak located at m/z = 863, corresponding to a molecular weight consistent with the addition of a single trimethylsilyl moiety to complex 1-V6IVO72− (Fig. S2). Collectively, these data suggest successful formation of the siloxide-functionalised assembly, 2-V6IVO6(OSiMe3)1−.
image file: d1sc02809j-s1.tif
Scheme 1 Synthesis of complex 2-V6IVO6(OSiMe3)1−.

To unambiguously identify the molecular structure of 2-V6IVO6(OSiMe3)1−, crystals suitable for analysis via single crystal X-ray diffraction (SCXRD) were grown from slow diffusion of diethyl ether into a concentrated solution of the product in dichloroethane (Fig. 2 and Table 1, see Table S1 for complete crystallographic data and refinement parameters). Refinement of structural data confirms installation of a trimethylsilyl cation at a single vanadyl moiety. Broadly speaking, bond metrics of the Lindqvist ion remain constant following attachment of the silyl moiety to the surface of the assembly. The only exception to this observation is the V1–O bond lengths of the site differentiated vanadium centre. Indeed, activation of the vanadyl moiety via formation of the siloxide ligand manifests in a significant reduction of the V1–O1 bond order, as indicated by the elongation of the V1–O1 bond (V1–O1 (2-V6IVO6(OSiMe3)1−) = 1.768(3) Å vs. V[double bond, length as m-dash]Ot (avg, 1-V6IVO72−) = 1.606(1) Å; Ot = terminal oxido moiety). The V1–O1 bond length is reminiscent of previously reported mononuclear vanadium(IV) complexes bearing siloxide ligands (1.761–1.772 Å).40,41 Additionally, activation of the terminal oxido ligand following silylium coordination results in significant truncation of the V1–Oc bond (V1–Oc (2-V6IVO6(OSiMe3)1− = 2.019(3) Å vs. V–Oc (avg. 1-V6IVO72−) = 2.311(6) Å; Oc = central μ6-oxido moiety). The V1–Oc bond distance resembles values reported previously upon defect formation in POV–alkoxide clusters (2.0666(14)–2.120(5) Å).37–39

image file: d1sc02809j-f2.tif
Fig. 2 Crystal structures of 1-V6IVO72− and 2-V6IVO6(OSiMe3)1− shown with 30% probability ellipsoids. Hydrogen atoms, counter ions, and solvent molecules have been removed for clarity.
Table 1 Selected bond lengths and angles of complexes 1-V6IVO72−, 2-V6IVO6(OSiMe3)1−, and 3-VIIIV5IV(OSiMe3)2−
Bond/angle 1-V6IVO72− , e distrib. = VIV6 2-V6IVO6(OSiMe3)1−, e distrib. = VIV6 3-VIIIV5IVO6(OSiMe3)2−, e distrib. = VIIIVIV5
a Bond distances obtained from previously reported structural analysis of complex 1-V6IVO72−.51
Vn–Ot (avg.) 1.606(1) Å (n = 1–3) 1.595(6) Å (n = 2–6) 1.6010(4) Å (n = 2–6)
Vn–Oc (avg.) 2.311(6) Å (n = 1–3) 2.342(33) Å (n = 2–5) 2.323(30) Å (n = 2–5)
V1–O1 1.768(3) Å 1.929(4) Å
V1–Oc 2.019(3) Å 2.138(3) Å
V1–O1–Si1 174.8(2)° (major, 78%) 176.3(3)° (major 78%)
160.5(6)° (minor, 22%) 151.1(8)° (minor 22%)

Surface functionalisation of polyoxometalates (POMs) via the addition of silylium cations has been previously reported. However, most studies have been rooted in the development of methodologies for the organofunctionalisation of polyoxotungstates, focused on the formation of bridging siloxide ligands at the nucleophilic terminal oxide ligands of lacunary assemblies.42–50 By contrast, the formation of terminal siloxide moieties at the surface of plenary POMs, as observed following addition of TMSOTf to 1-V6IVO72−, is quite rare. The single example is reported in the case of a niobium-functionalised polyoxotungstate ion, [Nb2W4O19]4−.51 In this work, addition of tert-butyldimethylsilyl triflouromethylsulfonate results in the silylation of a lone terminal NbV[double bond, length as m-dash]O bond to generate a siloxide ligand, similar to that observed in the case of complex 2-V6IVO6(OSiMe3)1−. The authors credit the selectivity of this reaction for the terminal Nb[double bond, length as m-dash]O ligand to the steric bulk of the selected silylium cation, noting that smaller electrophiles (e.g. protons, methyl cations, trimethylsilylium ions), bind preferentially to the more nucleophilic, bridging oxide positions. Here, we have eliminated the possibility of the engagement of μ2-oxido ligands by fully saturating these bridging positions with alkoxide moieties.

Crystallographic characterisation of 2-V6IVO6(OSiMe3)1− provides atomic-level insight for local structural changes that occur upon proton uptake in VO2. Analysis of the site-differentiated vanadium ion within the POV–alkoxide cluster reveals a slight changes in V–O bond distances along the z-axis of V1, which, through ligand field analysis, would translate to the stabilization of the dz2 orbital of V1. This result is notable, as previous work has indicated that changes in the electronic properties of VO2 upon proton uptake occur as a result of the modification of the energetics of the dz2 orbital of the doped ions.8,10,13,52

Implications of silylium coordination on the electronic structure of 1-V6IVO72−

Toward the goal of using reduced POV–alkoxide clusters to model changes in the local electronic properties of VO2 upon proton uptake, we next analysed the electronic structure of complex 2-V6IVO6(OSiMe3)1−. The siloxide-functionalised POV–alkoxide cluster crystallises with each vanadium ion in a unique position within the unit cell, enabling the assignment of the oxidation states of individual vanadium centres using bond valence sum (BVS) calculations (Table S2).53,54 All vanadium centres within the hexavanadate assembly were identified as tetravalent ions, resulting in an overall oxidation state distribution of VIV6 for complex 2-V6IVO6(OSiMe3)1−. This assignment is identical to that of the parent cluster, 1-V6IVO72− (VIV6).

Additional evidence for the retention of the oxidation state distribution upon silylium coordination was obtained via X-ray photoelectron spectroscopy (XPS, Fig. 3). This analytical technique provides valuable information about the oxidation states and chemical environment of elements present in a solid sample by measuring the binding energies of photo-ejected core electrons. The XPS data of complex 1-V6IVO72− exhibits a single V 2p3/2 peak with a binding energy of 515.3 eV, consistent with literature values for VIV.55 Similarly, in the case of complex 2-V6IVO6(OSiMe3)1−, a single V 2p3/2 peak is observed with a binding energy of 515.7 eV, suggesting this Lindqvist ion also possesses a hexavanadate core comprised of solely tetravalent vanadium ions. In contrast, the mixed valent monoanionic parent cluster [V6O7(OMe)12]1− exhibits two peaks in the V 2p3/2 region with binding energies of 515.7 eV and 517.5 eV that correspond to VIV and VV, respectively (Fig. S3). The assignment of oxidation states in these clusters was corroborated by comparison to their solid-state analogues (i.e. VO2); in vanadium oxides, the oxidation state of vanadium is estimated by the difference in binding energy between the O 1s and V 2p3/2 orbital levels.55,56 These energy differences in 1-V6IVO72− and 2-V6IVO6(OSiMe3)1− are 14.7 and 14.9 eV, respectively, which matches well with previous reports on VO2.57 It is worth noting that a small shift in binding energy is observed when comparing the spectra of complexes 1-V6IVO72− and 2-V6IVO6(OSiMe3)1−, which is likely due to the sensitivity of core electrons to changes in their local bonding environment.58 The coordination of a silylium ion to the surface of the cluster induces a small shift toward higher binding energy in 2-V6IVO6(OSiMe3)1−, indicating that silylium coordination reduces the electron density around the vanadium metal centres across the entire assembly.

image file: d1sc02809j-f3.tif
Fig. 3 X-ray photoelectron spectra of the V 2p3/2 region of (a) 2-VIV6O6(OSiMe3)1− and (b) 1-VIV6O72−.

Comparison of the electronic structures of the siloxide functionalised POV–alkoxide cluster, 2-V6IVO6(OSiMe3)1−, and its chloride-doped congener, [VIIIVIV4VVO6Cl(OEt)12]−1, allows for comparison of the local electronic environments surrounding surface-bound cationic and substitutional anionic dopants in VO2, respectively.24 Previously, we reported that substitution of a surface oxido ligand for a chloride moiety results in the redistribution of electron density across the cluster core; two adjacent VIV ions disproportionate to form a high-valent vanadyl centre (VV[double bond, length as m-dash]O) and a site differentiated VIII–Cl moiety.24 The change in electronic structure upon formation of the chloride-doped assembly was confirmed by XPS; after fitting the data collected for [VIIIVIV4VVO6Cl(OEt)12]1−, an oxidation state distribution of VIIIVIV4VV was confirmed. In contrast, spectroscopic analysis of complex 2-V6IVO6(OSiMe3)1− reveals perturbations of the electronic structure of the fully-oxygenated parent cluster, 1-V6IVO72−, but does not support the presence of either VIII or VV centres in the functionalised assembly. Instead, both assemblies are found to possess an isovalent oxidation state distribution of vanadium ions within the Lindqvist assembly (e.g. VIV6).

With appropriate caution, consideration of the differences in the changes in the electronic structures of complexes 2-V6IVO6(OSiMe3)1− and [VIIIVIV4VVO6Cl(OEt)12]1−, compared to that of their fully-oxygenated precursors, [VIV6O7(OR)12]2−, can provide insight into the consequences of the formation of surface-bound and substitutional dopants, respectively, in VO2. While both categories of dopants have been credited with injecting carrier density into the bulk material, it is clear from the local changes in electronic structure upon dopant formation in POV–alkoxide clusters that the mechanisms by which these changes in electronic structure take place are distinct. In the case of the chloride-doped assembly, formal reduction of the doped vanadium ion is observed, requiring redistribution of electron density across the cluster core.24 This results in the formation of a charge-separated species, where the chloride-functionalised vanadium centre is electronically decoupled from the remaining vanadyl ions of the assembly. In contrast, silylium coordination to the surface of the POV–alkoxide cluster reveals only small changes in the electronic structure of the low-valent vanadium oxide assembly; the Lindqvist cluster can tolerate the addition of the surface-bound dopant without major reorganization of electron density within the three-dimensional framework. This result suggests that changes in carrier density upon cation uptake in these materials are likely due to the change in reducibility of the lattice upon dopant incorporation. These considerations will be discussed in depth in connection with electrochemical analysis of 2-V6IVO6(OSiMe3)1− (vide infra).

Electrochemical consequences of silylium binding to 1-V6IVO72−

To further interrogate changes in the electronic properties of the POV–alkoxide cluster upon silylium coordination, complex 2-V6IVO6(OSiMe3)1− was characterised via cyclic voltammetry (CV). Analysis of the electrochemical profile of 2-V6IVO6(OSiMe3)1− in dichloromethane reveals a voltammogram containing four quasi-reversible redox events (E1/2 = +0.84, +0.20, −0.47, −1.08 V vs. Fc+/0; Fig. 4, Table 2). Comparing this to the CV of 1-V6IVO72−, one notes striking similarities in the overall trace; the CV of the parent cluster likewise possesses four, one-electron redox events at potentials slightly shifted from that of 2-V6IVO6(OSiMe3)1− (E1/2 = +0.84, +0.20, −0.44, −0.93 V vs. Fc+/059,60,63 The open circuit potential (OCP) of complex 2-V6IVO6(OSiMe3)1− was measured to be −0.78 V vs. Fc+/0, suggesting the isovalent VIV6 charge state lies between the two most reducing events. This observation indicates that silylium coordination stabilises an additional reduction event for the cluster core, corresponding to a VIV/III process not observed in the parent assembly, 1-V6IVO72−.
image file: d1sc02809j-f4.tif
Fig. 4 Cyclic voltammograms of complexes 1-VIV6O72− (blue) and 2-VIV6O6(OSiMe3)1− (green) collected in dichloromethane with 0.1 M [nBu4N][PF6] as the supporting electrolyte (scan rate = 20 mV s−1). Note that impurities observed in the CV of 2-V6O6(OSiMe3)1− result from oxidative instability of the siloxide-functionalised cluster.
Table 2 Electrochemical parameters of complexes 1-V6IVO72− and 2-V6IVO6(OSiMe3)1− in dichloromethane
Redox couple 1-V6O72− 2-V6O6(OSiMe3)1−
E 1/2 E 1/2
a All E1/2 values are reported on the basis of data collected in dichloromethane and referenced vs. Fc+/0.
[VIV2VV4]/[VIV3VV3] +0.84
[VIV3VV3]/[VIV4VV2] +0.20 +0.84
[VIV4VV2]/[VIV5VV1] −0.44 +0.20
[VIV5V1]/[VIV6] −0.93 −0.47
[VIV6]/[VIII1VIV5] −1.08

Reduction to VIII within a POV–alkoxide is a rare observation. Previous reports reveal that the formation of trivalent vanadium centres in cluster assemblies is possible upon formation of an oxygen defect; however, these ions typically remain electronically decoupled from other vanadyl sites in the core.24,25,31,37–39 Evidence for this phenomenon is observed upon chemical oxidation of the cluster; in all reported examples, oxidised forms of oxygen-deficient POV–alkoxide clusters retain a site-differentiated VIII centre at the defect site (i.e. the oxidation event occurs across vanadyl ions composing the remainder of the cluster core).24,25,37,38

Inspection of the CV of 2-V6IVO6(OSiMe3)1− reveals its instability under electroanalytical conditions. Anodic scans of the cluster produce a mixture of species, as indicated by the appearance of new electrochemical features located at E1/2 = −0.704, −0.155, and +0.519 V vs. Fc+/0. The redox potentials of the “impurities” closely resemble values reported previously for the electrochemical profile of [VIIIVIV5O6(OMe)12]1−.38 The observed formation of the oxygen-deficient POV–alkoxide cluster is consistent with disproportionation of 2-V6IVO6(OSiMe3)1−, resulting in an equimolar mixture of [VIIIVIV5O6(OMe)12]1− and [VIV5VVO7(OMe)12]1− and half an equivalent of hexamethyldisiloxane. This type of reactivity resembles that previously invoked for the putative hydroxide-functionalised POV–alkoxide cluster.31 Monitoring a dichloromethane solution of 2-V6IVO6(OSiMe3)1−via1H NMR spectroscopy confirms the suspected instability of the assembly; after 24 hours, full conversion of 2-V6IVO6(OSiMe3)1− to a 1[thin space (1/6-em)]:[thin space (1/6-em)]1 ratio of [VIIIVIV5O6(OMe)12]1− and [VIV5VVO7(OMe)12]1− is observed (Fig. S4). Notably, as one extends the range of the oxidative sweep past subsequent oxidation events, the current responses of redox processes corresponding to [V6O6(OMe)12]n and [V6O7(OMe)12]n increase in magnitude, which suggests that the concentrations of the products of disproportionation increase as the stability of the siloxide-functionalised cluster decreases under more oxidizing conditions (Fig. S5).

Modelling hydrogen-atom uptake in VO2(R); structural, electronic, and mechanistic insights

Interested in isolating a molecular model for H-doped VO2, we next sought to explore the reduction chemistry of complex 2-V6IVO6(OSiMe3)1−. As described above, the electrochemical profile of 2-V6IVO6(OSiMe3)1− suggests that binding of the silylium cation to the cluster surface stabilises the reduction of the vanadium oxide assembly. This suggests the possibility of the formation of a reduced siloxide-functionalised POV–alkoxide with an oxidation state distribution of VIIIVIV5. Chemical reduction of complex 2-V6IVO6(OSiMe3)1− with CoCp2 (E1/2 = −1.33 V vs. Fc+/0 in dichloromethane65) results in isolation of the reduced assembly, [V6O6(OSiMe3)(OMe)12]2− (3-VIIIV5IVO6(OSiMe3)2−; Scheme 2), in moderate yield (23%) as confirmed by 1H NMR spectroscopy (Fig. S6) and elemental analysis.
image file: d1sc02809j-s2.tif
Scheme 2 Synthesis of 3-VIIIV5IVO6(OSiMe3)2−.

Crystals of 3-VIIIV5IVO6(OSiMe3)2− suitable for SCXRD were grown from diffusion of diethyl ether into a concentrated solution of the complex in acetonitrile. Broadly, the structure of compound 3-VIIIV5IVO6(OSiMe3)2− resembles that of 2-V6IVO6(OSiMe3)1− (Fig. S7,Table 1, see Table S1 for complete crystallographic data and refinement parameters); the reduced siloxide-functionalised assembly adopts a pseudo-C4v point group. Reduction of the cluster core manifests in slight perturbations of the bond metrics of the vanadyl ions of the assembly. In complex 3-VIIIV5IVO6(OSiMe3)2−, the average V[double bond, length as m-dash]Ot distance of 1.610(4) Å is slightly elongated from that of 2-V6IVO6(OSiMe3)1− (1.595(6) Å), while the V–Oc distances of 3-VIIIV5IVO6(OSiMe3)2− are shorter than that of 2-V6IVO6(OSiMe3)1−. These variations in V[double bond, length as m-dash]Ot and V–Oc bond distances are inconsistent with structural changes that typically occur across the Lindqvist core following reduction. For example, in the case of the reduction of [V6O7(OMe)12]−1 to 1-V6IVO72−, the V[double bond, length as m-dash]Ot bond distances remain approximately constant (1.60(8) Å vs. 1.606(1) Å, respectively), while significant elongation of the V–Oc distances are observed (2.297(38) Å vs. 2.311(6) Å).53

More substantial perturbations of V–O bond distances are observed upon inspection of the bond metrics of the siloxide-functionalised vanadium centre (V1). Significant elongation of the V1–O1 bond of complex 3-VIIIV5IVO6(OSiMe3)2− (1.929(4) Å) is noted in comparison to that of 2-V6IVO6(OSiMe3)1− (1.768(3) Å). The V1–O1 bond length is slightly longer than VIII–OSiR3 bond distances reported in the literature (1.852–1.883 Å).41,66 This point can be rationalised by discrepancies in coordination geometry of the vanadium centre (pseudo-octahedral VIII ion in the case of 3-VIIIV5IVO6(OSiMe3)2−vs. pseudo-tetrahedral VIII centres in previous work). Taken together, these observations suggest that reduction of the cluster core might be localised to the site-differentiated vanadium centre.

Like 2-V6IVO6(OSiMe3)1−, complex 3-VIIIV5IVO6(OSiMe3)2− crystallises with each vanadium centre possessing a unique position within the unit cell. As such, oxidation state distribution assignments of vanadium centres within the cluster core can be determined through BVS calculations (Table S3). Analysis of the bond metrics for each vanadium centre indicates an oxidation state distribution of VIIIVIV5; the site-differentiated vanadium centre (V1) possesses V–O bond distances consistent with reduction to VIII, while the vanadyl ions composing the remainder of the cluster core retain their tetravalent oxidation states. BVS calculations support the hypothesis that reduction of the cluster core is largely localised to the site-differentiated vanadium ion, resulting in the formation of a VIII–OSiMe3 moiety embedded within the Lindqvist assembly. Spectroscopic analysis of H-doped VO2via XPS suggests that vanadium atoms bound to hydrogenated oxygen atoms also undergo a localised reduction from VIV to VIII.8,10,15

Although comparison of 2-V6IVO6(OSiMe3)1− and 3-VIIIV5IVO6(OSiMe3)2− is useful for gaining insight into the perturbations of the assembly upon addition of electron density, a more useful analysis from the perspective of modelling H-atom doping in VO2 lies in the comparison of the structural data of 3-VIIIV5IVO6(OSiMe3)2− to 1-V6IVO72−. In analogy to spectroscopic data reported for H-doped VO2, the general positions of vanadium ions within the Lindqvist lattice do not change substantially upon formation of 3-VIIIV5IVO6(OSiMe3)2−; V[double bond, length as m-dash]Ot and V–Oc lengths of vanadyl ions composing 1-V6IVO72− and 3-VIIIV5IVO6(OSiMe3)2− are approximately identical. However, the addition of a hydrogen atom surrogate to a V[double bond, length as m-dash]Ot bond has significant structural ramifications on the site-differentiated vanadium centre. Elongation of the V1–O1 bond is consistent with reduction of the bond order as a result of silyl doping. Analogous atomistic perturbations occur in computational modelling of H-doped VO2(R) as a result of decreased charge of the oxido moiety upon hydrogenation.14 This distortion ultimately results in reduction of electron correlation due to increased occupancy of V d-levels, suppressing the metal-to-insulator phase transition and promoting metallicity of the bulk material. To date, these subtle structural changes have not been emphasised in discussions of empirical data collected for VO2; however, the lattice expansion observed following hydrogen-atom doping indicates that the structural changes supported by ab initio calculations, and experimentally modelled in our system, are indeed operative in bulk VO2.14

Finally, it is useful to discuss the implications of the formation of complex 3-VIIIV5IVO6(OSiMe3)2− in the context of published reports detailing the mechanism of hydrogen-atom uptake in VO2. As described above, the oxidation state distribution of vanadium centres in complex 3-VIIIV5IVO6(OSiMe3)2− (VIIIVIV5) reveals the net addition of a H-atom surrogate results in reduction of a site-differentiated metal centre within the vanadium(IV) oxide assembly. This result is analogous to the changes in electronic structure following hydrogen-atom uptake in VO2.8,10,15 However, the two systems are distinct from one another in their proposed mechanism of H-atom doping. In the case of the formation of the siloxide-functionalised POV–alkoxide cluster, silylium coordination precedes reduction of the assembly. Recall that similar reduction to a POV–alkoxide cluster bearing an oxidation state distribution of VIIIVIV5 is not possible in the fully oxygenated assembly (e.g.1-V6IVO72−), underscoring the importance of silylium binding to the surface of the assembly prior to injection of electron density.

The mechanism of hydrogen-atom addition to complex 1-V6IVO72− stands in contrast to the suggested mechanism for hydrogen-atom doping in VO2 when using acid as the proton source, where injection of electron density is considered the necessary first step to generate H atom dopants in the lattice.8,15 The disparity in these two mechanisms, which follow diverging pathways along the thermochemical square schemes (see Fig. 5), suggests the possibility of a third route that could apply to both systems: cation(proton)-coupled electron transfer. It is worth noting that metal oxide materials have been widely reported to perform proton-coupled electron transfer, yielding stable, hydrogenated nanocrystalline scaffolds.16,67–70 However, similar mechanisms have not yet been invoked for hydrogen-atom uptake in VO2.

image file: d1sc02809j-f5.tif
Fig. 5 Cation-coupled electron transfer square scheme for the uptake of hydrogen atoms/silyl radicals into vanadium oxide assemblies.


Here, we present the investigation of proton- and hydrogen atom-doping in POV–alkoxide clusters using silylium ions as bulky proton surrogates. This strategy affords isolation of a series of siloxide-functionalised POV–alkoxide clusters, complexes 2-V6IVO6(OSiMe3)1− and 3-VIIIV5IVO6(OSiMe3)2−, that serve as models of H-doped VO2. The molecular nature of the discrete Lindqvist ion allows for the investigation of atomistic perturbations of V–O bond lengths as a function of silylium coordination. This investigation reveals significant elongation of the V–O bond length following cation binding to the surface of the assembly. Our findings resemble ab initio calculations performed for H-doped VO2 and provide, for the first time, structural evidence for the formal reduction of the V–O bond order following interstitial doping of the lattice. Silylium coordination stabilises the injection of charge density into the lattice, facilitating the formation of a VIII ion at the site-differentiated metal centre. Altogether, these observations qualify the siloxide-functionalised POV–alkoxide cluster as an intriguing model for structural and electronic changes that occur upon hydrogen-uptake in VO2.

Additionally, the synthesis and characterisation of the siloxide-functionalised POV–alkoxide clusters present a distinct entry in the atomically precise modelling of dopants using these types of clusters. The site-differentiated vanadium centre that forms upon the surface coordination of a silylium ion appears to remain embedded within the electronic structure of the Lindqvist assembly. This behaviour is in contrast to our previous studies describing the incorporation of substitutional dopants (e.g. transition metals, halides) within the hexavanadate core.24,25,61,62,64,71–73 In these examples, the heteroatom-functionalised POV–alkoxide clusters adopt electronic structures where the site-differentiated ion is decoupled from the electronic structure of the remaining vanadyl ions of the Lindqvist core. These findings serve to deepen our understanding of the electronic architecture of molecular metal oxides, as well as indicate that the doping motif is relevant to the physicochemical effects observed for a doped material.

Our results suggest alternative mechanisms for electron–proton co-doping of VO2. Electrochemical characterisation of the siloxide-functionalised assemblies indicate that silylium coordination is required prior to the reduction of the vanadium oxide assembly. This finding is in contrast to previously reported mechanisms of H-doping in VO2(R) that suggest that the injection of electron density precedes the intercalation of protons.8,15 While it is possible that the molecular vs. extended structure of the VO2 assemblies is responsible for the observed differences in mechanism, this study poses important questions regarding the relevance of proton-coupled electron transfer in structurally dynamic materials. Due to the implications of the presented work on the rationalization of the consequences of hydrogen-doping in VO2, ongoing investigations in our laboratory are focused on the reactivity of POV–alkoxide clusters with hydrogen atoms and their surrogates, with a goal of more precisely defining the thermochemistry of surface-mediated proton coupled electron transfer (i.e. electron–proton co-doping) in VO2.


General considerations

All manipulations were carried out in the absence of water and oxygen using standard Schlenk techniques, or in a UniLab MBraun inert atmosphere drybox under a dinitrogen atmosphere. All glassware was oven dried for a minimum of 4 hours and cooled in an evacuated antechamber prior to use in the drybox. Solvents were dried and deoxygenated on a Glass Contour System (Pure Process Technology, LLC) and stored over activated 3 Å molecular sieves purchased from Fisher Scientific prior to use. Trimethylsilyl trifluoromethylsulfonate (TMSOTf) and cobaltocene (CoCp2) were purchased from Sigma Aldrich and used as received. [nBu4N]2[VIV6O7(OMe)12] (1-V6IVO72−) was prepared according to published procedure.53

1H NMR spectra were recorded at 500 and 400 MHz on Bruker DPX-500 and Bruker DPX-400 MHz spectrometers locked on the signal of deuterated solvents. All chemical shifts were reported relative to the peak of residual H signal in deuterated solvents. CD3CN and CDCl3 were purchased from Cambridge Isotope Laboratories, degassed by three freeze–pump–thaw cycles, and stored over fully activated 3 Å molecular sieves. Mass spectrometry analyses were performed on an Advion ExpressionL Compact Mass Spectrometer equipped with an electrospray probe and an ion-trap mass analyser. Direct injection analysis was employed in all cases with a sample solution in acetonitrile.

Single crystals of 2-V6IVO6(OSiMe3)1− and 3-VIIIV5IVO6(OSiMe3)2− were mounted on the tip of a thin glass optical fiber (goniometer head) and mounted on a XtaLab Synergy-S Dualflex diffractometer equipped with a HyPix-6000HE HPC area detector for data collection at 100.00(10) K. The structures were solved using SHELXT-2018/2 (ref. 74) and refined using SHELXL-2018/3.75 Elemental analyses were performed on a PerkinElmer 2400 Series II Analyzer, at the CENTC Elemental Analysis Facility, University of Rochester.

The XPS measurements of complexes 1-V6IVO72− and 2-V6IVO6(OSiMe3)1− were recorded with a Kratos Axis Ultra DLD system equipped with a monochromatic Al Kα ( = 1486.6 eV) radiation. During the measurements, pressure in the main chamber was kept below 1 × 10−7 mbar. Charge compensation was carried out via a neutraliser with current: 7 × 10−6 A, charge balance: 5 eV and filament bias: 1.3 V. The X-ray gun was set to 10 mA emission. Sample preparation was performed under a dinitrogen atmosphere glovebox. The powders were dissolved in dry dichloromethane to obtain a concentrated solution. The solution was drop-casted on cleaned Si wafers, which were grounded to the sample bar by carbon tape. Survey scans from 0 to 1200 eV were carried out to identify the elements present in the sample. Binding energies were referenced to the C 1s peak arising from adventitious carbon with binding energy of 284.8 eV.24 The C1s, O1s, Si 2p and V 2p core levels were recorded with a pass energy of 80 eV. 3 scans for vanadium and 2 scans for non-metals were performed. Multiple sample points were used to ensure reproducibility. The V 2p region was increased to 510–540 eV to simultaneously measure the O 1s and V 2p signal in one energy window. XPS data analysis was performed with CasaXPS (Version 2.3.1). The Touggard function was used for background subtraction. The O 1s and V 2p XPS signals were fitted with Gaussian curves. Binding energies were obtained from the same peak fits.

Concentrations of active species for electrochemical analysis (vanadium oxide cluster) and [nBu4N][PF6] used were 1 and 100 mM, respectively, in dichloromethane. Prior to running electrochemical experiments, the supporting electrolyte was recrystallised three times from ethanol and stored under dynamic vacuum. CV measurements were carried out using a Bio-Logic SP 150 potentiostat/galvanostat and the EC-Lab software suite. Glassy carbon disks (3 mm, CH Instruments, USA) were used as working electrodes. Working electrodes were polished using a microcloth pad and 0.05 μm alumina powder. Potentials recorded during CV were measured relative to a nonaqueous Ag/Ag+ reference electrode with 10 mM AgNO3 and 100 mM [nBu4N][PF6] in acetonitrile (Bio-Logic). A platinum wire served as the counter electrode. All experiments were carried out inside a nitrogen-filled glovebox (MBraun, USA). All CV measurements were IR compensated at 85% with impedance taken at 100 kHz using the ZIR tool included with the EC-Lab software. All redox events were referenced against a ferrocenium/ferrocene (Fc+/0) redox couple.

Synthesis of [nBu4N][VIV6O6(OSiMe3)(OMe)12] (2-V6IVO6(OSiMe3)1−). In a glove box, a 20 mL scintillation vial was charged with 1-V6IVO72− (0.040 g, 0.030 mmol) and 4 mL acetonitrile. TMSOTf (0.007 g, 0.030 mmol) was dissolved in 2 mL acetonitrile in a separate 20 mL scintillation vial. Both solutions were frozen in a liquid nitrogen cold well. The frozen solutions were taken out of the cold well and TMSOTf was added dropwise with a glass pipette in three portions to the frozen slurry of 1-V6IVO72− as it thawed. The colour of the solution changed from teal to dark green upon the addition of TMSOTf. The reaction was stirred for an additional 10 minutes after the complete addition of TMSOTf. Residual solvent was subsequently removed under reduced pressure to give a green solid. The crude solid was triturated with small portions of pentane (4 mL) followed by washing with 2 mL of diethyl ether. Dichloromethane was then used to extract the product. After dichloromethane was removed under reduced pressure, the product, 2-V6IVO6(OSiMe3)1−, was isolated as a dark green solid in good yield (85%, 0.035 g, 0.027 mmol). 1H NMR (500 MHz, CD3CN): δ = 27.50, 24.12, 3.06, 1.57, 1.34, 0.96, 0.76, −9.79 ppm. Elemental analysis for C31H81NO19V6Si 0.5 CH2Cl2 (MW = 1148.17 g mol−1) calcd (%): C, 32.95; H, 7.20; N, 1.22. Found (%): C, 32.99; H, 7.21; N 0.89.
Synthesis of [nBu4N][CoCp2][VIIIVIV5O6(OSiMe3)(OMe)12] (3-VIIIV5IVO6(OSiMe3)2−). A 20 mL scintillation vial was charged with complex 2-V6IVO6(OSiMe3)1− (0.055 g, 0.050 mmol) and 4 mL of dichloromethane. In a separate vial, CoCp2 (0.009 g, 0.050 mmol) was dissolved in 2 mL of dichloromethane. Both solutions were frozen in a liquid nitrogen cold well. While thawing, the solution of CoCp2 was added in three parts, drop-wise, to the solution of complex 2-V6IVO6(OSiMe3)1−. Following the complete addition of the reductant, the reaction mixture was stirred for an additional 5 minutes. Solvent was removed under reduced pressure, leaving a dark coloured residue at the bottom of the vial. The product was extracted with cold tetrahydrofuran and filtered through a bed of Celite. After removing tetrahydrofuran under vacuum, the product, complex 3-VIIIV5IVO6(OSiMe3)2−, was isolated as a blue-green solid in moderate yield (23%, 0.015 g, 0.011 mmol). Crystals suitable for X-ray analysis were grown from vapour diffusion of diethyl ether into a concentrated solution of 3-VIIIV5IVO6(OSiMe3)2− in acetonitrile; refinement of structural data revealed two tetrabutylammonium ions, suggesting cation exchange during crystallization. 1H NMR (500 MHz, CD3CN): δ = 24.77, 22.75, 3.07, 1.58, 1.33, 0.96, −8.41 ppm. Elemental analysis for C47H117N2O19V6Si·CH2Cl2 (MW = 1348.17 g mol−1) calcd (%): C, 40.23; H, 8.37; N, 1.95. Found (%): C, 40.47; H, 8.07; N 1.82.

Data availability

All experimental data supporting the findings in this article have been uploaded as part of the ESI.

Author contributions

S. C., E. S., and E. M. M. conceived and planned the experiments. S. C. and E. S. performed the synthesis of complexes and collected all experimental data, except those involving X-ray crystallography and X-ray photoelectron spectroscopy. K. R. S.-L., M. T., and K. E. K. analyzed clusters via X-ray photoelectron spectroscopy. W. W. B. collected single crystal X-ray diffraction data and solved the crystal structures. All authors contributed to the writing of the manuscript.

Conflicts of interest

There are no conflicts to declare.


S. C., E. S., and E. M. M. gratefully acknowledge financial support from the National Science Foundation through grant CHE-165195. E. M. M. is also a recipient of a Cottrell Scholar Award from the Research Corporation for Science Advancement and a Sloan Fellowship, which have provided additional resources to support this research. K. R. S. -L., M. T., and K. E. K. acknowledge funding from the American Chemical Society's Petroleum Research Fund (60898-DN13) and the National Science Foundation (CHE-2044462). The authors also acknowledge the UR Nano facilities for technical support with X-ray Photoelectron Spectroscopy experiments.


  1. C. Wu, F. Feng and Y. Xie, Chem. Soc. Rev., 2013, 42, 5157–5183 RSC .
  2. W. Wang, B. Jiang, W. Xiong, H. Sun, Z. Lin, L. Hu, J. Tu, J. Hou, H. Zhu and S. Jiao, Sci. Rep., 2013, 3, 3383 CrossRef PubMed .
  3. C. Niu, Funct. Mater. Lett., 2018, 11, 1850068 CrossRef CAS .
  4. S. Han, Z. Zou and S. Huo, Mater. Sci. Eng., 2018, 382, 022059 Search PubMed .
  5. Y. Cui, Y. Ke, C. Liu, Z. Chen, N. Wang, L. Zhang, Y. Zhou, S. Wang, Y. Gao and Y. Long, Joule, 2018, 2, 1707–1746 CrossRef CAS .
  6. A. Riapanitra, Y. Asakura and S. Yin, Funct. Mater. Lett., 2020, 13, 1951008 CrossRef CAS .
  7. C. Xia, Z. Lin, Y. Zhou, C. Zhao, H. Liang, P. Rozier, Z. Wang and H. N. Alshareef, Adv. Mater., 2018, 30, 1803594 CrossRef PubMed .
  8. Y. Chen, Z. Wang, S. Chen, H. Ren, L. Wang, G. Zhang, Y. Lu, J. Jiang, C. Zou and Y. Luo, Nat. Commun., 2018, 9, 1–8 CrossRef PubMed .
  9. T. E. G. Alivio, D. G. Sellers, H. Asayesh-Ardakani, E. J. Braham, G. A. Horrocks, K. E. Pelcher, R. Villareal, L. Zuin, P. J. Shamberger, R. Arróyave, R. Shahbazian-Yassar and S. Banerjee, Chem. Mater., 2017, 29, 5401–5412 CrossRef CAS .
  10. H. Yoon, M. Choi, T.-W. Lim, H. Kwon, K. Ihm, J. K. Kim, S.-Y. Choi and J. Son, Nat. Mater., 2016, 15, 1113–1119 CrossRef CAS PubMed .
  11. J. Zhang, H. He and B. Pan, Phys. Chem. Chem. Phys., 2013, 15, 4687–4690 RSC .
  12. J. Chenavas, J. C. Joubert, J. J. Capponi and M. Marezio, J. Solid State Chem., 1973, 6, 1–15 CrossRef CAS .
  13. J. Wei, H. Ji, W. Guo, A. H. Nevidomskyy and D. Natelson, Nat. Nanotechnol., 2012, 7, 357–362 CrossRef CAS .
  14. Y. Filinchuk, N. A. Tumanov, V. Ban, H. Ji, J. Wei, M. W. Swift, A. H. Nevidomskyy and D. Natelson, J. Am. Chem. Soc., 2014, 136, 8100–8109 CrossRef CAS PubMed .
  15. B. Li, L. Xie, Z. Wang, S. Chen, H. Ren, Y. Chen, C. Wang, G. Zhang, J. Jiang and C. Zou, Angew. Chem., Int. Ed., 2019, 58, 13711–13716 CrossRef CAS PubMed .
  16. C. F. Wise and J. M. Mayer, J. Am. Chem. Soc., 2019, 141, 14971–14975 CrossRef CAS .
  17. M. T. Pope, A. Muller, Polyoxometalate Chemistry: From Topology via Self-Assembly to Applications, Berlin, 2001 Search PubMed .
  18. N. I. Gumerova and A. Rompel, Nat. Rev. Chem., 2018, 2, 0112 CrossRef CAS .
  19. X. Lopez, J. J. Carbó, C. Bo and J. M. Poblet, Chem. Soc. Rev., 2012, 41, 7537–7571 RSC .
  20. P. Gouzerh and A. Proust, Chem. Rev., 1998, 98, 77–112 CrossRef CAS PubMed .
  21. V. W. Day, W. G. Klemperer and D. J. Maltbie, J. Am. Chem. Soc., 1987, 109, 2991–3002 CrossRef CAS .
  22. Q. Chen, D. P. Goshorn, C. P. Scholes, X. L. Tan and J. Zubieta, J. Am. Chem. Soc., 1992, 114, 4667–4681 CrossRef CAS .
  23. K. Kastner, J. T. Margraf, T. Clark and C. Streb, Chem. - Eur. J., 2014, 20, 12269–12273 CrossRef CAS .
  24. B. E. Petel, R. L. Meyer, M. L. Maiola, W. W. Brennessel, A. M. Müller and E. M. Matson, J. Am. Chem. Soc., 2020, 142, 1049–1056 CrossRef CAS PubMed .
  25. M. L. Maiola, B. E. Petel, W. W. Brennessel and E. M. Matson, Dalton Trans., 2020, 49, 16184–16192 RSC .
  26. B. E. Petel and E. M. Matson, Chem. Commun., 2020, 56, 13477–13490 RSC .
  27. C. Wu, Z. Hu, W. Wang, J. Yang and Y. Xie, Chem. Commun., 2008, 3891–3893 RSC .
  28. L. Liu, F. Cao, T. Yao, Y. Xu, M. Zhou, B. Qu, B. Pan, C. Wu, S. Wei and Y. Xie, New J. Chem., 2012, 36, 619–625 RSC .
  29. S. R. Popuri, M. Miclau, A. Artemenko, C. Labrugere, A. Villesuzanne and M. Pollet, Inorg. Chem., 2013, 52, 4780–4785 CrossRef CAS PubMed .
  30. S. Lee, I. N. Ivanov, J. K. Keum and H. N. Lee, Sci. Rep., 2016, 6, 19621 CrossRef CAS PubMed .
  31. E. Schreiber, B. E. Petel and E. M. Matson, J. Am. Chem. Soc., 2020, 142, 9915–9919 CrossRef CAS PubMed .
  32. I. Fleming, Chem. Soc. Rev., 1981, 10, 83 RSC .
  33. S. Murata, M. Suzuki and R. Noyori, J. Am. Chem. Soc., 1980, 102, 3248–3249 CrossRef CAS .
  34. R. Noyori, S. Murata and M. Suzuki, Tetrahedron, 1981, 37, 3899–3910 CrossRef CAS .
  35. J. K. Pagano, J. M. Dorhout, R. Waterman, K. R. Czerwinski and J. L. Kiplinger, Chem. Commun., 2015, 51, 17379–17381 RSC .
  36. J. Chu, T. G. Carroll, G. Wu, J. Telser, R. Dobrovetsky and G. Ménard, J. Am. Chem. Soc., 2018, 140, 15375–15383 CrossRef CAS PubMed .
  37. B. E. Petel, W. W. Brennessel and E. M. Matson, J. Am. Chem. Soc., 2018, 140, 8424–8428 CrossRef CAS PubMed .
  38. B. E. Petel, A. A. Fertig, M. L. Maiola, W. W. Brennessel and E. M. Matson, Inorg. Chem., 2019, 58, 10462–10471 CrossRef CAS PubMed .
  39. B. E. Petel, R. L. Meyer, W. W. Brennessel and E. M. Matson, Chem. Sci., 2019, 10, 8035–8045 RSC .
  40. J. K. Money, K. Folting, J. C. Huffman, D. Collison, J. Temperley, F. E. Mabbs and G. Christou, Inorg. Chem., 1986, 25, 4583–4589 CrossRef CAS .
  41. M. Rost, H. Görls, W. Imhof, W. Seidel and K. Thiele, Z. Anorg. Allg. Chem., 1998, 624, 1994–2000 CrossRef CAS .
  42. W. H. Knoth, J. Am. Chem. Soc., 1979, 101, 759–760 CrossRef CAS .
  43. R. Carlisle Chambers, E. J. Osburn Atkinson, D. McAdams, E. J. Hayden and D. J. Ankeny Brown, Chem. Commun., 2003, 2456–2457 RSC .
  44. I. Bar-Nahum, H. Cohen and R. Neumann, Inorg. Chem., 2003, 42, 3677–3684 CrossRef CAS PubMed .
  45. A. Mazeaud, Y. Dromzee and R. Thouvenot, Inorg. Chem., 2000, 39, 4735–4740 CrossRef CAS PubMed .
  46. R. C. Schroden, C. F. Blanford, B. J. Melde, B. J. S. Johnson and A. Stein, Chem. Mater., 2001, 13, 1074–1081 CrossRef CAS .
  47. A. Mazeaud, N. Ammari, F. Robert and R. Thouvenot, Angew. Chem., Int. Ed., 1996, 35, 1961–1964 CrossRef CAS .
  48. V. Duffort, R. Thouvenot, C. Afonso, G. Izzet and A. Proust, Chem. Commun., 2009, 6062–6064 RSC .
  49. B. Matt, S. Renaudineau, L. M. Chamoreau, C. Afonso, G. Izzet and A. Proust, J. Org. Chem., 2011, 76, 3107–3112 CrossRef CAS .
  50. A. J. Kibler, C. Martín, J. M. Cameron, A. Rogalska, J. Dupont, D. A. Walsh and G. N. Newton, Eur. J. Org. Chem., 2019, 2019, 456–460 CrossRef CAS .
  51. V. W. Day, W. G. Klemperer and C. Schwartz, J. Am. Chem. Soc., 1987, 109, 6030–6044 CrossRef CAS .
  52. C. Wu, F. Feng, J. Feng, J. Dai, L. Peng, J. Zhao, J. Yang, C. Si, Z. Wu and Y. Xie, J. Am. Chem. Soc., 2011, 133, 13798–13801 CrossRef CAS PubMed .
  53. J. Spandl, C. Daniel, I. Brüdgam and H. Hartl, Angew. Chem., Int. Ed., 2003, 42, 1163–1166 CrossRef CAS PubMed .
  54. C. Aronica, G. Chastanet, E. Zueva, S. A. Borshch, J. M. Clemente-Juan and D. Luneau, J. Am. Chem. Soc., 2008, 130, 2365–2371 CrossRef CAS PubMed .
  55. J. Mendialdua, R. Casanova and Y. Barbaux, J. Electron Spectrosc. Relat. Phenom., 1995, 71, 249–261 CrossRef CAS .
  56. G. Silversmit, D. Depla, H. Poelman, G. B. Marin and R. De Gryse, J. Electron Spectrosc. Relat. Phenom., 2004, 135, 167–175 CrossRef CAS .
  57. J. Tian, F. Liu, C. Shen, H. Zhang, T. Yang, L. Bao, X. Wang, D. Liu, H. Li, X. Huang, J. Li, L. Chen and H. Gao, J. Mater. Res., 2007, 22, 1921–1926 CrossRef CAS .
  58. D. J. O'Connor, B. A. Sexton and R. S. C. Smart, Surface analysis methods in materials science, Springer, Berlin ; New York, 2nd edn, 2003 Search PubMed .
  59. C. Daniel and H. Hartl, J. Am. Chem. Soc., 2009, 131, 5101–5114 CrossRef CAS PubMed .
  60. C. Daniel and H. Hartl, J. Am. Chem. Soc., 2005, 127, 13978–13987 CrossRef CAS PubMed .
  61. F. Li, S. H. Carpenter, R. F. Higgins, M. G. Hitt, W. W. Brennessel, M. G. Ferrier, S. K. Cary, J. S. Lezama-Pacheco, J. T. Wright, B. W. Stein, M. P. Shores, M. L. Neidig, S. A. Kozimor and E. M. Matson, Inorg. Chem., 2017, 56, 7065–7080 CrossRef CAS PubMed .
  62. F. Li, R. L. Meyer, S. H. Carpenter, L. E. VanGelder, A. W. Nichols, C. W. Machan, M. L. Neidig and E. M. Matson, Chem. Sci., 2018, 9, 6379–6389 RSC .
  63. L. E. VanGelder, A. M. Kosswattaarachchi, P. L. Forrestel, T. R. Cook and E. M. Matson, Chem. Sci., 2018, 9, 1692–1699 RSC .
  64. L. E. VanGelder and E. M. Matson, J. Mater. Chem. A, 2018, 6, 13874–13882 RSC .
  65. N. G. Connelly and W. E. Geiger, Chem. Rev., 1996, 96, 877–910 CrossRef CAS PubMed .
  66. A. S. Veige, L. M. Slaughter, E. B. Lobkovsky, P. T. Wolczanski, N. Matsunaga, S. A. Decker and T. R. Cundari, Inorg. Chem., 2003, 42, 6204–6224 CrossRef CAS PubMed .
  67. J. N. Schrauben, R. Hayoun, C. N. Valdez, M. Braten, L. Fridley and J. M. Mayer, Science, 2012, 336, 1298–1301 CrossRef CAS PubMed .
  68. C. N. Valdez, M. F. Delley and J. M. Mayer, J. Am. Chem. Soc., 2018, 140, 8924–8933 CrossRef CAS PubMed .
  69. S. M. Laga, T. M. Townsend, A. R. O'Connor and J. M. Mayer, Inorg. Chem. Front., 2020, 7, 1386–1393 RSC .
  70. R. G. Agarwal, H.-J. Kim and J. M. Mayer, J. Am. Chem. Soc., 2021, 143, 2896–2907 CrossRef CAS PubMed .
  71. F. Li, L. E. VanGelder, W. W. Brennessel and E. M. Matson, Inorg. Chem., 2016, 55, 7332–7334 CrossRef CAS PubMed .
  72. R. L. Meyer, W. W. Brennessel and E. M. Matson, Polyhedron, 2018, 156, 303–311 CrossRef CAS .
  73. L. E. VanGelder, W. W. Brennessel and E. M. Matson, Dalton Trans., 2018, 47, 3698–3704 RSC .
  74. G. Sheldrick, Acta Crystallogr., Sect. A: Found. Adv., 2015, 71, 3–8 CrossRef PubMed .
  75. G. Sheldrick, Acta Crystallogr., Sect. C: Struct. Chem., 2015, 71, 3–8 Search PubMed .


Electronic supplementary information (ESI) available: 1H NMR, ESI-MS (−ve), cyclic voltammogram, and X-ray photoelectron spectroscopy of 2-VIV6O6(OSiMe3)1−, 1H NMR spectra of dichloromethane solution of 2-VIV6O6(OSiMe3)1− in time intervals, 1H NMR spectra of 1-V6IVO72− and 3-VIIIV5IVO6(OSiMe3)2−, crystallographic parameters, and bond valence sum calculations of 2-V6IVO6(OSiMe3)1− and 3-VIIIV5IVO6(OSiMe3)2−, X-ray crystal structures of 2-V6IVO6(OSiMe3)1− and 3-VIIIV5IVO6(OSiMe3)2−. CCDC 2058288 and 2058289. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/d1sc02809j
Authors contributed equally to this work.

This journal is © The Royal Society of Chemistry 2021