Density Functional Theory-Based Molecular Modeling: Verification of Decisive Roles of Van der Waals Aggregation of Triiodide Ions for Effective Electron Transfer in Wet-Type N3-Dye-Sensitized Solar Cells

Density functional theory-based molecular modeling (DFT/MM) validates that KI and I2 undergo exothermic van der Waals (vdW) aggregation in acetonitrile (AN) or in the presence of 4-tert-butylpyridine (TBP), forming potassium triiodide (KI3) and, further mutual vdW aggregation leads to the formation of (KI3)2 and AN, (KI3)2 and (AN)2 and (KI3)2 and TBP in the AN-based Dye sensitized solar cells (DSSC) electrolytes. All KI3 aggregates have a very low energy gap, 0.17 eV, 0.14 eV and 0.05 eV of lowest unoccupied molecular orbital (LUMO) + 1 and LUMO, respectively, verifying efficient electron diffusion in μm-thick DSSC electrolytes. Hydrogen-bonding aggregation of anatase TiO2 model, Ti9O18H and OH, with N3 (proton) dye is also validated by DFT/MM, and the energy structure verifies unidirectional electron flow from highest occupied molecular orbital (HOMO) on thiocyanide (SCN) groups of N3 dye to LUMO on the TiO2 model at the aggregates. Further, DFT/MM for the aggregation of K+I3− with N3 verifies the most exothermic formation of the aggregate of N3 (proton) and K+I3−. The UV-Vis spectra of N3 (proton) and K+I3− is consistent with reported incident photocurrent efficiency (IPCE) action spectra (λ = 450–800 nm) of N3-sensitized DSSC, verifying that the N3 dye of N3 (proton) and K+I3− becomes an effective sensitizer in the anode / TiO2 / N3 (proton) / KI/I2 / acetonitrile (AN) / cathode structured DSSC. The energy structure of LUMO and LUMO + 1 of the aggregates, Ti9O18H and OH and N3 (proton), N3 and K+I3−, (KI3)2 and AN and (KI3)2 and TBP verifies high IPCE photocurrent and effective electron diffusion via KI3-aggregates in the DSSC of Ti9O18H and OH and N3 (proton) and K+I3−.


Introduction
As reported in the preceding papers [1][2][3][4][5][6], quantum chemistry molecular modeling, i.e., density functional theory-based molecular modeling (DFT/MM) can be regarded as theory-based "experiments". DFT/MM can be carried out very quickly using high-end supercomputer-like personal computers. DFT/MM are applicable to molecular aggregates which are induced by van der Waals (vdW) force (including hydrogen bonding and Coulomb interactions). DFT/MM for the structuring of aggregates gives the heat of formation, dipole moment for understanding solubility in polar solvent like water and acetonitrile (AN) and the energy structures of aggregates, i.e., surface orbital electron energy structures [almost degenerate highest molecular orbitals [HOMO (n), n = 0-4)] and the lowest unoccupied molecular orbitals [LUMO + n, n = 0-2)]. Their molecular orbital configurations are evaluated as theory-based "experimental" data for verification and prediction of interacting molecules. Semiconducting property of the aggregates can be verified by the orbital energy gap between LUMO and LUMO + 1 [2,3]. We performed DFT/MM using the B3LYP exchange correlation functional and the 6-31G(d) basis set with Spartan'18 (Wavefunction, Inc. Irvine, CA). Prior to the DFT calculation for electronic properties, the geometries of the molecular aggregates were determined by the molecular mechanics employing the Merck molecular force field (MMFF).
Since the first report of its high-performance by O'Regan and Grätzel [7], the dye-sensitized solar cell (DSSC) using KI / I 2 / acetonitrile electrolyte has been a focus for a fundamental understanding of microscopic electronic phenomena of TiO 2 , dye molecules and electrolytes [8]. As for wet-type DSSC electrolytes of acetonitrile (AN), KI and I 2 , redox couple of I −/ I 3 − was introduced to explain electron diffusion by electron hopping and exchange mechanisms [9]. Recently, we reported that DFT/MM verifies the semiconducting properties for quasi-solid state I 5 − N(CH 3 ) 4 + -based electrolytes of DSSC [1].
In addition, we clarified that PbI 6 4− /MeNH 3 + -based perovskite solar cells have excellent electron diffusion among panchromatic PbI 6 4− /MeNH 3 + matrixes in the solar cells [2,3]. We have ascribed excellence of perovskite solar cells to PbI 6 4− /MeNH 3 + -induced vdW aggregation. Here, we corroborate that the iodine element-based vdW aggregations also play an important role in unidirectional electron diffusion in wet-type DSSC using KI/I 2 /acetonitrile electrolytes. Figure 1 shows the EPM of I 2 , KI, and AN.

