Open Access Article
Joan C. Borea,
S. M. Al Rafata,
Clara A. Hoffert
a,
Allison A. Ellisa,
Cameron Blakea,
Cristian Celis-Barros
b,
Wei-Yuan Chen
a,
Aliaksei Boika
a,
Briana R. Schrage
*b and
Christopher J. Ziegler
*a
aDepartment of Chemistry, University of Akron, Akron, Ohio 44325, USA
bRadioisotope Science and Technology Division, Oak Ridge National Laboratory, 1 Bethel Valley Rd., Oak Ridge, Tennessee 37830, USA
First published on 2nd June 2026
Isoindoline-based chelates, in particular bis(arylimino)isoindolines, have shown extensive metal binding chemistry. Although this chemistry has been explored for the middle and late transition metal ions, little work has been carried out on early transition metal complexes. In this article, we present the first examples of vanadium coordinated using four bis(arylimino)isoindolines, in which the aryl groups are pyrazole, indazole, benzimidazole, and pyridine (ligands 1–4, respectively). We isolated five complexes using vanadyl sulfate or vanadyl acetylacetonate as the vanadium source. In all cases, the ligands bound in a meridional mode, and for four of the complexes, the vanadium ion was observed in the V(V) oxidation state. Three of the ligands (1–3) formed VO2 complexes with vanadyl sulfate and vanadyl acetylacetonate, but the bis(pyridylimino)isoindoline (ligand 4) formed a V(V) oxosulfonato complex with the former starting material and a vanadyl V(IV) acetylacetonate with the latter starting material. All metal compounds were structurally elucidated by X-ray crystallographic methods, and we probed their electronic structures using DFT methods.
This metal has also been employed for energy storage applications, specifically as the most well-known and well-studied single-component redox flow battery material.14–18 The vanadium redox flow battery utilizes vanadium's ability to exist in four oxidation states, with V(II)/V(III) and V(IV)/V(V) redox pairs existing in the anode and cathode respectively. Finally, in radiochemistry, 48V could be used as a potential imaging or therapeutic isotope; this nuclide is a beta emitter with a 16-day half-life.19,20 Vanadium is also used as a redox surrogate for neptunium because they have similar redox potentials vs. standard H electrode for the VO2+/VO2+ and NpO22+/NpO2+ couples at 1.0 V and 1.15 V respectively.21
For several years, we have been investigating the coordination chemistry of isoindoline-based chelates.22–30 Isoindoline is one of the subunits of phthalocyanine, the aromatic synthetic dye used as a bulk colorant and a specialized material.31–33 Since the 1950s, isoindoline chelates such as bis(pyridylimino)isoindoline (BPI) have been studied for their metal binding chemistry,28,34–37 and their metal complexes have been explored for their biomimetic properties,38,39 as catalysts for organic reactions,40,41 and for their novel electronic properties.42,43 However, much of this chemistry has focused on the middle and late transition metal ions; notably, few investigations to date have studied early transition metals. In particular, other than a few reports on its chemistry with molybdenum,44,45 no examples exist of BPI or other bis(arylimino)isoindoline complexes with metals from groups 3 to 6.
In this article, we present the first study into the vanadium chemistry of a series of chelating isoindoline-based ligands. These compounds, shown in Fig. 1, incorporate pyrazole, indazole, benzimidazole, and pyridine as peripheral heterocycles. All four ligands react rapidly with vanadyl sulfate and vanadyl acetylacetonate to afford 1
:
1 heteroleptic complexes. The ligands behave as monoanionic, tridentate ligands that bind in a planar geometry. Notably, all products with one exception exhibit vanadium in the +5 oxidation state. Three of the ligands form VO2 complexes upon reaction with vanadyl sulfate and vanadyl acetylacetonate. The BPI ligand, in contrast, forms a VO(SO4) complex in the former case and a vanadyl, or V(IV), acetylacetonate in the latter case. All compounds have been fully characterized by X-ray crystallography, spectroscopic methods, and we have elucidated their electronic structures with calculations. These compounds all represent a new class of tridentate, nitrogenous ligand chemistry with vanadium.
The structures of the four bis(arylimino)isoindolines are shown in Fig. 1. The four heterocycles that we incorporated into these ligands included pyrazole (1), indazole (2), benzimidazole (3), and pyridine (4). All four compounds had been previously synthesized and can be readily produced using the Siegl method with anhydrous CaCl2 as a template and catalyst.46 Of these four, the structures of compounds 1, 2, and 4 were previously elucidated,22,47 and we were able to determine the structure of compound 3, which is shown in Fig. S11. The structure of 3 shares many of the attributes present in heterocycle-modified diiminoisoindolines. First, the compounds are all planar with three nitrogen atoms facing the core of the molecule, which is the metal binding site. In all cases, the central isoindoline nitrogen atom is protonated and stabilized by hydrogen bonding interactions with the flanking heterocycle unprotonated nitrogen atom positions. For the three compounds with ionizable heterocycles (1–3), the protonated nitrogen positions do not face the central core; in the case of 3, the ionizable protons reside on the external nitrogen positions. When these ligands bind metal ions, they most typically do so in a meridional fashion, but we observed in our prior work on the chemistry of rhenium that these ligands can be deformed to coordinate in a facial configuration.28 In the case of the Re(CO)3 complexes, the coordination mode is enforced by the rigid facial configuration of the carbonyl groups.
We reacted ligands 1–4 with two vanadyl starting materials: vanadyl sulfate and vanadyl acetylacetonate. All the reactions afforded crystalline solids that ranged in color from yellow to red-orange, and the chemistry reaction observed for these ligands is shown in Scheme 1. For ligands 1–3, we observed the formation of the same product using both the vanadyl sulfate and vanadyl acetylacetonate starting material: compounds 1VO2–3VO2 (Fig. 2). The structures of these three compounds are shown in Fig. 3. In each complex, a single equivalent of ligand binds to the vanadium center in a tridentate, planar geometry. The central isoindoline nitrogen atom is deprotonated, resulting in a monoanionic chelate. For 1VO2–3VO2, the remainder of the vanadium coordination sphere is occupied by two O atoms, which we could assign as oxo groups. The V–O distances range from approximately 1.61 to 1.63 Å, which clearly falls in the V–O double bond range. The V–N bonds range between approximately 2.05 and 2.12 Å, with the central isoindoline V–N bond longer than the flanking V–N bonds. The metals in these complexes exhibit distorted trigonal bipyramidal geometries, with the N–V–N bonds from the flanking heterocycles occupying the axial positions. The axial N–V–N bonds are smaller than the ideal 180° (approximately 150°–165°), and the equatorial O–V–N bonds are, on average, larger than 120° (approximately 125°), with smaller O–V–O angles (approximately 109°). Based on structural features, we could assign the oxidation states in these three complexes as V(V). Additionally, 1VO2–3VO2 exhibited readily interpretable diamagnetic 1H nuclear magnetic resonance (NMR) spectra.
![]() | ||
| Fig. 3 The structures of compounds 1VO2–3VO2 with 35% thermal ellipsoids. Hydrogen atoms on carbon positions have been omitted for clarity. | ||
For bis(pyridyl) ligand 4, we observed reactivity dependent on the starting material, which is shown in Scheme 1. When VO(SO4) was employed as the vanadium source, the reaction resulted in a V(V) complex, 4VO(SO4), with a sulfate bound at the metal center, as shown on the left side of Fig. 4. In this compound, rather than having a second vanadium oxo bond, a sulfate binds in a bidentate fashion with V–O bond lengths of approximately 1.52 and 1.49 Å. The vanadium oxo bond is slightly shorter than that seen in 1VO2–3VO2, with a distance of 1.59 Å. BPI binds in a planar tridentate fashion, with the central isoindoline deprotonated, as seen in 1VO2–3VO2. One notable difference in 4VO(SO4) is that the V–N bond to the isoindoline nitrogen is shorter (approximately 2.03 Å) than that which is seen with the flanking pyridine groups (approximately 2.13 Å). The diamagnetic character for 4VO(SO4) was confirmed by 1H NMR spectroscopy, which revealed standard ligand chemical shifts.
![]() | ||
| Fig. 4 Structures of (left) 4VO(SO4) and (right) 4VO(acac) with 35% thermal ellipsoids hydrogen atom positions on carbon atoms have been omitted for clarity. | ||
When the BPI chelate reacted with VO(acac)2 as the vanadium source, we observed the formation of a vanadyl complex (4VO(acac)), with the metal in the +4 oxidation state. The structure of this lower-valent complex is also shown in Fig. 4. With a monoanionic BPI, a vanadyl unit, and a bidentate acetylacetonate, this compound clearly has V(IV) character. When compared with the similarly structured but higher-valent 4VO(SO4), the vanadium oxo bond is slightly longer (approximately 1.61 Å; this may not be chemically significant) but the V–O bonds to the acetylacetonate are much longer than those seen in the sulfate in 4VO(SO4), measuring approximately 2.02 and 2.17 Å.
Next, the FTIR spectra of the vanadium complexes were acquired and all V(V) complexes show a band between 961–984 cm−1, attributed to the vanadyl VO2 unit. The V(IV) VO mode in the 4VO(acac) complex is found at 983 cm−1. Additionally, the ligand imine stretching vibration bands are found between 1622–1645 cm−1, these stretches are similarly reported in other bisaryliminoisoindoline ligands.22 Finally, in the 4VO complexes, the sulfate containing complex portrays S
O stretching at 1062 and 1034 cm−1 as well as a SO4 bending vibration at 607 cm−1. The acac complex of 4VO reveals carbonyl stretching bands at 1575 and 1548 cm−1 due to the bidentate carbonyl anion.
Spectroscopically, we found evidence that the oxidation state in 4VO(acac) is V(IV), but we also observed that this complex can be readily oxidized in solution to V(V). In the solid state and in dimethylformamide (DMF) glass at low temperature, 4VO(acac) exhibited an electron paramagnetic resonance (EPR) spectrum consistent with a V(IV) ion in a rhombic field with hyperfine coupling due to the 7/2 nuclear spin (Fig. S6). In contrast, when 4VO(acac) was dissolved in dimethyl sulfoxide (DMSO) or CDCl3 in air, we observed diamagnetic chemical shifts in its NMR spectrum. DMSO is slightly oxidizing in potential, so it is not surprising that 4VO(acac) is oxidized under these conditions.
Surprisingly, only a few meridional nitrogenous ligand complexes of VO2 have been structurally characterized. Several terpyridine complexes of VO2 have been elucidated, but unlike the compounds presented in this report, the (tpy)VO2 compounds are all cationic, with slightly longer V–N bonds than those seen in 1VO2–3VO2.48–52 Notably, the central V–N bond is shorter than those seen in the flanking V–N bonds. Most tridentate planar VO2 complexes reported in the literature are anionic and asymmetric. The symmetric anionic complexes include a deprotonated bis(2-pyridylcarbonyl)amine VO2 complex,53 with V–N bonds of approximately 2.08 Å for the flanking positions and approximately 2.07 Å for the central amide unit, as well as two expanded porphyrin systems, in which the VO2 unit is bound to a tripyrrane segment of the macrocycle.54,55 In the latter case, the bond lengths vary widely, as seen in the flanking V–N bonds that range from approximately 2.08 to 2.15 Å, and may result from macrocyclic steric factors. In the pyrrane compounds, the central V–N bond is the shortest of the three, measuring approximately 2.04 and 2.06 Å in the two compounds. The asymmetric anionic ligand systems belong to two types: aryl-substituted amide with a deprotonated amide N-bound56 and azole-modified, hydrazine-based ligands.57 In both types of systems, the V–N bond lengths vary significantly between the different nitrogen donors, but typically, the shortest V–N bond is observed with the formally deprotonated nitrogen position.
We can also compare the vanadium structures of these isoindoline ligands with porphyrin and phthalocyanine vanadyl complexes. In the recently elucidated V(IV) vanadyl etioporphyrin structure, the V–O bond measures ∼1.59 Å and the V–N bonds measure approximately 2.07 Å.58 In vanadyl tetraphenyl porphyrin these bonds are similar in length, ∼1.62 and ∼2.08 Å respectively.59 The isoindoline-based phthalocyanine macrocycle has a similar axial V–O bond (∼1.58 Å) but appreciably shorter V–N bonds (∼2.03 Å), similar to what we observe in our tridentate chelate systems.60
All vanadium complexes are colored and strongly absorb in the visible region. The ultraviolet (UV)–visible spectra of the VO2 and VO complexes are shown in Fig. 5. The spectra of the free base ligands 1–4 are shown in Fig. S7–S10 and exhibit π → π* transitions with a vibronic fine structure.22 Complexation of the ligands with vanadium resulted in a bathochromic shift in the spectra. In the benzimidazole derivative 3VO2, a red-shifted shoulder is shown in the spectrum. Extinction coefficients range between 1.1 × 104 and 2 × 104 M−1 cm−1 for these compounds.
![]() | ||
| Fig. 5 The UV–visible spectra of the bis(arylimino)isoindoline vanadium complexes in DMF solution. Top spectra show 1VO2–3VO2, and bottom spectra show 4VO(SO4) and 4VO(acac). | ||
Finally, in order to further understand the electronic structure of the vanadium adducts, we carried out density functional theory (DFT) calculations. The calculated HOMO/LUMO energies and gaps of all vanadium complexes are summarized in Table 1. Complexes 1VO2–3VO2 display broadly similar frontier orbital energies, with HOMOs between −6.042 and −5.759 eV and LUMOs clustered near −3.19 to −3.13 eV. Their HOMO–LUMO separations (ΔEH–L) fall in a narrow window (2.57–2.92 eV).
| Complex | HOMO | LUMO | ΔEH–L |
|---|---|---|---|
| 1VO2 | −6.042 | −3.125 | −2.917 |
| 2VO2 | −5.759 | −3.186 | −2.573 |
| 3VO2 | −5.955 | −3.147 | −2.808 |
| 4VO(SO4) | −6.326 | −5.575 | −0.751 |
| 4VO(acac) | −5.601 | −3.088 | −2.513 |
Within this set, 2VO2 shows the smallest gap (2.573 eV), mainly due to a higher-lying HOMO (−5.759 eV) compared to 1VO2 and 3VO2, while the LUMOs remain essentially constant across the three. This pattern suggests that substituent- or ligand-dependent effects modulate the donor character (and thus HOMO energy) more strongly than the acceptor manifold (LUMO), which appears largely governed by the vanadium oxo fragment and the conserved coordination environment.
A notably different behavior is observed for 4VO(SO4), which exhibits a dramatically reduced HOMO–LUMO gap (0.751 eV), driven primarily by a much lower LUMO energy (−5.575 eV). This indicates stabilization of low-lying acceptor orbitals relative to the other complexes and implies a substantially altered electronic structure for the sulfate-containing species. In contrast, 4VO(acac) resembles the 1VO2–3VO2 series more closely (ΔEH–L = 2.513 eV) possibly due to lower oxidation state of the vanadium, where the pyridine-based isoindoline ligand cannot stabilize the 3d orbitals enough as in 4VO(SO4) (Fig. 6). It is noteworthy that looking solely at the HOMO–LUMO gaps will not be able to qualitatively explain the UV-visible spectrum, especially in the case of 4VO(SO4). However, when the gaps between π → π* transitions as shown in Fig. 6, can qualitatively explain the similarities in their absorption spectra.
Next, the energy decomposition analysis (EDA) results (Table 2) provide insight into the balance of stabilizing and destabilizing contributions governing ligand binding. Across all complexes, the Pauli term is (as expected) strongly destabilizing (+164 to +184 kcal mol−1), while the interaction is stabilized by electrostatics (−263 to −296 kcal mol−1), orbital interactions (–157 to −200 kcal mol−1), and a smaller but non-negligible dispersion term (−8 to −16 kcal mol−1). The total interaction energies (ΔEint) follow: 4VO(SO4) (−327.8 kcal mol−1) > 1VO2 (−309.9 kcal mol−1) ≈ 3VO2 (–301.6 kcal mol−1) ≈ 2VO2 (–300.4 kcal mol−1) ≫ 4VO(acac) (–251.3 kcal mol−1).
| Complex | Pauli | Electrostatic | Orbital | Dispersion | Total |
|---|---|---|---|---|---|
| 1VO2 | 172.0 | −295.9 | −177.9 | −8.1 | −309.9 |
| 2VO2 | 164.1 | −279.2 | −177.0 | −8.4 | −300.4 |
| 3VO2 | 163.8 | −287.4 | −167.8 | −10.2 | −301.6 |
| 4VO(SO4) | 166.7 | −280.6 | −199.9 | −14.0 | −327.8 |
| 4VO(acac) | 183.5 | −262.6 | −156.8 | −15.5 | −251.3 |
The three dioxovanadate complexes display closely comparable total interaction energies (−300 to −310 kcal mol−1), indicating similar overall metal–ligand binding strengths under the chosen fragmentation. Their stabilization is dominated by a combination of electrostatic attraction (−279 to −296 kcal mol−1) and substantial covalent/orbital interaction (−168 to −178 kcal mol−1), with dispersion contributing modestly (−8 to −10 kcal mol−1). Small variations within the series track primarily with changes in electrostatics and orbital terms rather than dispersion, consistent with subtle electronic tuning by ligand substitution while preserving the same core VO2+ bonding framework.
4VO(SO4) shows the most stabilizing total interaction (−327.8 kcal mol−1). This also true for the orbital interaction term of the entire set (−199.9 kcal mol−1), along with the most stabilizing electrostatic term among the four complexes (−280.6 kcal mol−1) and enhanced dispersion (−14.0 kcal mol−1). The large orbital contribution points to particularly strong donor–acceptor (covalent) interactions between the isoindoline fragment and the [VO(SO4)]+ unit in this compound, consistent with a distinct electronic structure relative to the [VO2]+-containing complexes. This stronger orbital stabilization is also in line with the unusually low LUMO energy computed for 4VO(SO4), i.e., an electronically more accepting fragment. Though this is also complemented by the quasi-linear geometry of the Ooxo–V–OSO4 which stabilizes through orbital mixing with unoccupied V(V) 3d orbitals as shown in Fig. S16 and S17. Conversely, and as expected, the lower oxidation state in the 4VO(acac) complex causes the same isoindoline ligand show lower magnitudes for the attractive forces and a larger Pauli repulsion compared to 4VO(SO4). This ultimately causes a lower stabilization of the unoccupied 3d orbitals. The open-shell character of 4VO(acac) shows the unpaired 3d electron to be the most unstable electron in the system followed by the π-type molecular orbital seen in the other 4 isoindoline complexes.
NMR spectra were recorded on a 400 MHz spectrometer, and chemical shifts were given in parts per million relative to residual solvent resonances (1H NMR spectra). Infrared spectra were collected on a Thermo Scientific Nicolet iS5 that was equipped with an iD5 attenuated total reflectance detector. UV–visible spectra were recorded on a Shimadzu UV-2600 UV–visible spectrometer. Elemental analyses (C, H, N) were performed at Atlantic Microlab on a CHN elemental analyzer. Mass spectroscopy spectra were performed at UTA Research Facilities. EPR spectra were collected on a Bruker X-band ELEXSYS E-500 instrument at 130 K.
X-ray intensity data were measured on a Bruker D8 Venture diffractometer equipped with an Iμs 3.0 Mo X-ray source (λ = 0.71073 Å) and on a Bruker CCD-based diffractometer with dual Cu/Mo IμS microfocus optics (Cu Kα radiation, λ = 1.54178 Å, Mo Kα radiation, λ = 0.71073 Å). APEX4 software was used for data collection and unit cell determination. Crystals were mounted on a cryoloop using Paratone oil and placed under a stream of N2 at 100 K (Oxford Cryosystems) for 3 and 4VO(acac) and 300 K for all others. The detector was placed at a distance of 5.00 cm from the crystal. The data were corrected for absorption with the SADABS program. The structures were refined using the Bruker SHELXTL Software Package (Version 6.1)61 and were solved using direct methods until the final anisotropic full-matrix, least squares refinement of F2 converged. Disordered solvent in 3VO2 was squeezed out by PLATON.62
1VO2: yield: 0.107 g (97%). IR: cm−1 1625 (νC
N imine), 981 (νV
O). 1H NMR (300 MHz, d6-DMSO): 7.98 (s, 2H), 7.79 (s, 2H), 7.70 (s, 2H), 6.38 (s, 2H). HRMS (ESI-TOF, positive mode) m/z: calcd for C14H11N7O2V 360.0414, found 360.0406 [M + H]+.
2VO2: yield: 0.125 g (82%). IR: cm−1 1622 (νC
N imine), 961 (νV
O). 1H NMR (300 MHz, d6-DMSO): 13.46 (s, 2H), 8.14–8.12 (m, 4H), 7.76 (s, 2H), 7.58–7.55 (m, 4H), 7.34–7.31 (t, 2H). HRMS (ESI-TOF, positive mode) m/z: calcd for C22H15N7O2V 460.0727, found 460.0721 [M + H]+.
3VO2: yield: 0.052 g (37%). IR: cm−1 1645 (νC
N imine), (νV
O), 973 (νV
O). 1H NMR (300 MHz, d6-DMSO): 8.73 (s, 2H), 8.04 (s, 2H), 7.94–7.69 (t, 2H), 7.77 (s, 2H), 7.47–7.44 (d, 2H), 7.29 (m, 2H). Anal. Calcd for C22H19N7O4.5V (3VO2 + 2.5 H2O; 504.39): C, 52.39; H, 3.80; N, 19.44. Found: C, 52.38; H, 3.76; N, 19.82.
4VO(SO4): yield: 0.084 g (65%). IR: 1645 (νC
N imine), 984 (νV
O), 1062, 1034, 607 (νS–O). 1H NMR (300 MHz, CDCl3): 8.63–8.62 (d, 2H), 8.13 (s, 2H), 7.80–7.78 (t, 2H), 7.67 (s, 2H), 7.52 (m, 2H), 7.17–7.12 (t, 2H). Anal. Calcd for C18H12N5O5SV (4VO(SO4); 461.33): C, 46.86; H, 2.62; N, 15.18. Found: C, 47.08; H, 2.79; N, 15.27. HRMS (ESI-TOF, positive mode) m/z: calcd for C19H14N5O3V 411.05364, found 411.0528 [M-SO4 + formic acid (FA)]+.
4VO(acac): yield: 0.071 g (89%). IR: 1635 (νC
N imine), 1575, 1548 (νC
O AcAc), 953 (νV
O). 1H NMR (300 MHz, CDCl3): 8.63 (s, 2H), 8.09 (s, 2H), 7.79–7.77 (t, 2H), 7.67 (s, 2H), 7.47 (d, 2H), 7.12 (t, 2H), 5.50 (s, H on acac), 2.04 (2, H on acac –CH3). Anal. Calcd for C23H20N5O3.5V (4VO(acac) + 0.5 H2O; 473.38): C, 58.36; H, 4.26; N, 14.79. Found: C, 58.74; H, 4.18; N, 14.51. HRMS (ESI-TOF, positive mode) m/z: calcd for C19H14N5O3V 411.05364, found 411.0528 [M-acac + formic acid (FA)]+.
CCDC 2516084–2516089 contain the supplementary crystallographic data for this paper.65a–f
| This journal is © The Royal Society of Chemistry and the Centre National de la Recherche Scientifique 2026 |