Increasing Charge Carrier Mobility through Modifications of Terminal Groups of Y6: A Theoretical Study

The applications of non-fullerene acceptor Y6 with a new type of A1-DA2D-A1 framework and its derivatives have increased the power conversion efficiency (PCE) of organic solar cells (OSCs) up to 19%. Researchers have made various modifications of the donor unit, central/terminal acceptor unit, and side alkyl chains of Y6 to study the influences on the photovoltaic properties of OSCs based on them. However, up to now, the effect of changes of terminal acceptor parts of Y6 on the photovoltaic properties is not very clear. In the present work, we have designed four new acceptors—Y6-NO2, Y6-IN, Y6-ERHD, and Y6-CAO—with different terminal groups, which possess diverse electron-withdrawing ability. Computed results show that with the enhanced electron-withdrawing ability of the terminal group, the fundamental gaps become lower; thus, the wavelengths of the main absorption peaks of UV-Vis spectra red-shifts and total oscillator strength increase. Simultaneously, the electron mobility of Y6-NO2, Y6-IN, and Y6-CAO is about six, four, and four times faster than that of Y6, respectively. Overall, Y6-NO2 could be a potential NFA because of its longer intramolecular charge-transfer distance, stronger dipole moment, higher averaged ESP, enhanced spectrum, and faster electron mobility. This work provides a guideline for the future research on modification of Y6.


Introduction
During the past few years, the power conversion efficiency (PCE) of single-junction organic solar cells (OSCs) has been dramatically increased, and researchers have been creating new records. At present, the highest PCE of single-junction OSCs reaches 19% [1]. All these breakthroughs can be attributed to the emergence of the revolutionary Y6 (BTPTT-4F) nonfullerene acceptor (NFA), synthesized by Zou and co-workers [2]. With a new A 1 -DA 2 D-A 1 (A: acceptor group; D: donor group) molecular framework, Y6 consists of two strongly electron-withdrawing peripheral groups linked by a ladder-like multi-fused DAD group. Compared with the previous A-D-A type of NFAs such as ITIC and its derivatives, the introduced DA 2 D core in Y6 provides a stronger intramolecular electronic push-pull (IEPP) effect. This can further reduce the optical gap and enhance light harvesting in the nearinfrared (NIR) region [3]. So, the Y6-based OSC is superior to the ITIC-based one, although the latter's PCE is already up to 13.0% [4].
Until now, researchers have applied several strategies to make molecular modifications of Y6. The first one is to change the side alkyl chain of Y6. For instance, the side alkyl chains on the pyrrole ring of Y6 were substituted with 2-butyloctyl chains to synthesize the new acceptor, BTP-4F-12 [5], and a higher PCE of 16.4% (OSC: PBDB-TF/BTP-4F-12) was achieved, whereas the PCE of the prototype of the PBDB-TF/Y6 OSC is 15.3%. Another example is NFA L8-BO, which was synthesized by substituting the alkyl chain on the thiophene ring of Y6 with branched alkyl chains (2-butyloctyl) [6]; the PM6/L8-BO-based ternary OSCs had a high PCE of 18.32%. The second modification strategy is the change of the donor unit of Y6. 2 of 16 For example, by extending the π-conjugation length of Y6 by fusing one more thiophene onto one side of the thienothiophene unit, BP5T-4F was made [7], and the PCE of the OSC based on PM6/BP5T-4F reached 16.7%. Altering the central acceptor unit of Y6 is a third modification strategy. For example, the S atom on the central thiadiazole was replaced by an N atom, and an alkyl branch was attached to the N atom to make new molecule Y1-4F [8]; the PCE of the PBDB-TF/Y1-4F-based OSC was only 14.4%. Furthermore, the S atom on the central thiadiazole was replaced with an Se atom to obtain molecule Y6Se [9], and the OSC based on as-cast D18:Y6Se had a PCE of 17.7%, which is higher than that of the D18:Y6-based device (17.1%). The fourth modification strategy is to halogenate or expand the end groups of Y6. For instance, all fluorine atoms of Y6 were substituted with chlorine to obtain new NFA BTP-4Cl [10], and the PBDB-TF/BTP-4Cl-based OSC achieved a PCE of 16.5%, which was higher than the PCE (15.6%) of the PBDB-TF/Y6 OSC. Another example is BTP-FCl-FCl, in which there is one Cl substitution on each terminal group of Y6, and the PCE of the BTR-Cl/BPT-FCl-FCl OSC is 15.3% [11]. On the other hand, an asymmetric nonfullerene acceptor LL3, one of the Y6 derivatives, was designed and synthesized with one expanded norbornyl-modified 1,1-dicyanomethylene-3-indanone (CBIC) terminal group and a chlorinated 1,1-dicyanomethylene-3-indanone (IC-2Cl) terminal group, and the OSC based on LL3 had a PCE of 16.82%, which was higher than symmetric LL3-2Cl (14.4%) [12]. Another asymmetric NFA is BTP-S9, with two terminals (IC-2Cl and NC-2F (expanded ring)); the PM6:BTP-S9 OSC exhibited the efficiency of 17.56% [13]. As mentioned above, the modifications on alkyl chains, donor units, central acceptor units, and the end group of Y6 have been successfully conducted to investigate their effects on photovoltaic properties. Among them, the modification of end groups can change the stacking morphology and charge transport. The synthesis of modifying end groups is relative easy, and the stacking of terminal groups is of great significance [14]. So, it is particularly important. However, there is still a lack of research in this area, such as the substitution of terminal groups with non-halogen electron-withdrawing groups and the study of the substitution of entire end groups, which may have a significant impact on their photovoltaic performance.
Therefore, in the present work, we have designed and modeled a series of symmetric Y6 derivatives, i.e., Y6-NO 2 , Y6-IN, Y6-ERHD, and Y6-CAO, as shown in Figure 1. The electron-withdrawing ability of the terminal groups of these molecules is Y6-NO 2 > Y6 > Y6-IN > Y6-ERHD > Y6-CAO. Density functional theory (DFT) [15] and time-dependent density functional theory (TDDFT) [16] are used to study these molecules because they have proved to be highly appropriate and efficient for calculations of electronic structures and excited states of organic molecules with proper XC kernels [17]. The effect of the changes in terminal acceptor parts on the dipole moment, FMO energy, UV-Vis absorption spectra, electrostatic potential (ESP), exciton binding energy, electron/hole distribution, and charge carrier mobility have been systematically studied.
The purposes of this work are two-fold: (1) To investigate the effects of different terminal acceptor units on the photovoltaic properties of Y6.

