- 1Jiangsu Key Laboratory for Carbon-Based Functional Materials & Devices, Institute of Functional Nano & Soft Materials (FUNSOM), Soochow University, Suzhou, China
- 2King Mongkut's University of Technology Thonburi (KMUTT), Bangkok, Thailand
Dithienopicenocarbazole (DTPC), as the kernel module in A-D-A non-fullerene acceptors (NFA), has been reported for its ultra-narrow bandgap, high power conversion efficiency, and extremely low energy loss. To further improve the photovoltaic performance of DTPC-based acceptors, molecular engineering of end-capped groups could be an effective method according to previous research. In this article, a class of acceptors were designed via bringing terminal units with an enhanced electron-withdrawing ability to the DTPC central core. Their geometrical structures, frontier molecular orbitals, absorption spectrum, and intramolecular charge transfer and energy loss have been systematically investigated on the basis of density functional theory (DFT) and time-dependent density functional theory (TD-DFT) calculations. Surprisingly, NFA 4 highlights the dominance for its increased open circuit voltages while NFA 2, 7, and 8 exhibit great potential for their enhanced charge transfer and lower energy loss, corresponding to a higher short-circuit current density. Our results also manifest that proper modifications of the terminal acceptor with extensions of π-conjugation might bring improved outcomes for overall properties. Such a measure could become a feasible strategy for the synthesis of new acceptors, thereby facilitating the advancement of organic solar cells.
Introduction
Currently, organic solar cells (OSCs), recognized for their great potential as sustainable energy sources, have been attracting extensive attention from researchers (Brabec et al., 2001; Gao et al., 2018; Upama et al., 2020; Zhao et al., 2020). OSCs have displayed unparalleled properties like practically inexhaustible resource storage, flexible modification, rapid power payback time, and being free of pollution during synthesis and operation (Lewis, 2007; Cheng et al., 2009; Lin et al., 2012). Fullerene and its derivatives have primarily been employed as acceptors and their widespread use has made significant contributions to the improved performance of OSCs (Sonar et al., 2011; Zhang et al., 2015). Nevertheless, despite much effort being devoted toward the development of fullerene derivatives, there still remain several troublesome obstacles to be overcome, such as complex synthesis, week absorption abilities, poor light harvesting properties, and hardly tunable band gaps (Liang and Yu, 2010; Lin et al., 2012; Chen et al., 2015). In recent years, the above difficulties have given birth to non-fullerene acceptor materials, which offer several unique advantages, including low-cost facile synthesis, tunable flexibility of FMO energy levels, and superior capacity of light absorption in the visible region (Anthony, 2011; Ala'a et al., 2014; Li et al., 2015a; Lin et al., 2015). However, researchers are still persisting in advancing progress in the research of non-fullerene acceptors (NFA).
Among numerous NFAs, the A-D-A type acceptors have received considerable attention from researchers (Li et al., 2019). After abundant exploitation and rapid development, A-D-A acceptors have proven to be a wise design to date, and are likely to become potential alternatives to fullerenes (Kini et al., 2020). In order to optimize the photovoltaic properties of OSCs, indacenodithienothiophene (ITs) have been widely used as central “D” cores to fabricate some small molecule NFAs. A typical molecule is ITIC, which can be used in discussions on photovoltaic performances of different polymer donor materials (Shen et al., 2018). And within various polymers, the cell of PBDTS-TDZ:ITIC offered the PCE of 13.35%, which is the best record in ITIC-based devices, as is reported by Peng's group (Xu et al., 2018). Moreover, Hou's group developed a new acceptor IT-4F via adding the fluorine atom to the IC “A” units and achieved an increased efficiency of 13.7% (Zhao et al., 2017). Another promising central core is BTP, as Yuan et al. demonstrated a remarkable PCE of over 15% for the novel acceptor BTP-4F, which is structured by end capping the BTP central “D” core with the 2F-DCI “A” units (Yuan et al., 2019). Later, Cui et al. replaced the fluorine atom at the terminal of BTP-4F with the chlorine atom to form a new acceptor, namely BTP-4Cl, and obtained a high efficiency of 16.5%, which is the highest among reported OSCs (Cui et al., 2019).
Recently, Yao et al. have synthesized a novel non-fullerene acceptor consisting of a highly electron-rich core, Dithienopicenocarbazole (DTPC), and two terminal electron-withdrawing units, 2-(5,6-difluoro-3-oxo-2,3-dihydro-1H-inden-1-ylidene) malononitrile (DFIC), which is proven to obtain considerably narrow gaps of 1.21 eV, quite high PCE up to 10.21%, and exceedingly low energy loss for only 0.45 eV simultaneously (Yao et al., 2018). However, DTPC-based acceptors would seem to be defective if compared to acceptors based on the core of IT and BTP, which are able to perform PCE over 13%, and still have more room for further research on higher PCE. Fortunately, as with A-D-A type acceptors, the central “D” core and the terminal “A” units are of equal importance and both of them play an important role in the overall performance of NFAs (Dey, 2019). In other words, besides the central core, end-capped groups may also affect the overall properties of photovoltaic devices to a great degree, and proper collocation of core and end groups could bring significantly improved outcomes (Shen et al., 2018). For instance, by molecular engineering of one more malononitrile group at the terminal, the efficiency highlights a noteworthy increase from 4.11% of FBR to 8.4% of BAF-4CN (Holliday et al., 2015; Gupta et al., 2017). Additionally, previous research has suggested that the device performance of different NFAs would vary depending upon the types of end-capped units according to their electron-withdrawing capability (Singh and Prakash, 2019), and extending the π-conjugation and enhancing the electron-withdrawing ability of terminal groups may be effective and feasible to promote the photovoltaic properties (Song et al., 2020). So here we pay our attention to suitable molecular engineering of extended π-conjugation and enhanced electron-withdrawing ability at the terminal units to modulate photophysical performance of NFAs.
Herein, we express a keen interest in DTPC-based molecules and focus on expanding the electron-withdrawing ability of terminal units to investigate the photovoltaic properties. It could be reasonably assumed that such a modification may not only maintain the original forte of narrow band gaps and low energy loss benefitting from the electron-rich core, but also possess a higher open circuit voltage and short circuit current, corresponding to improved PCE. In this article, we take DTPC-DFIC for reference (NFA 1) to innovatively design a series of new non-fullerene acceptors by the addition of enhanced electron-withdrawing terminal units to a DTPC core, as is shown in Figure 1. By means of density functional theory (DFT) and time-dependent DFT (TD-DFT) computational methodology, geometrical structures, frontier molecular orbitals, and absorption spectrum are systematically investigated, then open circuit voltage, energetic driving force, and energy loss are reliably predicted. The simulation of intramolecular charge transfer, charge transport, and stability properties are also covered in this article. Exploration of the above physical parameters closely correlated with PCE permits us to fully understand how such measures influence the photovoltaic performance of the device, and to seek potential alternatives to previous DTPC-based acceptors.
Theoretical Methods and Computational Details
There is no doubt that side alkyl chains have a notable impact on molecular aggregation and geometrical structures, hence they were considered for the investigated NFAs during the calculations. And for the sake of obtaining as accurate results as possible, all tasks in this article were carried out with the 6-31G(d,p) basis set. In order to increase the reliability of calculated value, three DTF functionals with different HFexc were used to calculate the highest occupied molecular orbital (HOMO) level of NFA 1, including B3LYP, PBE0, and MPW1K (Becke, 1993; Stephens et al., 1994; Adamo and Barone, 1999; Lynch et al., 2000). As can be seen in Supplementary Table 1, the result (−5.36 eV) calculated by PBE0 function showed great accordance with experimental data (−5.31 eV), so DFT/PBE0 was selected to perform ground-state-related calculations (Yao et al., 2018). The cationic and anionic geometries were also optimized at the PBE0 level to further evaluate IPs, EAs, λe, and λh. The absorption spectrum of investigated NFAs were simulated by TD-DFT/CAM-B3LYP under the solvent condition of dichloromethane with polarizable continuum model (PCM), as CAM-B3LYP (19% HFexc at short-range and 65% HFexc at long-range) was proven to be a valid function that could exactly take the effect of charge transfer into account (Gross and Kohn, 1985; Yanai et al., 2004). The charge density difference (CDD), electron–hole distributions, and inter fragment charge transfer (IFCT) were implemented in Multiwfn 3.7 package to understand intramolecular electronic transition (Lu and Chen, 2012). Incidentally, the donor molecule PBT7-Th was simplified into repeated units (PBT7-Th)3 to save computational costs, as is shown in Supplementary Figure 1, and its geometrical structure was optimized at the B3LYP/6-31G(d,p) level and demonstrated in Supplementary Figure 2. All the calculations above were performed using the Gaussian 09 software package (Frisch et al., 2013).
Results and Discussion
Geometrical Structure
Without any structural constraints, all ground-state geometrical structures of NFAs are optimized and shown in Supplementary Figure 3. Taking NFA 1 and 4 as example, we show their optimized structure from the top view and side view in Figure 2. The dihedral angles and bond lengths between the DTPC core and end-capped acceptor are listed in Supplementary Table 2. Obviously, all A-π-A type NFAs based on the DTPC core are presented as quasi-planar conformations. All dihedral angles of investigated NFAs are extremely small, mostly close to zero. Taking molecular vibration when optimizing into consideration, the distortion could be overlooked, and excellent planarity is achieved between the core and acceptor. As is known, the core of DTPC plays a significant role in expanding conjugate planes of molecules when used for acceptors. Fortunately, for all A-D-A skeletal structures investigated, large planarity not only benefits from the DTPC core, but also relies on the coplanarity of the core and acceptor. Improvement of the planarity of NFAs can effectively narrow the energy gap and broaden the optical absorption range of NFAs, which also enhances the electron mobility of NFAs. After extending π-conjugation of terminal groups, the molecular structures still maintain superior conjugation and coplanarity, which could lead to narrow band gaps and excellent optical properties.
Bond length alternation (BLA) is a construct that could be applied to describe the conjugate degree of a molecule, which is defined as the average bond length difference between a single bond and double bond in a chain, as is calculated by the equation below (Meyers et al., 1994; Lee et al., 2001; Yang and Kertesz, 2006; Wang et al., 2014):
Table 1 shows the bond length of a, b, c, d, e, and BLA values for neutral, anion, and cation. Values of bond length range from 1.374 to 1.461 Å, which is in the middle of C–C bonds (1.54 Å) and C=C bonds (1.33 Å). Therefore, several bonds measured are not typical single or double bonds. BLA values of neutral are obviously larger than anion, but a little smaller than cation, indicating that all molecules are suitable for acceptors rather than donors. Due to good conjugation of electronic structure, NFAs prefer to be changed into an anionic state when delocalized. Further, larger BLA indicates larger differences between bonds and higher energy levels. The neutral BLA values of NFA 4 (0.0537 Å), 6 (0.0549 Å), and 7 (0.0523 Å) are a little higher than those of others, which could be conducive to raising their LUMO level. And the differences between neutral (0.0537 Å/0.0549 Å) and anionic state (0.0204 Å/0.0211 Å) of NFA 4 and 6 are slightly larger, which could make exciton delocalization better and further to make electron transport more favorable. Therefore, molecular engineering of terminal units will change the conjugate degree of the whole structure slightly, and some of these changes would have very positive results.
Frontier Molecular Orbital
It is essential to study the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) of any organic molecule, because frontier molecular orbital (FMO) exerts conclusive effects on open circuit voltages and plays an extremely important role in the property of excited-state and optical performance (Ku et al., 2011; Li, 2012; Janssen and Nelson, 2013). All FMO energy levels of investigated NFAs are depicted in Figure 3.
Clearly, all NFAs perform perfect adaptation to the donor PTB7-Th, ensuring minimal energy loss on electron transfer and as high an internal voltage as possible, and the introduction of an electron-withdrawing group changes their properties of energy levels greatly (Nielsen et al., 2015). Compared with NFA 1, NFA 2, 4, and 6–8 possess higher LUMO and HOMO energy levels, while NFA 3 and 5 possess narrower band gaps. Apparently, LUMO of NFA 3 is lower than that of 2 and LUMO of NFA 5 is lower than that of 7. By contrast, their band gaps show the opposite tendency. Therefore, we infer that the introduction of nitrogen, which is usually considered as an atom with strong electronegativity, may lower LUMO values, although it could lead to a decrease in the energy gap. With comparisons of NFA 7 and 8, it is reasonable to deduce that introducing fluorine and chlorine could slightly narrow the band gap, which is consistent with experimental results. Hence, for A-π-A skeletal NFAs, proper modifications of the terminal acceptor could be accompanied with relatively high LUMO, further effectively increasing the value of open circuit voltages. However, excessive π-conjugations at end-capped acceptors may not always be a good choice. Instead, it could conspicuously decrease the LUMO further to have a negative impact on open circuit voltages.
To more wholly understand the factors that cause the difference in FMO levels, all NFAs are divided into three fragments, a core of DTPC and two end-capped groups, to calculate respective contributions to HOMO and LUMO. The results are collected in Supplementary Tables 3, 4, and FMOs of all NFAs are illustrated in Figure 4.
Combining graphic analysis and numerical comparison, we conclude that for HOMO, DTPC-core takes hold observably because its percentage of contributions is more than 80% while other moieties only sum up to occupy ~15%. And outcomes present significant variations in how much the three fragments contribute to LUMO medially, as the percentages are all around 30–40%. So, we can reasonably consider that HOMO of NFAs mainly come from the DTPC core while LUMO of NFAs are distributed in the whole molecule. Additionally, from these FMO figures, we are able to observe that for all NFAs, both the HOMO orbital and the LUMO orbital present π and π* features, which will promote charge transfer. It is also verified that resulting from its large conjugate planar peculiarity, DTPC is a prominent selection to be utilized as the core of narrow-band gap NFAs due to its favorable abilities to withdraw electron and holes. Last but not least, with the exploration of FMO distribution, we testify that it is feasible to raise LUMO energy levels by adjusting suitable end-cap acceptors.
Open Circuit Voltages and Energetic Driving Force
As we know, maximizing open circuit voltages (VOC) properly is commonly applied to facilitate the efficiency of OSCs, thus VOC is an indispensable parameter to examine the performance of photovoltaic devices, which is usually calculated using the equation (Scharber et al., 2006; Rand et al., 2007; Credgington and Durrant, 2012):
where represents the HOMO energy level of the donor, while represents the LUMO energy level of acceptor, and ΔV is considered as exciton binding energy for charge separation. As for non-fullerene acceptors, it is empirically referred to as 0.5 eV (Pan et al., 2016; Yin et al., 2018).
Energetic driving force (ΔEL−L) is another parameter to estimate whether a molecule is ideal for application in the device or not. ΔEL−L is equal to the LUMO difference between donor and acceptor, which is intuitionally expressed by the formula (Faist et al., 2012):
With the assistance of FMO energy levels discussed above, we were able to obtain the two dominant factors of all NFAs, which are demonstrated in Table 2. Firstly, the calculated VOC of NFA 1 (0.83 eV) is quite close to the experimental measured value (0.77 eV), which attests that our consequences are reliable in accordance with actual research (Yao et al., 2018). Furthermore, most VOC values of NFAs studied are higher than that of NFA 1, except for NFA 3 and 5. These NFAs with prior VOC may have great potential for being good candidates. From Table 2, the values of ΔEL−L are all higher than 0.3 eV, which ensures an adequate driving force for separated charges to transfer between donor and acceptor (Zhan and Yao, 2016). However, it is not negligible that a superfluous driving force does not mean a better performance. Several previous pieces of research have substantiated that there will be excess energetic loss degenerating the efficiency of the device if the value is more than 0.8 eV (Li, 2012), so NFA 3 and 5 could not be regarded as advisable acceptors in view of their unsuitable driving force. To conclude, unfortunately on account of an inferior performance on open circuit voltages and energetic driving force, NFA 3 and 5 will be seldom mentioned below in this article. In addition, more π-conjugations at end-capped units does not mean high open circuit voltages, so in order to improve VOC, the balance should be found. And according to the results, the structural modification at terminal groups like NFA 2, 4, 6, and 7 may be admissible.
Table 2. Calculated open circuit voltages VOC (eV), energetic driving force ΔEL−L (eV) for investigated NFAs.
Absorption Spectra
As far as we know, in addition to open circuit voltage, another crucial parameter affecting the efficiency of OSCs is short-circuit current density (JSC), which is determined by photo-physical processes and photoexcitation properties to a considerable extent (Mühlbacher et al., 2006; Murphy and Frechet, 2007). In this case, it is a requisite for non-fullerene acceptors to possess wide absorption spectra matching well with that of donor molecules and solar spectrums (Nielsen et al., 2015; Sui et al., 2017). Herein, we further explored the absorbing ability of NFAs in order to seek the influence of substituting different π-conjugate units at the terminal group compared with parent molecules. Accordingly, we figure out the simulated absorption spectra in Figure 5 and summarized the excitation energy, maximum absorption peaks, largest oscillator strength, light harvesting efficiency, normalized integrated spectral area, and major configurations in Table 3. The normalized integrated spectral area is defined as the proportion of the integrated spectral area of other NFAs to that of NFA 1 in the whole absorption range, and light harvesting efficiency is calculated as ηλ = 1-10-f, where f corresponds to the highest oscillator strength (Song et al., 2020).
Table 3. The excitation energies (eV), associated wavelengths (nm), oscillator strengths, light harvesting efficiency, normalized integrated spectral area, and orbital contribution.
From Figure 4, the broad absorption of NFAs ranges from 300 to 900 nm and highly overlaps with the visible spectrum and near infra-red spectrum. Delightfully, all NFAs exhibit perfect matching for the donor due to wider absorption in the area of 600–800 nm, where PTB7-Th presents a weak absorption. Meanwhile, PTB7-Th makes up for the shortfall that NFAs have in deficient absorption in the region of 450–600 nm. Compared with the absorption wavelength of NFA 1, NFA 4, 7, and 8 show slight bathochromic shifts, while NFA 6 shows hypsochromic shifts of 14 nm. More importantly, the values of the normalized integrated spectral area show similar tendencies that most NFAs, other than NFA 6, have larger spectral areas which means an enhanced capacity for light absorbing. The oscillator strengths at the maximum absorption wavelengths are sorted in the order of 4 = 7 = 8 = 2 > 1 > 6, and, generally, higher oscillator strength is related to better absorption coefficients and transition probabilities. Hence, NFA 2, 4, 7, and 8 are proven to be hopeful for high-efficient acceptors in terms of their admirable absorbing ability The excitation energies are in the order of 6 > 1 = 4 > 7 > 8 > 2, excellently corresponding to the sequence of band gap value. In consequence, we deduce that the broader band gap of NFA 6 degenerates its optical performance, including absorption peak, spectral area, and oscillator strength. Inspection of Table 3 reveals that the strongest absorption peak is originated by S0→ S1 electronic transitions derived from HOMO to LUMO. Learning from the FMO discussion above, we conclude that HOMO mainly distributes in the core of the molecule while LUMO spreads in the whole molecule, and both exhibit π and π* orbital features. As a result, the transition nature of main absorption could be assigned to π-π** transition with local excitation and charge-transfer excitation.
Consequently, NFA 2, 4, 7, and 8 with intense absorption spectra and superior sunlight-absorbing ability tend to have enhanced short-circuit current density and be advantageous candidates, which could be due to adopting different terminal π-conjugation groups properly. It implies that suitable molecular engineering of end-capped units significantly affects the optical absorption of the designed acceptors and could be an effective strategy to enhance the short-circuit current density and further enlarge the photovoltaic energy conversion of the device. Unremarkable performance on molecular absorption and charge transfer ability of NFA 6 may be attributed to its special configuration, indicating that such a molecular π-conjugation at the end-cap group is inadaptable for this system, and NFA 6 would not be discussed in detail in this article.
Energy Loss
According to previous research, NFA 1 is considered as a non-fullerene acceptor with high power conversion efficiency and low energy loss. Evidently, energy loss (Eloss) is a vital factor to estimate the capacity of the device and it is encouraging to obtain an Eloss value as small as possible. For the sake of thorough analysis on the cell performance of designed acceptors, we calculate energy loss of investigated NFAs corresponding to the equation below (Veldman et al., 2009; Li et al., 2015b):
where VOC represents open circuit voltage (which has been discussed above) and represents a relatively smaller optical band gap of donor and acceptor. According to experimental results, NFA based on DTPC-core exhibits an extremely small band gap, thus we employed the optical gap of NFA and calculated it by the following formula:
where λabs,edge is equal to absorption edge wavelength, which is depicted in Figure 6, and parameters with respect to energy loss are collected in Table 4.
Compared with NFA 1, four designed molecules all possess smaller Eloss values owing to their lower and higher VOC. Indeed, Eloss of NFA 1 may be overrated as the calculated value (0.52 eV) is slightly higher than the experimental data (0.45 eV), indicating that other NFAs would demonstrate better performance in reality than in theoretical prediction (Yao et al., 2018). As a rule, any extra loss of energy is generally connected to thriftless driving force and inefficient charge transfer that affects the PCE of the device negatively. Hence, it may be promising to utilize the designed acceptors to improve the performance of OSCs due to the disappearance of limitation by large energy loss, and extended electron-withdrawing groups at the terminal may become favorable to decrease energy loss for non-fullerene molecules.
Intramolecular Charge Transfer
Photoexcitation always accompanies exciton dissociation to burst the restraints of Coulomb attraction (Lemaur et al., 2005; Clarke and Durrant, 2010). With the intention of characterizing exciton properties during the excitation process, we paid attention to the study of intramolecular charge transfer (Sui et al., 2019). Here we selected NFA 1 and 4 as examples and drew their hole and electron distribution for the S0→ S1 transitions in Figure 7A, then we partitioned it into three fragments and applied the inter fragment charge transfer (IFCT) method to calculate the actual amount of charge transfer, as is shown in Figures 7B,C. The results of other NFAs are available in Supplementary Figure 4, and concrete data associated with the excitation process is collected in Supplementary Table 5. From Figure 7A, it is visualized to discover that holes mainly concentrate on the center of the molecule while electrons localize along the whole molecule. To be exact, the terminal units are almost fully occupied by electrons. That is to say, during the process of photoexcitation, partial electrons are photoinduced and delocalized from the center of the molecule to end-capped group, leading to the generation of holes at the DTPC-core, revealing typical charge-transfer (CT) excitation feature. Meanwhile, some electrons are photoinduced locally at the center of the molecule, presenting localized excitation (LE) feature. Consequently, it is confirmed that photoexcitation of the designed molecule is displayed in the form of both charge-transfer excitation and localized excitation, or rather π-π** transition, which is in exact accordance with the analysis above. It can be accurate to verify the above-mentioned conclusion and further quantify the charge transfer from Figure 7C that fragments 2 (DTPC) provide nearly equivalent electrons for fragments 1 and fragments 3, ~0.18 electrons. This can be attributed to the A-D-A skeletal structures in which DTPC holds a strong electron-rich property, profiting from its large conjugate planes, while end-capped units are electron-withdrawing. As a result, we believe that electron-withdrawing substituents at the terminal group may significantly affect its interaction with DTPC-core in the excited states and change the possibility of the effective exciton dissociation.
Figure 7. Hole and electron distribution for the S0→ S1 transitions of NFA 1,4 (A), and heat maps associated with S0→ S1 transitions for hole–electron overlap in various fragment of NFA 1,4 (B,C).
Additionally, it is in our interest to deeply explore how the end-capped group affects the intramolecular charge transfer by charge density difference (CDD) maps. Firstly, we figured out CDD maps of NFA 1 and 4 between ground state and first excited state in Figure 8A, where the violet and turquoise colors stand for an increase and a decrease in electron density, respectively. Then we focused attention on the charge-transfer length (DCT) which represents the central distance between two regions of density decrease and increase, as is exhibited in Figure 8B. It is worth noting that, for the NFAs investigated, the barycenters of the positive electron density and negative electron density are both positioned at the core of DTPC owing to their central symmetric structure, making it difficult to describe the charge-transfer length exactly. Fortunately, there have been several previous researches selecting DCT of the half-molecule on behalf of the charge-transfer length, offering a novel approach to address the problem (Wang et al., 2015; Chen et al., 2017; Song et al., 2020). Utilizing their experience for reference, we chose half-molecules of NFA 1 and 4 as the subject and figured out their charge density differences map and charge-transfer length in Figures 8C,D. The CDD maps and DCT in half-molecules of other NFAs can be achieved in Supplementary Figure 5, and relevant values, including q, DCT, and variation in dipole moment Δμhalf, are listed in Table 5.
Figure 8. Charge density differences Δρ between ground state and first excited state and charge transfer distance DCT (red arrows) for NFA 1 and 4 (A,B) and half of NFA 1 and 4 (C,D).
Table 5. Transferred charge q(|e−|) for NFAs and qhalf(|e−|), charge-transfer distance Dhalf(Å), variation in dipole moments Δμhalf (D) for the transition S0→ S1 for half-NFAs.
There is no doubt that the designed acceptors all possess slightly larger values of q, DCT, and Δμhalf compared to NFA 1 from Table 5, proving the deduction that the capacity of charge transfer increases with the addition of electron-withdrawing ability at the terminal group. Particularly, profiting from enhanced intramolecular charge transfer, NFA 2, 7, and 8 may prominently overshadow other NFAs in the aspect of photovoltaic performance. To sum up, with additional photocurrent contributions, NFA 2, 7, and 8 may be good at the charge-transfer properties upon photoexcitation, thereby increasing the short circuit current, and introduction of π-conjugated units at the terminal group could become a considerable way to improve the JSC of a device.
Charge Transport and Stability Properties
For the purpose of estimating charge transport and stability properties, we summarize electron affinities (EAs), ionization potentials (IPs), reorganization energies of electron (λe) and hole (λh), and absolute hardness (η) in Table 6.
According to Marcus theory, the rates of charge transfer are greatly determined by λe and λh, which could be calculated by the equation below (Marcus, 1993; Köse et al., 2007):
where () is the energy of cation (anion) molecule that is calculated based on the optimized neutral molecule, and | | are the energy of the cation/anion/neutral that are calculated based on their optimized structure, respectively. Similarly, | is energy of the neutral molecule that is calculated based on the optimized cation (anion) molecule.
In most cases, absolute hardness serves as a prediction of stability properties, which can be predicted by the following formula (Cheung and Troisi, 2010):
It is widely known that lower ionization potentials and higher electron affinities contribute to better efficiency of electron injection and exciton separation (Liu et al., 2010; Gui et al., 2012). Compared to NFA 1, other NFAs possess lower values of ionization potentials and NFA 8 possesses higher values of electron affinity, revealing that NFA 8 may do well in the transport ability of holes and electrons. Synthesizing the consequence of electron and hole, it is manifest that other NFAs possess lower reorganization energies than NFA 1, which means less barriers in the charge transfer and higher rates of electron and hole transfer. From Table 6, it is common for the designed acceptors to show slightly smaller values of stability than that of NFA 1. Taking steric hindrances into consideration, we suppose that there is little difference in stability properties between NFA1 and designed acceptors. Overall, it suggests that the introduction of different terminal units is favorable for the enhancement of hole–electron binding and molecular ionization, and would not cause significant variation in the stability of the acceptors.
Conclusion
In this article, we design a class of novel non-fullerene acceptors based on dithienopicenocarbazole as the kernel module theoretically, and investigate the relevant qualitative parameters associated with JSC and VOC of these proposed molecules systematically. NFA 4 may be a good candidate due to its outstanding VOC and obviously increased LUMO energy levels, as well as its slightly intense light-absorption ability. Though NFA 2, 7, and 8 do not highlight their value of VOC like NFA 4, they exhibit superior optoelectronic properties of strong charge transfer and small energy loss, which corresponds to higher JSC, and could still be optimal acceptors with preferable photovoltaic performance. Therefore, we draw a conclusion that proper molecular modifications with the extension of electron-withdrawing units at the terminal group could make a conducive influence for OPV properties together with increased open circuit voltages, improved sunlight-absorbing ability, enhanced charge transfer, and decreased energy loss. It should also be noted that theoretical calculation is of enormous importance to molecule design to avoid unnecessary attempts. Finally, we believe this work could provide helpful guidelines for better and faster designing of DTPC-based photovoltaic molecules, and could shed light on the study of high-performance non-fullerene acceptors materials.
Data Availability Statement
All datasets generated for this study are included in the article/Supplementary Material.
Author Contributions
YL and HW conceived the idea and initiated this project. HW and JF wrote the manuscript. YL, NR, and HW contributed to fruitful discussions and supervision of the project. All authors discussed the results and commented on the manuscript.
Funding
This work was supported by the National Key R&D Program of China (Grant Nos. 2017YFB0701600 and 2017YFA0204800), National Natural Science Foundation of China (Grant Nos. 51761145013, 21703145, and 21673149), the 111 Project, Joint International Research Laboratory of Carbon-Based Functional Materials and Devices, the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD), and Collaborative Innovation Center of Suzhou Nano Science & Technology.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2020.580252/full#supplementary-material
References
Adamo, C., and Barone, V. (1999). Toward reliable density functional methods without adjustable parameters: the PBE0 model. J. Chem. Phys. 110, 6158–6170. doi: 10.1063/1.478522
Ala'a, F. E., Sun, J.-P., Hill, I. G., and Welch, G. C. (2014). Recent advances of non-fullerene, small molecular acceptors for solution processed bulk heterojunction solar cells. J.Mater.Chem. A 2, 1201–1213. doi: 10.1039/C3TA14236A
Anthony, J. E. (2011). Small-molecule, nonfullerene acceptors for polymer bulk heterojunction organic photovoltaics. Chem. Mater. 23, 583–590. doi: 10.1021/cm1023019
Becke, A. D. (1993). A new mixing of hartree–fock and local density-functional theories. J. Chem. Phys. 98, 1372–1377. doi: 10.1063/1.464304
Brabec, C. J., Cravino, A., Meissner, D., Sariciftci, N. S., Fromherz, T., Rispens, M. T., et al. (2001). Origin of the open circuit voltage of plastic solar cells. Adv. Funct. Mater. 11, 374–380. doi: 10.1002/1616-3028(200110)11:5<374::AID-ADFM374>3.0.CO;2-W
Chen, F. Y., Zhang, W. X., Tian, T. J., Bai, B. L., Wang, H. T., and Li, M. (2017). Substituent effect on intramolecular charge transfer of symmetric methoxy-substituted bi-1, 3, 4-oxadiazole derivatives. J. Phys. Chem. A 121, 8399–8407. doi: 10.1021/acs.jpca.7b08845
Chen, W. Q., Yang, X., Long, G. K., Wan, X. J., Chen, Y. S., and Zhang, Q. C. (2015). A perylene diimide (PDI)-based small molecule with tetrahedral configuration as a non-fullerene acceptor for organic solar cells. J.Mater.Chem.C 3, 4698–4705. doi: 10.1039/C5TC00865D
Cheng, Y. J., Yang, S. H., and Hsu, C. S. (2009). Synthesis of conjugated polymers for organic solar cell applications. Chem. Rev. 109, 5868–5923. doi: 10.1021/cr900182s
Cheung, D. L., and Troisi, A. (2010). Theoretical study of the organic photovoltaic electron acceptor PCBM: Morphology, electronic structure, and charge localization. J. Phys. Chem. C 114, 20479–20488. doi: 10.1021/jp1049167
Clarke, T. M., and Durrant, J. R. (2010). Charge photogeneration in organic solar cells. Chem. Rev. 110, 6736–6767. doi: 10.1021/cr900271s
Credgington, D., and Durrant, J. R. (2012). Insights from transient optoelectronic analyses on the open-circuit voltage of organic solar cells. J. Phys. Chem. Lett. 3, 1465–1478. doi: 10.1021/jz300293q
Cui, Y., Yao, H. F., Zhang, J. Q., Zhang, T., Wang, Y. M., Hong, L., et al. (2019). Over 16% efficiency organic photovoltaic cells enabled by a chlorinated acceptor with increased open-circuit voltages. Nat. Commun. 10:2515. doi: 10.1038/s41467-019-10351-5
Dey, S. (2019). Recent progress in molecular design of fused ring electron acceptors for organic solar cells. Small 15:1900134. doi: 10.1002/smll.201900134
Faist, M. A., Kirchartz, T., Gong, W., Ashraf, R. S., McCulloch, I., de Mello, J. C., et al. (2012). Competition between the charge transfer state and the singlet states of donor or acceptor limiting the efficiency in polymer: fullerene solar cells. J. Am. Chem. Soc. 134, 685–692. doi: 10.1021/ja210029w
Frisch, M., Trucks, G., Schlegel, H., Scuseria, G., Robb, M., Cheeseman, J., et al. (2013). Gaussian 09, Revision D.01. Wallingford, CT: Gaussian Inc.
Gao, W., Liu, T., Ming, R. J., Luo, Z., Wu, K. L., Zhang, L., et al. (2018). Near-infrared small molecule acceptor enabled high-performance nonfullerene polymer solar cells with over 13% efficiency. Adv. Funct. Mater. 28:1803128. doi: 10.1002/adfm.201803128
Gross, E., and Kohn, W. (1985). Local density-functional theory of frequency-dependent linear response. Phys. Rev. Lett. 55:2850. doi: 10.1103/PhysRevLett.55.2850
Gui, K., Mutkins, K., Schwenn, P. E., Krueger, K. B., Pivrikas, A., Wolfer, P., et al. (2012). A flexible n-type organic semiconductor for optoelectronics. J. Mater. Chem. 22, 1800–1806. doi: 10.1039/C1JM14089B
Gupta, V., Bagui, A., and Singh, S. P. (2017). Molecular engineering of highly efficient small molecule nonfullerene acceptor for organic solar cells. Adv. Funct. Mater. 27:1603820. doi: 10.1002/adfm.201603820
Holliday, S., Ashraf, R. S., Nielsen, C. B., Kirkus, M., Röhr, J. A., Tan, C.-H., et al. (2015). A rhodanine flanked nonfullerene acceptor for solution-processed organic photovoltaics. J. Am. Chem. Soc. 137, 898–904. doi: 10.1021/ja5110602
Janssen, R. A., and Nelson, J. (2013). Factors limiting device efficiency in organic photovoltaics. Adv. Mater. 25, 1847–1858. doi: 10.1002/adma.201202873
Kini, G. P., Jeon, S. J., and Moon, D. K. (2020). Design principles and synergistic effects of chlorination on a conjugated backbone for efficient organic photovoltaics: a critical review. Adv. Mater. 32:1906175. doi: 10.1002/adma.201906175
Köse, M. E., Mitchell, W. J., Kopidakis, N., Chang, C. H., Shaheen, S. E., Kim, K., et al. (2007). Theoretical studies on conjugated phenyl-cored thiophene dendrimers for photovoltaic applications. J. Am. Chem. Soc. 129, 14257–14270. doi: 10.1021/ja073455y
Ku, J., Lansac, Y., and Jang, Y. H. (2011). Time-dependent density functional theory study on benzothiadiazole-based low-band-gap fused-ring copolymers for organic solar cell applications. J. Phys. Chem. C 115, 21508–21516. doi: 10.1021/jp2062207
Lee, W. H., Lee, H., Kim, J. A., Choi, J. H., Cho, M., Jeon, S. J., et al. (2001). Two-photon absorption and nonlinear optical properties of octupolar molecules. J. Am. Chem. Soc. 123, 10658–10667. doi: 10.1021/ja004226d
Lemaur, V., Steel, M., Beljonne, D., Brédas, J.-L., and Cornil, J. (2005). Photoinduced charge generation and recombination dynamics in model donor/acceptor pairs for organic solar cell applications: a full quantum-chemical treatment. J. Am. Chem. Soc. 127, 6077–6086. doi: 10.1021/ja042390l
Lewis, N. S. (2007). Toward cost-effective solar energy use. Science 315, 798–801. doi: 10.1126/science.1137014
Li, C., Fu, H. T., Xia, T., and Sun, Y. M. (2019). Asymmetric nonfullerene small molecule acceptors for organic solar cells. Adv. Energy Mater. 9:1900999. doi: 10.1002/aenm.201900999
Li, H. Y., Hwang, Y. J., Courtright, B. A., Eberle, F. N., Subramaniyan, S., and Jenekhe, S. A. (2015a). Fine-tuning the 3D structure of nonfullerene electron acceptors toward high-performance polymer solar cells. Adv. Mater. 27, 3266–3272. doi: 10.1002/adma.201500577
Li, W. W., Hendriks, K. H., Furlan, A., Wienk, M. M., and Janssen, R. A. (2015b). High quantum efficiencies in polymer solar cells at energy losses below 0.6 eV. J. Am. Chem. Soc. 137, 2231–2234. doi: 10.1021/ja5131897
Li, Y. F. (2012). Molecular design of photovoltaic materials for polymer solar cells: toward suitable electronic energy levels and broad absorption. Acc. Chem. Res. 45, 723–733. doi: 10.1021/ar2002446
Liang, Y. Y., and Yu, L. P. (2010). A new class of semiconducting polymers for bulk heterojunction solar cells with exceptionally high performance. Acc. Chem. Res. 43, 1227–1236. doi: 10.1021/ar1000296
Lin, Y. Z., Li, Y. F., and Zhan, X. W. (2012). Small molecule semiconductors for high-efficiency organic photovoltaics. Chem. Soc. Rev. 41, 4245–4272. doi: 10.1039/c2cs15313k
Lin, Y. Z., Wang, J. Y., Zhang, Z. G., Bai, H. T., Li, Y. F., Zhu, D. B., et al. (2015). An electron acceptor challenging fullerenes for efficient polymer solar cells. Adv. Mater. 27, 1170–1174. doi: 10.1002/adma.201404317
Liu, C.-C., Mao, S.-W., and Kuo, M.-Y. (2010). Cyanated pentaceno [2, 3-c] chalcogenophenes for potential application in air-stable ambipolar organic thin-film transistors. J. Phys. Chem. C 114, 22316–22321. doi: 10.1021/jp1099464
Lu, T., and Chen, F. W. (2012). Multiwfn: a multifunctional wavefunction analyzer. J. Comput. Chem. 33, 580–592. doi: 10.1002/jcc.22885
Lynch, B. J., Fast, P. L., Harris, M., and Truhlar, D. G. (2000). Adiabatic connection for kinetics. J. Phys. Chem. A 104, 4811–4815. doi: 10.1021/jp000497z
Marcus, R. A. (1993). Electron transfer reactions in chemistry. Theory and experiment. Rev. Mod. Phys. 65:599. doi: 10.1103/RevModPhys.65.599
Meyers, F., Marder, S., Pierce, B., and Bredas, J. (1994). Electric field modulated nonlinear optical properties of donor-acceptor polyenes: sum-over-states investigation of the relationship between molecular polarizabilities (alpha., beta., and. gamma.) and bond length alternation. J. Am. Chem. Soc. 116, 10703–10714. doi: 10.1021/ja00102a040
Mühlbacher, D., Scharber, M., Morana, M., Zhu, Z., Waller, D., Gaudiana, R., et al. (2006). High photovoltaic performance of a low-bandgap polymer. Adv. Mater. 18, 2884–2889. doi: 10.1002/adma.200600160
Murphy, A. R., and Frechet, J. M. (2007). Organic semiconducting oligomers for use in thin film transistors. Chem. Rev. 107, 1066–1096. doi: 10.1021/cr0501386
Nielsen, C. B., Holliday, S., Chen, H.-Y., Cryer, S. J., and McCulloch, I. (2015). Non-fullerene electron acceptors for use in organic solar cells. Acc. Chem. Res. 48, 2803–2812. doi: 10.1021/acs.accounts.5b00199
Pan, Q. Q., Li, S. B., Wu, Y., Sun, G. Y., Geng, Y., and Su, Z. M. (2016). A comparative study of a fluorene-based non-fullerene electron acceptor and PC61BM in an organic solar cell at a quantum chemical level. RSC Adv. 6, 81164–81173. doi: 10.1039/C6RA08364A
Rand, B. P., Burk, D. P., and Forrest, S. R. (2007). Offset energies at organic semiconductor heterojunctions and their influence on the open-circuit voltage of thin-film solar cells. Phys. Rev. B 75:115327. doi: 10.1103/PhysRevB.75.115327
Scharber, M. C., Mühlbacher, D., Koppe, M., Denk, P., Waldauf, C., Heeger, A. J., et al. (2006). Design rules for donors in bulk-heterojunction solar cells—towards 10% energy-conversion efficiency. Adv. Mater. 18, 789–794. doi: 10.1002/adma.200501717
Shen, F. G., Xu, J. Z., Li, X. M., and Zhan, C. L. (2018). Nonfullerene small-molecule acceptors with perpendicular side-chains for fullerene-free solar cells. J.Mater.Chem.A 6, 15433–15455. doi: 10.1039/C8TA04718A
Singh, and Prakash, S. (2019). Impact of end groups on the performance of non-fullerene acceptors for organic solar cell applications. J. Mater. Chem. A 7, 22701–22729. doi: 10.1039/C9TA08620J
Sonar, P., Lim, J. P. F., and Chan, K. L. (2011). Organic non-fullerene acceptors for organic photovoltaics. Energy Environ. Sci. 4, 1558–1574. doi: 10.1039/c0ee00668h
Song, C. P., Wang, Z. X., Li, J., Chen, Y. N., Zhao, F., and Zhang, H. Y. (2020). Extension of π-conjugation and enhancement of electron-withdrawing ability at terminal indenedione for A-π-D-π-A small molecules for application in organic solar cells. Org. Electron. 81:105679. doi: 10.1016/j.orgel.2020.105679
Stephens, P. J., Devlin, F., Chabalowski, C., and Frisch, M. J. (1994). Ab initio calculation of vibrational absorption and circular dichroism spectra using density functional force fields. J. Phys. Chem. 98, 11623–11627. doi: 10.1021/j100096a001
Sui, M.-Y., Li, M.-Y., Sun, G.-Y., and Su, Z.-M. (2019). An updated strategy for designing non-fullerene acceptors by the lowest singlet and triplet states excitation: influence of periodical substitution from O, S, and Se to Te for BAE derivatives. J. Phys. Chem. C 123, 11397–11405. doi: 10.1021/acs.jpcc.9b00498
Sui, M. Y., Pan, Q. Q., Yin, H., Sun, G. Y., Geng, Y., and Su, Z. M. (2017). Design of hexabenzocoronene derivatives as non-fullerene acceptors in organic photovoltaics by bridging dimers and modulating structural twists. Solar RRL 1:1700060. doi: 10.1002/solr.201700060
Upama, M. B., Mahmud, M. A., Conibeer, G., and Uddin, A. (2020). Trendsetters in high-efficiency organic solar cells: toward 20% power conversion efficiency. Solar RRL 4:1900342. doi: 10.1002/solr.201900342
Veldman, D., Meskers, S. C., and Janssen, R. A. (2009). The energy of charge-transfer states in electron donor–acceptor blends: insight into the energy losses in organic solar cells. Adv. Funct. Mater. 19, 1939–1948. doi: 10.1002/adfm.200900090
Wang, G. Y., Kan, Y. H., Geng, Y., Duan, Y. A., Wang, L., Wu, H. Q., et al. (2014). Theoretical study on the transport properties of oligothiophene–diketopyrrolopyrrole derivatives: quinoidal versus aromatic. Theor. Chem. Acc. 133:1453. doi: 10.1007/s00214-014-1453-1
Wang, H. T., Liu, H. M., Bai, F. Q., Qu, S. N., Jia, X. S., Ran, X., et al. (2015). Theoretical and experimental study on intramolecular charge-transfer in symmetric bi-1, 3, 4-oxadiazole derivatives. J. Photochem. Photobiol. A Chem. 312, 20–27. doi: 10.1016/j.jphotochem.2015.07.006
Xu, X. P., Yu, T., Bi, Z. Z., Ma, W., Li, Y., and Peng, Q. (2018). Realizing over 13% efficiency in green-solvent-processed nonfullerene organic solar cells enabled by 1, 3, 4-thiadiazole-based wide-bandgap copolymers. Adv. Mater. 30:03973. doi: 10.1002/adma.201703973
Yanai, T., Tew, D. P., and Handy, N. C. (2004). A new hybrid exchange–correlation functional using the Coulomb-attenuating method (CAM-B3LYP). Chem. Phys. Lett. 393, 51–57. doi: 10.1016/j.cplett.2004.06.011
Yang, S. J., and Kertesz, M. (2006). Bond length alternation and energy band gap of polyyne. J. Phys. Chem. A 110, 9771–9774. doi: 10.1021/jp062701+
Yao, Z. Y., Liao, X. F., Gao, K., Lin, F., Xu, X. B., Shi, X. L., et al. (2018). Dithienopicenocarbazole-based acceptors for efficient organic solar cells with optoelectronic response over 1000 nm and an extremely low energy loss. J. Am. Chem. Soc. 140, 2054–2057. doi: 10.1021/jacs.7b13239
Yin, H., Sui, M. Y., Pan, Q. Q., Sun, G. Y., and Geng, Y. (2018). A series of bowl-shaped PDI dimers designed for organic photovoltaic cells through engineering N-annulated bridge towards potential alternatives of PDI bridged dimer acceptors. Dyes Pigm. 148, 394–404. doi: 10.1016/j.dyepig.2017.09.032
Yuan, J., Zhang, Y. Q., Zhou, L. Y., Zhang, G. C., Yip, H.-L., Lau, T.-K., et al. (2019). Single-junction organic solar cell with over 15% efficiency using fused-ring acceptor with electron-deficient core. Joule 3, 1140–1151. doi: 10.1016/j.joule.2019.01.004
Zhan, C. L., and Yao, J. N. (2016). More than conformational “twisting” or “coplanarity”: molecular strategies for designing high-efficiency nonfullerene organic solar cells. Chem. Mater. 28, 1948–1964. doi: 10.1021/acs.chemmater.5b04339
Zhang, Q., Kan, B., Liu, F., Long, G. K., Wan, X. J., Chen, X. Q., et al. (2015). Small-molecule solar cells with efficiency over 9%. Nat. Photonics 9, 35–41. doi: 10.1038/nphoton.2014.269
Zhao, C. C., Wang, J. X., Jiao, J. Q., Huang, L. J., and Tang, J. G. (2020). Recent advances of polymer acceptors for high-performance organic solar cells. J.Mater.Chem. C 8, 28–43. doi: 10.1039/C9TC05567C
Keywords: organic solar cell, molecular engineering, density functional theory, dithienopicenocarbazole, non-fullerene acceptors (NFAs)
Citation: Feng J, Wang H, Rujisamphan N and Li Y (2020) Theoretical Design of Dithienopicenocarbazole-Based Molecules by Molecular Engineering of Terminal Units Toward Promising Non-fullerene Acceptors. Front. Chem. 8:580252. doi: 10.3389/fchem.2020.580252
Received: 05 July 2020; Accepted: 28 September 2020;
Published: 05 November 2020.
Edited by:
Nawee Kungwan, Chiang Mai University, ThailandReviewed by:
Patchanita Thamyongkit, Chulalongkorn University, ThailandKen Sakata, Toho University, Japan
Copyright © 2020 Feng, Wang, Rujisamphan and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Youyong Li, eXlsaSYjeDAwMDQwO3N1ZGEuZWR1LmNu
†These authors have contributed equally to this work