Electronic and geometric structure e ﬀ ects on one-electron oxidation of ﬁ rst-row transition metals in the same ligand framework †

Developing new transition metal catalysts requires understanding of how both metal and ligand properties determine reactivity. Since metal complexes bearing ligands of the Py5 family (2,6-bis-[(2-pyridyl)methyl] pyridine) have been employed in many ﬁ elds in the past 20 years, we set out here to understand their redox properties by studying a series of base metal ions (M = Mn, Fe, Co, and Ni) within the Py5OH (pyri-dine-2,6-diylbis[di-(pyridin-2-yl)methanol]) variant. Both reduced (M II ) and the one-electron oxidized (M III ) species were carefully characterized using a combination of X-ray crystallography, X-ray absorption spectroscopy, cyclic voltammetry, and density-functional theory calculations. The observed metal – ligand interactions and electrochemical properties do not always follow consistent trends along the periodic table. We demonstrate that this observation cannot be explained by only considering orbital and geometric relaxation, and that spin multiplicity changes needed to be included into the DFT calculations to reproduce and understand these trends. In addition, exchange reactions of the sixth ligand coordinated to the metal, were analysed. Finally, by including published data of the extensively characterised Py5OMe (pyridine-2,6-diylbis[di-(pyridin-2-yl)methoxymethane])complexes, the special characteristics of the less common Py5OH ligand were extracted. This comparison highlights the non-innocent e ﬀ ect of the distal OH functionalization on the geometry, and consequently on the electronic structure of the metal complexes. Together, this gives a complete analysis of metal and ligand degrees of freedom for these base metal complexes, while also providing general insights into how to control electrochemical processes of transition metal complexes.


Introduction
The optimization of multi-step redox catalysis performed by molecular complexes requires a fundamental understanding of the factors that determine performance. Redox potentials of individual one-electron processes, the coupling of electron and proton transfer and the binding affinities of reactants and products are often intertwined properties that only a systematic experimental and theoretical study can disentangle. 1 As the properties of coordination complexes are determined by the metalligand interactions, such analysis should describe their influence in a systematic manner. 2 Presently available studies where the metal site is varied within the same coordination sphere mainly focus on a phenomenological description without developing deep insight into the reasons for different reactivities. 3,4 Here we begin to fill this gap by a detailed study of one-electron redox reactions in a series of first-row transition metal complexes with the same pentapyridyl ligand framework (Fig. 1). This fundamental approach makes it possible to focus initially on the role of the metal in a well-defined reaction step, followed by a separate analysis of the ligand dimension.
The family of pentapyridyl ligands stemming from Py5, (2,6-bis-[(2-pyridyl)methyl]pyridine) has been extensively adopted in the past two decades to complex different first-row transition metals. 5 Feringa et al. reported the first example of Fe II and Mn II complexes with the methoxy-substituted ligand (Py5OMe), highlighting how the apical coordination site is particularly prone to ligand substitution and therefore can reversibly bind possible substrate candidates for catalytic reactions. 6 Simultaneously, Stack et al. studied catalytic applications of the same complexes for the oxidation of organic substrates as they mimic the cofactor of lipoxygenases. [7][8][9] This team also reported structural parameters, magnetic susceptibility, electrochemical behaviour, and optical properties of the series of [M(Py5OMe)Cl] + (M = Mn, Fe, Co, Ni, Cu and Zn) complexes, but only in their reduced M II oxidation state. 4 During the following two decades, Py5 complexes have also been applied as anti-tumour agents, 10,11 redox mediators for dye-sensitized solar cells, 12,13 and materials with tunable magnetic properties. [14][15][16] In the field of the artificial photosynthesis, Py5-type ligands have been reported to be a suitable scaffold for synthesizing both proton reduction [17][18][19][20] and water oxidation catalysts. [21][22][23][24] With the idea of creating a starting point for further functionalization, our group introduced the Py5OH ligand and employed it in homogeneous water oxidation catalysis using either Co II or Fe II as the metal centre. 25,26 Our recent study on the [Fe II (Py5OH)Cl] + complex showed that, unlike [Fe II (Py5OMe)Cl] + , it undergoes a spin-crossover from a high spin (HS) to low spin (LS) configuration, even in the presence of a weak-field ligand like Cl. 15 This infers that the peripheral R-groups of the ligand can affect the electronic structure of the metal centre.
Using the Py5OH ligand framework ( Fig. 1; R = OH) we synthesized here a series of metal complexes with the general formula [M II (Py5OH)Cl](PF 6 ) (abbreviated as [M II -Cl]), where M = Mn, Fe, Co and Ni. Their one-electron redox potentials were then determined by cyclic voltammetry, and the geometric and electronic structures of both the reduced [M II -Cl] and oxidized [M III -Cl] species were investigated in powder and/or dissolved form by employing a combination of singlecrystal X-ray diffraction (XRD), synchrotron X-ray absorption spectroscopy (XAS) and density-functional theory (DFT) calculations.
The effect of the ligand sphere on the redox potential was then studied by exchange of the apical chloride ligand. This is especially interesting with regard to catalytic reactions as its exchangeability is likely important for substrate activation. To facilitate this, we also prepared the Cl-free [Fe II (Py5OH)MeOH] (ClO 4 ) 2 complex, which upon dissolving in dimethylformamide or acetonitrile resulted in two additional set of complexes: [Fe II -DMF] and [Fe II -MeCN]. Finally, the electrochemical properties of the [M II -Cl] complexes were compared to those reported previously for [M II (Py5OMe)Cl] + , abbreviated [M II -Cl] OMe . Together, this work represents a consistent set of modifications for the pentapyridyl complexes in Fig. 1, varying metal (M), apical ligand (X) and pentapyridyl substituent (R). Combined with the careful experimental and theoretical analysis of a well-defined one-electron redox event, this provides fundamental insight into the electrochemical processes of transition metal complexes. and [Co II -Cl] from the present study showed good agreement with those published previously. 25,26 In all the complexes, the transition metal was embedded in a pocket composed of the five pyridines with a sixth-axial chloride ligand. The metal environment had approximate local C 4v symmetry with the Fe-Cl vector along the C 4 -proper axis. Fig. 3 shows that all the metal-nitrogen bond distances (average of the equatorial M-N eq. , blue; axial M-N axial , red) decreased following the order of the periodic table, for M-N axial from 2.249 Å in [Mn II -Cl] to 2.065 Å in [Ni II -Cl]. This was not the case for the metal-chlorine bond distance: although the Mn-Cl bond (2.458 Å) was longer than the Fe-Cl bond (2.419 Å), it increased again for the subsequent Co    15 To get closer to the conditions of the CV experiments, the XAS parameters were also obtained after dissolving the four samples in acetonitrile with electrolyte. The XAS spectra of the [Mn II -Cl], [Co II -Cl] and [Ni II -Cl] complexes remained unchanged as compared to the spectra recorded for the powder samples (see Fig. S3-S6 †). By contrast, large differences in the XAS spectra were observed for [Fe II -Cl] under the two conditions. A reasonable XAS fit of the dissolved Fe-complex at 20 K could be obtained with three components: (i) 40% of the chloride-bound LS form that was also observed for the powder, (ii) 40% of a HS form with long metal-ligand bonds, and (iii) 20% of a LS structure where the chloride ligand had exchanged with a solvent molecule. This was supported by a good EXAFS fit at 150 K where only two components, 80% of (ii) and 20% (iii), were necessary (see Fig. S8 and S9 †). Dissolving the sample in an electrolyte solution thus led to partial ligand exchange and incomplete SCO down to 20 K.
The electronic and geometric structures of the [M II -Cl] complexes were calculated for all spin multiplicities consistent with 3d n configurations (n = 5-8). Correct spin-state energetics are in general challenging to calculate with DFT due to strong functional dependence. [27][28][29][30] The results are especially sensitive to the amount of HF exchange, but as the functional dependence varies with the type of bonding it is difficult to get accurate results for a wide range of complexes. 29 To address this challenge, we used the SCO in [Fe II -Cl] as a reference point for selecting the functional. 15 For [Fe II -Cl] the B3LYP* functional (15% HF exchange) gave good results. It favoured the quintet over the singlet by 5.3 kcal mol −1 at room temperature, see Fig. 4. This value was 2.2 kcal mol −1 higher than what can be expected from the determined spin-transition temperature of the powder sample at 80 K. 15 As the solvated chloride-coordinated sample had remaining HS species down to 20 K, the B3LYP* calculations with 15% HF exchange should thus represent the complex in solution with relatively high accuracy. Calculations of the other complexes showed that they all favour high-spin states, even more so than iron, see Fig. 4 (blue bars) and Table S2. † At RT, the spin-state energetics for [Co II -Cl] was very close to that of [Fe II -Cl], but compared to [Fe II -Cl] the LS form of the Co-complex was less well-stabilized at low temperatures, so that the calculations did not predict this to be a spin-crossover complex, see Table S5. † The corresponding DFT geometric structures were, in general, in good agreement with XRD data, see Fig. 3 and Table 1. The Fe-N eq. distances had absolute deviations of less than 0.02 Å, which meant that trends in bond distances were Fig. 3 Comparison of distances from single-crystal X-ray diffraction collected at 150 K (full triangle, solid line) and calculated distances using DFT; (half-filled triangle, dotted line). M-Eq. plane, metal displacement from the equatorial plane, which is defined by the four nitrogen atoms: N2, N2i, N3 and N3i (see Fig. 2). also well reproduced. The agreement is poorer for the axial bonds, with too short M-Cl bonds (−0.1 Å on average) and much too long for M-N axial bonds (+0.2 Å on average). This led to calculated metal centre positions that were further displaced from the equatorial plane than measured by XRD (by 0.1 Å). However, the structural trends were still very well reproduced, with all M-N distances (M-N eq. and M-N axial ) decreasing following the order in the periodic table, while the M-Cl distances showed the above-described deviation for Co and Ni.
The choice of functional and basis set size did not affect the calculated distances by more than 0.02 Å, see Table S3, † and the very long Fe-N axial bond is present also with the local BP86 functional. Despite the deviations of the axial bond distances, the HS structures still give the best agreement with experiment. Structures with lower spin multiplicities all have significantly shorter Fe-N eq. bonds, which led to an underestimation by 0.2 Å, with no significant improvement in the Fe-N axial distance, see Table S2. †