Dipole Moments and Electrostatic Potential of the Modified Molecules
Dipole moment has important effects on molecular packing and charge separation in thin film [18]. As shown in Table 1, Y6-NO2, Y6, Y6-IN, Y6-ERHD, and Y6-CAO possess dipole moments of 15.56, 0.82, 4.16, 6.31, and 2.96 Debye, respectively. We therefore can conclude that the modification of the terminal groups has a great influence on their dipole moments. In particular, the strongest electron-withdrawing ability of the terminal groups of Y6-NO2 causes the longest distance between the positive and negative charge centers, thus Y6-NO2 has the largest dipole moment. Y6-IN, Y6-ERHD, and Y6-CAO have moderate dipole moments (4.16, 6.31, and 2.96 Debye, respectively), which are all larger than that of Y6. This should benefit crystallinity in thin film [19]. One good example is the BTP-4Cl, synthesized by Hou and co-workers [10]. It has a dipole moment of 1.44 Debye, which is slightly larger than that of Y6 (0.64 Debye). So, the crystallinity of BTP-4Cl is better than that of Y6, and the PCE of the OSC with BTP-4Cl as an NFA and PM6 as a donor rises to 16.1% [10]. On the contrary, another opposite example is BTPT-4F NFA (A1-DA2D-A1 type), which was synthesized by Fan et al. [7]. The PCE of a single-junction P2F-EHp/BTPT-4F OSC is just 1.09%. There exist two reasons for this sharp reduction: One is that the BTPT-4F molecule does not possess two C11H23 side alkyl chains and the other is that too large a dipole moment (12.34 Debye) results in strong aggregation and phase separation. Y6-NO2, possessing the largest dipole moment (15.56 Debye), may have strong aggregation and phase separation. In summary, on account of the similar structures and the same side alkyl chains, Y6-IN, Y6-ERHD, and Y6-CAO, with moderately larger dipole moments, may possess better performance than Y6 in OSCs.

