Next Article in Journal
Discovery, Development, and Patent Trends on Molnupiravir: A Prospective Oral Treatment for COVID-19
Previous Article in Journal
Chloropyridinyl Esters of Nonsteroidal Anti-Inflammatory Agents and Related Derivatives as Potent SARS-CoV-2 3CL Protease Inhibitors
Previous Article in Special Issue
Structural Relevance of Intramolecular H-Bonding in Ortho-Hydroxyaryl Schiff Bases: The Case of 3-(5-bromo-2-hydroxybenzylideneamino) Phenol
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

The Contribution of Density Functional Theory to the Atomistic Knowledge of Electrochromic Processes

Dipartimento di Chimica e Tecnologie Chimiche, Università della Calabria, 87036 Rende, CS, Italy
*
Author to whom correspondence should be addressed.
Molecules 2021, 26(19), 5793; https://doi.org/10.3390/molecules26195793
Submission received: 13 July 2021 / Revised: 8 September 2021 / Accepted: 19 September 2021 / Published: 24 September 2021
(This article belongs to the Special Issue Fundamentals and Applications in Quantum Chemistry)

Abstract

:
In this review, we provide a brief overview of the contribution that computational studies can offer to the elucidation of the electronic mechanisms responsible for the electrochromism phenomenon, through the use of the density functional theory (DFT) and its time-dependent formulation (TDDFT). Although computational studies on electrochromic systems are not as numerous as those for other physico-chemical processes, we will show their reliability and ability to predict structures, excitation energies, and redox potentials. The results confirm that these methods not only help in the interpretation of experimental data but can also be used for the rational design of molecules with interesting electrochromic properties to be initiated for synthesis and experimental characterization.

1. Introduction

Electrochromism refers to the phenomenon in which a chemical system changes color following the application of an electric field of suitable value, returning to the original color once the applied field is eliminated. The change in the absorption band is due to the triggering of a redox reaction induced by the electric field. Compounds possessing more than two redox states are subject to different color changes as the applied electric field reaches the given redox potential. The whole process is described in Scheme 1.
Although the term electrochromic originally referred to a color change from a transparent (absorption in the UV region) to a colored species, today it is extended to all shifts of the absorption band in the entire electromagnetic spectrum caused by the application of an electric field.
Since Platt’s pioneering work in 1961 [1], electrochromism has been the subject of a large number of studies for its potential applications in different fields of modern technologies [2,3,4,5,6,7]. To date, electrochromic materials are mainly used for the manufacture of smart and intelligent windows that allow one to control the flow of heat and light [8,9,10], but show promising applications in glare reduction [11], mirroring [12], low-consumption display panels [13,14], transistors and diodes [15,16,17], catalysts [18], solar cells [19,20], camouflage materials [21], memory [22,23] and energy storage [24,25] devices, molecular machines [26], electrodes for batteries and supercapacitors [27,28,29,30,31], and antibacterial agents [32,33,34].
The most-studied systems are the viologen (1,1′-disubstituted-4,4′-bipyridinium salts) and its derivatives, for which many experimental studies have shown intrinsic properties that make them excellent candidates for a series of technological applications [35,36,37,38,39,40].
In recent years, however, the investigation has expanded to other promising organic and inorganic systems including those containing boron [41,42,43], quinodimethane derivatives [44,45], metal-phtalocyanine [46], metal oxides-based film [7,47,48], oligoaniline [49], antracenes [50,51], azacenes [52], polymethine and polyimides [53,54], and some organometallic complexes [55,56].
From a theoretical point of view, the studies are scarcer and have concerned some viologen-like [57], dithiolodithiole and thiophene [58], thiophene−pyrrole-containing rings [59], borepins [60], poly-dioxythiophene [61], and methyl ketone [62] systems. It is also worth mentioning certain investigations in which experimental data are coupled to theoretical data, demonstrating how the two approaches can be combined to provide useful insights on these interesting processes [36,39,53,56,63,64,65,66,67,68,69,70,71].
Despite this amount of work, accurate information at the atomistic level of electrochromic processes is somehow lacking, and modern theoretical methods can help to shed light on these aspects. Among these, those based on the Density Functional Theory (DFT) have proved to be very useful since they are able to provide accurate results even for medium to large systems.
In this short review, we summarize the contribution of DFT methods in improving the current knowledge of the electrochromism phenomenon at the atomistic level. Particular emphasis will be devoted to those studies in which they computed and analyzed properties often difficult to obtain experimentally. Moreover, the outcomes of such investigations are not only useful to explain the material properties but also to build a rational design strategy to synthesize new molecules with given properties.
The computed properties we are interested in include structural determination, absorption behaviors and their origins, redox potentials, and electronic features.
The structures of the considered systems are reported in Scheme 2.

2. Results and Discussion

2.1. Structural Determination