Electrochemistry
In cyclic voltammetry (CV) experiments, sweeping from the equilibrium potential towards the positive direction, reversible oxidation processes were observed for all [M II -Cl] complexes (Fig. 5). This was tentatively assigned to a metal-centred oneelectron oxidation from [M II -Cl] → [M III -Cl]. The experimentally determined order of the reduction potentials of the redox couple [M III -Cl]/[M II -Cl] was: Co (0.08 V) < Fe (0.33 V) < Mn (0.58 V) < Ni (1.17 V), which deviates from the trend of increasing potentials along the periodic table that may be expected on the basis of increasing nuclear charge density. The vertical lines in Fig. 5 indicate the calculated redox potentials, which will be discussed in detail after describing the characterization of the oxidized complexes. In the case of [Fe II -Cl], a second reversible redox wave at E 1/2 = 0.83 V was observed, which had an overlapping shoulder feature at E p = 0.72 V in the oxidative scan (marked with * in Fig. 5). This shoulder feature was assigned to the quasi-reversible response of the redox-active couple Cl 2 /Cl − by the observed increase of this feature after addition of tetrabutylammonium chloride (TBACl) to [Fe II -Cl] (Fig. S10 †) and by a separate voltammogram of TBACl (Fig. S11 †). The second reversible redox wave was likely not an additional one-electron oxidation of Fe III to Fe IV , since there was a significant difference in current intensity between the two waves. Instead, the higher redox potential feature likely arose from a secondary species that, based on the cathodic peak analysis, made up about 25% of the total compound. Table 1 Metal-ligand distances (Å) derived from single-crystal X-ray diffraction (XRD), synchrotron X-ray absorption (EXAFS), and DFT calculations. EXAFS data were obtained at 20 K from complexes dissolved in acetonitrile with electrolyte, while XRD data were collected at 150 K.  Given the observation of the irreversible Cl 2 /Cl − redox wave, the most likely candidate was the complex where the Cl − ligand was replaced by a solvent molecule, which was also observed in the EXAFS data. This possibility will be analysed in detail in a later section.
We note that due to the difference in temperature between the experiments, the LS 1 [Fe II -Cl] species observed with EXAFS at 20 K was not present in the room temperature CV experiment. Thus, we estimate that the room temperature composition was 75% HS 5 [Fe II -Cl] and 25% 1 [Fe II -MeCN], in good agreement with the CV peak analysis.