Dipole Moments and Electrostatic Potential of the Modified Molecules
Dipole moment has important effects on molecular packing and charge separation in thin film [18]. As shown in Table 1, Y6-NO 2 , Y6, Y6-IN, Y6-ERHD, and Y6-CAO possess dipole moments of 15.56, 0.82, 4.16, 6.31, and 2.96 Debye, respectively. We therefore can conclude that the modification of the terminal groups has a great influence on their dipole moments. In particular, the strongest electron-withdrawing ability of the terminal groups of Y6-NO 2 causes the longest distance between the positive and negative charge centers, thus Y6-NO 2 has the largest dipole moment. Y6-IN, Y6-ERHD, and Y6-CAO have moderate dipole moments (4.16, 6.31, and 2.96 Debye, respectively), which are all larger than that of Y6. This should benefit crystallinity in thin film [19]. One good example is the BTP-4Cl, synthesized by Hou and co-workers [10]. It has a dipole moment of 1.44 Debye, which is slightly larger than that of Y6 (0.64 Debye). So, the crystallinity of BTP-4Cl is better than that of Y6, and the PCE of the OSC with BTP-4Cl as an NFA and PM6 as a donor rises to 16.1% [10]. On the contrary, another opposite example is BTPT-4F NFA (A 1 -DA 2 D-A 1 type), which was synthesized by Fan et al. [7]. The PCE of a single-junction P2F-EHp/BTPT-4F OSC is just 1.09%. There exist two reasons for this sharp reduction: One is that the BTPT-4F molecule does not possess two C 11 H 23 side alkyl chains and the other is that too large a dipole moment (12.34 Debye) results in strong aggregation and phase separation. Y6-NO 2 , possessing the largest dipole moment (15.56 Debye), may have strong aggregation and phase separation. In summary, on account of the similar structures and the same side alkyl chains, Y6-IN, Y6-ERHD, and Y6-CAO, with moderately larger dipole moments, may possess better performance than Y6 in OSCs. Furthermore, the electrostatic potential (ESP) maps of Y6-NO 2 , Y6, Y6-IN, Y6-ERHD, and Y6-CAO monomers are presented in Figure 2. Owing to the A 1 -DA 2 D-A 1 type of molecular frame, the negatively charged regions of these molecules are mainly located at the terminal group and the central benzothiadiazole (BT) unit, while the thiophene rings and side alkyl chains on the skeleton are positive. From Y6-IN to Y6 to Y6-NO 2 , with the increasing of the electron-withdrawing ability of the terminal groups, the averaged ESP of them increases from 64 meV to 457 meV. Because of the significant differences of the structures of terminal groups of Y6-ERHD and Y6-CAO (changed conjugated skeleton), both the averaged ESPs of two molecules are smaller than that of Y6. Generally, the electron acceptance ability is worse with smaller averaged ESP, so Y6-IN, Y6-ERHD, and Y6-CAO are not better acceptors from this perspective. Furthermore, the electrostatic potential (ESP) maps of Y6-NO2, Y6, Y6-IN, Y6-ERHD, and Y6-CAO monomers are presented in Figure 2. Owing to the A1-DA2D-A1 type of molecular frame, the negatively charged regions of these molecules are mainly located at the terminal group and the central benzothiadiazole (BT) unit, while the thiophene rings and side alkyl chains on the skeleton are positive. From Y6-IN to Y6 to Y6-NO2, with the increasing of the electron-withdrawing ability of the terminal groups, the averaged ESP of them increases from 64 meV to 457 meV. Because of the significant differences of the structures of terminal groups of Y6-ERHD and Y6-CAO (changed conjugated skeleton), both the averaged ESPs of two molecules are smaller than that of Y6. Generally, the electron acceptance ability is worse with smaller averaged ESP, so Y6-IN, Y6-ERHD, and Y6-CAO are not better acceptors from this perspective.

