The structural, electronic, and optical properties of PBDTTBT are comprehensively studied by density functional theory to rationalize the experimentally observed properties. Periodic boundary conditions method is employed to simulate the polymer block and calculate effective charge mass from the band structure calculation to describe charge transport properties. Moreover, both time-dependent density functional theory and a set of multidimensional visualization techniques are used to characterize the exciton dissociation ability in the PCBM: PBDTTBT interface. These theoretical methods and calculation techniques not only promote deep understanding of the connection between chemical structures and optical and electronic properties of the donor-acceptor system but also can be used to rationally design a novel donor-acceptor system. Based on the same calculated methods as PBDTTBT, four copolymers PBDTTTP, PBDTTTO, PBDTTTPD, and PBDTTFPD are designed to study their potentials as donors in polymer BHJ. The results indicate that PBDTTBT’s well conjugation benefits its good stability, and its wide and strong absorption spectra in the range of visible light, appropriate FMO levels, well charge transport, and favorable exciton dissociation lead to its photovoltaic performance. Furthermore, through comparing the four designed polymers with PBDTTBT, we conclude that the four designed polymers have stronger exciton dissociation ability and larger open-circuit voltage and external quantum efficiencies. Consequently, the four designed copolymers are promising candidates for polymer BHJ solar cells.
China Scholarship CouncilNatural Science Foundation of Yangtze Normal University2103XJQN007Ministry of EducationZ2015143Research Project from Chongqing Committee of EducationKJ150120201. Introduction
Solar cells are one critical technology for solving world energy needs. Traditional solar cells made from high-purity silicon have been commercialized, but their applications are finite because of high cost and weight. Polymer photovoltaic cells with a bulk heterojunction (BHJ) active layer have attracted attention due to their potentially low cost, lightweight, flexible, and easy manufacturing [1–4]. In the BHJ structure, the active layer that is sandwiched between anode and cathode constitutes of electron-donating conjugated polymer and electron acceptor, and its properties are the most determining factors in the whole performance of polymer BHJ [5].
In the past few decades, many progresses have been made in polymer BHJ which have reached power conversion efficiencies (PCE) of 11% [6]. The limiting factors of PCE are the stabilities, ability of harvesting the photon flux from the sun, effective exciton separation, and carrier mobilities. In order to improve the device performances, one can develop new device architectures [7, 8], synthesize new polymer donors [9, 10] and new electron acceptors [11–14], or work on both ends.
In recent years, developing low bandgap polymers has provided an alternative approach for achieving high PCEs resulting from their important roles in photoelectric conversion process: capturing solar photons from the sun, exciton generation and separation, carrier injection, hole transport, and affecting the size of the open-circuit voltage (Voc) and short-circuit current (JSC), thus affecting PCE (η=VOC×JSC×FF/Pin). In order to develop low-bandgap conjugated polymers, the most powerful strategy is to incorporate electron-rich donor segments and electron-deficient acceptor segments into the polymer backbone. Due to the push-pull interaction, efficient internal charge transfer (ICT) can take place from the donor (D) to the acceptor (A) upon photoexcitation, thus leading to a new absorption band at longer wavelengths [15] and producing appropriate molecular energy levels for high Voc, good charge transport, and high JSC, thus improving PCE [16].
In recent years, some benzo[1,2-b:4,5-b′]dithiophene-containing polymers have been applied to polymer BHJ and field-effect transistor (FET) [17–22] on account of its electron-rich structures. The relatively large and planar molecular structures of thiophene (TH) derivatives not only keep a high mobility but also help to promote cofacial π-π stacking. Given the symmetrical planar structure and high mobility of benzo[1,2-b:4,5-b′]dithiophene (BDT), Huo and his coworkers synthesized a polymer based on BDT, that is, PBDTTBT [23]. PBDTTBT was well soluble in common organic solvents due to its strong octyl chains and presented an excellent thermal stability with a decomposition temperature of 337°C without the protection of an inert atmosphere. The group fabricated a polymer solar cell device based on PBDTTBT donor and PC70BM (a fullerene derivative) acceptor in 2010. The structure of the device is ITO/PEDOT-PSS/polymer: PC70BM (1 : 2, w/w)/Ca/Al, and its Voc, JSC, fill factor (FF), and PCE are 0.92 V, 10.70 mA·cm−2, 57.5%, and 5.66%, respectively [23].
In order to rationalize the experimentally observed properties of known materials and predict those of unknown ones, theoretical investigations on the structure features, charge transport properties, and exciton dissociation abilities of PBDTTBT are indispensable. In the modeling process, we replace the octyls with butyls to reduce the computation cost. Analogously, we replace benzo[c][1,2,5]thiadiazole with four pyridazine derivatives and designed four polybenzo[1,2-b:4,5-b′]dithiophene derivatives, PBDTTTP, PBDTTTO, PBDTTTPD, and PBDTTFPD and study their potentials as donors in polymer BHJ based on PC70BM. Several parameters of determining the performance of solar cells, electronic and structural properties, open-circuit voltage (VOC), charge transfer properties, exciton dissociation abilities, and theoretical PCE have been investigated.
This paper’s distribution is as follows. In Section 2, we describe the computational methods to obtain electronic and structural properties. In Section 3, we detailedly describe the influences on the performance of solar cells by comparing PBDTTBT with designed copolymers, and the conclusion is in Section 4.
2. Calculation Methods
Density functional theory (DFT) has been broadly used to investigate the properties of organic compounds because its high accuracy is reasonable with the ab initio method and less computational time cost, and B3LYP, a hybrid functional, is widely used in calculating organic systems [24–27]. In this paper, DFT and time-dependent DFT (TD-DFT) [28] are employed to gain the qualitative properties of all the compounds at b3lyp/6-31g∗ level [29, 30]. All the optimized structures are the global minima on the potential energy surface. Periodic boundary conditions (PBC) [31] method is used to optimize block polymers at B3LYP/6-31G(d) and B3PW91/6-31G(d) levels [32]. All the calculations are implemented by Gaussian 09 package [33]. As to the copolymers, density topological analyses are examined by atom in molecule (AIM) analyses [34]. Moreover, the nucleus-independent chemical shift (NICS) [35, 36] is also calculated for copolymers at the B3LYP/6-31G∗ level, and NICS defined as the negative of the magnetic shielding at a ring critical point (RCP) is obtained from the AIM analyses. Furthermore, the bonding characteristics were also investigated by natural bond orbital (NBO) theory [37–40] and AOMix program packages [41].
3. Results and Discussion3.1. Geometry Optimization
The drafts of the studied copolymers are depicted in Figure 1, and the optimized structures of parent molecules of copolymers by B3LYP/6-31G(d) along with the abbreviations of each segments are shown in Figure 2. The selected bond lengths and bond angles of these copolymers are also given in Figure 1; herein, the C1-C2 (in PBDTTBT) or C2-C3 (in other copolymers) bond is defined as the central bond that connects the donor and acceptor. All the studied copolymers have the same central bonds (1.45 Å), and their dihedral angles are smaller than 37°, which suggest that all the polymers are rigid backbones. The nitrogen- (N-) hydrogen (H) or sulfur (S)/oxygen- (O-) nitrogen (N) interactions forming stable six or five-member rings reduce the dihedral angles and keep the molecular coplanarity, thus benefitting to the rigidity of the copolymers. Besides, the densities of bond critical points (BCPs) and RCPs as well as the distances between the two interactional atoms are also presented in Figure 1. It is worth noting that the angles of S1-C4-C5-C6 and C7-C8-C9-S2 in the four designed copolymers are smaller than C4-C5-C6-C7 and C8-C9-C10-C11 angles in PBDTTBT, which leads to better coplanarity than PBDTTBT.
Structural parameters of copolymers and locations for calculated NICS.
The stereograph of optimized parent molecules of copolymers by DFT//B3LYP/6-31G(d).
In order to acquire the charge population analysis on central bonds in studied polymers, we calculate the characters of bond critical points (BCPs) in central bonds and list the data in Table 1S in Supplementary Materials. The data indicate that the two atoms in central bonds are relatively accumulated due to sharing interactions.
Furthermore, the HOMO and LUMO diagrams of monomers in Figure 1S not only can qualitatively illustrate the electronic cloud distribution but also reflect the electron-donating and electron-accepting segments. We can clearly see from Figure 1S that the electronic cloud distribution of HOMO of all the monomers localizes near the polymeric axis and mainly localizes on BDT segments in the four designed copolymers, while that of LUMO localizes near the TBT, TTP, TTO, TTPD, and TFPD segments (the segments in rectangle in Figure 2). The molecular orbital diagrams illustrate that the BDT segments contributing to HOMOs in the five copolymers are as electron-donating segments, and the left segments (TBT, TTP, TTO, TTPD, and TFPD) mainly contributing to LUMOs are as electron-accepting segments. Moreover, in order to quantitatively view the components of the HOMOs and LUMOs, the density of state (DOS) and partial DOS (PDOS) of monomers are calculated and given in Figure 3. As shown in Figure 3, TBT segment contributes to both HOMO and LUMO; nevertheless, BDT mainly contributes to HOMO to some degree. The DOS diagrams for other polymers show that the HOMOs are mainly contributed by BDT segments, and the LUMOs are almost entirely from TBT, TTP, TTO, TTPD, and TFPD segments.
DOS and PDOS of monomers.
In a word, BDT as a donor mainly contributes to HOMO and the left segments mainly contribute to LUMO, and the push-pull interactions in D-A in the four copolymers are formed.
3.2. Conjugational Properties
Generally speaking, favorable structural stability of molecule originates from its favorable conjugated properties to some degree. In order to understand the structural stability of the studied polymers, we investigate their conjugation properties.
NICS is comprehensively used to express the aromaticity of molecules because it can clearly and simply monitor the condition of ring current. Aromatic systems have pretty negative NICS values, antiaromatic systems have strongly positive NICS values, and nonaromatic cyclic systems should have NICS values close to zero [35, 42–44]. The NICS values for the repeated units of polymers are calculated and listed in Table 2S, and the positions of all the rings in the molecules are shown in Figure 1.
In Table 1S, the ring NICS values along the polymeric axis (a, b, and c) are more positive than the individual TH (NICS, −13.9 at B3LYP/6-31G(d) level) and benzene (NICS, −9.7 at B3LYP/6-31G(d) level), which originates from the electrons delocalizing to the whole molecule, thus, forming conjugated systems. As to the same rings, the NICS values of center rings are more negative than the terminal rings, such as the NICS values in a2 and a3 rings are more negative than those in a1 and a4, which illustrate that the ring currents mainly accumulate in the molecular center. In addition, the electrons of the molecules delocalize from donors toward acceptors along the polymeric axis because there are push-pull interactions in the systems; therefore, a2, a3, a4, and c rings are full of ring currents but a1 and b rings lack ring currents. Further, inside rings, that is, d rings, where the NICS values are most negative in respective systems, since electron-rich sulfur and nitrogen atoms in the three systems provide large ring currents. Nevertheless, the carbon atoms in PBDTTTPD and PBDTTFPD are relatively electron-deficient in comparison of nitrogen atoms, which results in small ring currents, thus leading to more positive NICS values in d rings.
Through analyzing the ring currents in RCPs, we find PBDTTTP and PBDTTTO have the same conjugation as PBDTTBT and are relatively stable copolymers.
3.3. Absorption Spectra
As shown in Figure 4(a), the experimental absorption spectrum of PBDTTBT has three absorption bands from 300 to 700 nm both in chloroform and as a solid film, and the main absorption peak of PBDTTBT is located at approximately 581 nm in solution [23]. The simulated absorption spectrum of PBDTTBT monomer at the TD-DFT/B3LYP 6-31G(d) level is depicted in Figure 4(a) for comparison, and the effect of the solvent (chloroform) within polarizable continuum model (PCM) [44] is taken into account during the calculation. The main absorption peak of calculated PBDTTBT monomer is located at 607 nm. The agreements between measured and simulated spectra are well in overall spectral shape, though their main absorption peaks are different. Furthermore, the calculated electronic transitions, absorption wavelengths (λ), oscillator strengths (f), and main configuration of PBDTTBT monomer were shown in Table 1. Herein, the S1 ← S0 electronic transition mainly comes from HOMO to LUMO (absorption process). In the process of S1 ← S0 electronic transition, charges transfer from the conjugation section of homocyclic rings along the polymeric axis to TBT section as shown in the HOMO and LUMO diagrams of PBDTTBT in Figure 1S. In comparison of the charge difference densities of the main electronic excited states of PBDTTBT in Figure 2S, we find that the charge transfer takes place from BDT to TBT, whereas the alkyl groups hardly participate in charge transfer.
(a) Experimental absorption spectra of PBDTTBT as a film and in chloroform solution [23] and simulated absorption spectra of PBDTTBT monomer (the smaller one) in chloroform solution. (b) Absorption spectra of monomers of the four designed copolymers in chloroform solution.
Maximum absorption wavelengths (λ), oscillator strengths (f), the lowest excitation energies, and main configuration of monomers.
Monomer
λ (nm)
f
Eex (eV)
Main configuration
BDTTBT
607
0.81
2.04
HOMO → LUMO (99%)
BDTTTP
706
0.61
1.76
HOMO → LUMO (99%)
BDTTTO
718
0.63
1.73
HOMO → LUMO (100%)
BDTTTPD
538
0.96
2.30
HOMO → L + 1 (97%)
BDTTFPD
540
0.85
2.30
HOMO → L + 1 (83%)
At the same calculation level, the simulated absorption spectra of the four designed copolymers are given in Figure 4(b). As shown in Figure 4(b), the absorption spectra of monomers for BDTTTP and BDTTTO have two absorption bands in the range of 300–800 nm, especially the absorption bands (Q band) near to 700 nm, thus benefitting to harvesting influx photons from the sun, and the monomers of BDTTTPD and BDTTFPD have wide absorption bands in the range of 400–600 nm. The above results suggest that the parent molecules of designed copolymers have broad absorption in the visible region.
Moreover, the detailed maximum absorption wavelengths (λ), oscillator strengths (f), the lowest excitation energy (Eex), and the main configuration of monomers are given in Table 1. In Table 1, the maximum absorption peaks of BDTTTP and BDTTTO belonging to Q bands are red-shifted compared with BDTTBT and the Eex values decrease, thus benefitting electronic transition, whereas, these of BDTTTPD and BDTTFPD are blue-shifted. Generally speaking, the larger oscillator strengths (f), the larger electronic transition probability. As shown in Table 1, the oscillator strengths (f) of all the studied monomers are larger than 0.61, and that of BDTTTPD is the largest one (0.96) of all the f. Combining Figure 4 and Table 1, we find that the main configurations with the maximum wavelength of the parent molecules belong to single electron transitions and originate from HOMO to LUMO mainly assigned to π → π∗ transition; as shown in Figure 1S, the electrons transfer from donor (BDT) to acceptors (TBT, TTP, TTO, TTPD, and TFPD).
In comparison to the experimental and theoretical absorption spectra of BDTTBT, we find that the PBDTTBT has wide and strong absorption spectra. Moreover, the designed copolymers have wide and strong absorption, and they have smaller transition energy in the visible region. Therefore, these designed copolymers also have the favorable ability of harvesting the flux photons from the sun.
3.4. Frontier Molecular Orbitals
The properties of frontier molecular orbitals (FMOs) of polymers seriously affect stable and photovoltaic properties. In order to harvest the maximum of the photon flux from the sun and obtain high short-circuit current (JSC), the bandgap (Eg) of the polymers should lie between 1.3 and 1.9 eV [45]. Further, its HOMO energy level should be between −5.2 and− 5.8 eV if the donor can keep stable toward oxidation from air, meanwhile, its LUMO level should be between −3.7 and− 4.0 eV as shown in Figure 5. Rather, the open-circuit voltage (VOC) of PSC is eventually confined by the difference between the HOMO of the donor and the LUMO of the acceptor [12, 46]. It is useful to research the molecular FMOs because the relative levels of the occupied and virtual orbitals can provide reasonable qualitative indications for processes of exciton generation and dissociation.
Photoelectric conversion process diagram and energy levels of polymers. (exp represents the values of the experimental measures).
In the present work, we calculate the LUMO and HOMO energy levels of the five copolymers by exploring several DFT methods. Since only the values at B3LYP/6-31G(d) and B3PW91/6-31G(d) levels get close to the experimental values of PBDTTBT (the detailed calculated results are given in Table 3S in Supplementary materials); therefore, we only list the LUMO and HOMO energy levels based on the above methods in Figure 5.
In Figure 5, we can see that the value at the B3LYP/6-31G(d) level (−5.01 eV) is higher than the one at B3PW91/6-31G(d) level (−5.11 eV), and the two calculated values are higher than the experimental one (−5.31 eV). The difference between experimental and calculated values is about 0.2–0.3 eV, and the different environments between experiment and calculation should be responsible for the difference. As to the designed copolymers, the HOMO levels from B3LYP/6-31G(d) level are higher than the corresponding ones from B3PW91/6-31G(d) level as the PBDTTBT; moreover, all the HOMO levels of designed copolymers are lower than the one in PBDTTBT at the same calculated level; therefore, we predict that the experimental HOMO levels of the designed copolymers are lower than the experimental one of PBDTTBT (−5.31 eV), and the designed copolymers have stronger antioxidation properties than PBDTTBT.
In the formula of power conversion efficiency: η=VOC×JSC×FF/Pin, the three parameters, VOC, JSC, and FF, determine the solar cell performance directly. The parameter, open-circuit voltage (VOC) formed in the process of carrier transport, is popularly used to estimate the maximum PCE [9]. As shown in Figure 5, the experimental VOC (exp) is gained when the JSC is equal to zero in J-V curves. There are two models to describe theoretical VOC: one is the metal-insulator-metal (MIM) model [47], the other model is the HOMOD-LUMOA offset model [9, 48]. Moreover, Lo et al. suggest that Voc is described by a combined MIM model with HOMOD-LUMOA offset model; in that model, they find that when the work function deviation (Δϕelectrodes) of ITO from Al electrode is in the range− 3 and 0 eV, Voc enlarges linearly with Δϕelectrodes as prescribed by the MIM model. Outside this range Voc depends on HOMOD-LUMOA offset model [49]. In the work, we assume that Δϕelectrodes is outside the range −3 and 0 eV; thus, we employed the HOMOD-LUMOA offset model, and the Voc of a conjugated polymer-PC70BM solar cell can be estimated by [9]
(1)Voc=1eEDonorHOMO−EPCBMLUMO−0.3V,where e is the elementary charges, EPCBMLUMO is equal to −4.3 eV (PC70BM), and the 0.3 V is an empirical factor to offset the exciton binding energy [50]. In comparison with the experimental and simulated VOC values of PBDTTBT, the simulated VOC is lower than the measured one by 0.51 (at the B3LYP/6-31G∗ level) and 0.41 eV (at the B3PW91/6-31G∗ level), which primarily originates from the HOMO level difference between simulation and experiment. For the sake of comparing the VOC size of designed polymers with PBDTTBT, we also simulate the VOC values of designed polymers by HOMOD-LUMOA offset model. As shown in Figure 5, the Voc of the designed copolymers are larger than the one of PBDTTBT based on either B3LYP or B3PW91 functional resulting from the deeper HOMO levels of designed copolymers. Therefore, the four designed copolymers are promising to improve the Voc relative to PBDTTBT in the experiment.
3.5. Charge Transfer Properties
On the basis of band-like theory, the bandwidth (BW) and electron effective mass (m∗) are beneficial parameters for predicting the hole and electron-transporting ability of polymers [51–53]. The effective mass of carrier at the band edge representing mobility was gained as the square of ħ multiplied by the reciprocal of the curvature from Ek with k, and the formulation is defined as
(2)1m∗=1ℏ∂2Ek∂k2.
The kinetic model of mobility (μ) is given by the following formula:
(3)μ=eTm∗.
The BW and m∗ data are given in Table 2. According to the band theory, the wider the BW, the smaller the effective mass, and the larger the kinetic model of mobility [54]. The data in Table 2 investigate that PBDTTBT has wide valence band (0.25 eV) and small m∗ (−4.45 × 104). By using the same methods as PBDTTBT, we calculate the bandwidths and effective masses for designed polymers and list the data in Table 2.
Bandwidth (eV) and effective mass (m∗).
Polymer
Conduction band
Valance band
BW
m∗(×104)
BW
m∗(×104)
PBDTTBT
0.08
9.45
0.25
−4.45
PBDTTTP
0.13
5.82
0.38
−2.91
PBDTTTO
0.16
5.04
0.34
−2.91
PBDTTTPD
0.03
12.6
0.33
−2.91
PBDTTFPD
0.03
12.6
0.30
−3.29
The data in Table 2 investigate that all the designed polymers have wider valence bands (PBDTTTP: 0.38 eV, PBDTTTO: 0.34 eV, PBDTTTPD: 0.33 eV, and PBDTTFPD: 0.30 eV) and smaller m∗(PBDTTTP: −2.91 × 104, PBDTTTO: −2.91 × 104, PBDTTTPD: −2.91 × 104, and PBDTTFPD: −3.29 × 104) in comparison of PBDTTBT. Moreover, PBDTTTP and PBDTTTO also have wider conduction bands and smaller m∗ than PBDTTBT; however, PBDTTTPD and PBDTTFPD have narrower conduction bandwidths and valence bandwidths through comparing PBDTTBT. The above investigations indicate that all the designed polymers have better hole transport properties than PBDTTBT; further, PBDTTTP and PBDTTTO have better electron transport properties than PBDTTBT through conduction band properties.
3.6. Exciton Dissociation
After exciton formation, excitons transport to the donor-acceptor interface. In the interface, exciton dissociation competes with possible recombination. The rates of exciton dissociation and charge recombination are evaluated using the Marcus theory [55]:
(4)K=4π3h2λkBTVDA2exp−ΔG+λ24λkBT,where λ is the reorganization energy, VDA is the electronic coupling between donor and acceptor, ΔG is the free energy change for the electron transfer reaction, kB is the Boltzmann constant, h is the Planck’s constant, and T is the temperature (we set T = 298 K in our calculations). In the exciton dissociation and charge recombination, ΔG = ΔGCT and ΔG = ΔGCR, respectively. The electronic coupling can be calculated through employing the generalized Mulliken-Hush (GMH) model [56]. In the work, we calculate the ratio of exciton dissociation rate (KCT) and charge recombination rate (KCR); therefore, the parameters λ, ΔGCT, and ΔGCR are mainly studied.
The reorganization energy λ is composed of inner and outer reorganization energy. The inner originates from the change in equilibrium geometry of the donor (D) and acceptor (A) sites consecutive to the gain or loss of electronic charge upon electron transfer. The outer arises from electronic and nuclear polarization/relaxation of the surrounding medium. The inner is composed of two terms
(5)λ=λ1A+λ2D,(6)λ1A=EA−−EA,(7)λ2D=ED−ED+,
Here, EA− and EA are the energies of the neutral acceptor A at the anionic geometry and optimal ground-state geometry, respectively, and ED and ED+ are, accordingly, the energies of the radical cation at the neutral geometry and optimal cation geometry. In the condition of corrole-fullerene dyads in nonpolar solvent, the overall reorganization energy is 0.5 eV. Therefore, we assume here a value of 0.5 eV for the overall reorganization energy in our calculations. All the reorganization energies of researched polymers are given in Table 3. The reorganization energy of PBDTTBT, PBDTTTP, PBDTTTO, and PBDTTFPD is 0.69 eV except for 0.67 eV of PBDTTTPD. The ΔGCR can be estimated with
(8)ΔGCR=EIPD−EEAA,where EIPD and EEAA are the ionization potential of the donor and electron affinity of the acceptor, respectively. The calculated ΔGCR values are successively −1.76 eV, −1.92 eV, −2.10 eV, −2.09 eV, −2.15 eV for PBDTTBT, PBDTTTP, PBDTTTO, PBDTTTPD, PBDTTFPD shown in Table 3. As shown in Table 3, we also estimated the ΔGCT values in accordance with Rehm-Weller equation [57]
(9)ΔGCT=−ΔGCR−ΔE0−0−Ebwhere ΔE0−0 is the energy of the lowest excited state of free-base donor and Eb is the exciton binding energy. The calculated ΔGCT are depicted in Table 3, and they are −0.31 eV (PBDTTBT), −0.14 eV (PBDTTTP), −0.29 eV (PBDTTTO), −0.21 eV (PBDTTTPD), and −0.15 eV (PBDTTFPD). According to the estimated values of λ, ΔGCR, and ΔGCT, we calculate the ratio of exciton dissociation rate and charge recombination rate (KCT : KCR), and the data of KCT : KCR listed in Table 3 indicate that PBDTTBT has large KCT : KCR (1.62 × 106). Moreover, the designed polymers PBDTTTP, PBDTTTO, PBDTTPD, and PBDTTFPD have larger KCT : KCR values (2.36 × 107, 2.13 × 1011, 1.90 × 1011, and 2.70 × 1011, resp.) than PBDTTBT, which suggest that the designed polymers have better exciton abilities than PBDTTBT.
Reorganization (λ), exciton dissociation energy (ΔGCT), charge recombination energy (ΔGCR), the ratio of exciton dissociation rate and charge recombination rate (KCT : KCR), and external quantum efficiency (EQE).
Polymer
PBDTTBT
PBDTTTP
PBDTTTO
PBDTTTPD
PBDTTFPD
λ (eV)
0.69
0.69
0.69
0.67
0.69
ΔGCR (eV)
−1.76
−1.92
−2.10
−2.09
−2.15
ΔGCT (eV)
−0.31
−0.14
−0.29
−0.21
−0.15
KCT : KCR
1.62 × 106
2.36 × 107
2.13 × 1011
1.90 × 1011
2.70 × 1011
EQE
1.3%
4.1%
3.8%
3.7%
3.6%
3.7. External Quantum Efficiency
According to the studies of Scharber et al. in 2006 [9], the external quantum efficiency (EQE) of solar cell is related to the bandgap energy and frontier molecular orbital energy levels under the assumptions that any contribution to the short-circuit current from photons absorbed by the fullerene is neglected and the FF is set to 65%. The calculated EQE data of all the studied polymers are listed in Table 3. From Table 3, we can see that the EQE values of all the designed polymers are larger than 3.5% (especially PBDTTTP: 4.1%) and larger than that of PBDTTBT (1.3%).
In a word, the designed polymers, especially PBDTTTP, have better photovoltaic properties than PBDTBT; therefore, the designed polymers are the promising candidates for polymer BHJ.
4. Conclusions
In this contribution, density functional theory (DFT) and time-dependent DFT (TD-DFT) calculations have been employed to model PBDTTBT to rationalize the relationship between the experimentally observed properties and the structural features in polymer BHJ. Besides, according to PBDTTBT, we design four benzodithiophene-like copolymers, PBDTTTP, PBDTTTO, PBDTTTPD, and PBDTTFPD. In order to investigate their potentials as donors in polymer BHJ based on PC70BM, several parameters of determining the performance of solar cells, structural properties, absorption spectra, frontier molecular orbitals, charge transfer properties, and external quantum efficiency have been researched. PBDTTBT’s well conjugation benefits its good stability, and its well photovoltaic performance mainly results from its strong absorption spectra in the range of visible light, appropriate FMO levels, well charge transport, and favorable exciton dissociation. Comparing the four designed copolymers with PBDTTBT, we conclude that the four designed copolymers have stronger exciton dissociation ability and larger open-circuit voltage and external quantum efficiencies. Consequently, the four designed copolymers are promising candidates for polymer BHJ solar cells.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the Research Project from Chongqing Committee of Education (KJ15012020), “Chunhui Plan” Cooperating Project from the Ministry of Education (Z2015143), Natural Science Foundation of Yangtze Normal University (Grant no.2103XJQN007), and Chinese Scholarship Council. The authors are grateful to the anonymous referees for their suggestions.
Supplementary Materials
Table 1S: BCP properties, Wiberg bond index, and electronic configurations for bonding orbitals of central bonds. Table 2S: NICS for polymers at RCPs. Table 3S: HOMO and LUMO energy levels and bandgaps from the experiment and calculation for PBDTTBT. Figure 1S: HOMO (H) and LUMO (L) diagrams for monomers and dimers. Figure 2S: charge difference densities of the main electronic excited states of PBDTTBT. Figure 3S: orbital interaction diagrams (B3LYP/6-31G ∗) for dimers.
KrebsF. C.JorgensenM.NorrmanK.HagemannO.AlstrupJ.NielsenT. D.FyenboJ.LarsenK.KristensenJ.A complete process for production of flexible large area polymer solar cells entirely using screen printing—first public demonstration200993442244110.1016/j.solmat.2008.12.0012-s2.0-60449088983KrebsF. C.Polymer solar cell modules prepared using roll-to-roll methods: knife-over-edge coating, slot-die coating and screen printing200993446547510.1016/j.solmat.2008.12.0122-s2.0-60449105587KrebsF. C.GevorgyanS. A.AlstrupJ.A roll-to-roll process to flexible polymer solar cells: model studies, manufacture and operational stability studies200919305442545110.1039/b823001c2-s2.0-67649388867KrebsF. C.All solution roll-to-roll processed polymer solar cells free from indium-tin-oxide and vacuum coating steps200910576176810.1016/j.orgel.2009.03.0092-s2.0-67649173107ThompsonB. C.FréchetJ. M. J.Polymer–fullerene composite solar cells2008471587710.1002/anie.2007025062-s2.0-3754904353018041798LiangC.WangY.LiD.JiX.ZhangF.HeZ.Modeling and simulation of bulk heterojunction polymer solar cells20141274678610.1016/j.solmat.2014.04.0092-s2.0-84899824214HadipourA.de BoerB.WildemanJ.KooistraF. B.HummelenJ. C.TurbiezM. G. R.WienkM. M.JanssenR. A. J.BlomP. W. M.Solution-processed organic tandem solar cells200616141897190310.1002/adfm.2006001382-s2.0-33749189209KimJ. Y.LeeK.CoatesN. E.MosesD.NguyenT. Q.DanteM.HeegerA. J.Efficient tandem polymer solar cells fabricated by all-solution processing2007317583522222510.1126/science.11417112-s2.0-3444750494817626879ScharberM. C.MühlbacherD.KoppeM.DenkP.WaldaufC.HeegerA. J.BrabecC. J.Design rules for donors in bulk-heterojunction solar cells—towards 10% energy-conversion efficiency200618678979410.1002/adma.2005017172-s2.0-33645392241BundgaardE.KrebsF. C.Low band gap polymers for organic photovoltaics2007911195498510.1016/j.solmat.2007.01.0152-s2.0-34247359132WienkM. M.KroonJ. M.VerheesW. J. H.KnolJ.HummelenJ. C.van HalP. A.JanssenR. A. J.Efficient methano[70]fullerene/MDMO-PPV bulk heterojunction photovoltaic cells200342293371337510.1002/anie.2003516472-s2.0-004253048312888961BrabecC. J.CravinoA.MeissnerD.SariciftciN. S.FromherzT.RispensM. T.SanchezL.HummelenJ. C.Origin of the open circuit voltage of plastic solar cells200111537438010.1002/1616-3028(200110)11:5<374::AID-ADFM374>3.0.CO;2-WBackerS. A.SivulaK.KavulakD. F.FrechetJ. M. J.High efficiency organic photovoltaics incorporating a new family of soluble fullerene derivatives200719122927292910.1021/cm070893v2-s2.0-34347211016KooistraF. B.KnolJ.KastenbergF.PopescuL. M.VerheesW. J. H.KroonJ. M.HummelenJ. C.Increasing the open circuit voltage of bulk-heterojunction solar cells by raising the LUMO level of the acceptor20079455155410.1021/ol062666p2-s2.0-3384779725417253699CampbellA. J.BradleyD. D. C.AntoniadisH.Dispersive electron transport in an electroluminescent polyfluorene copolymer measured by the current integration time-of-flight method200179142133213510.1063/1.14061432-s2.0-0035473546ChengY. J.YangS. H.HsuC. S.Synthesis of conjugated polymers for organic solar cell applications2009109115868592310.1021/cr900182s2-s2.0-7324913888819785455HouJ.ChenH.-Y.ZhangS.ChenR. I.YangY.WuY.LiG.Synthesis of a low band gap polymer and Its application in highly efficient polymer solar cells200913143155861558710.1021/ja90649752-s2.0-7035065607719821569ChenH.-Y.HouJ.ZhangS.LiangY.YangG.YangY.YuL.WuY.LiG.Polymer solar cells with enhanced open-circuit voltage and efficiency200931164965310.1038/nphoton.2009.1922-s2.0-70350774564PanH.LiY.WuY.LiuP.OngB. S.ZhuS.XuG.Low-temperature, solution-processed, high-mobility polymer semiconductors for thin-film transistors2007129144112411310.1021/ja067879o2-s2.0-3424711194117362006PanH.WuY.LiY.LiuP.OngB. S.ZhuS.XuG.Benzodithiophene copolymer—a low-temperature, solution-processed high-performance semiconductor for thin-film transistors200717173574357910.1002/adfm.2006012422-s2.0-36849000763HouJ.ParkM.-H.ZhangS.ChenL.-M.LiJ.-H.YangY.Bandgap and molecular energy level control of conjugated polymer photovoltaic materials based on benzo[1,2-b:4,5-b′]dithiophene200841166012601810.1021/ma800820r2-s2.0-51549089043HuoL.HouJ.ChenH.-Y.ZhangS.JiangY.ChenT. L.YangY.Bandgap and molecular level control of the low-bandgap polymers based on 3,6-Dithiophen-2-yl-2,5-dihydropyrrolo[3,4-c]pyrrole-1,4-dione toward highly efficient polymer solar cells200942176564657110.1021/ma90129722-s2.0-69549103024HuoL.HouJ.ZhangS.ChenH.-Y.YangY.A polybenzo[1,2-b:4,5-b′]dithiophene derivative with deep HOMO level and its application in high-performance polymer solar cells20104981500150310.1002/anie.2009069342-s2.0-7664910142020101663LimaI. T.SousaL.FreitasR. d. S.Ribeiro JuniorL. A.de Sousa JúniorR. T.de Silva FilhoD. A.A DFT study of a set of natural dyes for organic electronics2017231234310.1007/s00894-017-3512-82-s2.0-8503513720529147781MohrT.AroulmojiV.RavindranR. S.MüllerM.RanjithaS.RajarajanG.AnbarasanP. M.DFT and TD-DFT study on geometries, electronic structures and electronic absorption of some metal free dye sensitizers for dye sensitized solar cells20151351066107310.1016/j.saa.2014.07.0942-s2.0-84907362406LiY.PulleritsT.ZhaoM.SunM.Theoretical characterization of the PC60BM:PDDTT Model for an organic solar cell201111544218652187310.1021/jp20406962-s2.0-80455129711XieX. H.ShenW.HeR. X.LiM.A density functional study of furofuran polymers as potential materials for polymer solar cells201334102995300410.5012/bkcs.2013.34.10.29952-s2.0-84889246847ParrG.YangW.1989Oxford, NY, USAUniversity PressBeckeA. D.Density-functional thermochemistry. III. The role of exact exchange19939875648565210.1063/1.4649132-s2.0-0000189651LeeC.YangW.ParrR. G.Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density198837278578910.1103/PhysRevB.37.7852-s2.0-0345491105KudinK. N.ScuseriaG. E.Linear-scaling density-functional theory with Gaussian orbitals and periodic boundary conditions: efficient evaluation of energy and forces via the fast multipole method20006124164401645310.1103/PhysRevB.61.164402-s2.0-0000154314PerdewJ. P.ZiescheP.EschrigH.1991BerlinAkademie VerlagFrischM. J.TrucksG. W.SchlegelH. B.ScuseriaG. E.RobbM. A.CheesemanJ. R.ScalmaniG.BaroneV.MennucciB.PeterssonG. A.NakatsujiH.CaricatoM.LiX.HratchianH. P.IzmaylovA. F.BloinoJ.ZhengG.SonnenbergJ. L.HadaM.EharaM.ToyotaK.FukudaR.HasegawaJ.IshidaM.NakajimaT.HondaY.KitaoO.NakaiH.VrevenT.MontgomeryJ. A.Jr.PeraltaJ. E.OgliaroF.BearparkM.HeydJ. J.BrothersE.KudinK. N.StaroverovV. N.KobayashiR.NormandJ.RaghavachariK.RendellA.BurantJ. C.IyengarS. S.TomasiJ.CossiM.RegaN.MillamJ. M.KleneM.KnoxJ. E.CrossJ. B.BakkenV.AdamoC.JaramilloJ.GompertsR.StratmannR. E.YazyevO.AustinA. J.CammiR.PomelliC.OchterskiJ. W.MartinR. L.MorokumaK.ZakrzewskiV. G.VothG. A.SalvadorP.DannenbergJ. J.DapprichS.DanielsA. D.FarkasÖ.ForesmanJ. B.OrtizJ. V.CioslowskiJ.FoxD. J.2009Wallingford, CT, USAGaussian, Inc.BaderR. F. W.199022Oxford, UKOxford University PressSchleyerP. v. R.MaerkerC.DransfeldA.JiaoH.van Eikema HommesN. J. R.Nucleus-independent chemical shifts: a simple and efficient aromaticity probe1996118266317631810.1021/ja960582d2-s2.0-001119049728872872ShenW.LiM.HeR.ZhangJ.LeiW.The electronic and structural properties of nonclassical bicyclic thiophene: monomer, oligomer and polymer200748133912391810.1016/j.polymer.2007.04.0452-s2.0-34249776373CarpenterJ. E.WeinholdF.Analysis of the geometry of the hydroxymethyl radical by the “different hybrids for different spins” natural bond orbital procedure1988169416210.1016/0166-1280(88)80248-32-s2.0-45449123708ReedA. E.CurtissL. A.WeinholdF.Intermolecular interactions from a natural bond orbital, donor-acceptor viewpoint198888689992610.1021/cr00088a0052-s2.0-0011083499FosterJ. P.WeinholdF.Natural hybrid orbitals1980102247211721810.1021/ja00544a0072-s2.0-33744614883ReedA. E.WeinstockR. B.WeinholdF.Natural population analysis198583273574610.1063/1.4494862-s2.0-36148995600GorelskyS. I.LeverA. B. P.Electronic structure and spectra of ruthenium diimine complexes by density functional theory and INDO/S. Comparison of the two methods20016351-218719610.1016/S0022-328X(01)01079-82-s2.0-0001694717FuY.ShenW.LiM.Theoretical analysis on the electronic structures and properties of PPV fused with. Electron-withdrawing unit: monomer, oligomer and polymer200849102614262010.1016/j.polymer.2008.03.0512-s2.0-43049103038von Ragué SchleyerP.JiaoH.van Eikema HommesN. J. R.MalkinV. G.MalkinaO. L.An evaluation of the aromaticity of inorganic rings: refined evidence from magnetic properties199711951126691267010.1021/ja9719135Rague von SchleyerP.JiaoH.What is aromaticity?199668220921810.1351/pac199668020209KroonR.LenesM.HummelenJ. C.BlomP. W. M.de BoerB.Small bandgap polymers for organic solar cells (polymer material development in the last 5 years)200848353158210.1080/155837208022318332-s2.0-49449095269KosterL. J. A.MihailetchiV. D.RamakerR.BlomP. W. M.Light intensity dependence of open-circuit voltage of polymer:fullerene solar cells20058612123509/1123509/3MihailetchiV. D.BlomP. W. M.HummelenJ. C.RispensM. T.Cathode dependence of the open-circuit voltage of polymer: fullerene bulk heterojunction solar cells200394106849685410.1063/1.16206832-s2.0-0345377522BauschlicherC. W.Jr.LawsonJ. W.Current-voltage curves for molecular junctions: the effect of substituents200775, article 11540610.1103/physrevb.75.1154062-s2.0-33847786828LoM. F.NgT. W.LiuT. Z.RoyV. A. L.LaiS. L.FungM. K.LeeC. S.LeeS. T.Limits of open circuit voltage in organic photovoltaic devices2010961111330310.1063/1.33603362-s2.0-77949713634BredasJ. L.BeljonneD.CoropceanuV.CornilJ.Charge-transfer and energy-transfer processes in π-conjugated oligomers and polymers: a molecular picture2004104114971500410.1021/cr040084k2-s2.0-1024422914815535639HaddonR. C.SiegristT.FlemingR. M.BridenbaughP. M.LaudiseR. A.Band structures of organic thin-film transistor materials19955101719172410.1039/jm9950501719CornilJ.CalbertJ. P.BrédasJ. L.Electronic structure of the pentacene single crystal: relation to transport properties200112361250125110.1021/ja005700i2-s2.0-003585741811456688ChengY. C.SilbeyR. J.da Silva FilhoD. A.CalbertJ. P.CornilJ.BrédasJ. L.Three-dimensional band structure and bandlike mobility in oligoacene single crystals: a theoretical investigation200311883764377410.1063/1.15390902-s2.0-0037460607ZhuJ. G.ZhengW. X.ZhengJ. G.SunX. S.WangH. T.2005Beijing, ChinaScience PressMarcusR. A.Electron transfer reactions in chemistry: theory and experiment (Nobel lecture)19933281111112110.1002/anie.1993111132-s2.0-33748216903VoityukA. A.Estimation of electronic coupling in π-stacked donor-bridge-acceptor systems: correction of the two-state model20061246, article 06450510.1063/1.21662322-s2.0-32644488347KavarnosG. J.TurroN. J.Photosensitization by reversible electron transfer: theories, experimental evidence, and examples198686240144910.1021/cr00072a0052-s2.0-0039972823