DFT-Based Electrostatic Potential Map (EPM) and Van der Waals Aggregation
Energies 2020, 13, x FOR PEER REVIEW 2 of 15 orbital electron energy structures [almost degenerate highest molecular orbitals [HOMO (n), n = 0-4)] and the lowest unoccupied molecular orbitals [LUMO + n, n = 0-2)]. Their molecular orbital configurations are evaluated as theory-based "experimental" data for verification and prediction of interacting molecules. Semiconducting property of the aggregates can be verified by the orbital energy gap between LUMO and LUMO + 1 [2,3]. We performed DFT/MM using the B3LYP exchange correlation functional and the 6-31G(d) basis set with Spartan'18 (Wavefunction, Inc. Irvine, CA). Prior to the DFT calculation for electronic properties, the geometries of the molecular aggregates were determined by the molecular mechanics employing the Merck molecular force field (MMFF).
Since the first report of its high-performance by O'Regan and Grätzel [7], the dye-sensitized solar cell (DSSC) using KI / I2 / acetonitrile electrolyte has been a focus for a fundamental understanding of microscopic electronic phenomena of TiO2, dye molecules and electrolytes [8]. As for wet-type DSSC electrolytes of acetonitrile (AN), KI and I2, redox couple of I −/ I3 − was introduced to explain electron diffusion by electron hopping and exchange mechanisms [9]. Recently, we reported that DFT/MM verifies the semiconducting properties for quasi-solid state I5 − N(CH3)4 + -based electrolytes of DSSC [1]. In addition, we clarified that PbI6 4− /MeNH3 + -based perovskite solar cells have excellent electron diffusion among panchromatic PbI6 4− /MeNH3 + matrixes in the solar cells [2,3]. We have ascribed excellence of perovskite solar cells to PbI6 4− /MeNH3 + -induced vdW aggregation. Here, we corroborate that the iodine element-based vdW aggregations also play an important role in unidirectional electron diffusion in wet-type DSSC using KI / I2 / acetonitrile electrolytes. In the colorful electrostatic potential map (ESPM), the property range (kJ/mol, eV) indicates potential energy of surface electron density of molecules and of vdW aggregates. Blue surface color is particularly electron poor and subject to electrophilic attack (electron poor vdW surface), and red surface color is particularly electron rich, subject to nucleophilic attack (electron rich vdW surface). As quantitative structured activity, volume (Å 3 ) of molecules is determined by the CPK model, reflecting surface electron number of molecules and aggregates. It may be added that the blue region corresponds to the LUMO rich region, and the red region to the HOMO rich region.

DFT-Based Electrostatic Potential Map (EPM) and Van der Waals Aggregation
Pale blue EPM of I2 indicates its electrophilic nature and the top surface is particularly electrophilic. Colorful EPM of KI and AN indicate the presence of both electrophilic and nucleophilic In the colorful electrostatic potential map (ESPM), the property range (kJ/mol, eV) indicates potential energy of surface electron density of molecules and of vdW aggregates. Blue surface color is particularly electron poor and subject to electrophilic attack (electron poor vdW surface), and red surface color is particularly electron rich, subject to nucleophilic attack (electron rich vdW surface). As quantitative structured activity, volume (Å 3 ) of molecules is determined by the CPK model, reflecting surface electron number of molecules and aggregates. It may be added that the blue region corresponds to the LUMO rich region, and the red region to the HOMO rich region.
Pale blue EPM of I 2 indicates its electrophilic nature and the top surface is particularly electrophilic. Colorful EPM of KI and AN indicate the presence of both electrophilic and nucleophilic surface, reflecting large dipole, 10.68 debye and 3.81 debye, respectively. It is worth noting that the volumes of I 2 and Energies 2020, 13, 3027 3 of 14 KI are much larger than that of AN. Van der Waals aggregation occurs between electrophilic electron surface and nucleophilic one.