FMO and Gap Energy
The HOMO, LUMO, and gap energy are shown in Figure 3. As the electronwithdrawing ability of the terminal acceptor unit decreases, the energy gaps of these acceptors generally becomes higher (except for Y6-ERHD) and the energy of the LUMO increases. This may lead to the rise of VOC. Y6-NO2 has the smallest energy gap and the lowest LUMO due to the strongest electron-withdrawing ability of -NO2. In addition, the HOMO energy generally decreases towards the stronger electron-withdrawing ability of the terminal units. This may not be favorable for hole transfer.

FMO and Gap Energy
The HOMO, LUMO, and gap energy are shown in Figure 3. As the electron-withdrawing ability of the terminal acceptor unit decreases, the energy gaps of these acceptors generally becomes higher (except for Y6-ERHD) and the energy of the LUMO increases. This may lead to the rise of V OC . Y6-NO 2 has the smallest energy gap and the lowest LUMO due to the strongest electron-withdrawing ability of -NO 2 . In addition, the HOMO energy generally decreases towards the stronger electron-withdrawing ability of the terminal units. This may not be favorable for hole transfer.

UV-Vis Spectra
As can be seen from Figure 4, the simulated UV-Vis spectra of Y6-NO 2 , Y6, Y6-IN, Y6-ERHD, and Y6-CAO in film are provided. Table S1 shows the results of the calculated excited states of these acceptors. For Y6-NO 2 , Y6, Y6-IN, Y6-ERHD, and Y6-CAO, the wavelength of the maximum absorption peaks is at 784.17, 741.48, 646.09, 645.41, and 601.25 nm, respectively, and the total oscillator strength in the visible region is 3.77, 3.70, 3.34, 3.22, and 2.77, correspondingly. Compared with Y6, the wavelength of the absorption peak of Y6-NO 2 is red-shifted by 43 nm and its total oscillator strength is increased by 0.07. With the rise of electron-withdrawing ability (Y6-CAO < Y6-ERHD < Y6-IN < Y6 < Y6-NO 2 ), the oscillator strength increases and the UV-Vis spectrum is red-shifted (see Figures 4 and 5), which is consistent with the reduction of the energy gaps. The wavelength of absorption peaks and oscillator strength of Y6-NO 2 are larger than those of Y6, owing to the strongest electron-withdrawing ability.

UV-Vis Spectra
As can be seen from Figure 4, the simulated UV-Vis spectra of Y6-NO2, Y6, Y6-IN, Y6-ERHD, and Y6-CAO in film are provided. Table S1 shows the results of the calculated excited states of these acceptors. For Y6-NO2, Y6, Y6-IN, Y6-ERHD, and Y6-CAO, the wavelength of the maximum absorption peaks is at 784.17, 741.48, 646.09, 645.41, and 601.25 nm, respectively, and the total oscillator strength in the visible region is 3.77, 3.70, 3.34, 3.22, and 2.77, correspondingly. Compared with Y6, the wavelength of the absorption peak of Y6-NO2 is red-shifted by 43 nm and its total oscillator strength is increased by 0.07. With the rise of electron-withdrawing ability (Y6-CAO < Y6-ERHD < Y6-IN < Y6 < Y6-NO2), the oscillator strength increases and the UV-Vis spectrum is red-shifted (see Figures 4 and 5), which is consistent with the reduction of the energy gaps. The wavelength of absorption peaks and oscillator strength of Y6-NO2 are larger than those of Y6, owing to the strongest electron-withdrawing ability.