Structures of oxidized complexes
To understand the redox potentials, also the products of the one-electron oxidation need to be analysed. Thus, the metal complexes were electrochemically oxidized in acetonitrile and analysed with XAS. The K-edge positions of the oxidized materials displayed a shift to higher energy compatible with a complete one-electron metal oxidation, see  Table 1). These structural changes observed by EXAFS spectroscopy are compared below with distances obtained by DFT calculations.
Oxidation of the [Mn II -Cl] complex shortened the average Mn-N distance by about 0.11 Å and the Mn-Cl distance by 0.21 Å (Table 1). A direct EXAFS fit of the oxidized form requires a large Debye-Waller factor when a single Mn-N shell is used (0.085 Å, compared to around 0.065 Å for the other complexes, see Table S4 †), indicating a large Jahn-Teller distortion. DFT calculations of the oxidized Mn complex give the lowest energy for a quintet 5 [Mn III -Cl] structure ( Fig. 4), with the triplet 3 [Mn III -Cl] 6.7 kcal mol −1 higher in energy (Table S2 † Table S2 †) as compared to the quintet, which gives better agreement with EXAFS distances (Fig. 6). However, the fit to the EXAFS spectrum was not significantly improved ( Fig. S12 †). Instead, an alternative quintet structure 5 [Mn III -Cl]′ with a shorter Mn-Cl bond of 2.26 Å (Table 1 and Fig. 6) and elongated bonds in the equatorial plane, lies only 0.6 kcal mol −1 above the most stable quintet and gives a much better fit to EXAFS (Fig. S12 †). Considering the small energy difference and the good structural match, this alternative quintet structure will also be considered in further comparisons and analyses.
For [Co III -Cl], EXAFS fitting suggests that upon oxidation the average Co-N distance shortens from 2.14 Å to 1.97 Å, and the average Co-Cl distance from 2.40 Å to 2.20 Å. The majority of reported Co III complexes exist in a LS singlet configuration, 31 and the distances found here for [Co III -Cl] were also consistent with such an assignment. DFT calculations also prefer 1 [Co III -Cl] over [Co III -Cl] is also much higher in energy, see Table S2. † The calculations reproduce the short metal-ligand bonds of the EXAFS structure, although they are overestimated by approx. 0.04 Å.
In contrast to the Mn and Co complexes, the oxidation of [Fe II -Cl] to [Fe III -Cl] did not lead to a similar shortening of the bonds and connected EXAFS changes. While the EXAFS spectrum of [Fe II -Cl] required three components to be fitted at 20 K, the EXAFS data obtained with the oxidized complex, [Fe III -Cl], could be simulated well with a single species that had a short Fe-Cl distance (2.22 Å) and short average Fe-N distance (1.99 Å), see Table 1. This suggests the oxidation to a LS 2 [Fe III -Cl] structure. Unlike the reduced complexes, there was thus no equilibrium between different spin multiplicities or axial ligands (Cl − vs. solvent) after oxidation. Energetically, DFT calculations favour the low-spin over the high-spin 6 [Fe III -Cl] structure by 1.1 kcal mol −1 already at 298 K, and this tendency increases with decreasing temperature (Table S5 †). The calculated average Fe-N bonds (2.04 Å) in the 2 [Fe III -Cl] structure are 0.05 Å longer than predicted by experiment (Fig. 6), but still a much better fit than the 2.17 Å of the average Fe-N bond distance obtained for the 6

Simulated redox potentials
Redox potentials of the [M III -Cl]/[M II -Cl] couples were simulated using the optimized DFT structures of reduced and oxidized complexes. Calculations of absolute redox potentials of transition metal complexes often show significant errors, up to 0.5 eV. [32][33][34][35][36][37] Highly charged complexes are in general challenging due to large absolute solvation effects. 38 Techniques to reduce errors include explicit modelling of the solvent, 39 and multi-step calculations to reach a neutral reference complex. 40 Here, the main target was to obtain relative potentials between structurally similar complexes, rather than absolute values, and thus a standard approach was selected. The results of our DFT calculations are shown as vertical lines in the CVs plotted in Fig. 5, and are also reported in Table S6. † In agreement with our experiments, the lowest redox potential was found for [Co II -Cl]. The theoretical value of 0.18 V, is +0.10 V higher than the experimental value (0.08 V). This deviation is well within the expected accuracy of the calculations, and is likely due to error cancellation. For [Fe II -Cl], the calculated value was 0.37 V, which is in even better agreement with the experiment (0.33 V). Experimentally, [Mn II -Cl] had a higher redox potential (0.58 V) than both [Fe II -Cl] and [Co II -Cl]. Calculations with the lowest energy structure, 5 [Mn III -Cl], gave a value of 0.34 V (−0.24 V vs. exp.), which placed it at the same level as that of the iron complex. Using the alternative 5 [Mn III -Cl]′ structure would instead give a redox potential of 0.37 V, which is a minor improvement compared to the experiment. Finally, the calculated redox potential of [Ni II -Cl] was 1.19 V, much higher than the other complexes, and in excellent agreement with the experiment (1.17 V).