Verification of KI 3 Formation by Van der Waals Aggregation of I 2 and KI in Electrolytes
We speculate that the minimum number of molecules which may form vdW aggregates in AN should be three, and one of them is AN. Figure 2 shows ESPM of equilibrium geometry of vdW aggregates of (I 2 ) 2 AN, (KI) 2 AN, and acetonitrile (AN) 3 . When compared to ESPM of I 2 , and KI and AN in Figure 1, all of trimeric aggregates result from vdW aggregation between the electron-poor electrophilic surface and the electro-rich nucleophilic surface.
Energies 2020, 13, x FOR PEER REVIEW 3 of 15 surface, reflecting large dipole, 10.68 debye and 3.81 debye, respectively. It is worth noting that the volumes of I2 and KI are much larger than that of AN. Van der Waals aggregation occurs between electrophilic electron surface and nucleophilic one.

Verification of KI3 Formation by Van der Waals Aggregation of I2 and KI in Electrolytes.
We speculate that the minimum number of molecules which may form vdW aggregates in AN should be three, and one of them is AN. Figure 2 shows ESPM of equilibrium geometry of vdW aggregates of (I2)2AN, (KI)2AN, and acetonitrile (AN)3. When compared to ESPM of I2, and KI and AN in Figure 1, all of trimeric aggregates result from vdW aggregation between the electron-poor electrophilic surface and the electro-rich nucleophilic surface. All vdW aggregates form exothermically, suggesting stability in AN. The most stable vdW aggregate of (KI)2AN (heat of formation is −49.6 kcal/mol) gives a colorful ESPM, predicting further vdW aggregation with KI and I2 in AN. Figure 3 shows ESPM of I2 and KI and AN, KI3 and I2 and KI, and KI3 and (AN)2. All vdW aggregates give large heat of formation and, interestingly, I2 aggregates with KI in aggregate of I2 and KI and AN. Geometry of I3 − gives angle of 172° comparable to that of geometry I3 − in KI3 and (AN)2. DFT/MM-determined UV-Vis spectrum comparison (max: 447 nm and 414 nm) verifies in situ formation of KI3 from I2 and KI in AN. In situ formation of KI3 from I2 and KI is also confirmed in vdW aggregate of KI3 and I2 and KI. All vdW aggregates form exothermically, suggesting stability in AN. The most stable vdW aggregate of (KI) 2 AN (heat of formation is −49.6 kcal/mol) gives a colorful ESPM, predicting further vdW aggregation with KI and I 2 in AN. Figure 3 shows ESPM of I 2 and KI and AN, KI 3