UV-Vis Spectra
As can be seen from Figure 4, the simulated UV-Vis spectra of Y6-NO2, Y6, Y6-IN, Y6-ERHD, and Y6-CAO in film are provided. Table S1 shows the results of the calculated excited states of these acceptors. For Y6-NO2, Y6, Y6-IN, Y6-ERHD, and Y6-CAO, the wavelength of the maximum absorption peaks is at 784.17, 741.48, 646.09, 645.41, and 601.25 nm, respectively, and the total oscillator strength in the visible region is 3.77, 3.70, 3.34, 3.22, and 2.77, correspondingly. Compared with Y6, the wavelength of the absorption peak of Y6-NO2 is red-shifted by 43 nm and its total oscillator strength is increased by 0.07. With the rise of electron-withdrawing ability (Y6-CAO < Y6-ERHD < Y6-IN < Y6 < Y6-NO2), the oscillator strength increases and the UV-Vis spectrum is red-shifted (see Figures 4 and 5), which is consistent with the reduction of the energy gaps. The wavelength of absorption peaks and oscillator strength of Y6-NO2 are larger than those of Y6, owing to the strongest electron-withdrawing ability.

Exciton Binding Energy
The exciton binding energy is the energy cost to separate excitons into free holes and electrons in organic solar cells. The small exciton binding energy should be beneficial to charge separation, Jsc, and PCE [20]. We provide the calculated exciton binding energy of Y6-NO 2 , Y6, Y6-IN, Y6-ERHD, and Y6-CAO, obtained with ωB97X/6-31+G(d) in film, as shown in Figure 6. The exciton binding energy is calculated with Equations (5)-(7). The values of Y6-NO 2 , Y6, Y6-IN, Y6-ERHD and Y6-CAO are 0.455, 0.454, 0.483, 0.477, and 0.501 eV, respectively. Y6 possesses the smallest exciton binding energy, and the E b of Y6-NO 2 is similar to that of Y6.

Reorganization Energy of Electron Transfer
Reorganization energy (λ) is an important parameter influencing charge mobility, and it contains both inner and external parts. The inner reorganization energy (λ i ) is the energy for the change of geometric structure in the process of charge transport. The inner reorganization energy of Y6-NO 2 , Y6, Y6-IN, Y6-ERHD, and Y6-CAO was calculated with

Exciton Binding Energy
The exciton binding energy is the energy cost to separate excitons into free holes and electrons in organic solar cells. The small exciton binding energy should be beneficial to charge separation, Jsc, and PCE [20]. We provide the calculated exciton binding energy of Y6-NO2, Y6, Y6-IN, Y6-ERHD, and Y6-CAO, obtained with ωB97X/6-31+G(d) in film, as shown in Figure 6. The exciton binding energy is calculated with Equations (5)

Reorganization Energy of Electron Transfer
Reorganization energy (λ) is an important parameter influencing charge mobility, and it contains both inner and external parts. The inner reorganization energy (λi) is the energy for the change of geometric structure in the process of charge transport. The inner reorganization energy of Y6-NO2, Y6, Y6-IN, Y6-ERHD, and Y6-CAO was calculated with Equation (8)

Exciton Binding Energy
The exciton binding energy is the energy cost to separate excitons into free holes and electrons in organic solar cells. The small exciton binding energy should be beneficial to charge separation, Jsc, and PCE [20]. We provide the calculated exciton binding energy of Y6-NO2, Y6, Y6-IN, Y6-ERHD, and Y6-CAO, obtained with ωB97X/6-31+G(d) in film, as shown in Figure 6. The exciton binding energy is calculated with Equations (5)

Reorganization Energy of Electron Transfer
Reorganization energy (λ) is an important parameter influencing charge mobility, and it contains both inner and external parts. The inner reorganization energy (λi) is the energy for the change of geometric structure in the process of charge transport. The inner reorganization energy of Y6-NO2, Y6, Y6-IN, Y6-ERHD, and Y6-CAO was calculated with Equation  The external reorganization energy (λo) means the energy required for the nuclear relaxation and electron polarization of surrounding molecules during charge transfer. The results are calculated according to Equation (9), as shown in Table 2. Y6 and Y6-CAO have the smallest and biggest total reorganization energy, respectively. New molecules with modified terminal groups all have a bigger total reorganization energy than Y6. The external reorganization energy (λ o ) means the energy required for the nuclear relaxation and electron polarization of surrounding molecules during charge transfer. The results are calculated according to Equation (9), as shown in Table 2. Y6 and Y6-CAO have the smallest and biggest total reorganization energy, respectively. New molecules with modified terminal groups all have a bigger total reorganization energy than Y6.