Exchange reactions at the sixth ligand position
For [Fe II -Cl], the changes in EXAFS data when dissolved in acetonitrile, together with the presence of a second oxidation event in the CV, suggested that the Cl − ligand can be exchanged. Therefore, chloride-free iron complexes, where the sixth ligand position was occupied by a solvent molecule, were prepared (Scheme 1). Crystals suitable for XRD diffraction were obtained from a DMF solution, and the crystal structure of [Fe II -DMF] (Fig. S2 †) showed that the oxygen of the amide group in DMF coordinated to the iron. For this material, the XRD data showed that the average Fe-N bond length was, at 150 K, shorter by approximately 0.2 Å than in the [Fe II -Cl] complex discussed above, indicating singlet configuration (Table 1). 41 DFT calculations of the similar [Fe II -MeCN] complex gave very good agreement with the experimental (XRD) distances, with a maximum deviation of 0.02 Å.
Preparing the chloride-free complex in methanol yielded a solid material where MeOH is in the apical coordination site. When this product was dissolved in acetonitrile, methanol was replaced by an acetonitrile molecule ([Fe II -MeCN]), as observed in a HR-MS conducted at low ionization energy (Fig. S13 †). Taken together, these data suggest that the solvent molecule at the apical coordination site was weakly bound and could be easily exchanged. The EXAFS data of [Fe II -MeCN] shown in Fig. 7 could be fitted using a short metal-nitrogen bond length of 1.98 Å (Table 1), consistent with a LS Fe II centre. The LS configuration was further supported by a 1 H-NMR spectrum of [Fe II -MeCN] that showed the presence of a diamagnetic metal complex (Fig. S14 †). DFT calculations also favoured the LS form (by 3.2 kcal mol −1 ) at 298 K and gave metal-ligand distances within 0.02 Å of the EXAFS analysis.
In contrast to the cyclic voltammogram of [Fe II -Cl] in acetonitrile that showed two reversible redox waves (E 1/2 = 0.33 V and 0.83 V), the CV of [Fe II -MeCN] showed only one reversible redox wave with a potential that matched the position of the second wave of the CV labelled [Fe II -Cl] (Fig. 5 and S15 †). interchanged easily in solution. An electrochemical contribution to this ligand exchange can be excluded as we observed the same shift in the equilibrium by monitoring the Cl-titration by UV-Vis spectroscopy in pure acetonitrile (Fig. S16 †). The UV-Vis data obtained in presence of TBAPF 6 revealed that the supporting electrolyte partially shifted this equilibrium in favour of the [Fe II -MeCN] species (Fig. S17 †). This effect was probably due to the increased ionic strength of the medium, which stabilized free chloride in the solvent. This explains, why even after addition of two equivalents of Cl − a mixture of [Fe II -Cl] and [Fe II -MeCN] was still observed ( Fig. 8 and S16a †). In addition, the Cl 2 /Cl − redox couple was also observed at E p = 0.72 V in the oxidative scan in Fig. 8.
In contrast to the [M II -Cl] complexes, the oxidation of the [Fe II -MeCN] complex occurred without significant structural changes ( Fig. 7 and Table 1), although the oxidation was witnessed by the edge shift in the XANES region (Fig. S7 †). The near-constant metal-ligand distances (Table 1)

Orbital energetics
The above experiments and DFT analysis established that the M II to M III oxidations can be assigned to clean one-electron oxidations with the expected metal-ligand bond shortenings for the Mn-and Co complexes, while a more complex situation was identified for the Fe-complex, where the apical position was exchangeable for [Fe II -Cl] and different spin states can be stabilized depending on temperature, while only one structure was observed after oxidation to Fe III . To rationalize these results and the observed spin multiplicities within simple models, we discuss the electronic structure of our complexes using a molecular orbital framework.
The introducing an apical ligand (X) to a hypothetical complex with six equivalent ligands (O h point group) decreases Thus. a C 4v symmetry is a fair assignment of the local metal environment. In reality, the symmetry deviates from C 4v to a certain extent since there is no authentic C 4 rotational axis giving the Py5OH ligand. This deviation, however, does not change the principles in the following section.
The calculated electronic structure of the Py5 complexes are consistent with a previous analysis that described the pyridine rings as being predominantly σ-donor and π-acceptor ligands, 5 while Cl − as σ and π donor. In a C 4v symmetry, the orbital degeneracy from the well-known cubic environment is lifted. Consequently, the metal 3d-dominated t 2g and e g orbitals in a virtual O h symmetry transform to b 2 + e and a 1 + b 1 respectively, see Fig. 9. The lowest b 2 level is related to a d xy orbital (equatorial pyridines π-bonding), above which is two-fold degenerated e level related to d xz and d yz orbitals (pyridine π-bonding and Cl − π-antibonding). These three orbitals will be labelled as π-type orbitals. Higher in the energy are a 1 and b 1 levels, both σ-type orbitals related to d z 2 (Cl − /pyridine σ-antibonding) and d x 2 −y 2 (pyridine σ-antibonding), respectively. The ordering of these two levels will be discussed in more detail below.
Going from left to right in the periodic table lowers the energy level of the metal orbitals because the additional nuclear charge is only partially screened by the extra electron. This lowering is even more pronounced when increasing the formal oxidation state. For ligand-donor bonding, lowering the metal level leads to smaller energy differences between metal and ligand orbitals, and thus a stronger interaction. For ligand-acceptor binding, it instead leads to a larger energy difference between metal and empty ligand energy levels and weaker interactions. In general, the energy difference between the π-levels and the antibonding σ-levels increases, which corresponds to a larger ligand-field splitting. These general principles can now be used to rationalize trends in geometric structure, spin-state energetics, and electrochemical behaviour.  Table S4. † The shoulder at 0.72 V marked with (*) was assigned to the oxidation of Cl − . The reductive feature at −0.4 V marked with ( ‡) was tentatively assigned to a further reaction between Cl 2 and the Fe complex.