Van der Waals Aggregation of Potassium Triiodide (KI3) and 4-Tertbyutylpyridine (TB) Effect on DSSC Electrolytes.
Interestingly, DFT/MM of vdW aggregate of KI3 and I2 and KI reveals that the energy gap of ELUMO + 1 and ELUMO is 0.19 eV, predicting that the aggregate is semiconducting once it is formed. With this prediction in mind, DFT/MM is extended to vdW aggregates of (KI3)2 and AN, and of (KI3)2 and (AN)2 to clarify the reason why the electrolyte of KI /I2 /AN become semiconducting. In addition, since it is well-known that presence of 4-tert-butylpyridine (TBP) in the electrolytes improves open circuit voltage (Voc) of DSSC. Open circuit voltage depends on electron leak-free unidirectional photocurrent [8,10]. DFT/MM is also extended to vdW aggregate of (KI3)2 and TBP. Figure 4 shows the EPSM of equilibrium geometry of vdW aggregates of (KI3)2 and AN, (KI3)2 and (AN)2, (KI3)2 and TBP. The heat of formation values verify exclusive formation of vdW dimer of KI3 in wet type DSSC electrolytes. UV-Vis spectra support the orange color of the electrolytes. Energy gaps of LUMO + 1 and LUMO are all less than 0.2 eV, and their configurations explain unidirectional electron flow from LUMO + 1 to LUMO. Interestingly, the vdW aggregate of (KI3)2 and TBP gives the smallest energy gap of 0.05 eV.
In order to validate conductivity of electron-running state of the vdW aggregates of KI3 -based electrolytes, EPM and energy structures of electron injected [(KI3)2 and AN] .− and [(KI3)2 and TBP] .− are examined ( Figure 5). Reddish ESPM indicates nucleophilic surface state, both negative values of heat of their electron acceptance verify their electron acceptability, and hopping between them is favorable in the electrolytes. As for surface energy structures, Ea (HOMO) is identical with electron singly occupied molecular orbital (SOMO). Accepted electron (spin density) is located in bLUMO. We understand that spin may run on bLUMO at electron running process. The energy gap of bLUMO and SOMO (0.46 eV and 0.35 eV, respectively), is small enough for spin (unpaired electron) to run under irradiance of solar energy (h).

Van der Waals Aggregation of Potassium Triiodide (KI 3 ) and 4-Tertbyutylpyridine (TB) Effect on DSSC Electrolytes
Interestingly, DFT/MM of vdW aggregate of KI 3 and I 2 and KI reveals that the energy gap of ELUMO + 1 and ELUMO is 0.19 eV, predicting that the aggregate is semiconducting once it is formed. With this prediction in mind, DFT/MM is extended to vdW aggregates of (KI 3 ) 2 and AN, and of (KI 3 ) 2 and (AN) 2 to clarify the reason why the electrolyte of KI /I 2 /AN become semiconducting. In addition, since it is well-known that presence of 4-tert-butylpyridine (TBP) in the electrolytes improves open circuit voltage (Voc) of DSSC. Open circuit voltage depends on electron leak-free unidirectional photocurrent [8,10]. DFT/MM is also extended to vdW aggregate of (KI 3 ) 2 and TBP. Figure 4 shows the EPSM of equilibrium geometry of vdW aggregates of (KI 3 ) 2 and AN, (KI 3 ) 2 and (AN) 2 , (KI 3 ) 2 and TBP. The heat of formation values verify exclusive formation of vdW dimer of KI 3 in wet type DSSC electrolytes. UV-Vis spectra support the orange color of the electrolytes. Energy gaps of LUMO + 1 and LUMO are all less than 0.2 eV, and their configurations explain unidirectional electron flow from LUMO + 1 to LUMO. Interestingly, the vdW aggregate of (KI 3 ) 2 and TBP gives the smallest energy gap of 0.05 eV.
In order to validate conductivity of electron-running state of the vdW aggregates of KI 3 -based electrolytes, EPM and energy structures of electron injected [(KI 3 ) 2 and AN] .− and [(KI 3 ) 2 and TBP] .− are examined ( Figure 5). Reddish ESPM indicates nucleophilic surface state, both negative values of heat of their electron acceptance verify their electron acceptability, and hopping between them is favorable in the electrolytes. As for surface energy structures, Ea (HOMO) is identical with electron singly occupied molecular orbital (SOMO). Accepted electron (spin density) is located in bLUMO. We understand that spin may run on bLUMO at electron running process. The energy gap of bLUMO and SOMO (0.46 eV and 0.35 eV, respectively), is small enough for spin (unpaired electron) to run under irradiance of solar energy (hν).