Electron/Hole Distribution of Monomers
The electron/hole distributions of the first LE states of five molecules are shown in Figure 8. The intramolecular charge-transfer amounts (q CT ) of five molecules are different (0.414-0.576 e), indicating that changing the terminal acceptor unit affects q CT . The intramolecular charge transition distance (D CT ) of Y6-NO 2 , Y6, Y6-IN, Y6-ERHD, and Y6-CAO is 3.83, 1.47, 1.34, 1.29, and 0.91 Å, respectively, demonstrating that D CT rises with the enhanced electron-withdrawing ability of terminal groups. So, the intermolecular charge transfer in Y6-NO 2 may be easier than that in Y6.

Quadrupole Moment along π-π Stacking Direction, ∆LUMO and ∆HOMO
Y6-NO 2 , Y6 and Y6-ERHD and Y6-CAO have a positive Q π, whereas Y6-IN has a negative Q π as shown in Figure 9. Those of Y6-NO 2 (177.55 ea 0 2 ), Y6 (100.30 ea 0 2 ) and Y6-ERHD (33.09 ea 0 2 ) are big. Except for Y6-CAO (4.36 ea 0 2 ), the ∆LUMO between monomers and corresponding dimers with M configurations is quite different from that between monomers and corresponding dimers with S configurations. Y6-NO 2 , Y6 and Y6-ERHD show a large positive Q π, which can account for the large energetic shift with different molecular orientations, as observed in M and S configurations [21]. However, there is no law for ∆HOMO.
Y6-NO2, Y6 and Y6-ERHD and Y6-CAO have a positive Qπ, whereas Y6-IN has a negative Qπ as shown in Figure 9. Those of Y6-NO2 (177.55 ea0 2 ), Y6 (100.30 ea0 2 ) and Y6-ERHD (33.09 ea0 2 ) are big. Except for Y6-CAO (4.36 ea0 2 ), the ΔLUMO between monomers and corresponding dimers with M configurations is quite different from that between monomers and corresponding dimers with S configurations. Y6-NO2, Y6 and Y6-ERHD show a large positive Qπ, which can account for the large energetic shift with different molecular orientations, as observed in M and S configurations [21]. However, there is no law for ΔHOMO.

Binding Energy of Dimers and Intermolecular Distance of Charge Transfer
In this section, the binding energies and the geometry center distances (D GC ) of Y6-NO 2 , Y6, Y6-IN, Y6-ERHD, and Y6-CAO dimers with M and S configurations (see Figure 10) are presented, as shown in Figure 11 and Tables S2 and S3. We note that all the side alkyl chains are not contained when measuring the geometry center distances since they do not have noticeable contributions to the FMOs. Figure 11 shows that the Y6-IN and Y6-CAO dimers have larger binding energies, suggesting that the Y6-IN and Y6-CAO molecules may be more tightly packed in the dimer. In addition, Y6-NO 2 , Y6-IN, Y6-ERHD, and Y6-CAO M-configuration dimers have larger binding energy, showing that they may be more stable than S-configuration dimers, except for Y6. On the contrary, for Y6, the S-configuration dimer is more stable, and this is consistent with the experimental results [22]. Clearly, the change of the terminal group makes the stable configuration of the dimer transform from the S-type to the M-type. Figure 12 presents the stable dimer configurations of these molecules. We shall focus on these dimers for charge transfer and mobility. We also calculated the dihedral angles of the side chains of the dimers with both M and S configurations to understand the effect of side chains on molecular stacking in dimers (SI Figures S1 and S2). However, we could not find any law regarding the dihedral angles of side chains. Almost in all the dimers with both M and S configurations, the orientations of the side chains near the end group stacking of the dimer tend to be horizontal, and the orientations of the side chains far from the end group stacking are arbitrary [23].