Trends in geometric and electronic structure for different metals
Starting with the geometric structures of the reduced [M II -Cl] complexes, there is a clear trend of decreasing metal-ligand bond distances from Mn to Ni, see Fig. 3. This is consistent with a decrease in ionic radius with higher nuclear charge Z. However, the change in ionic radius from Mn to Ni is only 0.03 Å, while the average M-N distances decrease by 0.15 Å. At the same time, the M-Cl distance decreases when going from [Mn II -Cl] to [Fe II -Cl], but then increases slightly for later elements. As all reduced complexes are predicted to be HS, the number of electrons in the π-type orbitals increases from three to six along with the periodic table, while the number of electrons in σ-type orbitals is constant at two. Going from 6 [Mn II -Cl], which has one unpaired electron in each 3d orbital, to 5 [Fe II -Cl] adds a spin-down electron in the lowest orbital d xy .
This orbital is weakly bonding with respect to the equatorial Fe-N and the effect of an extra electron is a strengthening of the Fe-N bonds and a slight decrease in metal-ligand bond lengths. At the same time, going from Mn to Fe lowers the energy of the metal 3d orbitals, which leads to stronger σbonding and a reduction in all bond distances. As both effects work in the same direction, all bond distances contract.
The results are more interesting when going from Fe to Co. In 4 [Co II -Cl] the extra electron goes into the d xz /d yz level, which is bonding with respect to the pyridines and antibonding with respect to Cl − . The added electron leads to shorter M-N bonds and longer M-Cl bonds. At the same time, increased donor bonding should have an overall contracting effect. Still, for the metal-chlorine bond, the anti-bonding effect dominates, and the bond length increases. Finally, going from 4 [Co II -Cl] to 3 [Ni II -Cl] puts another electron into the d xz /d yz orbitals, which leads to similar changes as seen when going from Fe to Co. The molecular orbital analysis can thus fully explain the bond length trends for the reduced complexes.
The metal-ligand distances of the oxidized complexes show a complex behaviour upon exchange of the metal (Fig. 6). The computational analysis of the [Mn III -Cl] compound is complicated because of the deviations between the experiment and the lowest-energy 5 [Mn III -Cl] (π 3 σ 1 ) structure. The short Mn-Cl bond from EXAFS (2.22 Å) suggests that an electron is removed from the anti-bonding σ-type d z 2 orbital upon oxidation, but the calculations instead give for 5 [Mn III -Cl] a long Mn-Cl and shorter Mn-N eq. bonds, see Table S2, † consistent with an electron taken from the d x 2 −y 2 orbital. The alternative 5 [Mn III -Cl]′ structure, which with shorter axial and longer Mn-N eq. bonds, gives a better EXAFS match (Fig. 6), is more consistent with an empty d z 2 orbital. It is thus possible that the calculations remove the electron from the wrong σ orbitals in this oxidation process.
As the octahedral HS Mn III complex shows a Jahn-Teller distortion, the difference between the two calculated structures can also be viewed as a change in distortion axis.
As the calculations consistently overestimate the axial Mn-N bond distances, it is possible that they artificially favour the 5 [Mn III -Cl] structure with a Jahn-Teller-like elongation in this direction. Overestimations of the Mn III JT distortion have also been previously observed in DFT calculations. 42 For 2 [Fe III -Cl] theory and experiment largely agree, see Fig. 6. The 2 [Fe III -Cl] complex has a LS (π 5 σ 0 ) electron configuration and compared to 5 [Mn III -Cl] this means extra electrons in d xy and one of the d xz , d yz orbitals while removing the final σ electron. This leads to short axial as well as equatorial bonds. Going further from 2 [Fe III -Cl] to LS 1 [Co III -Cl] (π 6 σ 0 ) adds an electron to d xz /d yz , and as these are antibonding with respect to Cl − , together with the general contraction with larger Z there should be small changes in the M-Cl distance while the M-N distances decrease when going from Fe to Co (Fig. 6). Incidentally, the extra electron is here added to the same orbital level as for the comparison of the reduced iron and cobalt complexes. Finally, the trend when going from 1 [Co III -Cl] to 2 [Ni III -Cl] (π 6 σ 1 ) can only be derived from calculations as the oxidized nickel complex could not be obtained. They show a significant increase in M-Cl bond length and a slight increase in the average M-N bond length. The increase in bond distances is expected as the added electron is placed in an σ antibonding orbital. As the bond distances increase more along the C 4 -axis, the calculations suggest that the electron ends up in the d z 2 rather than the d x 2 −y 2 orbital. LS Ni III complex is expected to show a similar Jahn-Teller-like distortion as HS Mn III , and the long M-Cl distance is similar to that of the 5 [Mn III -Cl] structure. However, as seen from the [Mn III -Cl] calculations, it is possible that the structure instead would have a distortion along the equatorial M-N plane.
The spin-state energetics are governed by two major factors, the ligand-field strengths and the differences in exchange stabilization between electron configurations. If the former is sufficiently large, the benefit of being at the lower orbital levels competes with and overcomes the loss of exchange stabilization and the electrons fill lower b 2 and e orbitals, resulting in LS configurations. The exception is Ni II where both spin configurations have two electrons in the upper two levels. In general, both ligand-field strengths and the exchange interactions increase along the periodic table. However, with changing d-electron count, the difference in the number of ligandfield excitations and exchange interactions change between LS and HS states. The resulting spin-state energetics in Fig. 4, therefore, do not follow any linear trend but are rather consistent with the spectrochemical series of metals. Taking the reduced complexes as an example, [Mn II -Cl], it has a relatively weak ligand field and the sextet has all electrons with the same spin orientation, leading to a highly stable configuration as seen in Fig. 4. Going to d 6 Fe II and d 7 Co III , all electrons can no longer have the same spin orientation, and the ligand-field splitting increases, resulting in the lower relative stability of the HS configurations. As mentioned above, for d 8 [Ni II -Cl] b 2 and e orbitals (b 2 2 e 4 ) are filled in all configurations, and the high spin is favoured because it allows for a 1 and b 1 to be half-filled with electrons in the same spinorientation.