In order to show the accuracy of DFT methods in determining geometric parameters with accuracy, we report the calculated structures for quadruply benzannulated nonplanar borepin (1) [60] and for 4-(6-(4-Metylphenyl)-[1,2]dithiolo[4,3-c][1,2]dithiol-3-yl)benzonitrile (2) [58], using three different exchange and correlation functionals (see Scheme 2 and Table 1).
The values of the main geometrical parameters reported in Table 1 clearly show very good agreement with the previously reported X-ray crystallographic data [71].
In particular, for diphenyl- [1,2] dithiol [4,3-c] [1,2] we obtain C-C, C-S, and S-S bond length values that are very similar for the three functionals used and in the range of the experimental errors. The same agreement is found for the C-C-C, C-C-S, and C-S-S valence angles for which the error, with respect to the experimental counterpart [71], is about one degree.
Concerning the dihedral angles φ1 and φ2, a higher deviation from the experiment is noted for all the employed functionals. This can be explained by considering that, in the solid phase, the value of this parameter is influenced by the other molecules present in the crystalline structure, while our result is related to a single molecule.
For borepin molecules, we note that the maximum deviation from the experimental B-C and C-C bond lengths [41] is found in the case of the M06 functional, which does, however, provide values in the range of experimental error (Table 1).
Theoretical values obtained for the C-B-C valence angles and the C-C-C-C torsional angle obtained, using the three functionals, are in excellent agreement with the experimental X-ray results [41].
As previously mentioned, the experimental structures relative to the species with different oxidation states obtained upon the application of the electric field are often missing and difficult to measure. The use of calculations based on quantum mechanics can be useful to overcome this lack of information. In Figure 1, we report the main geometrical parameters obtained for viologen 3 as a function of the oxidation state.
The bond between the carbon atom of the CH3 group and the nitrogen of the pyridine ring (C1-N) undergoes a decrease, going from +2 to +1 radical to neutral charge states while the N-C2 increases following the consecutive addition of electrons. The C2-C3 and C4-C5 bond lengths involved in the conjugation of the ring decrease.
A more significant structural change is evidenced by analyzing the inter-ring dihedral angle that assumes a value of 39° in the case of the more-stable 32+ oxidation state but becomes planar in the case of 3+1 (0.13°) and 3 (0.02°) states.
These two examples account well for the importance of determining the structural variations induced by the change in the oxidation state on electrochromic systems, achievable by using reliable theoretical tools.

2.2. Excitation Energies

To assess the performance of DFT-based methods, we report the computed excitation energies for 3 [57] and squaraine 4 [70] in their different oxidation states by using the B3LYP, M06, and wB97XD exchange correlation functionals. The range-separated hybrid wB97XD functional has been chosen since it also takes into account weak interactions. The results are reported in Table 2 together with the corresponding experimental UV-Vis values [70,72].
For system 3, the maximum deviation with respect to the experiment [72] (113 nm) is found in the case of the wB97XD functional that also shows an average error of about 57 nm considering all the oxidation states. The other two functionals used give similar results with average errors of 24 and 20 nm for M06 and B3LYP, respectively. The higher average error was obtained with the functional WB97XD, although it falls within an acceptable error range (0.1–0.3 eV), which could be due to the fact that at short range, it contains 22% of Hartree-Fock exchange and in the intermediate region is smoothly described by an error function parameter of 0.2.
In discussing the spectral behavior following the increase of the applied electric field, we note that they are in close agreement with the experiment. In fact, the dication 32+ species has a measured maximum absorption band of 337 that is well reproduced by all the employed functionals.
The addition of an electron generating the 3+ cation radical slightly red-shifts this band and two new bands in the Vis/NIR region appear (447 and 1100 nm). Both M06 and B3LYP functionals reproduce these spectral changes well. A further injection of an electron causes the 3 species to change the spectral feature that now presents two bands at 558 and 712 nm. Again, the obtained excitation energies for this oxidation state give results very close to the experimental counterparts. The analysis of the excitation behavior allows one to individuate the origin of the transitions. Table 2 shows the absorption at 368 nm (B3LYP) is mainly HOMO-LUMO (H→L) in nature (70 %), while it becomes βHOMO-1→LUMO in the radical cation. The other two peaks present in the Vis/NIR spectrum are due to the βH→L and αH→L transitions, respectively. Likewise, the computed excitation energies for the 3 neutral molecule allows one to assign the two recorded absorptions to the transitions HOMO→ LUMO + 3 and H→L, respectively (see Table 2).
The squaraines containing thiophene 4 are interesting since their structures contain donor–acceptor–donor moieties, which make them useful for a variety of applications. Molecule 4 is also important because it also has, in addition to the neutral one, positive (+1, +2) and negative (−1) oxidation states that are generated in Vis/NIR spectroelectrochemistry experiments. Table 2 shows how for these systems, the excitation energies obtained are also close to the experimental counterparts especially when the M06 and B3LYP functionals are used. A major deviation is found in the neutral species, but it is due to the possible intermolecular aggregation as demonstrated recently [70]. The average error obtained with these methods (about 25 nm) confirms their reliability in the treatment of electronic properties of complex systems. The analysis of the HOMO and LUMO orbitals involved in the transitions (Figure S1) allows one to establish that the band in the 450–500 nm range can be attributed to a charge-transfer between olate donors to the carbon atoms of the central cyclobutenediylium fragment and the thienyl moiety, for the 42+ species [70], while the main transition computed for molecule 4 is a π−π* band, as also demonstrated by the computed density difference plot with the first excited state considered (Figure S2).

2.3. Redox Potentials