Electrostatic Potential Map (ESPM) of TiO2
Anatase Nanoparticle and Conductivity. Figure 6 shows the ESPM of anatase TiO2 model of Ti9O18H and OH proposed by Jono and Yamashita (Jono-Yamashita model) [11]. The validity of the Jono-Yamashita model for TiO2 nanoparticle (NP) has been confirmed elsewhere [12]. Colorful ESPM verifies that three-dimensional vdW aggregates may form TiO2 thin films and result in giving porous structures when sintered at high temperature.
To understand the nature of electron acceptability and electron-transfer in TiO2, electron running molecular orbital configurations of LUMO and LUMO + 1 are shown as the side view (a) and the view from the top (b) of the TiO2 model. The energy gap of LUMO + 1 and LUMO is 0.43 eV, which verifies and predicts that TiO2 model and their vdW aggregates become semiconducting when LUMO and LUMO + 1 accept photoelectron.
(a)  Figure 6 shows the ESPM of anatase TiO 2 model of Ti 9 O 18 H and OH proposed by Jono and Yamashita (Jono-Yamashita model) [11]. The validity of the Jono-Yamashita model for TiO 2 nanoparticle (NP) has been confirmed elsewhere [12]. Colorful ESPM verifies that three-dimensional vdW aggregates may form TiO 2 thin films and result in giving porous structures when sintered at high temperature.

Electrostatic Potential Map (ESPM) of TiO2
Anatase Nanoparticle and Conductivity. Figure 6 shows the ESPM of anatase TiO2 model of Ti9O18H and OH proposed by Jono and Yamashita (Jono-Yamashita model) [11]. The validity of the Jono-Yamashita model for TiO2 nanoparticle (NP) has been confirmed elsewhere [12]. Colorful ESPM verifies that three-dimensional vdW aggregates may form TiO2 thin films and result in giving porous structures when sintered at high temperature.
To understand the nature of electron acceptability and electron-transfer in TiO2, electron running molecular orbital configurations of LUMO and LUMO + 1 are shown as the side view (a) and the view from the top (b) of the TiO2 model. The energy gap of LUMO + 1 and LUMO is 0.43 eV, which verifies and predicts that TiO2 model and their vdW aggregates become semiconducting when LUMO and LUMO + 1 accept photoelectron. (a)

Electrostatic Potential Map (ESPM) of Sensitizing N3 Dyes and Perspective of Van der Waals
Aggregation. Figure 7 shows ESPM of two kinds of N3 (tetramethylammonium) (a homolog of N719) and N3 (proton) with UV-Vis spectra and electron energy structures. Colorful ESPM suggests that they are likely not to aggregate each other, but orange-colored SCN group sites may be vdW-aggregated by electrophilic potassium triiodide (K + I3 − ). Each energy gap of LUMO + 1 and LUMO, 0.36 eV and 0.10 eV predicts that both N3 dyes are conductive.
Recently we reported that DFT/MM predict that DFT/MM-determined UV-Vis absorption spectra are strongly affected by vdW aggregation with other molecules, which was exemplified by UV-Vis spectrum of benzene solution [3]. Figure 7 reveals that DFT/MM-based UV-Vis spectrum of N3 (tetramethylammonium) is quite comparable to the incident photocurrent efficiency (IPCE) action spectra of N3-sensitized DSSC [13]. However, UV-Vis spectrum of N3 (proton) gives 380 nm red shifted one with decreased strength of absorption. DFT/MM-determined three-dimensional molecular structures verify that the difference comes from a change of coordination structure of ligands, bipyridyl and SCN ligands around Ru(II) ion (see Section 3).     Recently we reported that DFT/MM predict that DFT/MM-determined UV-Vis absorption spectra are strongly affected by vdW aggregation with other molecules, which was exemplified by UV-Vis spectrum of benzene solution [3]. (tetramethylammonium) is quite comparable to the incident photocurrent efficiency (IPCE) action spectra of N3-sensitized DSSC [13]. However, UV-Vis spectrum of N3 (proton) gives 380 nm red shifted one with decreased strength of absorption. DFT/MM-determined three-dimensional molecular structures verify that the difference comes from a change of coordination structure of ligands, bipyridyl and SCN ligands around Ru(II) ion (see Section 3).