Trends in redox potential for different metals
An important target of this study has been to understand the factors that determine the redox potentials of different base metals in the Py5OH complexes. Starting from a basic atomic picture, the redox potential increases when going from left to right in the periodic table as the energy of the 3d levels decrease with increasing nuclear charge. For the Py5OH complexes, the electrochemical experiments instead gave the trend Co < Fe < Mn < Ni. To analyse the effect of the ligand environment, we consider electronic, geometric, and spin degrees of freedom. From Fig. 5 it is clear that the redox potentials of all complexes except [Mn II -Cl] are well described by the calculations, not only in terms of relative but also absolute potentials. The main discrepancy is that manganese is too easy to oxidize according to the calculations. However, the deviation, −0.24 eV, is within the expected accuracy of the method, and the calculations reproduce the overall experimental trend well enough to be used to rationalize the trends in redox potential.
The oxidation process will be analysed in three schematic steps. The first step is the removal of an electron from the highest occupied molecular orbital, which in all cases is the σtype d x 2 −y 2 orbital, and the subsequent orbital relaxation. This is followed by geometry relaxation to the minimum energy geometry of the oxidized species. Finally, in cases where a more stable spin multiplicity existed, spin change was allowed.
The description of [Mn II -Cl] oxidation is straightforward: removing an electron is equivalent to going from 6 [Mn II -Cl] (π 3 σ 2 ) to 5 [Mn III -Cl] (π 3 σ 1 ). Including only orbital relaxation would lead to a hypothetical redox potential of 1.11 V, see Fig. 10. This is followed by geometric relaxation of the oxidized complex to its minimum energy, which leads to a drop in the potential to 0.34 V. As the quintet is also the lowest energy of [Mn III -Cl], this is also the final calculated redox potential. Using instead the alternative 5 [Mn III -Cl]′ structure, which would be consistent with the removal of a d z 2 electron, would give a very similar analysis, although with a slightly higher final potential (0.37 V). The relative order of d z 2 and d x 2 −y 2 might thus not be correctly predicted, but this does not seem to have a significant effect on the analysis.
For [Fe II -Cl], the oxidation is a 5 [Fe II -Cl] (π 4 σ 2 ) to 2 [Fe III -Cl] (π 5 σ 0 ) transition. Taking an electron from a σ-type orbital leads initially to a 4 [Fe III -Cl] π 4 σ 1 configuration at 1.20 V, which then drops to 0.73 V upon geometry relaxation, see Fig. 10. However, for iron, there is an additional change in spin multiplicity from a quartet to a doublet, which gives the final calculated potential of 0.37 V. The analysis for the 4 [Co II -Cl] (π 5 σ 2 ) to 1 [Co III -Cl] (π 6 σ 0 ) transition is similar to the case for iron. The oxidation to 3 [Co III -Cl] (π 5 σ 1 ) gives potentials of 1.48 and 0.98 V after orbital and geometric relaxation, respectively. For [Co III -Cl], the singlet is much more stable than the triplet, which leads to a drop in the potential down to 0.19 V. Finally, for [Ni II -Cl] the electron is also taken from a σ-type orbital as it goes from 3 [Ni II -Cl] (π 6 σ 2 ) to 2 [Ni III -Cl] (π 6 σ 1 ). As expected, the removal of a σ electron is more difficult in [Ni II -Cl] than in the other complexes, leading to a potential of 1.19 V after geometric relaxation. As there is no further change in spin multiplicity, this is also the final redox potential, which is then much higher than for any other complex in the series. From the trends in Fig. 10 it is thus clear that if only orbital and geometry relaxation would occur, the potentials of all our Py5OH complexes would follow the order in the periodic table, i.e. show a higher redox potentials with higher Z. It is thus the difference in energy stabilization from potential changes in spin multiplicity that determines the final order of the redox potentials. The calculations can thus be used to rationalize the changes in electrochemical properties between different base metals.

Exchange reactions at the sixth ligand position
Exchange reactions at the sixth ligand position are likely required for the activation of substrates in catalytic reactions involving the [M II -Cl] complexes. From the experimental data, it is clear that [Fe II -Cl] can release Cl − in the electrolyte acetonitrile solution. To get insight into the process, the reaction energy for chloride substitution by the solvent acetonitrile was calculated for the four metal complexes using DFT. The calculations show a trend in Cl − binding energies consistent with the distance analysis (Fig. 11), where binding in [Fe II -Cl] is stronger than in [Mn II -Cl] due to the better overlap between metal and ligand donor orbitals. The Fe-Cl bond is also stronger than in [Co II -Cl], because the latter has additional electrons in the anti-bonding d xz /d yz orbitals as described above. Finally, the binding energy to [Ni II -Cl] is the highest in the series, likely due to further metal-ligand mixing with higher Z.
The ligand exchange equilibrium is instead determined by the energy costs for solvating the Cl − ion and the binding energy of the acetonitrile ligand, which binds stronger to Fe and Ni compared to Mn and Co. Compared to the chloride complex, acetonitrile favours the singlet state by 8.5 kcal mol −1 , which makes [Fe II -MeCN] a low-spin complex.
However, the rest of the [M II -MeCN] complexes remain HS, even though the energy difference is only 1.8 kcal mol −1 for [Co II -MeCN]. That means all complexes except iron have electrons in the anti-bonding d z 2 orbitals. The reason for stronger bonding in Ni compared to Co could be a combination of improved overlap and the fact that the extra electron now appears in a d xz /d yz orbital that is bonding between metal and MeCN.
Taken together, the calculations show that the exchange energy, which is the difference between the two binding energies, is close to zero for [Fe II -Cl] and [Ni II -Cl], while it is slightly positive for [Mn II -Cl] and [Co II -Cl], see Fig. 11. For [Mn II -Cl], [Fe II -Cl] and [Co II -Cl] this reproduces the trend observed in the experiments, as we observed the exchange of the apical ligand only for the [Fe II -Cl] complex. By contrast, for [Ni II -Cl] the thermodynamically expected ligand exchange was not observed, which may indicate a kinetic limitation of this process. Interestingly, for the related [Ni II (Py5Me)Cl] + complex the chloride was recently reported to exchange with the solvent. 23 Calculations of [Fe III -Cl] show that in the oxidized state, the corresponding exchange reaction for a neutral solvent molecule is unfavourable, which is consistent with the EXAFS observations of [Fe III -Cl] described above. The difference in calculated redox potential between [Fe II -Cl] and [Fe II -MeCN] can also be attributed to the difference in charge, with the anionic ligand favouring lower redox potentials.