The possibility to obtain reliable redox potentials at the theoretical level allows one to examine the energetics of the redox processes of a chemical species. Furthermore, it is useful to set up the spectroelectrochemical and cyclovoltametric experiments as we have indications on the potentials of the electric fields to apply.
Theoretically, the redox potentials can be obtained following the Nernst equation
E = −Δ(ΔGs)redox/nF
where the Faraday constant (F) and the number of electrons transferred in the process (n) are present. The free energy change, Δ(ΔGs)redox, is computed on the basis of the Born–Haber thermodynamic cycle [73]. For a given system (A), the Δ(ΔGs)redox can be obtained by computing the solvation energies (Δ Gs) for the neutral and cationic state by using the following equation:
Δ(ΔGsolv)redox = ΔGgredox + ΔGs(Red) − ΔGs(Ox)
where ΔGgredox, ΔGs(Red), and ΔGs(Ox) are the change of standard Gibbs free energy for the gas phase reaction and solvation energies of reduced and oxidized species, respectively.
In order to compare the results with the experimental data, the computed redox potentials are provided relative to that of ferrocene (E(Fc+/Fc).
In Table 3, we report the computed reduction potentials (V) for a series of electrochromic materials obtained by using the M06 exchange correlation potential, that in our experience is the most suitable for these determinations [57,60]. Comparison with the relative experimental data [41,43,72,74,75] is, generally, quite good.
A more detailed comparison between computed and experimental data reveals that the difference varies from 0 to 0.43 and the average error turns out to be 0.15 and 0.7 V for the first and second reduction potentials, respectively.

2.4. Electronic Properties

During the redox process, radical species can be obtained. As an example, in viologen-like systems, starting from the 32+, after the first reduction step, the 3+ radical species is formed. It is well known that the stability of the radicals strongly depends on their capability to delocalize the unpaired electron: Higher delocalization means great stability. Theoretical methods allow one to compute the spin distribution over the entire set of molecules as reported in Figure 2 for the 3+ radical.
The examination of this figure indicates that the unpaired electron is not only localized on the electron-rich nitrogen, but it is distributed on all the conjugated pyridine rings and on the C-C connecting bond with resulting system stabilization.
Furthermore, the computation of the molecular electrostatic potential allows one to follow the charge flow during the reduction steps. The visualization of what happens after the electron addition for the 32+ + e→3+. + e→3 reduction behavior is given in Figure 3.
The analysis indicates how the positive charge, although reduced, continues to be distributed over the entire molecular surface upon the first reduction step. The second reduction process is characterized by different charge distribution with negative charges concentrated in the aromatic rings and between them.
The knowledge of the frontier molecular orbitals and their composition can give further information on the electronic changes occurring during the redox processes. In Figure 4, we report the HOMO–LUMO energies and the orbital composition for the three oxidation states of 3. It is interesting to note that both HOMO and LUMO orbitals are significantly destabilized in going from the +2 to 0 oxidation states. The HOMO–LUMO energy gap changes during the reduction process, assuming the highest value in 32+ (5.27 eV) and the lowest in the neutral (3) oxidation state. The changes in the molecular orbital composition evidence the effect of the addition of electrons in going from 32+ to 3, especially for the HOMO orbital.

3. Materials and Methods

Full geometry optimizations in solvent media have been performed using the 6–31G* basis set in conjunction with the hybrid B3LYP [76,77], meta-hybrid M06 [78], and the range-separated hybrid wB97XD [79,80] exchange-correlation functionals. After the geometry optimization step, frequency calculations have been performed, at the same level of theory, to verify if the obtained structures are real minima. The unrestricted Kohn–Sham formalism has been used in the case of radical species. On these species, the analysis of the 〈S2〉 value (in all cases, reported at about 0.77) ensures the lack of spin contamination. The time-dependent density functional linear response theory (TDDFT) has been used to obtain the excitation energies starting from the previous optimization structures and employing the 6–31+G* basis set. Solvent effects have been taken into account by using the non-equilibrium polarizable continuum model [81]. All the computations have been performed with the Gaussian 09 code [82].

4. Conclusions

In this contribution, we have summarized the great potential of density functional theory in explaining the atomistic processes involved in the electrochromism of some important molecular materials. The main benefits of the theoretical computational investigation can be briefly summarized:
A reliable prediction of the molecular structures with errors comparable to that coming from solid-state X-ray measurements. In addition, DFT computations allow one to determine the structures of systems with unstable oxidation states that are, consequently, difficult to measure experimentally in laboratory conditions;
The reproduction of the excitation energies in all the molecular oxidation states is close to the experimental counterparts, with an error of about 0.2–0.3 eV. It is possible to assign the transitions and analyze the orbital origin and possible charge transfer processes;
The redox potentials can be easily obtained with results very close to the experimental data;
Spin charge distribution, molecular electrostatic potentials, charge distribution, and orbital picture substantially contribute to the interpretation of experimental data and shed light on their atomistic mechanisms.
We believe that the performances highlighted herein are sufficiently positive to encourage future theoretical studies with the aim of deepening knowledge on electrochromism in organometallic systems with multiple oxidation states, ionic materials, and molecules inserted in a polymeric or liquid crystal matrix, and at the same time, to rationally design new molecules with particular electrochromic and chemical properties to be suggested for syntheses and future applications.

Supplementary Materials

The following are available online, Figure S1: Plots of the molecular orbitals involved in the transitions for 4 and 42+, Figure S2: Computed density difference plot for molecule 4 with the first excited state considered.

Author Contributions

B.C.D.S., M.E.A., T.M., N.R., and M.T. have equally contributed to the work. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The University of Calabria is gratefully acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Platt, J.R. Electrochromism, a possible change in color producible in dyes by an electric field. J. Chem. Phys. 1961, 34, 862–864. [Google Scholar] [CrossRef]
  2. Granqvist, C.G. Oxide electrochromics: An introduction to devices and materials. Sol. Energy Mater. Sol. Cells 2012, 99, 1–13. [Google Scholar] [CrossRef]
  3. Granqvist, C.-G. Electrochromic materials: Out of a niche. Nat. Mater. 2006, 5, 89−90. [Google Scholar] [CrossRef] [PubMed]
  4. Mortimer, R.J. Electrochromic materials. Chem. Soc. Rev. 1997, 26, 147–156. [Google Scholar] [CrossRef]
  5. Mortimer, R.J. Electrochromic Materials. Annu. Rev. Mater. Res. 2011, 41, 241–268. [Google Scholar] [CrossRef]
  6. Deb, S.K. Opportunities and challenges in science and technology of WO3 for electrochromic and related applications. Sol. Energy Mater. Sol. Cells 2008, 92, 245–258. [Google Scholar] [CrossRef]
  7. Cai, G.; Wang, J.; Lee, P.S. Next-Generation Multifunctional Electrochromic Devices. Acc. Chem. Res. 2016, 49, 469–1476. [Google Scholar] [CrossRef]
  8. Runnerstrom, E.L.; Llordés, A.; Lounis, S.D.; Milliron, D.J. Nanostructured electrochromic smart windows: Traditional materials and nir-selective plasmonic nanocrystals. Chem. Commun. 2014, 50, 10555–10572. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Granqvist, C.G. Electrochromics for smart windows: Oxide-based thin films and devices. Thin Solid Films 2014, 564, 1–38. [Google Scholar] [CrossRef]
  10. Lampert, C.M. Chromogenic smart materials. Mater. Today 2004, 7, 28–35. [Google Scholar] [CrossRef]
  11. Mortimer, R.J. Switching colors with electricity: Electrochromic materials can be used in glare reduction, energy conservation and chameleonic fabrics. Am. Sci. 2013, 101, 38–46. [Google Scholar] [CrossRef]
  12. Varaprasad, D.V.; Habibi, H.; McCabe, I.A.; Lynam, N.R. Electrochromic Mirrors and Devices. European Patent EP0758929B1, 6 August 1997. [Google Scholar]
  13. Mortimer, R.J.; Dyer, A.L.; Reynolds, J.R. Electrochromic organic and polymeric materials for display applications. Displays 2006, 27, 2–18. [Google Scholar] [CrossRef] [Green Version]
  14. Beneduci, A.; Cospito, S.; La Deda, M.; Veltri, L.; Chidichimo, G. Electrofluorochromism in π-conjugated ionic liquid crystals. Nat. Commun. 2014, 5, 3105. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Wadhwa, K.; Nuryyeva, S.; Fahrenbach, A.C.; Elhabiri, M.; Platas-Iglesias, C.; Trabolsi, A. Intramolecular redox-induced dimerization in a viologen dendrimer. J. Mater. Chem. C 2013, 1, 2302–2307. [Google Scholar] [CrossRef]
  16. Nguyen, Q.V.; Martin, P.; Frath, D.; Della Rocca, M.L.; Lafolet, F.; Bellinck, S.; Lafarge, P.; Lacroix, J.C. Highly Efficient Long-Range Electron Transport in a Viologen-Based Molecular Junction. J. Am. Chem. Soc. 2018, 140, 10131–10134. [Google Scholar] [CrossRef]
  17. Yu, W.J.; Liao, L.; Chae, S.H.; Lee, Y.H.; Duan, X.F. Toward Tunable Band Gap and Tunable Dirac Point in Bilayer Graphene with Molecular Doping. Nano Lett. 2011, 11, 4759–4763. [Google Scholar] [CrossRef] [Green Version]
  18. Wu, K.; Rodriguez-Cordoba, W.E.; Liu, Z.; Zhu, H.; Lian, T. Beyond Band Alignment: Hole Localization Driven Formation of Three Spatially Separated Long-Lived Exciton States in CdSe/CdS Nanorods. ACS Nano 2013, 7, 7173–7185. [Google Scholar] [CrossRef]
  19. Maufroy, A.; Favereau, L.; Anne, F.B.; Pellegrin, Y.; Blart, E.; Hissler, M.; Jacquemin, D.; Odobel, F. Diketopyrrolopyrrole-zinc porphyrin, a tuned panchromatic association for dye-sensitized solar cells. J. Mater. Chem. A 2013, 1, 7572–7575. [Google Scholar]
  20. Suzuka, M.; Hara, S.; Sekiguchi, T.; Oyaizu, K.; Nishide, H. Polyviologen as the charge-storage electrode of an aqueous electrolyte- and organic-based dye-sensitized solar cell. Polymer 2015, 68, 353–357. [Google Scholar] [CrossRef]
  21. Beaupre, S.; Breton, A.C.; Dumas, J.; Leclerc, M. Multicolored electrochromic cells base on poly(2,7-carbazole) derivatives for adaptive camouflage. Chem. Mater. 2009, 21, 1504–1513. [Google Scholar] [CrossRef]
  22. Sonmez, G.; Sonmez, H.B. Polymeric electrochromics for data storage. J. Mater. Chem. 2006, 16, 2473–2477. [Google Scholar] [CrossRef]
  23. Liu, G.; Chen, Y.; Gao, S.; Zhang, B.; Li, R.W.; Zhuang, X.D. Recent Advances in Resistive Switching Materials and Devices: From Memories to Memristors. Eng. Sci. 2018, 4, 4–43. [Google Scholar]
  24. Yang, P.; Sun, P.; Mai, W. Electrochromic energy storage devices. Mater. Today 2016, 19, 394–402. [Google Scholar] [CrossRef]
  25. Tong, Z.; Tian, Y.; Zhang, H.; Li, X.; Ji, J.; Qu, H. Recent advances in multifunctional electrochromic energy storage devices and photoelectrochromic devices. Sci. China Chem. 2017, 60, 13–37. [Google Scholar] [CrossRef]
  26. Cantekin, S.; Markvoort, A.J.; Elemans, J.A.A.W.; Rowan, A.E.; Nolte, R.J.M. Allosterically Controlled Threading of Polymers through Macrocyclic Dimers. J. Am. Chem. Soc. 2015, 137, 3915–3923. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Janoschka, T.; Martin, N.; Hager, M.D.; Schubert, U.S. An Aqueous Redox-Flow Battery with High Capacity and Power: The TEMPTMA/MV System. Angew. Chem. Int. Ed. 2016, 55, 14427–14430. [Google Scholar] [CrossRef] [PubMed]
  28. Hu, B.; DeBruler, C.; Rhodes, Z.; Liu, T.L. Long-Cycling Aqueous Organic Redox Flow Battery (AORFB) toward Sustainable and Safe Energy Storage. J. Am. Chem. Soc. 2017, 139, 1207–1214. [Google Scholar] [CrossRef]
  29. Luo, J.; Hu, B.; Debruler, C.; Liu, T.L. A π-Conjugation Extended Viologen as a Two-Electron Storage Anolyte for Total Organic Aqueous Redox Flow Batteries. Angew. Chem. Int. Ed. 2018, 57, 231–235. [Google Scholar] [CrossRef]
  30. De Bruler, C.; Hu, B.; Moss, J.; Luo, J.; Liu, T.L. A Sulfonate-Functionalized Viologen Enabling Neutral Cation Exchange, Aqueous Organic Redox Flow Batteries toward Renewable Energy Storage. ACS Energy Lett. 2018, 3, 663–668. [Google Scholar] [CrossRef]
  31. Wang, L.; Sun, S.; He, Y.; He, N.; Zhang, F.; Yao, Y.; Zhang, B.; Zhuang, X.; Chen, Y. Viologen-bridged polyaniline based multifunctional heterofilms for all-solid-state supercapacitors and memory devices. Eur. Polym. J. 2018, 98, 125–136. [Google Scholar] [CrossRef]
  32. Yang, X.; Liu, L.; Yang, W. Direct preparation of monodisperse core-shell microspheres with surface antibacterial property by using bicationic viologen surfmer. Polymer 2012, 53, 2190–2196. [Google Scholar] [CrossRef]
  33. Yuan, S.J.; Xu, F.J.; Kang, E.T.; Pehkonen, S.O. Modification of Surface-Oxidized Copper Alloy by Coupling of Viologens for Inhibiting Microbiologically Influenced Corrosion. J. Electrochem. Soc. 2007, 154, C645–C657. [Google Scholar] [CrossRef]
  34. Asaftei, S.; Huskens, D.; Schols, D. HIV-1 X4 Activities of Polycationic “Viologen” Based Dendrimers by Interaction with the Chemokine Receptor CXCR4: Study of Structure–Activity Relationship. J. Med. Chem. 2012, 55, 10405–10413. [Google Scholar] [CrossRef] [PubMed]
  35. Shah, K.W.; Wang, S.-X.; Soo, D.X.Y.; Xu, J. Viologen-based electrochromic materials: From small molecules, polymers and composites to their applications. Polymers 2019, 11, 1839. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Woodward, A.N.; Kolesar, J.M.; Hall, S.R.; Saleh, N.-A.; Jones, D.S.; Walter, M.G. Thiazolothiazole fluorophores exhibiting strong fluorescence and viologen-like reversible electrochromism. J. Am. Chem. Soc. 2017, 139, 8467–8473. [Google Scholar] [CrossRef]
  37. Wu, Y.; Zhou, J.; Phelan, B.T.; Mauck, C.M.; Stoddart, J.F.; Young, R.M.; Wasielewski, M.R. Probing distance dependent charge-transfer character in excimers of extended viologen cyclophanes using femtosecond vibrational spectroscopy. J. Am. Chem. Soc. 2017, 139, 14265–14276. [Google Scholar] [CrossRef]
  38. Ho, K.-C.; Lu, H.-C.; Yu, H.-F. Viologens-based electrochromic materials and devices. In Electrochromic Smart Materials—Fabrication and Applications; Xu, J., Chua, M.H., Shah, K.W., Eds.; Royal Society of Chemistry: Cambridge, UK, 2018; pp. 372–405. [Google Scholar]
  39. Bridges, C.R.; Borys, A.M.; Beland, V.A.; Gaffen, J.R.; Baumgartner, T. Phosphoryl- and phosphonium-bridged viologens as stable two- and three-electron acceptors for organic electrodes. Chem. Sci. 2020, 11, 10483–10487. [Google Scholar] [CrossRef] [PubMed]
  40. Ding, J.; Zheng, C.; Wang, L.; Lu, C.; Zhang, B.; Chen, Y.; Li, M.; Zhai, G.; Zhuang, X. Viologen-inspired functional materials: Synthetic strategies and applications. J. Mater. Chem. A 2019, 7, 23337–23360. [Google Scholar] [CrossRef]
  41. Schickedanz, K.; Radtke, J.; Bolte, M.; Lerner, H.-W.; Wagner, M. Facile Route to Quadruply Annulated Borepins. J. Am. Chem. Soc. 2017, 139, 2842–2851. [Google Scholar] [CrossRef] [Green Version]
  42. Messersmith, R.E.; Tovar, J.D. Assessment of the aromaticity of borepin rings by spectroscopic, crystallographic and computational methods: A historical overview. J. Phys. Org. Chem. 2015, 28, 378–387. [Google Scholar] [CrossRef]
  43. Messersmith, R.E.; Yadav, S.; Siegler, M.A.; Ottosson, H.; Tovar, J.D. Benzo [b] thiophene Fusion Enhances Local Borepin Aromaticity in Polycyclic Heteroaromatic Compounds. J. Org. Chem. 2017, 82, 13440–13448. [Google Scholar] [CrossRef] [PubMed]
  44. Dong, S.; Gopalakrishna, T.Y.; Han, Y.; Phan, H.; Tao, T.; Ni, Y.; Liu, G.; Chi, C. Extended Bis(anthraoxa)quinodimethanes with Nine and Ten Consecutively Fused Six-Membered Rings: Neutral Diradicaloids and Charged Diradical Dianions/Dications. J. Am. Chem. Soc. 2019, 141, 62–66. [Google Scholar] [CrossRef] [PubMed]
  45. Dong, S.; Herng, T.S.; Gopalakrishna, T.Y.; Phan, H.; Lim, Z.L.; Hu, P.; Webster, R.D.; Ding, J.; Chi, C. Extended Bis(benzothia)-quinodimethanes and Their Dications: From Singlet Diradicaloids to Isoelectronic Structures of Long Acenes. Angew. Chem. Int. Ed. 2016, 55, 9316. [Google Scholar] [CrossRef]
  46. Hunt, C.; Peterson, M.; Anderson, C.; Chang, T.; Wu, G.; Scheiner, S.; Ménard, G. Switchable Aromaticity in an Isostructural Mn Phthalocyanine Series Isolated in Five Separate Redox States. J. Am. Chem. Soc. 2019, 141, 2604–2613. [Google Scholar] [CrossRef] [Green Version]
  47. Wei, Y.X.; Zhou, J.L.; Zheng, J.M.; Xu, C.Y. Improved stability of electrochromic devices using Ti-doped V2O5 film. Electrochim. Acta 2015, 166, 277–284. [Google Scholar] [CrossRef]
  48. Wu, W.; Wang, M.; Ma, J.; Cao, Y.; Deng, Y. Electrochromic metal oxides: Recent progress and prospect. Adv. Electron. Mater. 2018, 4, 1800185. [Google Scholar] [CrossRef]
  49. Gong, C.; Yang, J.; Cao, L.; Gao, Y. Self-Doped Oligoaniline Electrochromic Devices: Fabrication and Effect of the Oligoaniline Molecular Architecture. ChemElectroChem 2017, 4, 521–532. [Google Scholar] [CrossRef]
  50. Cao, L.; Gong, C.; Yang, J. Conducting tetraaniline derivatives with fast switching time, enhanced contrast and coloration efficiency. Electrochim. Acta 2016, 192, 422–430. [Google Scholar] [CrossRef]
  51. Tsuneyasu, S.; Ichikawa, T.; Nakamura, K.; Kobayashi, N. Electrochemical Stability of Diphenylanthracene and Its Effect on Alternating-Current-Driven Blue-Ligh Electrochemiluminescence Properties. ChemElectroChem 2017, 4, 1731–1735. [Google Scholar] [CrossRef] [Green Version]
  52. Ji, L.; Friedrich, A.; Krummenacher, I.; Eichhorn, A.; Braunschweig, H.; Moos, M.; Hahn, S.; Geyer, F.L.; Tverskoy, O.; Han, J.; et al. Preparation, Properties, and Structures of the Radical Anions and Dianions of Azapentacenes. J. Am. Chem. Soc. 2017, 139, 15968–15976. [Google Scholar] [CrossRef]
  53. Asai, K.; Fukazawa, A.; Yamaguchi, S. A Near-Infrared Dye That Undergoes Multiple Interconversions through Acid–Base Equilibrium and Reversible Redox Processes. Angew. Chem. Int. Ed. 2017, 56, 6848–6852. [Google Scholar] [CrossRef]
  54. Hsiao, S.H.; Chen, Y.Z. Electroactive and ambipolar electrochromic polyimides from arylene diimides with triphenylamine N-substituents. Dye. Pigment. 2017, 144, 173–183. [Google Scholar] [CrossRef]
  55. Taylor, J.O.; Pižl, M.; Kloz, M.; Rebarz, M.; McCusker, C.E.; McCusker, J.K.; Záliš, S.; Hartl, F.; Vlček, A. Optical and Infrared Spectroelectrochemical Studies of CNSubstituted Bipyridyl Complexes of Ruthenium (II). Inorg. Chem. 2021, 60, 3514–3523. [Google Scholar] [CrossRef]
  56. Waldie, K.M.; Ramakrishnan, S.; Kim, S.-K.; Maclaren, J.K.; Chidsey, C.E.D.; Waymouth, R.M. Multielectron Transfer at Cobalt: Influence of the Phenylazopyridine Ligand. J. Am. Chem. Soc. 2017, 139, 4540–4550. [Google Scholar] [CrossRef] [PubMed]
  57. De Simone, B.C.; Marino, T.; Russo, N. TDDFT investigation on methylviologen, 3,7-diazabenzophosphole, and helical helquat electrochromic systems. Theor. Chem. Acc. 2016, 135, 118. [Google Scholar] [CrossRef]
  58. De Simone, B.C.; Mazzone, G.; Marino, T.; Russo, N.; Toscano, M. Time-Dependent Density Functional Computations of the Spectrochemical Properties of Dithiolodithiole and Thiophene Electrochromic Systems. Materials 2017, 10, 981. [Google Scholar] [CrossRef] [Green Version]
  59. De Simone, B.C.; Marino, T.; Prejanò, M.; Russo, N. Can fused thiophene−pyrrole-containing rings act as possible new electrochromic dyes? A computational prediction. Theor. Chem. Acc. 2016, 135, 238. [Google Scholar] [CrossRef]
  60. De Simone, B.C.; Mazzone, G.; Marino, M.; Russo, N.; Toscano, M. On the Electrochromic Properties of Borepins: A Computational Prediction. ACS Omega 2018, 3, 9556–9563. [Google Scholar] [CrossRef]
  61. Wheeler, D.L.; Rainwater, L.E.; Green, A.R.; Tomlinson, A.L. Modeling electrochromic poly-dioxythiophene-containing materials through TDDFT. Phys. Chem. Chem. Phys. 2017, 19, 20251–20258. [Google Scholar] [CrossRef]
  62. Alberto, M.E.; De Simone, B.C.; Marino, T.; Russo, N.; Toscano, M. Photophysical properties of methyl ketone based multi-responsive electrochromic materials: A theoretical investigation. J. Mol. Liq. 2021, 338, 116576. [Google Scholar] [CrossRef]
  63. Cappello, C.; Therien, D.A.B.; Staroverov, V.N.; Lagugn-Labarthet, F.; Gilroy, J.B. Optoelectronic, Aggregation, and Redox Properties of Double-Rotor Boron Difluoride Hydrazone Dyes. Chem. Eur. J. 2019, 25, 5994–6006. [Google Scholar] [CrossRef] [Green Version]
  64. Okamoto, Y.; Tanioka, M.; Muranaka, A.; Miyamoto, K.; Aoyama, T.; Ouyang, X.; Kamino, S.; Sawada, D.; Uchiyama, M. Stable Thiele’s Hydrocarbon Derivatives Exhibiting Near-Infrared Absorption/Emission and Two-Step Electrochromism. J. Am. Chem. Soc. 2018, 140, 17857–17861. [Google Scholar] [CrossRef]
  65. Higashino, T.; Kumagai, A.; Sakaki, S.; Imahori, H. Reversible p-system switching of thiophene-fused thiahexaphyrins by solvent and oxidation/reduction. Chem. Sci. 2018, 9, 7528. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Ye, Q.; Chang, J.; Shi, X.; Dai, G.; Zhang, W.; Huang, K.-W.; Chi, C. Stable 7,14-Disubstituted-5,12-Dithiapentacenes with Quinoidal Conjugation. Org. Lett. 2014, 16, 3966–3969. [Google Scholar] [CrossRef]
  67. Quinton, C.; Alain-Rizzo, V.; Dumas-Verdes, C.; Miomandre, F.; Clavier, G.; Audebert, P. Redox- and Protonation-Induced Fluorescence Switch in a New Triphenylamine with Six Stable Active or Non-Active Forms. Chem. Eur. J. 2015, 21, 2230–2240. [Google Scholar] [CrossRef] [PubMed]
  68. Alberto, M.E.; De Simone, B.C.; Cospito, S.; Imbardelli, D.; Veltri, L.; Chidichimo, G.; Russo, N. Experimental and theoretical characterization of a new synthesized extended viologen. Chem. Phys. Lett. 2012, 552, 141–145. [Google Scholar] [CrossRef]
  69. De Simone, B.C.; Quartarolo, A.D.; Cospito, S.; Veltri, L.; Chidichimo, G.; Russo, N. Theoretical and experimental investigation on the near-infrared and UV−Vis spectral regions of a newly synthesized triarylamine electrochromic system. Theor. Chem. Acc. 2012, 131, 1225–1233. [Google Scholar] [CrossRef]
  70. Maltese, V.; Cospito, S.; Beneduci, A.; De Simone, B.C.; Russo, N.; Chidichimo, G.; Janssen, R.A.J. Electro-optical Properties of Neutral and Radical Ion Thienosquaraines. Chem. Eur. J. 2016, 22, 10179–10186. [Google Scholar] [CrossRef]
  71. Schipper, D.J.; Moh, L.C.H.; Muller, P.; Swager, T.M. Dithiolodithiole as a Building Block for Conjugated Materials. Angew. Chem. Int. Ed. 2014, 53, 5847–5851. [Google Scholar] [CrossRef]
  72. Bockman, T.M.; Kochi, J.K. Isolation and oxidation-reduction of methylviologen cation radicals. Novel disproportionation in charge-transfer salts by x-ray crystallography. J. Org. Chem. 1990, 55, 4127–4135. [Google Scholar] [CrossRef]
  73. Born, M. Die Elektronenaffinität der Halogenatome. Verhandl. Deut. Physik. Ges. 1919, 21, 679–685. [Google Scholar]
  74. Durben, S.; Baumgartner, T. 3,7-Diazadibenzophosphole Oxide: A Phosphorus-Bridged Viologen Analogue with Significantly Lowered Reduction Threshold. Angew. Chem. Int. Ed. 2011, 50, 7948–7952. [Google Scholar] [CrossRef]
  75. Pospíšil, L.; Bednárová, L.; Štepánek, P.; Slavícek, P.; Vávra, J.; Hromadová, M.; Dlouhá, H.; Tarábek, J.; Teplý, F. Intense Chiroptical Switching in a Dicationic Helicene-Like Derivative: Exploration of a Viologen-Type Redox Manifold of a Non-Racemic Helquat. J. Am. Chem. Soc. 2014, 136, 10826–10829. [Google Scholar] [CrossRef] [PubMed]
  76. Becke, A.D. Density-Functional Thermochemistry. III. The Role of Exact Exchange. J. Chem. Phys. 1993, 98, 5648–5652. [Google Scholar] [CrossRef] [Green Version]
  77. Lee, C.; Yang, W.; Parr, R.G. Development of the Colle-Salvetti Correlation-Energy Formula into a Functional of the Electron Density. Phys. Rev. B Condens. Matter Mater. Phys. 1988, 37, 785–789. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Zhao, Y.; Truhlar, D.G. The M06 Suite of Density Functionals for Main Group Thermochemistry, Thermochemical Kinetics, Noncovalent Interactions, Excited States, and Transition Elements: Two New Functionals and Systematic Testing of Four M06-Class Functionals and 12 Other Functionals. Theor. Chem. Acc. 2008, 120, 215–241. [Google Scholar]
  79. Chai, J.-D.; Head-Gordon, M. Long-range corrected hybrid density functionals with damped atom–atom dispersion corrections. Phys. Chem. Chem. Phys. 2008, 10, 6615–6620. [Google Scholar] [CrossRef] [Green Version]
  80. Chai, J.D.; Head-Gordon, M. Systematic optimization of long-range corrected hybrid density functionals. J. Chem. Phys. 2008, 128, 84106–84155. [Google Scholar] [CrossRef] [PubMed]
  81. Cossi, M.; Barone, V. Solvent effect on vertical electronic transitions by the polarizable continuum model. J. Chem. Phys. 2000, 112, 2427–2435. [Google Scholar] [CrossRef]
  82. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; et al. Gaussian 09; Revision, D.01; Gaussian, Inc.: Wallingford, CT, USA, 2009. [Google Scholar]
Scheme 1. Schematic representation of a generic electrochromic process.
Scheme 1. Schematic representation of a generic electrochromic process.
Molecules 26 05793 sch001
Scheme 2. Structures of the considered systems. 1: Quadruply benzannulated nonplanar borepin; 2: 4-(6-(4-Metylphenyl)-[1,2]dithiolo[4,3-c][1,2]dithiol-3-yl)benzonitrile; 3: Metylviologen; 4: Squaraine; 5: 5mesityl-5H-diacetonafto[5,6-bc:5′, 6′ef]borpin; 6: dibenzo[b]thiophene-fused borepin; 7: 3,7-diazabenzophosphole.
Scheme 2. Structures of the considered systems. 1: Quadruply benzannulated nonplanar borepin; 2: 4-(6-(4-Metylphenyl)-[1,2]dithiolo[4,3-c][1,2]dithiol-3-yl)benzonitrile; 3: Metylviologen; 4: Squaraine; 5: 5mesityl-5H-diacetonafto[5,6-bc:5′, 6′ef]borpin; 6: dibenzo[b]thiophene-fused borepin; 7: 3,7-diazabenzophosphole.
Molecules 26 05793 sch002
Figure 1. B3LYP/6–31+G* bond length variations during the oxidation process of 3.
Figure 1. B3LYP/6–31+G* bond length variations during the oxidation process of 3.
Molecules 26 05793 g001
Figure 2. Electron spin distribution of 3+ computed at the B3LYP/6–31+G* level of theory.
Figure 2. Electron spin distribution of 3+ computed at the B3LYP/6–31+G* level of theory.
Molecules 26 05793 g002
Figure 3. B3LYP/6–31+G* molecular electrostatic potential for the different oxidation states of 3 species.
Figure 3. B3LYP/6–31+G* molecular electrostatic potential for the different oxidation states of 3 species.
Molecules 26 05793 g003
Figure 4. Frontiers molecular orbital values and orbital shape for the different oxidation states of 3.
Figure 4. Frontiers molecular orbital values and orbital shape for the different oxidation states of 3.
Molecules 26 05793 g004
Table 1. (a) Main geometrical parameters for molecules 1 and 2, computed by using different exchange-correlation functionals by using the 6- 31+G* basis sets. The optimizations have been performed employing the Polarizable Continuum Model considering the chloroform and dichloromethane solvent for 1 and 2, respectively. (b) Numeration of the considered parameters.
Table 1. (a) Main geometrical parameters for molecules 1 and 2, computed by using different exchange-correlation functionals by using the 6- 31+G* basis sets. The optimizations have been performed employing the Polarizable Continuum Model considering the chloroform and dichloromethane solvent for 1 and 2, respectively. (b) Numeration of the considered parameters.
a)X-RayB3LYPM06PBE0
1
B-C11.5631.5641.5591.560
B-C191.5951.5951.5851.590
C9-C101.5021.4971.4901.490
C1BC18127.3127.1127.3127.3
C1BC19116.4116.4116.3116.3
C8C9C10C1138.939.239.239.0
2
C1-C21.3651.3651.3611.363
C2-C31.4441.4471.4411.441
C2-S31.7571.7831.7711.764
S3-S42.0842.1282.1192.099
C1-C2-S3126.5125.7126.0126.3
C3-C2-S3113.4113.4113.6113.5
C2-S3-S495.895.195.195.4
C2-C3-C4120.1120.7120.4120.1
φ1154.6140.6142.2142.9
φ226.1−34.7−33.6−32.7
Molecules 26 05793 i001
Table 2. Excitation energies (nm) and oscillator strength of three and four molecules in different oxidation states computed in dichlorometane solvent by using different exchange-correlation potentials and the 6–31+G* basis set. The % of the transition is that obtained at the B3LYP level. Experimental data are taken from refs [70,72].
Table 2. Excitation energies (nm) and oscillator strength of three and four molecules in different oxidation states computed in dichlorometane solvent by using different exchange-correlation potentials and the 6–31+G* basis set. The % of the transition is that obtained at the B3LYP level. Experimental data are taken from refs [70,72].
SystemB3LYPM06wB97XDTransition (%)λexp
32+368 (0.23)366 (0.24)313 (0.30)H→L (70)337 [72]
3+399 (0.10)
433 (0.15)
1152 (0.11)
399 (0.13)
440 (0.12)
1169 (0.12)
355 (0.16)
403 (0.15)
1213 (0.15)
βH-1→L (84)
βH→L (75)
H→L (99)
400 [72], 447 [72]
1050–1100 [72]
3553 (0.1)
721 (0.2)
557 (0.1)
698 (0.2)
455 (0.0)
599 (0.3)
H→L+3 (57)
H→L (68)
558 [72], 712 [72]
42+513 (1.39)504 (1.49)450 (1.74)H→L (69)450 [70]
4+575 (1.61)
842 (0.00)
581 (1.62)
825 (0.00)
565 (1.64)
733 (0.05)
βH→L (71)
αH→L (74)
570 [70]
840 [70]
4572 (1.67)575 (1.64)585 (1.70)H→L (71)660 [70]
4494 (1.30)
709 (0.07)
487 (1.12)
690 (0.08)
444 (1.22)
643 (0.14)
αH→L (83)
βH→L (84)
660 [70]
Table 3. Reduction potentials (V) for studied compounds computed by using the M06/6–31+G* level of theory. The used theoretical and experimental solvents are 5 and 6 in tetrahydrofuran; 3 and 7 in dichloromethane. a from refs [41,72,74,75].
Table 3. Reduction potentials (V) for studied compounds computed by using the M06/6–31+G* level of theory. The used theoretical and experimental solvents are 5 and 6 in tetrahydrofuran; 3 and 7 in dichloromethane. a from refs [41,72,74,75].
MoleculePotentialTheory (M06)Experimental aΔtheo-exp
5EIred−1.60−1.490.11
EIIred−2.00−1.840.16
6EIred−2.24−2.140.10
3EIred−1.18−1.090.09
EIIred−1.42−1.520.00
7EIred−0.74−0.510.23
EIIred−1.23−1.000.23
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

De Simone, B.C.; Alberto, M.E.; Marino, T.; Russo, N.; Toscano, M. The Contribution of Density Functional Theory to the Atomistic Knowledge of Electrochromic Processes. Molecules 2021, 26, 5793. https://doi.org/10.3390/molecules26195793

AMA Style

De Simone BC, Alberto ME, Marino T, Russo N, Toscano M. The Contribution of Density Functional Theory to the Atomistic Knowledge of Electrochromic Processes. Molecules. 2021; 26(19):5793. https://doi.org/10.3390/molecules26195793

Chicago/Turabian Style

De Simone, Bruna Clara, Marta Erminia Alberto, Tiziana Marino, Nino Russo, and Marirosa Toscano. 2021. "The Contribution of Density Functional Theory to the Atomistic Knowledge of Electrochromic Processes" Molecules 26, no. 19: 5793. https://doi.org/10.3390/molecules26195793

APA Style

De Simone, B. C., Alberto, M. E., Marino, T., Russo, N., & Toscano, M. (2021). The Contribution of Density Functional Theory to the Atomistic Knowledge of Electrochromic Processes. Molecules, 26(19), 5793. https://doi.org/10.3390/molecules26195793

Article Metrics

Back to TopTop