Verification of Van der Waals Aggregation of N3 Dye with Anatase TiO 2 Model
Here, let us consider vdW aggregation of anatase TiO 2 and N3 dye molecule, which may allow us to gain insights into the electronic natures of the TiO 2 /N3 (proton) dye interface in DSSC. The problem of dye-anchoring on anatase TiO 2 had been argued several times in the past, e.g., ester formation with bipyridyl carboxylic acid ligand (bpy-dca), coordination of the bipyridyl ligand to Ti(IV), and hydrogen bonding with the bipyridyl ligand [14][15][16]. DFT/MM now verifies that N3 dye can be anchored by van der Waals force of hydrogen bonding between the hydroxyl group on TiO 2 and the carboxyl group of N3 dye (Figure 8). Heat of formation of −13.1 kcal/mol is not large enough to form strong vdW aggregation. We failed to obtain UV-Vis absorption spectra of the aggregate. However, surface electron energy structures, EHOMO and ELUMO and their configurations verify that photoexcitation of N3 dye induces electron injection from HOMO's on SCN ligands to LUMO's on anatase TiO 2 .

Verification of Van der Waals Aggregation of N3 Dye with Anatase TiO2 Model
Here, let us consider vdW aggregation of anatase TiO2 and N3 dye molecule, which may allow us to gain insights into the electronic natures of the TiO2 /N3 (proton) dye interface in DSSC. The problem of dye-anchoring on anatase TiO2 had been argued several times in the past, e.g., ester formation with bipyridyl carboxylic acid ligand (bpy-dca), coordination of the bipyridyl ligand to Ti(IV), and hydrogen bonding with the bipyridyl ligand [14][15][16]. DFT/MM now verifies that N3 dye can be anchored by van der Waals force of hydrogen bonding between the hydroxyl group on TiO2 and the carboxyl group of N3 dye (Figure 8). Heat of formation of −13.1 kcal/mol is not large enough to form strong vdW aggregation. We failed to obtain UV-Vis absorption spectra of the aggregate. However, surface electron energy structures, EHOMO and ELUMO and their configurations verify that photoexcitation of N3 dye induces electron injection from HOMO's on SCN ligands to LUMO's on anatase TiO2.
The most stabilized vdW aggregate of N3 (proton) and K + I3 − can be regarded as a new type of triiodide aggregated ruthenium dye. DFT/MM assigns the absorption maximum to molecular orbital components. All of them are regarded as results of photoexcitation from HOMO's on SCN ligands to LUMO's on bipyridyl ligands ( Figure 10).
TBP predict further vdW aggregation. Especially, vdW aggregation between N3 and K + I3 − , and (KI3)2 and AN or (KI3)2 and TBP is apparent in view of ESPM. The above-mentioned high Voc = 0.6-0.77 mV and the fill factors 0.65-0.72 of DSSC are validated as follows: Photo-formed electron runs through degenerate LUMO and LUMO + 1. Each LUMO + 1 and LUMO configuration visualize effective electron transfer in DSSC devices by following comparable potential energy of the KI3-assisted aggregates ( Figure 11).
Photoelectron on irradiated Ti9O18H and OH and N3 (proton) has a potential from −3.22 eV to −3.53 eV and transfer to LUMO + 1 and LUMO with potential from −4.20 eV to −4.63 eV. Photoelectron may have potential energy of 0.98-1.1 eV on anode to cathode. In general, clear-cut Voc becomes measurable by effective electron transfer via electrolytes of (KI3)2 and AN, and of (KI3)2 and TBP. The measured Voc = 0.6-0.77 mV is less than theoretical ones of 0.98-1.1 eV, and means that photoelectron leakage is not suppressed enough on N3-anchored TiO2 NP electrodes. Improved conductivity of (KI3)2 and TBP (energy gap: 0.05 eV) is consistent with enhancement of conversion efficiency of solar cells. The IPCE action spectrum ranges are 350-800 nm; at 350 nm (about 50%), at 500-630 nm (70-80%) and at 800 nm (about 5%) [13]. The excellent IPCE action spectra for N3-sensitized DSSC is now validated theoretically due to effective visible light-absorption of N3 and K + I − dyes on porous TiO2 films. Figure 12 shows a comparison of the UV-Vis spectra of N3 (proton), N3 and K + I3 − , and N3 (tetramethylammonium) as a homologue of N719. UV-Vis allowed transition occurs between highest occupied molecular orbital (HOMO), HOMO − 1, and HOMO − 2 on SCN ligands to lowest unoccupied molecular orbital (LUMO), LUMO + 1, and LUMO + 2 on bpy-dca ligands. In addition, when the molecular structures are compared, the coordination structures of N3 and K + I3 − and N3 (tetramethylammonium) are quite comparable, explaining that the respective 280 nm and 380 nm blue shift (bathochromic shift) comes from ligand-coordination change induced by vdW and The above-mentioned high V oc = 0.6-0.77 mV and the fill factors 0.65-0.72 of DSSC are validated as follows: Photo-formed electron runs through degenerate LUMO and LUMO + 1. Each LUMO + 1 and LUMO configuration visualize effective electron transfer in DSSC devices by following comparable potential energy of the KI 3 -assisted aggregates ( Figure 11).
Photoelectron on irradiated Ti 9 O 18 H and OH and N3 (proton) has a potential from −3.22 eV to −3.53 eV and transfer to LUMO + 1 and LUMO with potential from −4.20 eV to −4.63 eV. Photoelectron may have potential energy of 0.98-1.1 eV on anode to cathode. In general, clear-cut V oc becomes measurable by effective electron transfer via electrolytes of (KI 3 ) 2 and AN, and of (KI 3 ) 2 and TBP. The measured V oc = 0.6-0.77 mV is less than theoretical ones of 0.98-1.1 eV, and means that photoelectron leakage is not suppressed enough on N3-anchored TiO 2 NP electrodes. Improved conductivity of (KI 3 ) 2 and TBP (energy gap: 0.05 eV) is consistent with enhancement of conversion efficiency of solar cells.
The IPCE action spectrum ranges are 350-800 nm; at 350 nm (about 50%), at 500-630 nm (70-80%) and at 800 nm (about 5%) [13]. The excellent IPCE action spectra for N3-sensitized DSSC is now validated theoretically due to effective visible light-absorption of N3 and K + I − dyes on porous TiO 2 films. Figure 12 shows a comparison of the UV-Vis spectra of N3 (proton), N3 and K + I 3 − , and N3 (tetramethylammonium) as a homologue of N719. UV-Vis allowed transition occurs between highest occupied molecular orbital (HOMO), HOMO − 1, and HOMO − 2 on SCN ligands to lowest unoccupied molecular orbital (LUMO), LUMO + 1, and LUMO + 2 on bpy-dca ligands. In addition, when the molecular structures are compared, the coordination structures of N3 and K + I 3 − and N3 (tetramethylammonium) are quite comparable, explaining that the respective 280 nm and 380 nm blue shift (bathochromic shift) comes from ligand-coordination change induced by vdW and Coulomb interaction of N3 with I 3 − and tetramethylammonium. The hyperchromic effect of increase of strength (ε) and the bathochromic effect (the red shift) in absorption spectrum of N3 (tetramethylammonium) validates that N719 is a respectable sensitizer of DSSC.
Energies 2020, 13, x FOR PEER REVIEW 12 of 15 Coulomb interaction of N3 with I3 − and tetramethylammonium. The hyperchromic effect of increase of strength (ε) and the bathochromic effect (the red shift) in absorption spectrum of N3 (tetramethylammonium) validates that N719 is a respectable sensitizer of DSSC.