Effects of Py5 modifications
The properties of the complexes can be further modified by synthetic changes to the Py5 ligand. Here the focus is on a comparison of the hydroxyl-substituted Py5OH and previously published results of the methoxy-substituted Py5OMe ligands. The steric hindrance of the methoxy groups in Py5OMe-ligated complexes induces a tilting of the axial pyridine ring relative to the equatorial plane and a nonsymmetrical pyridine coordination in the equatorial plane. 15 In the case of Py5OH, the perpendicular pyridine is aligned with the axial axis and a less distorted octahedral geometry is observed. We previously calculated that the distortion in Py5OMe leads to a 2.9 kcal mol −1 destabilization of the LS state, 15  Electrochemically, the Py5OMe complexes have redox potentials that are 0.08 V higher than the Py5OH complexes, see Table S6. † The exception is Ni, where the reported redox potential of the Py5OMe complex is 0.78 V lower. 4 Based on the good agreement between calculated and measured half potentials for the [Ni III -Cl]/[Ni II -Cl] redox couple in this study, we assume that the reported values for [Ni II (Py5OMe)Cl] + are not reflecting the same redox couple.
To analyse the differences between the two versions of the ligand, the redox potential of the [Fe II (Py5OMe)Cl] + complex was calculated using the same protocol as for [Fe II -Cl]. With this ligand, the ferric complex is also LS, although by a very small margin at 298 K (0.7 kcal mol −1 ). The calculated redox potential is 0.11 V higher than for the [Fe II -Cl] complex, in good agreement with the experimental difference of 0.07 V. Breaking down the process into the different steps shows that the main contribution is the increased energy required to remove the σ electron (+0.18 V), while geometric relaxation and changes in spin multiplicity work in the opposite direc- tion ([Fe II -Cl] OMe in Fig. 10). In the calculations, the electron is taken from the antibonding d x 2 −y 2 orbital, and the structural distortion with the Py5OMe ligand lowers the d x 2 −y 2 level because of poorer overlap. This effect is also visible from the longer equatorial bond distances, by 0.08 Å, in [Fe II -Cl] OMe . If the observed differences in redox potential could be directly translated into differences in barrier heights, a 0.1 eV difference corresponds to a hundred-fold effect on reaction rates. This shows the potential of modifying Py5 type ligands for tuning the reactivity of these complexes.

Conclusions
Herein we reported the first comprehensive analysis of the effects of substituting the central base metal, the auxiliary ligand, and distal modification to the ligand framework have been on the redox potential of Py5 coordinated complexes. This analysis was based on a careful structural characterization of both reduced and oxidized species. The trends in structure, spin-state energies, redox potentials, and ligand-exchange reactivity were then rationalized in a molecular orbital picture, supported by DFT calculations. The deviations from a continuous increase in redox potential with increasing nuclear charge were explained by significant stabilization from changes in spin multiplicity for the Fe and especially Co complexes, leading to large stabilizations of the oxidized states and correspondingly lower redox potentials. Moreover, we demonstrated that the exchange dynamics of the apical ligand, which is observed for Fe in both EXAFS and CV data, can be explained by the favourable binding of the solvent molecule, rather than a weak Fe-Cl interaction. Finally, the analysis quantifies the role of Py5 ligand modifications far from the metal site, which had, in line with previously reported spin-state energetics, 4,15 also measurable effects on the electrochemical behaviour. Thereby, this study shows that even a simple one-electron oxidation reaction can entail unexpected complexity and tunability. Thus, developing efficient and stable multi-electron catalysis with base metals requires a detailed study of all steps of the catalytic cycle. We also demonstrated that DFT calculations could be employed to understand, in most cases, the trends observed. We thus expect that this systematic approach will, in the long run, facilitate deriving design principles allowing the development of efficient and stable catalysts via predictive calculations.

Synthesis
All starting reagents were obtained from commercial sources and used as received. All glassware were cleaned and dried overnight at 120°C. The synthesis of the metal complexes was conducted under dry Argon atmosphere. The synthesized samples were stored in air without observing any degradation. The complexes were fully characterized to verify the purity of the desired product by means of 1 H-NMR, FT-IR, UV-VIS, HR-MS and elemental analysis as described in the ESI. † The oxidised metal complexes were obtained by electrolysis in acetonitrile solution (1 mM) as described in the next section.
Caution: Perchlorate salts are potentially explosive and should be handled with care.
Py5OH was synthesized following a procedure introduced by us recently. 15 To a dry THF solution (40 ml) of 2,6-dibromopyridine (0.65 g, 2.72 mmol) in a 100 ml 3-neck round bottom flask, an excess of Mg (0.25 g, 102 mmol) was added. After sonicating the mixture for 20 minutes at 35°C using an ultrasonic bath (45 kHz frequency, model USC300TH, VWR Collection) most of the starting solid Mg was dissolved and the transparent solution turned dark. A THF solution (20 ml) of di (2-pyridyl)ketone (1.00 g, 5.44 mmol) was added dropwise to the Grignard reagent with subsequent formation of a white precipitate. The mixture was stirred for 48 hours followed by the addition of 10% HCl (30 mL). The organic solvent was evaporated and the aqueous solution was washed with CH 2 Cl 2 (2 × 50 mL). Neutralization of the aqueous solution with saturated Na 2 CO 3 -solution was followed by extraction with CH 2 Cl 2 (3 × 100 mL). The CH 2 Cl 2 solvent was evaporated and the Scheme 1 Synthetic route for the ligand and complexes investigated in this work. Elemental analysis. The X-ray spectroscopy and electrochemical data of the [M II -Cl] and [Fe II -MeOH] complexes were collected using the powder samples obtained by the synthesis procedures described below, i.e. without any recrystallization step. Below we report the elemental analyses for these powder samples. A fraction of each powder sample was utilized for obtaining the crystals for the single crystal X-ray diffraction structures reported for each complex. The deviations of up to 1.4% from the calculated values indicate a small level of impurities in the powder samples, which we assign to non-redox active and non-target metal containing species, as no signals were found in the EXAFS, XANES and electrochemical data that could be assigned to impurities.  The UV-Vis spectrum shows a strong absorption below 300 nm with a shoulder at 320 nm (ε = 0.75 × 10 3 M −1 cm −1 ), see Fig. S15. † Solid FT-IR (KBr) of the complex shows the same vibration modes as the ligand with a blueshift of 4 cm −1 and the characteristic P-F stretching at 842 cm −1 from the PF 6 − ion. 43

X-ray diffraction
Single crystals suitable for measurements were obtained for all four [M II -Cl] complexes and the [Fe II -DMF] complex by slow diffusion of toluene into a dimethylformamide (DMF) solution of the complex. All the measurements were performed using graphite-monochromatized Mo Kα radiation at 150 K using a Bruker D8 APEX-II equipped with a CCD camera. Data reduction was performed with SAINT. 44 Absorption corrections for the area detector were performed using SADABS. 45,46 The structure was solved by direct methods and refined by fullmatrix least-squares techniques against F 2 using all data (SHELX). 45,46 All non-hydrogen atoms were refined with anisotropic displacement parameters if not stated otherwise. Hydrogen atoms were constrained in geometric positions to their parent atoms. The structures were refined using the OLEX2 suite of programs. 47 See Table S1 † for the detailed parameters.

