The Halogenation Effects of Electron Acceptor ITIC for Organic Photovoltaic Nano-Heterojunctions

Molecular engineering plays a critical role in the development of electron donor and acceptor materials for improving power conversion efficiency (PCE) of organic photovoltaics (OPVs). The halogenated acceptor materials in OPVs have shown high PCE. Here, to investigate the halogenation mechanism and the effects on OPV performances, based on the density functional theory calculations with the optimally tuned screened range-separated hybrid functional and the consideration of solid polarization effects, we addressed the halogenation effects of acceptor ITIC, which were modeled by bis-substituted ITIC with halogen and coded as IT-2X (X = F, Cl, Br), and PBDB-T:ITIC, PBDB-T:IT-2X (X = F, Cl, Br) complexes on their geometries, electronic structures, excitations, electrostatic potentials, and the rate constants of charge transfer, exciton dissociation (ED), and charge recombination processes at the heterojunction interface. The results indicated that halogenation of ITIC slightly affects molecular geometric structures, energy levels, optical absorption spectra, exciton binding energies, and excitation properties. However, the halogenation of ITIC significantly enlarges the electrostatic potential difference between the electron acceptor and donor PBDB-T with the order from fluorination and chlorination to bromination. The halogenation also increases the transferred charges of CT states for the complexes. Meanwhile, the halogenation effects on CT energies and electron process rates depend on different haloid elements. No matter which kinds of haloid elements were introduced in the halogenation of acceptors, the ED is always efficient in these OPV devices. This work provides an understanding of the halogenation mechanism, and is also conducive to the designing of novel materials with the aid of the halogenation strategy.


