Spectroelectrochemical and Theoretical Study of [Si(ttpy)2](PF6)4: A Potential Polychromatic Electrochromic Dye

Complexes consisting of earth-abundant main group metals such as silicon with polypyridine ligands are of interest for a variety of optical and electronic applications including as electrochromic colorants. Previous spectroelectrochemical studies with tris(2,2′-bipyridyl)silicon(IV) hexafluorophosphate, [Si(bpy)3](PF6)4, demonstrated an ability to control the color saturation of the potential electrochromic dye, with the intensity of the dye’s green color increasing as the charge state sequentially reduces from 4+ to 1+. In this study, the synthesis of bis(4′-(4-tolyl)-2,2′:6′,2″-terpyridine)silicon(IV) hexafluorophosphate, [Si(ttpy)2](PF6)4, is reported along with electrochemical and spectroelectrochemical analyses. Computational modeling (density functional theory) is used to further elucidate the electrochromic properties of previously reported Si(bpy)3n+ species and the new Si(ttpy)2n+ species. While the homoleptic tris(bidentate)silicon(IV) complexes are attractive as electrochromic dyes for tunable color saturation, the bis(tridentate)silicon(IV) complexes are attractive as polychromatic electrochromic dyes.


Introduction
Silicon, an earth-abundant element with a strong tendency to maintain the tetravalent oxidation state, provides attractive attributes to coordination complexes for electronic and optical applications. Polypyridine ligands can bind Si(IV) to form robust and inert hexacoordinate silicon complexes capable of reversible redox chemistry with multiple stable oxidation states [1]. For example, the tris(2,2 -bipyridyl)silicon(IV) complex cation, [Si(bpy) 3 ] 4+ was shown via cyclic voltammetry to possess three chemically reversible oneelectron reductions at reduction potentials significantly less negative than other bipyridine complexes due to the tetravalent nature of the silicon ion [2]. Computational modeling of [Si(bpy)] n+ complexes by Wieghardt and coworkers indicated silicon remains tetravalent in the reduced species with all the reductions being ligand localized [3]. Their modeling showed Si(bpy) 3 1+ , for example, contains three ligand-localized SOMOs (≤3% Si character), and Mulliken spin density population analysis places one unpaired spin on each ligand and <0.01 on the silicon atom. It therefore has the predicted electronic structure best described as [Si IV (bpy˙) 3 ] + . Furthermore, using density functional theory (DFT) with broken-symmetry they were also able to show that the previously reported black, air-sensitive compound [Si(bpy) 3 ] 0 is best described as [Si IV (bpy • ) 2 (bpy 2-)] 0 , and it possesses a singlet diradical ground state with a coupling constant J calcd of −490 cm -1 , which is remarkably close to the J obs of −374 cm -1 [4,5].
We previously demonstrated that polypyridylsilicon(IV) complexes are attractive candidates for cathodic colorants, which are an integral component of dual-active electrochromic windows and displays [6]. In dual-active electrochromics, an electrochromic Molecules 2022, 27, 8521 2 of 10 dye attached to an anode is reinforced by a complimentary cathodic electrochromic dye (cathodic colorant) [7][8][9][10]. The cathodic colorant should be colorless (or nearly colorless) in its oxidized or bleached state and intensely colored upon reduction. The N,N -dialkyl and N,N -diaryl-4,4 -bipyridinium salts are perhaps the best known and widely used cathodic colorants, since they possess robust and reliable spectroelectrochemistry, particularly for blue or violet coloration [11][12][13][14][15]. Polypyridyl transition metal complexes have also been evaluated as promising electrochromic dyes [16], but the metal to ligand charge transfer bands (MLCT) tend to preclude the desired colorless or bleached state. By contrast polypyridylsilicon(IV) complexes possess several advantages: (1) they exhibit multiple stable oxidation states with well-behaved chemically reversible electrochemistry, (2) they are often colorless in the fully oxidized (+4) state, (3) they have very low reduction potentials attributed to the tetravalent silicon center, which could enable low-voltage devices and applications, (4) their reduced ligand based transitions appear in the visible portion of the spectrum, and (5) they can be synthesized easily from earth-abundant elements.
In this manuscript we report the spectroelectrochemical examination of [Si(ttpy) 2 ](PF 6 ) 4 . We also provide an experimental and theoretical comparison of the optoelectronic properties of the bis(terpyridyl)silicon(IV) motif with the previously reported tris(bipyridyl)silicon(IV)motif and discuss their relative merits for electrochromic applications.
We previously demonstrated that polypyridylsilicon(IV) complexes are attractive candidates for cathodic colorants, which are an integral component of dual-active electrochromic windows and displays [6]. In dual-active electrochromics, an electrochromic dye attached to an anode is reinforced by a complimentary cathodic electrochromic dye (cathodic colorant) [7][8][9][10]. The cathodic colorant should be colorless (or nearly colorless) in its oxidized or bleached state and intensely colored upon reduction. The N,N′-dialkyl and N,N′-diaryl-4,4′-bipyridinium salts are perhaps the best known and widely used cathodic colorants, since they possess robust and reliable spectroelectrochemistry, particularly for blue or violet coloration [11][12][13][14][15]. Polypyridyl transition metal complexes have also been evaluated as promising electrochromic dyes [16], but the metal to ligand charge transfer bands (MLCT) tend to preclude the desired colorless or bleached state. By contrast polypyridylsilicon(IV) complexes possess several advantages: (1) they exhibit multiple stable oxidation states with well-behaved chemically reversible electrochemistry, (2) they are often colorless in the fully oxidized (+4) state, (3) they have very low reduction potentials attributed to the tetravalent silicon center, which could enable low-voltage devices and applications, (4) their reduced ligand based transitions appear in the visible portion of the spectrum, and (5) they can be synthesized easily from earth-abundant elements.
In this manuscript we report the spectroelectrochemical examination of [Si(ttpy)2](PF6)4. We also provide an experimental and theoretical comparison of the optoelectronic properties of the bis(terpyridyl)silicon(IV) motif with the previously reported tris(bipyridyl)silicon(IV)motif and discuss their relative merits for electrochromic applications.
[Si(ttpy)2]I4 was converted to the hexafluorophosphate salt, [Si(ttpy)2](PF6)4 by metathesis with ammonium hexafluorophosphate (Scheme 1).  Cyclic voltammagrams of [Si(ttpy) 2 ](PF 6 ) 4 in acetonitrile with tetrabutylammonium hexafluorophosphate (TBAPF 6 ) as electrolyte indicated multiple stable redox states, with four chemically reversible, single electron reduction waves (Figure 1 top). The ∆E p (∆E p = E c − E a ) for the waves ranged from 76 to 86 mV at a scan rate of 200 mV/s (Table 1). For comparison, previously reported values for the first four reduction waves of [Si(bpy) 3 ](PF 6 ) 4 and [Si(terpy) 2 ](PF 6 ) 4 , where terpy is 2,2 :6 ,2 -terpyridine, are also included [2] The first reduction, E o 1/2 (4+/3+) = −0.403 V, is comparable to the previously reported first reduction potentials of [Si(terpy) 2 ](PF 6 ) 4 and [Si(bpy) 3 ](PF 6 ) 4. The low voltage required for reduction of polypyridylsilicon(IV) species is advantageous for electrochromic applications in that it allows for lower operational voltage and power demands and allows matching to low voltage anodic colorants. The difference between sequential reduction potentials is sufficient to exceed 99.9% mole fractions for the 4+, 2+, and neutral species. The cluded [2] The first reduction, E / (4+/3+) = −0.403 V, is comparable to the previously reported first reduction potentials of [Si(terpy)2](PF6)4 and [Si(bpy)3](PF6)4. The low voltage required for reduction of polypyridylsilicon(IV) species is advantageous for electrochromic applications in that it allows for lower operational voltage and power demands and allows matching to low voltage anodic colorants. The difference between sequential reduction potentials is sufficient to exceed 99.9% mole fractions for the 4+, 2+, and neutral species. The 3+ and 1+ species can achieve a maximum mole fraction of 96.4% and 99.1%, respectively (Figure 1 bottom). The UV-vis spectrum of [Si(ttpy)2](PF6)4 in acetonitrile as a function of applied potential was obtained using an argon flushed spectroelectrochemical cell with a 1.0 mm path length that consisted of a printed gold honeycomb working electrode, a gold counter electrode, and a Ag/AgCl reference electrode standardized to the ferrocene/ferrocenium couple ( Figure 2). Acetonitrile was selected as the solvent due to the low solubility of [Si(ttpy)2](PF6)4 in less polar solvents and the wide UV-vis and electrochemical windows acetonitrile provides. The sample was dissolved in acetonitrile with TBAPF6 (0.1M) as the electrolyte. Reversible transformations were observed with clean isosbestic points for each a Previously reported data from [2].
The UV-vis spectrum of [Si(ttpy) 2 ](PF 6 ) 4 in acetonitrile as a function of applied potential was obtained using an argon flushed spectroelectrochemical cell with a 1.0 mm path length that consisted of a printed gold honeycomb working electrode, a gold counter electrode, and a Ag/AgCl reference electrode standardized to the ferrocene/ferrocenium couple ( Figure 2). Acetonitrile was selected as the solvent due to the low solubility of [Si(ttpy) 2 ](PF 6 ) 4 in less polar solvents and the wide UV-vis and electrochemical windows acetonitrile provides. The sample was dissolved in acetonitrile with TBAPF 6 (0.1M) as the electrolyte. Reversible transformations were observed with clean isosbestic points for each of the subsequent reductions from the 4+ state to the neutral species. The 4+ species could be recovered upon oxidation. The fully oxidized [Si(ttpy) 2 ](PF 6 ) 4 has a λ max = 408 nm for its lowest energy peak with a low-energy tail extending out to approximately 460 nm. Consequently, [Si(ttpy) 2 ](PF 6 ) 4 has a slight yellowish appearance and not the preferred colorless appearance of [Si(bpy) 3 ](PF 6 ) 4 . Photos of [Si(ttpy) 2 ](PF 6 ) 4 in its five separate redox states surrounded by the honeycomb electrode are provided in Figure 3 along with previously obtained photos of [Si(bpy) 3 ](PF 6 ) 4 in its four chemically reversible redox states. Isolated spectra of [Si(ttpy) 2 ](PF 6 ) 4 and [Si(bpy) 3 ](PF 6 ) 4 for each oxidation state measured at its maximum mole fraction are included in Figure 4 along with theoretical modeling results, which are discussed later.       Figure 4a were previously reported in [6].
The net effect for electrochromic application is a voltage-tunable color intensity as seen in Figure 3 with increasing reduction leading to a more intense green color. By comparison, the [Si(ttpy) 2 ] 4+ is seen to proceed through multiple color transitions starting from a yellow species (4+), to a cyan (3+ and 2+) and finally a magenta color (1+ and neutral). The UV-vis spectrum (Figure 4b.) shows that the first reduction of [Si(ttpy) 2 ] 4+ introduces a broad peak centered around 650 nm, which we attribute to formation of a tolyl-terpyridine radical anion and a species best represented as [Si IV (ttpy)(ttpy • )] 3+ . The second electron appears to generate a second tolyl-terpyridine radical anion [Si IV (ttpy • ) 2 ] 2+ with a doubling of the broad band at 650 nm. The third reduction adds a new peak around 500 nm in addition to the 650 nm peak, which are assigned to the ttpy dianion and ttpy radical anion, respectively, for [Si IV (ttpy 2-) (ttpy • )] 1+ . The fourth reduction reduces the second ttpy radical anion, producing [Si IV (ttpy 2-) 2 ] 0 consistent with the doubling of the ttpy dianion peak at 500 nm and a disappearance of the radical anion peak at 650 nm.
This explanation is consistent with a ligand-localized, "high-spin" approximationelectrons occupy separate localized ligands until forced to occupy an already reduced ligand. Due to the non-involvement of the tetravalent silicon center, we summarize the species as follows.
[Si IV (bpy) 3  To test this theory, DFT and TD-DFT calculations were performed using the B3LYP functional with the 6-31G* basis set. For this study, spin-restricted Kohn-Sham (RKS) solutions were obtained for closed-shell (S = 0) species and unrestricted Kohn-Sham (UKS) solutions were obtained for open shell species. Although, Wieghardt demonstrated the broken symmetry (BS) DFT approach provided spin-coupled solutions for reduced Si(bpy) 3 n+ species in the gas phase that were in better agreement with experimental observations, we limited our analysis in this study to the RKS and UKS solutions due to computational cost and difficulties with integrating BS-DFT with solvation and TD-DFT.
Each species was optimized and confirmed to be a local minimum (no negativefrequency vibrational modes). An optimized structure of [Si(ttpy) 2 ] 4+ , S = 0, with solvation (acetonitrile) along with the HOMO and LUMO of the species is shown in Figure 5. The HOMO and HOMO-1 were nearly degenerate, ∆E < 0.0001 Hartrees, with the HOMO-1 localized on the opposite ligand. Similarly, the LUMO and LUMO+1 were also nearly degenerate with ∆E < 0.0001 Hartrees. The results of all calculations are summarized in Table 2. The formatted checkpoint files for all 30 calculations are available for download in the Supplementary Materials. For both species the high-spin states were generally found to be more stable. A notable exception was for [Si(ttpy) 2 ] 1+ in the gas state, for which the low-spin state (S = 0.5) was found to be slightly more stable than the high-spin state (S = 1.5). The situation was reversed for [Si(ttpy) 2 ] 1+ when solvation effects were included though. In addition, the closed shell, low-spin (S = 0) solution is the most stable for neutral [Si(ttpy) 2 ].
Molecules 2022, 27, 8521 7 of 10 2. The formatted checkpoint files for all 30 calculations are available for download in the Supplementary Materials. For both species the high-spin states were generally found to be more stable. A notable exception was for [Si(ttpy)2] 1+ in the gas state, for which the lowspin state (S = 0.5) was found to be slightly more stable than the high-spin state (S = 1.5). The situation was reversed for [Si(ttpy)2] 1+ when solvation effects were included though. In addition, the closed shell, low-spin (S = 0) solution is the most stable for neutral [Si(ttpy)2]. Simulated UV-vis structures were constructed from the RKS and UKS ground state structures using TD-DFT calculations (500 transitions, solvent=acetonitrile). Figure 4c,d shows the results using the high-spin ground states where available and the lower row (e and f) shows the results from using the low-spin ground states where available. Note that for [Si(ttpy)2], only the low-spin (S = 0, closed shell) ground state was included in TD-DFT calculations.
For [Si(bpy)3] n+ , TD-DFT (B3LYP/6-31G*, with solvation) provides qualitative agreement between the simulated and experimental UV-vis spectra of the respective complexes when the high-spin ground state is used. Specifically, in Figure 4c, both the high energy peak around 390 nm and the broad absorbance from 600-1000 nm that are associated with the bipyridine radical anion are evident (although slightly shifted to higher energies) and they are seen to increase in intensity with increasing reductions of the overall charge. Further optimization of functionals and basis sets should improve the absolute energy alignment of these transitions. On the other hand, TD-DFT from the low-spin ground state (4e) is unsuccessful in reproducing the experimental features.
The distinction is a little less clear with Si(ttpy)2. The spectra of the three unambiguous species [Si(ttpy)2] 4+ , [Si(ttpy)2] 3+ , and [Si(ttpy)2] 0 nicely correlate with the observed transitions in the UV-vis (although shifted in energy). With respect to [Si(ttpy)2] 2+ and [Si(ttpy)2] 1+ , both simulated spectra appear to show peaks consistent with the experimentally observed peaks at 500 nm and 650 nm. The inability to clearly distinguish a better model between these two solutions, suggests spin-coupling may be considerable in these species and a broken symmetry BS-DFT approach may be required for optimal reproduction of the spectroelectrochemical properties for these two states.  Simulated UV-vis structures were constructed from the RKS and UKS ground state structures using TD-DFT calculations (500 transitions, solvent = acetonitrile). Figure 4c,d shows the results using the high-spin ground states where available and the lower row (e and f) shows the results from using the low-spin ground states where available. Note that for [Si(ttpy) 2 ], only the low-spin (S = 0, closed shell) ground state was included in TD-DFT calculations.
For [Si(bpy) 3 ] n+ , TD-DFT (B3LYP/6-31G*, with solvation) provides qualitative agreement between the simulated and experimental UV-vis spectra of the respective complexes when the high-spin ground state is used. Specifically, in Figure 4c, both the high energy peak around 390 nm and the broad absorbance from 600-1000 nm that are associated with the bipyridine radical anion are evident (although slightly shifted to higher energies) and they are seen to increase in intensity with increasing reductions of the overall charge. Further optimization of functionals and basis sets should improve the absolute energy alignment of these transitions. On the other hand, TD-DFT from the low-spin ground state (4e) is unsuccessful in reproducing the experimental features.
The distinction is a little less clear with Si(ttpy) 2 . The spectra of the three unambiguous species [Si(ttpy) 2 ] 4+ , [Si(ttpy) 2 ] 3+ , and [Si(ttpy) 2 ] 0 nicely correlate with the observed transitions in the UV-vis (although shifted in energy). With respect to [Si(ttpy) 2 ] 2+ and [Si(ttpy) 2 ] 1+ , both simulated spectra appear to show peaks consistent with the experimentally observed peaks at 500 nm and 650 nm. The inability to clearly distinguish a better model between these two solutions, suggests spin-coupling may be considerable in these species and a broken symmetry BS-DFT approach may be required for optimal reproduction of the spectroelectrochemical properties for these two states.

Materials and Methods
All quantum chemical calculations were performed using density functional theory (DFT) realized in Gaussian 16 Rev. C.01 program packages [17][18][19][20]. We performed calculations applying Becke three parameter hybrid functional [21] which uses the non-local correlation provided by the Lee, Yang, and Parr expression-B3LYP [22,23]. As a basis set we have employed standard basis set 6-31G* based on a gaussian type of functions. When applicable, solvation was included using the scrf = (solvent = acetonitrile) command in Gaussian 16, which employs a Polarizable Continuum Model (PCM) [24]. The combination of functional and basis set was selected based on the success of previous attempts to model the electronic spectra of hexacoordinate silicon compounds using time-dependent DFT approach. Simulated UV-vis spectra were constructed by transforming calculated UV-vis transitions to Gaussian curves with 0.2 eV FWHM. For all calculated structures optimized structures were obtained. Frequencies calculations demonstrated absence of negative values indicating that all obtained structures correspond to the global minimum.
All reagents were used as received without further purification. SiI 4 and 2-picoline were all stored and used under a nitrogen atmosphere in the glovebox. Acetonitrile was HPLC-grade solvent from Sigma AldrichNMR spectra were acquired using a JEOL 300 MHz NMR spectrometer and referenced to tetramethylsilane (TMS).
Cyclic voltammetry was performed using a Princeton Applied Research model 273A potentiostat/galvanostat employing a conventional three-electrode setup consisting of a platinum disk working electrode, a Ag/AgCl reference electrode, and a platinum wire auxiliary electrode. Positive feedback iR compensation was routinely used. Acetonitrile with 0.100 M TBAPF 6 was used as the solvent and supporting electrolyte and ferrocene was used as an internal standard. Ferrocene was used as an internal standard with all potentials reported as V versus Fc/Fc+. Plots of I pf vs. V 1/2 for all reduction waves are included in the (Supplemental Materials Figures S6-S9).
The spectroelectrochemical cell was purchased from Pine Research Instrumentation and consisted of a covered quartz cuvette with a 1.0 mm pathlength with a printed gold or platinum honeycomb working electrode, a matching (gold or platinum) counter electrode, and a Ag/AgCl reference electrode standardized to the ferrocene/ferrocenium couple ( Figure S10). Samples were prepared in acetonitrile with 0.100 M TBAPF 6 and degassed with argon that had been bubbled through acetonitrile. Potentials were then applied to the solution with a Princeton Applied Research model 173 potentiostat/galvanostat and analyzed with an Agilent 8453 Diode Array Spectrometer. UV-Vis spectra were taken in 25 mV increments with at least 45 s of equilibrium time between applied voltage values. 4 -p-tolyl-2,2 :6 ,2"-terpyridine (ttpy): The ttpy ligand was synthesized following the procedure of Liu and coworkers [25]. 100 mL of ethanol was stirred in a 250 mL RBF equipped with a reflux condenser. 4 Figure S1.
[Si(ttpy) 2 ]I 4 : Ttpy ligand (0.244 g, 0.76 mmol) was combined with 0.16 g (0.30 mmol) of SiI4 and 3.00 mL of pyridine in a heavy walled glass pressure vessel equipped with a stir bar inside a nitrogen-filled glovebox. The vessel was transferred to an oil bath heated to 125 • C for 3.5 h. Following the reaction, the vessel was opened and 1 mL of methanol was added to quench unreacted residual iodosilanes. The contents of the vessel were transferred to a centrifuge tube, and the contents were centrifuged to afford a dark precipitate. The supernatant was decanted and the precipitate was rinsed twice with chloroform and twice with diethyl ether before drying overnight in a vacuum oven. This produced 0.302 g of [Si(ttpy) 2 ]I 4 which was used in the next step without futher purification. 1 Figure S2.
[Si(ttpy) 2 ](PF 6 ) 4 : 0.050 g (0.042 mmol) of [Si(ttpy) 2 ]I 4 was dissolved in 40 mL of water in a centrifuge tube. Addition of 5 mL (large excess) of saturated NH 4 PF 6 in water caused immediate precipitation of a yellow-colored salt which was isolated by centrifugation and rinsed with deionized water two more times. 1

Conclusions
Polypyridylsilicon(IV) complexes continue to hold promise as cathodic colorants in electrochromic applications. Due to the non-involvement of the tetravalent silicon center in the redox activity of the complex, ligand localized reductions dominate the observed spectroelectrochemical properties. Furthermore, a structure-property relationship is evident with homoleptic tris-bidentate complexes such as Si(bpy) 3 n+ providing a motif for voltage tunable variable color-saturation cathodic colorants and homoleptic bis-tridentate complexes such as Si(ttpy) 2 n+ providing a motif for voltage tunable polychromatic cathodic colorants. TD-DFT (B3LYP/6-31G*, with solvation) has been shown to accurately simulate all of the features of the UV-vis spectra of Si(bpy) 3 n+ throughout its positive charge states, although further benchmarking of the functionals and basis sets to better match the energy levels with observed energy levels is warranted. On the other hand, TD-DFT (B3LYP/6-31G*, with solvation) failed to accurately model some of the open-shell states of Si(ttpy) 2 n+ suggesting a broken symmetry DFT approach allowing for radical coupling may be necessary to model these states. Continued experimental and theoretical development of polypyridylsilicon(IV) complexes is still needed to facilitate broader electrochromic applications.

Informed Consent Statement: Not applicable.
Data Availability Statement: All data are fully available and included in the manuscript.