Charge Carrier Mobility
The charge-transfer distance in these dimers, shown in Table 3, is estimated by using the geometry center distance. In general, a short charge-transfer distance is conducive for enhancing electronic coupling but may not be beneficial for carrier mobility (see Equation (13)). On the other hand, electronic coupling represents the overlap between the FMOs of adjacent molecules and is critical for charge carrier mobility. The one-dimensional charge carrier mobility is calculated by assuming that the charge carrier can only jump to its nearest neighbor and its continuous jumps are not correlated [24]. The charge-transfer rate constants, electronic coupling, charge carrier mobility, and other related parameters of stable dimers are displayed in Table 3. We also provide these data for other dimer configurations in SI Table S4.       To further understand the reason for the differences in electronic coupling and electron mobility among these dimers, we measured the geometric center distances between the end groups in these dimers, as shown in Figure 13. The end group geometric center distances of Y6-NO 2 (3.77 Å) and Y6-IN (3.70 Å) are shorter than that of Y6 (6.89 Å), resulting in stronger electron coupling. On the other hand, although the geometric center distance of Y6-ERHD (3.90 Å) is much shorter than that of Y6 (6.89 Å), its terminal conjugated ring is too small when they stack with each other, leading to very low electronic coupling and a low charge-transfer rate constant. For Y6-CAO, the flexible end groups make the stacks less rigid, leading to high electronic coupling and electron mobility. Clearly, Y6-NO 2 possesses the fastest electron mobility because it has strong electronic coupling and the longest charge-transfer distance simultaneously. Based on these results, we can conclude that the change of terminal group has a great influence on the electron mobility of these Y6 derivatives.

Methods
All molecules were constructed with GaussView version 5.0 [25]. All calculations were calculated by using Guassian09 Rev E.01 software package [26]. In consideration of environmental effects in the thin film of OSCs, we adopted an integral equation formalism variant-polarizable continuum model (IEF-PCM), and the dielectric constant was set to 3.0 [27,28]. The Y6-NO2, Y6-IN, Y6-ERHD, and Y6-CAO molecules were initially optimized with the combination of the B3LYP [29] hybrid density functional and 6-31G(d) basis set, which has proven to be reliable for the geometric optimization of organic molecules [30,31]. This combination also was used for the calculation of the inner reorganization energy. If not mentioned otherwise, the 6-31+G(d) basis set was default for all other calculations since it was more accurate for excited state calculations and weak