Conclusion
Experimental and theoretical UV-Vis spectral analyses of sensitizing dyes on TiO2 NPs were not successful enough so far [17][18][19][20]. As we reported [3,5], Spartan (Q-Chem) (Wavefunction, Inc. Irvin. CA, USA) is a comprehensive ab-initio quantum-chemistry software for accurate predictions of molecular structures, reactivities, and vibrational, electronic and NMR spectra. We are very proud of the verification power of density functional theory-based molecular modeling (DFT/MM) using Spartan (involving Q-Chem), and now demonstrate that vdW aggregation power of triiodide ion contributes to solar light harvesting of N3 and effective electron diffusion in in situ formed triiodide ion-based electrolytes. The strong tendency of vdW aggregation via I3 − at interfaces of DSSC facilitates DSSC fabrication with high performance of fill factor (ff) and respectable open circuit voltage (Voc).
As a perspective conclusion, DFT/MM is a useful tool for verifying and predicting conductivity of electronic molecular aggregates, IPCE action spectra and theoretical Voc of solar cells, as mentioned in the DFT/MM of perovskite solar cells [3].
During editing period of the article, we succeeded to obtain DFT/MM-based UV-Vis spectrum of van der Waals aggregate of N3(proton) with two acetonitrile (AN) molecules (see Supplementary).