Introduction
Organic photovoltaics (OPVs) have received extensive attention in recent years due to their light weight, low cost, semitransparency, flexibility, etc. The design and development of high-performance OPV devices are of great significance to solving the energy crisis [1][2][3]. Up to now, the power conversion efficiency (PCE) record of OPVs was over 18% [4][5][6][7]. Most of the electron processes for OPV power conversion happen in the active layer which is composed of both electron donor and acceptor materials. As the core of the OPV device, the changing active layer materials can tune photovoltaic performances [8,9]. In the past few years, the design and development of new electron donor and acceptor materials for OPV active layer have promoted the rapid improvement of PCE. Fullerene and its derivatives were usually used as electron acceptor materials due to their high electron affinity and electron mobility. However, their shortcomings of fullerenes limit their extensive application in OPV, including for complicated manufacturing processes, poor absorption in the visible and infrared region, and scare energy level adjustability. Non-fullerenes acceptors (NFAs) have many advantages over fullerenes, such as tunable energy levels and high molar absorption coefficients, resulting in a situation where NFAs have gradually replaced fullerenes as electron acceptors in OPV.
Backbone modification, side-chain engineering, and end-group substitution are common strategies for designing new electron donor and acceptor materials of OPV. Among these methods, halogen substitution has been shown to be effective for improving OPV performance. A fluorine atom is a type of halogen atom with the strongest electronegativity (3.98) and the smallest size (the van der Waals radius of 147 pm), which can effectively adjust the π-electron properties without causing unexpected steric hindrance. Experiments also proved that some chlorinated and brominated systems showed similar or stronger capability for adjusting π-conjugated electronic properties [10][11][12], which were caused by the larger dipole moments of C-Cl and C-Br bonds than that of a C-F bond [13]. Therefore, it is extremely important to research the influence of halogenation on materials' physical properties and to explore the relationship between physical properties and photovoltaic performance.
Developing donor or acceptor materials by halogenation improved OPV performance [14,15]. Pei et al. found that the chlorinated backbone of a isoindigo-based polymer donor and acceptor (PC 71 BM) showed an appropriate degree of phase separation, reduced film crystallinity, and improved OPV performance (which was ascribed to the donor molecular stacking orientation [16]). Forrest et al. synthesized the BT-CIC by chlorinating BT-IC, and found that the chlorination caused the improved planarity of molecular geometry and significant red-shift of absorption spectrum (about 60 nm) with the 1.33 eV bandgap, resulting in an external quantum efficiency (EQE) of 75% in the wavelength region of 630-850 nm and of more than 11% of PCE [17]. Yan et al. reported that the C8-ITCC-Cl NFA which is chlorinated thienyls in C8-ITCC increased PCE from 10.8% to 12.7%, and found the chlorination of NFA induced the stronger noncovalent interactions which is favorable to improve the layered π-π stacking [18]. This highlighted the potential of halogenated thiophene end groups in NFA to improve PCE. Yasuda et al. designed a new type of NFA BDT-ID-X (X = F, Cl, Br, I), and found that subtle halogenations at the NFA end group can significantly affect the photoelectronic properties and film morphology [12]. Among them, Br and I substitutions can drastically reduce the surface free energies of donor/acceptor interfaces, forming fine bulk heterojunction (BHJ) for OPV. Hou et al. investigated the fluorinated electron donor PBDB-TF, fluorinated NFA BTP-4F, and chlorinated BTP-4Cl, and found that the halogenation can increase open-circuit voltage (V OC ) and PCE (15.6% and 16.5%, respectively) by reducing energy losses [19]. The corresponding energy losses are 0.60 eV and 0.53 eV, and the non-radiative energy losses are 0.25 eV and 0.21 eV, respectively, which are lower than those of most reported OPVs. Cao et al. designed the NFA BDTBO-4F and BDTBO-4Cl through halogenation strategy and adopted PM6 as electron donor for fabricating OPV. Hole transfer at donor/acceptor interfaces was still efficient though the highest occupied molecular orbital (HOMO) offsets between the donor and acceptor are smaller than 0.1 eV, and the PM6:BDTBO-4F device exhibited a low charge recombination (CR) rate, high charge mobility, and suitable phase separation [20]. Also, halogenation at different substitution positions can cause changes in π-π stacking distance, supermolecular interaction, miscibility, and steric hindrance, which affect the aggregation behavior in blends [21][22][23][24][25][26][27].
In addition to experimental studies, many theoretical works were also conducted to understand the halogenation effects on OPV material properties [28,29]. Bredas et al. compared the electronic properties of polymers PBT4T-2OD and PffBT4T-2OD for understanding the correlation between solution temperature-dependent aggregation and solid-state stacking, and found that fluorination can enhance the planarity of the donor polymer backbone and increase the interaction among polymers, promoting molecule packing order for higher PCE [30]. Koehler et al. found that fluorination tends to reduce exciton binding energy (E exb ), improve exciton dissociation (ED) efficiency, promote carrier mobility, restrict bimolecular CR, and reduce the resistance of bulk heterojunction [31]. Hou et al. reported that halogenation generates a larger difference of molecular surface average electrostatic potential (ESP) between the electron donor and acceptor [32]. Therefore, it induces a stronger intermolecular electric field at donor/acceptor interface which can increase ED efficiency and then enlarge the short-circuit current density (J sc ) [33]. Based on the designed compounds, Mahmood et al. found that the fluorination of the end-groups increased the dipole moment difference between the ground state and the first excited state with the reduced E exb , and also enhanced electronic coupling for charge transport [34]. Our previous work indicated that electron donor halogenation can remarkably influence energy level alignment at heterojunction interfaces, and the improved PCE induced by donor halogenation can be mainly ascribed to the increased charge transfer (CT) excitation energies and suppression of CR processes [35].
The ITIC is a superstar of NFA and has been received extensive research attention since it inaugurated a new era of OPV [31,36,37]. The A-D-A type molecular architecture of ITIC, containing large planar conjugated backbone and strong electron-withdrawing end groups, promotes the π-electron delocalization and intramolecular charge transfer/transport. The study of the fluorinated and chlorinated ITIC by Hou et al. found that chlorination has a stronger ability to stabilize molecular energy levels and broaden the absorption spectrum than fluorination, and the position of halogen atoms can affect the crystallinity of halogenated ITIC [38,39]. Among the OPV devices based on bis-substituted halogenation (F, Cl, Br) of ITIC as NFAs and PBDB-T as electron donor, the fluorinated system generated the largest J sc , while the brominated system exhibited the highest FF (about 0.71) with a small amount of agglomeration area in the blended film, leading to the reduction of EQE and J sc [40].
Overall, why halogenation of active layer materials usually improve OPV performance can be qualitatively attributed to two aspects. On the one hand, halogenation tunes electronic structures and excitation properties, including stabilizing molecular energy levels, enhancement of light absorption, increasing EQE and, thus J sc , and reducing energy losses and then increasing V oc [41]. On the other hand, halogenation improves aggregation, packing order, crystallization and film morphologies due to the strengthened intramolecular or intermolecular interaction, achieving effective ED and charge transfer/transport [24,40,[42][43][44]. However, the fundamental understanding of halogenation mechanisms, including the effects on electron processes at interface and OPV work principles, still remains uninvestigated. Here, to investigate acceptor halogenation effects, we selected PBDB-T:ITIC and PBDB-T:IT-2X (X = F, Cl, Br) as the model systems whose molecular structures are shown in Figure 1, and analyzed the molecular geometries, electronic structures, optical properties, excitation properties, ESP and the kinetics of electron processes at the heterojunction interfaces.

Computational Methods
The molecular structures of PBDB-T, ITIC and IT-2X (X = F, Cl, Br) were optimized at the level of ωB97XD/6-31G(d,p) [45,46]. Two repeat units were considered as the molecular structure of polymer donor PBDB-T in calculations, since the calculated frontier molecular orbitals (MOs) of polymer donor D18 whose molecular structure is similar to that of PBDB-T indicated that the frontier MOs are mainly localized on two repeat units [47]. Based upon the optimized molecular structures of PBDB-T, ITIC and IT-2X (X = F, Cl, Br), the PBDB-T:ITIC and PBDB-T:IT-2X (X = F, Cl, Br) heterojunction interfaces were modeled by using their complexes, in which the face-on configuration was considered due to the experimental results and the promoted π-π stacking [48]. To reduce the computational cost, the alkyl chains (-C 2 H 5 , -C 4 H 9 in the polymer donor and -C 6 H 13 in the acceptors) in the complexes were substituted with -CH 3 , since the side chain only affects the active layer morphology and does not affect other properties [49]; this simplification has been widely used in other works [31,50,51]. The constructed PBDB-T:ITIC and PBDB-T:IT-2X (X = F, Cl, Br) complexes were further optimized using ωB97XD/6-31G(d,p) methods. The excitation properties of donor and acceptor molecules, as well as the PBDB-T:ITIC and PBDB-T:IT-2X (X = F, Cl, Br) complexes, were calculated using time-dependent density functional theory (TDDFT) methods. The solid environment of OPV devices was simulated by self-consistent reaction field (SCRF) method with polarization continuum model (PCM) [52], requiring the static and dynamic dielectric constant (ε s and ε d ) to describe the polarization effects of solid environment. The ε s and ε d of models are 3.5 and 3.3, respectively.
TDDFT results heavily depend on the selected DFT functional. The long range corrected hybrid functional with the optimally tuned range separation parameter ω exhibited excellent performance in calculations of CT excitation. However, the optimally tuned range separation parameter in solid phase for long range corrected hybrid functional is quite small. Zheng et al. presented the optimally tuned screened range-separated hybrid (OT-SRSH) functional method, and gave the reasonable physical interpretation for optimally tuned range separation parameter [53]. In this work, the OT-SRSH functional method was applied with LC-PBE functional and 6-31G(d,p) basis sets for excitations, ESP and electronic coupling calculations. The following methods were adopted for optimally tuning ω [35,54,55]: are the HOMO level to the neutral donor, acceptor at neutral state and acceptor at anion state, respectively. The E ω D + , E ω A + , E ω A − are the single-point energy of the donor at cation state, acceptor at cation state, and acceptor at anion state, respectively. The E ω D and E ω A are the neutral state energy of the donor and acceptor, respectively. The optimally tuned ω is 0.155, 0.155, 0.167, 0.144, 0.155 Bohr −1 for PBDB-T, ITIC and IT-2X (X = F, Cl, Br), respectively.
In order to investigate the effect of acceptor halogenation on the electronic processes, the rate of CT, CR, and ED were calculated using Marcus' theory. The rate constant K can be calculated as following formula, where λ is the reorganization energy, including the internal-reorganization energy λ i and the external-reorganization energy λ ext , V represents the electronic coupling between the initial state and the final state, ∆G is the free energy change, K B is Boltzmann's constant, h is Planck's constant, and T is temperature (300 K was adopted). The λ i were estimated from the energies of donor and acceptor molecules [55,56]. The two-sphere model was adopted to calculate the λ ext [57]. The electronic couplings of CR processes were calculated using the Generalized Mulliken Hush (GMH) model [58,59], and the electronic couplings of CT and ED processes were calculated using two-state model [60]. The DFT and TDDFT calculations in this work were conducted using Gaussian 09 software [61]. The Multiwfn software was applied for the quantitative analysis of ESP, charge density difference (CDD), the transferred charges (∆q) and CT distance (∆d) [62][63][64].  Figure S1 in Supporting Material (SM), where the atomic serial numbers are also labeled. Table 1 lists the bond lengths, bond angles, and dihedral angles of the optimized PBDB-T, ITIC and IT-2X (X = F, Cl, Br) molecules, which can describe the distances and orientations among different moieties. More detailed geometric data are listed in Table S1 in SM. The dihedral angles between T 1 and BDD/BDT moieties in PBDB-T indicate the non-planar main chain and significant torsion among these fragments, which result from the steric hindrance of alkyl in side chain and dione in BDD moiety. For ITIC and IT-2X (X = F, Cl, Br), the DCI and DID are planar fragments, cyano groups are coplanar with DCI, and small torsion angles between DCI and DID (3.8 • , −7.6 • , −7.8 • , 3.2 • in ITIC and IT-2X (X = F, Cl, Br), respectively) indicate their quasi-planar structures, meaning that introducing halogen atom into DCI fragments at end sites cannot significantly affect the planarity of backbone. The molecular backbone planarity can promote electron delocalization and enhance intramolecular CT. Moreover, the intramolecular hydrogen bonds between the O in DCI and the adjacent H in DID enhance the backbone planarity of ITIC and IT-2X (X = F, Cl, Br). The average C-X (X = H, F, Cl, Br) bond lengths are 1.09, 1.34, 1.74, and 1.89 Å, respectively. The increase of C-X bond lengths is caused by the size of the halogen atoms. Hence, considering the end sites halogenation and the variations of torsion angles between DCI and DID fragments, the halogenation of ITIC cannot significantly reduce backbone planarity and increase steric hindrance. However, due to the strong electronegativity of halogen atom, halogenation can induce stronger inter-molecular hydrogen bonds, which are favorable for molecular packing order, and then for efficient charge transport/transfer processes.  Table 2 shows the energies of HOMO, the lowest unoccupied molecular orbital (LUMO) and other MOs, as well as the HOMO-LUMO gap (H-L gap ) of PBDB-T, ITIC, and IT-2X (X = F, Cl, Br) molecules. The data indicate that the fluorination, chlorination, and bromination downshift the HOMO about 0.03, 0.04, and 0.05 eV, respectively, and reduce LUMO about 0.02, 0.08, and 0.08 eV, respectively, corresponding to slight variation of H-L gap , agreeing experimental tendency of HOMO and LUMO variations by halogenation that were estimated by using electrochemical cyclic voltammetry measurements [40]. Also, halogenation of ITIC cannot effectively modify the energy differences between HOMO and HOMO-1, as well as between LUMO and LUMO+1, which were introduced as descriptors in machine-learning study of OPV materials [65]. Compared with ITIC, the downshift of HOMO and LUMO energies of IT-2X (X = F, Cl, Br) can be understood from the electronegativity of halogen atoms [66,67]. The strong electronegativity of halogen atoms causes more positive charges populated on the phenyl ring that bonded with halogen atoms, resulting in more electron localization and polarization effects [31]. Meanwhile, strong electronegativity of the halogen atom also changes the orbital hybridization of carbon atoms that bind with halogen atoms. The LUMO energy offset between donor and acceptor is usually considered as driving force for ED, while the energy offset between the LUMO of acceptor and the HOMO of donor is regarded as the driving force for CR. In terms of the data listed in Table 2, the halogenation of ITIC increase ED driving force and reduce CR driving force.   It was observed that the electron donor PBDB-T film exhibited good complementary absorption with those of acceptor ITIC and IT-2X (X = F, Cl, Br) [40]. PBDB-T has a good absorption characteristic in the region of 450-650 nm with absorption peaks at 615.5 nm and 567.5 nm, respectively. The absorption ranges of ITIC and IT-2X (X = F, Cl, Br) are between 600-815 nm. The complementary absorption of donor and acceptor in organic heterojunction are favorable to improve J sc by enhancing photon harvest [68]. The simulated absorption spectrum presented in Figure 2 well reproduced experimental results, including absorption profile, the wavelength at absorption maxima (λ max ) and relative shift among IT-2X (X = F, Cl, Br). For instance, the λ max of IT-2F is blue shifted 13.5 nm relative to those of IT-2Cl and IT-2Br, agreeing well with that of experimental 14 nm excellently [40].  Table 3 lists the main transition configurations, excited states characters (ESC), oscillator strengths (f ), excitation energies and wavelengths of the selected excited states. More detailed data of PBDB-T, ITIC and IT-2X (X = F, Cl, Br) are given in Table S2 in SM. The related MOs are presented in Figure 3, and more MOs are shown in Figure S2 in SM. The HOMO of PBDB-T is mainly contributed by BDT and thienyl fragments, whereas the LUMO is mainly composed of BDD, thienyl and BDT fragments. Hence, in terms of transition configurations and MOs, the main character of PBDB-T's S1 state is hybridized of CT and LE that is localized in donor fragment BDT (DLE). The PBDB-T's other excited states that are listed in Table 3 are CT or hybridizations of CT and DLE. Considering the f of PBDB-T's S1 state, the BDT and thienyl fragments are effective moieties in determining its light harvest efficiency. For ITIC and IT-2X (X = F, Cl, Br), their HOMOs are mainly contributed by DID fragment, while their LUMOs extend to DCI fragments. According to S1 transition configurations of ITIC and IT-2X (X = F, Cl, Br), the overlap between HOMO and LUMO means the LE of DID, which belongs to DLE. However, the relocation between HOMO and LUMO suggests CT excitations. Hence, the S1 states of ITIC and IT-2X (X = F, Cl, Br) hybridize the CT and DLE. The other excited states listed in Table 3 are also hybrid excitations. Since LE usually corresponds to large f that determine light harvest efficiency, DID is an effective chromophore for light harvest of ITIC and IT-2X (X = F, Cl, Br). CT excitations are also expected for increasing current density, efficient ED, and charge transfer/transport. The hybrid excitations combine the merits of CT and LE. Furthermore, halogenation doesn't significantly influence on excitation properties (transition configurations and ESC, etc.) of IT-2X (X = F, Cl, Br) except the slight change of excitation energies. The S1 excitation energies indicate that fluorination slightly increases excitation energy, while the chlorination and bromination slight reduce excitation energies. This can be understood from the transition configurations and halogenation effects on frontier MO energies. Moreover, there are some quasi-degenerated excited states (see Table  S2 in SM), for instance, S11 and S12, S13, and S14 of IT-2X (X = F, Cl, Br). The CT of these excited states can be enhanced due to coherence among quasi-degenerated states.   The CDD, CT distances ∆d and the transferred charges ∆q for the selected lowlying excited states of ITIC and IT-2X (X = F, Cl, Br) are shown in Figure S3 in SM. The interpenetrating of electron density increment/decrement regions and small ∆q support these excited states are hybrid excitations of CT and LE. Moreover, the CDD of ITIC and IT-2X (X = F, Cl, Br) suggest that the DCI and DID fragments are electron acceptor and donors, respectively, and the phenyl moieties in side chains also play the role of electron donor. In addition, the A-D-A type structure and symmetry of ITIC and IT-2X (X = F, Cl, Br) result into small ∆d of the excited states presented in Figure S3. Figure 4 shows the ESP of ITIC, IT-2X (X = F, Cl, Br), and PBDB-T molecules. For PBDB-T, the low ESP is mainly contributed by backbone, and the high ESP is dominated by side chains. For ITIC and IT-2X (X = F, Cl, Br), the high ESP are mainly contributed by alkyl in side-chains, and the low ESP mainly concentrate on O, S, and N atoms in backbone. Moreover, halogenation reduces the ESP at substitution sites. The average ESP values are −48.6, 63.8, 85.3, 105.9 and 106.4 Kcal/mol for PBDB-T, ITIC, IT-2F, IT-2Cl, and IT-2Br, respectively. Since the increased average ESP value of electron acceptor can increase the electron accepting ability [69], the halogenation of ITIC promotes their electron accepting ability, which is favorable for enhancing CT between donor-acceptor complexes [32]. Furthermore, halogenation of ITIC also generate the enlarged ESP differences between donor and acceptors, and further induce stronger interface electric field that can promote the ED and CT at the donor-acceptor interface [33]. Moreover, the enlarged ESP difference between donor and acceptors may also suppressing non-radiation recombination (which was regarded as one of the important energetic losses route in OPV) [70].

PBDB-T:ITIC and PBDB-T:IT-2X (X = F, Cl, Br) Properties
The optimized structures of PBDB-T:ITIC and PBDB-T:IT-2X (X = F, Cl, Br) complexes are shown in Figure 5, and the selected geometric parameters are given in Table 4. Compared with those of the corresponding geometric parameters of the ITIC and IT-2X (X = F, Cl, Br) molecules, the changes in bond lengths and angles are smaller than 0.001 Å and 1 • for the corresponding electron acceptors in complexes, respectively. The data in Table 4 indicate that the non-bond intermolecular interaction reduce the backbone planarity of ITIC and IT-2X (X = F, Cl, Br) in complexes. The distance between the mass centers of donor and acceptor in PBDB-T:ITIC and PBDB-T:IT-2X (X = F, Cl, Br) complexes are 4.32, 3.03, 5.86, and 3.19 Å, respectively. The maximum dihedral angle variations in complexes locate at the torsion between the donor DID and acceptor DCI fragments in ITIC and IT-2X (X = F, Cl, Br), and the corresponding values are 11 • , 19 • , 9 • , 13 • , respectively. Similarly, relative to those of isolated PBDB-T molecules, the changes of bond lengths and angles in complexes are very tiny. Due to the absence of intramolecular hydrogen bonds in PBDB-T, the PBDB-T's backbone is more flexible than ITIC and IT-2X (X = F, Cl, Br), corresponding to several remarkable dihedral angle changes of PBDB-T in complexes. The intermolecular non-bonding interactions can be evaluated from the binding energies (E b ) of the complexes. The calculated results show that halogenation increases E b (see Table 5), indicating that halogenation enhances the non-bonding intermolecular interaction between electron donor and acceptor, which is beneficial to improve the morphology regioregularity that was approved by relatively smooth surfaces of PBDB-T:IT-2X (X = F, Br) blended films [40]. The data of quadrupole moments of ITIC and halogenated systems (see Table S5 in SM) indicate halogenation results into an increase of the quadrupole, and then contribute to enhance quadrupole interaction between the donor and acceptor molecules in complexes [71]. Moreover, the large E b may enhance interface stability and improve resistance to thermal/light stress [72].   Table 5 also gives the HOMO and LUMO energies, as well as the H-L gap values of the PBDB-T:ITIC and PBDB-T:IT-2X (X = F, Cl, Br) complexes. Similar to that of ITIC and IT-2X (X = F, Cl, Br) molecules, the fluorination slightly elevates the LUMO of complex with broadening H-L gap , whereas the chlorination and bromination induce down-shift of LUMO of complexes, corresponding to H-L gap reductions. Among these complexes, the H-L gap of PBDB-T:IT-2Cl complex is the smallest. According to Shockley-Queisser theory, the lowering H-L gap tends to reduce energy loss [73]. Hence, chlorination and bromination may be helpful to suppress energy loss.   The absorption spectra of complexes that were presented in Figure 2b indicate that the halogenation not only induce the red-shift about 30-65 nm, but also enhance absorption in longer wavelength regions. The enhanced optical absorption is expected for prominent OPV materials since it can improve the utilization of photon energy from solar radiation.

Dihedral Angles PBDB-T:ITIC PBDB-T:IT-2F PBDB-T:IT-2Cl PBDB-T:IT-2Br
Organic materials with a smaller exciton binding energy (E exb ) usually show better charge separation efficiency, which impacts the free charge carrier generation in photon-toelectricity conversion of OPV. The calculated E exb of PBDB-T:ITIC and PBDB-T:IT-2X (X = F, Cl, Br) complexes are 0.63, 0.59, 0.62, and 0.60 eV, respectively, meaning that halogenation reduces E exb . The small E exb is favorable to reduce interfacial energy offsets toward efficient ED for free charge carrier generation in OPV [74][75][76]. In terms of the related parameters for calculating E exb listed in Table S6, the main reason for E exb reduction is halogenation increases electron affinity. Table 6 lists the excitation energies, corresponding oscillator strengths, excitation wavelengths, excited state characteristics and the main transition configurations of PBDB-T:ITIC and PBDB-T:IT-2X (X = F, Cl, Br) complexes. Details of more excited states are listed in Table S4 in SM. Figure 6 presents the selected MOs of PBDB-T:ITIC and PBDB-T:IT-2X (X = F, Cl, Br) complexes. More MOs are showing in Figure S4 in SM. According to transition configurations and MOs, the excited states' characters were assigned. The S1 states of all the complexes in this work are CT excitations, and the energies of local excitation in acceptor (ALE) and in donor (DLE) are higher than that of CT states, indicating energetic permission for ED. Comparing the CT state energy of PBDB-T:ITIC complex, the fluorination and bromination increase CT energy 0.11 and 0.02 eV, respectively, while the chlorination reduces CT energy 0.13 eV. Combined with the empirical formula of V OC [77], the large CT excitation energy is beneficial to improve the V OC . Meanwhile, the halogenation effects on CT excitation energies of these complexes are similar to that on H-L gap , which can be understood from transition configurations and weak intermolecular interaction between PBDB-T and acceptor in complexes. In addition, couple of excited states with higher energy of these complexes exhibit hybridization of CT and local excitation (see Table S4 in SM), which are favorable to charge carrier generation by hot ED processes [78]. Table 6. Electronic transition energies (in eV), corresponding oscillator strengths (f ), excitation wavelengths (in nm), excited states characters (ESC), and main transition configurations with coefficients larger than 10% for the selected excited states of PBDB-T:ITIC and PBDB-T:IT-2X (X = F, Cl, Br) complexes. (LC-PBE/6-31G**, the ω unit is Bohr −1 ; α = 0.200, β = 0.086; ε s = 3.5, ε d = 3.3).

States
Main  The ∆q and ∆d of the PBDB-T:ITIC and PBDB-T:IT-2X (X = F, Cl, Br) complexes are also calculated to quantitatively describe the CT and LE of the complexes. The CDD, ∆q and ∆d of the selected excited states are given in Figure 7 (more details to see Figure S5 in SM). Apparently, the ∆q labeled in Figure 7 and Figure S4 support the assignments of excited state characters. The halogenation of the complexes increases ∆q of CT excitations. The ∆d for most of excited CT states are in range of typical intermolecular distances between donor and acceptor in organic photovoltaic heterojunctions. Furthermore, the order of ∆q related to S1 states of the halogenated complexes is same as that of experimental J SC [40]. The ∆q was demonstrated as one of the critical factors to determine J SC [79].  Table 7 lists the calculated λ (including λ i and λ ext ), ∆G and V of the CT, CR and ED processes, which are required for calculating rate constants. Tables S6-S7 in SM provide the detailed parameters for calculating V. Since the halogenation site is at the end group of ITIC, the halogenation influences on λ are very tiny for each process (<0.01 eV). Meanwhile, the halogenation effects on ∆G are smaller than 0.1 eV for each process. As to V, apart from chlorination of ED process, the halogenation reduce V for each process. Compared with the rates of PBDB-T:ITIC complex, the fluorination increase CT rate and decrease CR rate. However, both the chlorination and bromination reduce CT rates and increase ED rates. For ED processes, the fluorination and chlorination decrease ED rates, while the bromination increases ED rate. In real OPV devices, CT, CR, and ED processes are competitive, resulting into ultimate performance. Considering the CT and ED rates are more than 13 orders larger than CR rates, meaning the efficient ED in these OPV devices. In the experiment, the ratio of the photocurrent density (J ph ) and its saturation (J sat ) is used to evaluate exciton dissociation efficiency and charge collection [40]. Charge collection is mainly determined by cathode/anode layers and their interfaces. Due to the same OPV architecture and experimental condition, the J sat must be same. Hence, the J ph /J sat ratios expose the ED efficiency. The J ph /J sat ratios of PBDB-T:ITIC, PBDB-T:IT-2F, PBDB-T:IT-2Cl and PBDB-T:IT-2Br were 93.3%, 92.2%, 90.6%, and 94.8%, respectively [40], agreeing well with the tendency of the calculated ED rates. Table 7. Calculated reorganization energies λ (in eV), free energy changes ∆G (in eV), electron coupling V (in eV) and rate constants K (in s −1 ) of the PBDB-T:ITIC and PBDB-T:IT-2X (X = F, Cl, Br) complexes for the three electron processes.

Conclusions
In this work, based on the quantum chemical calculations and the selected model systems which include donor PBDB-T, acceptors ITIC and IT-2X (X = F, Cl, Br), the halogenation mechanism and its effects on photovoltaic performances were studied. The main conclusions are as followings: Halogenation of electron acceptor ITIC in OPV heterojunction slightly affect the molecular geometric structures, MO energy levels, optical absorption spectra, exciton binding energies, and excitation properties. However, the halogenation of ITIC significantly enlarge the ESP difference between electron acceptor and donor PBDB-T with the order from fluorination and chlorination to bromination that can induce an electric field at the heterojunction interface, and the halogenation can also increase the transferred charges of CT states for the complexes, corresponding to the increase of J sc . Meanwhile, the halogenation effects on CT energies and electron processes rates depend on different elements. No matter which kinds of elements were introduced in the halogenation of the acceptor, the ED is always efficient in these OPV devices since the ED rates are more than 13-fold larger than CR rates. Hence, the main halogenation effects on OPV performances can be ascribed to its influences on ESP, transferred charges, and electron processes rates.
This work provides help in understanding the mechanism of halogenation influences on OPV performances of acceptor material properties, and is also conducive to the design of new high-performance acceptor materials with the aid of halogenation strategy.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.