X-ray absorption spectroscopy
X-ray absorption spectroscopy measurements were performed at the KMC-3 beamline at the BESSY II synchrotron facility in Berlin, Germany. The incident X-ray energy was scanned through the Mn, Fe, Co, and Ni K-edge regions using a silicon (111) double-crystal monochromator. Measurements were performed with samples positioned at 45°with respect to the incident beam in a liquid helium-cooled cryostat (Oxford Instruments). Samples in powder form were mixed with boron nitride in a weight ratio of 1 : 30. Solution samples of the complexes (1 mM) were prepared in dry acetonitrile. A 10 μm Fe, Co, or Ni foil (Goodfellow Cambridge Limited) positioned behind the sample served as an energy calibration standard. For Mn, a thin layer of KMnO 4 was used. Energy calibration was done by assigning the position of the maximum of the first derivative of the absorption of the Fe, Co, and Ni foils to energies 7112 eV, 7709 eV, and 8333 eV, respectively. The position of the KMnO 4 pre-edge was set to 6543.3 eV. Kα fluorescence signals from the samples were recorded with a 13-element silicon drift detector (RaySpec) positioned perpendicular to the incident beam. For each sample, 8-16 scans were taken in the energy range from 120 eV before the corresponding K-edge up to 750 eV after the edge. Each scan was collected on a new sample spot, to avoid possible radiation damage; additionally, three consecutive scans at the same sample spot confirmed that there was no observable radiation damage on the time scale of the XAS measurement. E 0 used for EXAFS extraction was 6539 eV (Mn), 7115 eV (Fe), 7710 eV (Co) and 8334 eV (Ni). EXAFS simulations were done with the FEFF 9.0 software 48 (using settings NLEG 6, CRITERIA 12 5, RPATH 7, SCF 7 1 30 0.05). The amplitude reduction factor S 0 2 was 0.90 for Mn, Co and Ni, and 0.85 for Fe. Least-squares fitting of k 3 -weighted EXAFS data by changing of interatomic distances for the first four single-scattering shells and Debye-Waller factors for all shells (with all multiple-scattering shells having the same Debye-Waller factor) was done with in-house software (SimXLite). The k 3 -weighted EXAFS data were fitted between k = 1.6 Å −1 and 13 Å −1 . Fit parameter errors were determined as described previously. 49

Bulk electrolysis
The solution samples of the reduced complexes (1 mM) in 0.1 M TBAPF 6 acetonitrile were oxidized employing a custommade continuous flow electrosynthesis cell. 50 The applied potential was selected by recording a slow (10 mV s −1 ) cyclic voltammogram in steady condition. The extent of electrolysis was monitored by recording the current response (typically around 0.2 mA at 0.05 ml min −1 flow rate) and changes in the UV-Vis spectrum to provide an estimation of the percentage of the conversion. The oxidized sample solution was collected and frozen in liquid nitrogen for XAS measurements.

Cyclic voltammetry
All solvents, chemicals, glassware, and electrochemical tools were carefully dried. A pseudo reference electrode that consisted of a AgCl coated Ag wire immersed in a 0.1 M TBAPF 6 acetonitrile solution contained in a porous glass tipped tube was employed. The electrochemical stability of the Ag pseudo reference electrode was tested over 2 days by recording the cyclic voltammetry of ferrocene. A drift of 0.014 V was observed in the redox potential (E 1/2 (Fc + /Fc) = 0.47 V vs. Ag pseudo reference, in MeCN). All potentials quoted in this study are relative to the Fc + /Fc couple. The glassy carbon working electrode (3 mm diameter) was polished with alumina particles (0.05 µm) immediately prior use. The counter electrode was a platinum rod polished with sandpaper before use. Unless stated otherwise, the following parameters were used to record cyclic voltammetry on our samples: scan rate: 100 mV s −1 , step potential 0.002 V. The sample concentration was 0.5 mM with tetrabutylammonium hexafluorophosphate supporting electrolyte (100 mM) in acetonitrile.

Density-functional theory
Density-functional theory (DFT) calculations were performed using Gaussian09 E.01. 51 Initial geometry optimization for all complexes was performed using the B3LYP functional with the Lanl2DZ basis set for transition metals, and 6-31G( p,d) for all other atoms. Thermal contributions were calculated from the subsequent Hessian calculations. Structures were then reoptimized using the B3LYP-D3 functional with the larger Lanl2TZ (f ) basis set for the transition metals, and 6-311+G(2df,2pd) for all other atoms. All optimizations were performed in acetonitrile solvent using the default Polarizable Continuum Model (PCM) in Gaussian09. Final energies were calculated with the B3LYP* functional, 52 using the D3 parameters from B3LYP. 53 This functional was previously shown to give good results for spin-state energetics in the spin-crossover complex [Fe II -Cl]. 15 Final energy calculations are done using the Solvation Model based on Density (SMD). This model gives good absolute solvation energies, which is important for the calculation of half potentials. 54  including the apical ligand X. For reactions with differences in concentration between components, a configurational entropy factor has been added according to Boltzmann's formula, S = k B ln W, where k B is Boltzmann's constant and W is the number of solvent molecules per complex. Calculated redox potentials are obtained using the reaction [M II -X] ⇌ [M III -X] + e − . The energy of the solvated electron was calculated using 4.28 V for the absolute potential of the standard hydrogen electrode (SHE). 55 The choice of reference value affects the absolute potentials, but not the comparison between complexes. Calculated values are given in V vs. Fc + /Fc for comparison with experimental results, with 0.40 V for the E 1/2 of Fc + /Fc vs. the SHE. 56 The analysis of the different contributions to the redox potential was made by separating the calculations of the oxidized species into three steps. First, the energy for removing a spin-up σ electron was calculated using the structure of the reduced complex (orbital contribution). Second, the structure was optimized using the same spin multiplicity (geometric). Finally, the energy of the most stable spin multiplicity was calculated (spin).

Conflicts of interest
There are no conflicts to declare.