Methods
All molecules were constructed with GaussView version 5.0 [25]. All calculations were calculated by using Guassian09 Rev E.01 software package [26]. In consideration of environmental effects in the thin film of OSCs, we adopted an integral equation formalism variant-polarizable continuum model (IEF-PCM), and the dielectric constant was set to 3.0 [27,28]. The Y6-NO 2 , Y6-IN, Y6-ERHD, and Y6-CAO molecules were initially optimized with the combination of the B3LYP [29] hybrid density functional and 6-31G(d) basis set, which has proven to be reliable for the geometric optimization of organic molecules [30,31]. This combination also was used for the calculation of the inner reorganization energy. If not mentioned otherwise, the 6-31+G(d) basis set was default for all other calculations since it was more accurate for excited state calculations and weak interactions [32]. We also selected the range separated hybrid (RSH) density functional ωB97X with tunable parameter ω [33] because ωB97X has proven to be accurate for Y6 and its derivatives [34]. By tuning the parameter ω with IEF-PCM, the gained absolute values of HOMO and LUMO energies are almost equal to those of EA and IP. We used the following equations to optimize ω [33]: where E ω (N) is the energy of the neutral system; ε ω HOMO (N) is the energy of the HOMO of the neutral system; E ω (N + 1) is the energy of the anion; and E ω (N − 1) is the energy of the cation. The optimized ω of Y6, Y6-NO 2 , Y6-IN, Y6-ERHD and Y6-CAO is 0.016, 0.013, 0.015, 0.014, and 0.016, respectively. Multiwfn 3.8 software was adopted to simulate absorption spectra with the output of TDDFT calculations [35]. The full width at half maximum (FWHM) of absorption peaks was default (0.6667 eV). For the wavelength of simulated absorption peaks, the following formula was used: where λ ave represents the wavelength of the absorption peak; f n means the oscillator strength of the nth excited state; and λ n donates the wavelength of the nth excited state [34]. Exciton binding energy (E b ) in thin films is a key factor determining charge transfer and separation. We used the formulas as follows [36]: where E b donates exciton binding energy; E opt g means the energy of the optical gap; E fund g is the energy of the fundamental gap; N represents a neutral system; N + 1 is a negative; and N − 1 is a positive. The inner reorganization energy is the energy required for the change of nuclear geometric structure during charge transfer. The following equation was used to calculate the inner reorganization energy [37,38]: where λ i, e means the inner reorganization energy during electron transfer; E 0 (A − ) is the total energy of the neutral acceptor in the optimized geometry of the anion; E 0 A 0 represents the total energy of the neutral acceptor in the optimized geometry of the neutral molecule; E − A 0 denotes the energy of the negative acceptor in the optimized geometry of the neutral molecule; and E − (A − ) is energy of the negative acceptor in the optimized geometry of the anion.
Taking into account environmental impact, external reorganization energy (λ o ) is defined as the energy required for the molecular nuclear relaxation and electronic polarization of the surrounding molecules of acceptor dimers during the charge-transfer process. The formula for calculating the external reorganization energy (λ o ) is [39,40] is: where ∆q is the transferred element charge; ε 0 is the vacuum permittivity (ε 0 = 8.85 × 10 −12 F/m); and ε opt and ε s represent the optical and static permittivity of the solvents, respectively, which were set to 2.0 and 3.0, correspondingly. In addition, R D and R A donate the effective radius of the donor and acceptor, respectively. In these dimers, R D is equal to R A , and r DA was set to the effective radius of the dimers. Basis set superposition error (BSSE) correction was applied to correct the binding energy of dimers. The binding energy was obtained by the following equation [41]: where E AB is the total energy of the dimer; E A represents the energy of monomer A; E B denotes the energy of monomer B; and E BSSE means the BSSE energy. The semi-classical Marcus formula is written as follows [39,40]: where k B is the Boltzmann constant; V denotes the electronic coupling, T represents the absolute temperature and is set to 298.15 K; λ means the reorganization energy, and in this paper, it is the sum of the inner reorganization energy and external reorganization energy; is the reduced Planck constant; and the free energy difference between initial and final states was set to zero.
The electronic coupling was defined by the following equation [42]: where H i f means the non-diagonal element of the two-state Hamiltonian matrix; H ii /H f f denotes the Hamiltonian of the initial/final state of charge transfer, respectively; and S i f is the overlap between initial and final states. The overall electronic coupling was calculated at CAM-B3LYP/6-31G(d) level of theory [43]. The carrier mobility could be computed by the Einstein-Smoluchowski formula [44]: where e is the elementary charge; D means the diffusion coefficient of the charge carrier; for the one-dimensional case, D = r 2 k/2 (supposing that the charge carrier can only jump to its nearest neighbor sites and its continuous jumps are uncorrelated) [24]; k denotes the chargetransfer rate constant; and r is the charge-transfer distance between two adjacent molecules, which is estimated by the distance between geometry centers of neighboring molecules. In order to simulate molecular stacking, we modeled the dimer configuration of these designed molecules, which is simple but very useful. Moreover, we assumed that the charge carrier can only jump to the nearest location, which is reasonable [45]. So, only dimer configurations that contain two neighboring monomers were considered in this work. There are different configurations for dimers, but experimental results have shown that the terminal-to-terminal stacking of the dimers of Y6 and its derivatives are dominant [43].