Conclusions
Experimental and theoretical UV-Vis spectral analyses of sensitizing dyes on TiO 2 NPs were not successful enough so far [17][18][19][20]. As we reported [3,5], Spartan (Q-Chem) (Wavefunction, Inc. Irvin. CA, USA) is a comprehensive ab-initio quantum-chemistry software for accurate predictions of molecular structures, reactivities, and vibrational, electronic and NMR spectra. We are very proud of the verification power of density functional theory-based molecular modeling (DFT/MM) using Spartan (involving Q-Chem), and now demonstrate that vdW aggregation power of triiodide ion contributes to solar light harvesting of N3 and effective electron diffusion in in situ formed triiodide ion-based electrolytes. The strong tendency of vdW aggregation via I 3 − at interfaces of DSSC facilitates DSSC fabrication with high performance of fill factor (ff ) and respectable open circuit voltage (V oc ). As a perspective conclusion, DFT/MM is a useful tool for verifying and predicting conductivity of electronic molecular aggregates, IPCE action spectra and theoretical V oc of solar cells, as mentioned in the DFT/MM of perovskite solar cells [3].
During editing period of the article, we succeeded to obtain DFT/MM-based UV-Vis spectrum of van der Waals aggregate of N3(proton) with two acetonitrile (AN) molecules (see Supplementary).
It is now verified that the aggregation of N3 with AN, i.e. solvation of AN to N3, also affect energy structures of N3, giving four λmax: 692 nm (0.0008), 823 nm (0.0242), 875 nm (0.0164), 875 nm (0.0164). Compared to λmax of N3 (proton)K+I 3 − , we conclude that the K + I 3 − aggregation to N3 is much more dramatic, causing favorable energy structure change of N3 (proton) followed by coordination geometry change.

Materials and Methods
The DFT-based MM with the SPARTAN'18 (Wavefunction, Inc. Irving, CA, USA) program code can be conducted by using high-end personal computers, deriving equilibrium geometries of van der Waals (vdW) aggregates such as I 2 and KI, and their aggregates with a Ru complex (N3 dye) and TiO 2 NP, along with the isosurface of the electrostatic potential based on the electronic charge distribution, thus helping us to gain insights into the nature of the iodine-based vdW bonds. Similar to the previous works [1][2][3][4][5][6], the DFT-based MM with B3LYP/6-31G(d) was employed to obtain the equilibrium geometry and electronic structure of vdW aggregates. Prior to the DFT calculation, the geometries of the molecular aggregates were determined by the molecular mechanics with the Merck molecular force field (MMFF) [21], so that the intermolecular vdW interaction and the electrostatic interaction was taken into account. The following calculation at the B3LYP/6-31G(d) level results in equilibrium geometry, van der Waals bonding is correctly described.
To simulate UV-Vis spectra revealing allowed electronic transitions with wavelength and intensity, we used the time-dependent DFT within the linear response regime [22,23]. The surface electron density and electrostatic potential map were used to elucidate molecular aggregation with vdW and electrostatic Coulomb interaction without covalent bonding. For a review of the quantum mechanical methods used in SPARTAN, see reference [24][25][26][27].