Theoretical calculations. The chemical structures of the polymer donor PM6, the acceptor BO-4Cl and the additive MT are depicted in Fig. 1a. MT possesses a strong dipole moment of 1.73 Debye (Fig. 1b) due to the high atomic electronegativity of the oxygen atom, promoting the appearance of prismatic crystals (Figure S1). This characteristic facilitates robust intermolecular interactions with the adjacent molecules. Hence, we conducted density functional theory (DFT) calculations to determine the electrostatic potential (ESP) distribution of PM6, BO-4Cl and MT. The singlet point energy calculations were performed using the latest version of ORCA quantum chemistry software20 (version 5.0.4) with PBE0 functional and def2-SVP basis set. DFT-D3 with BJ-damping was applied to account for weak interactions and improve the accuracy of the calculations. As illustrated in Fig. 1c, the conjugated polymer backbones of PM6 exhibit significant negative charges on their surfaces. Conversely, most of the surface of BO-4Cl carries positive charges, except for the negatively charged regions located on the cyano groups. Due to the polarization effect, the oxygen atoms on the hydroxyl groups naturally carry negative charges. The opposite polarity of charge distribution indicates that the solid additive primarily interacts with the acceptor.
To confirm the presence of intermolecular interactions, we optimized the geometric structures of PM6, BO-4Cl and MT (Figure S2) using the GFN2-xTB method.21 Subsequently, we constructed assembled geometric structures of the BO-4Cl:MT combinations. Considering the charge characteristics and steric hindrance, MT exhibits a propensity to interact with the core units (referred to as type1) and the end groups (referred to as type2) of BO-4Cl (Fig. 1d). The binding energies between BO-4Cl and MT are calculated using the formula: \({{E}_{binding}={E}_{complex}-(E}_{BO-4Cl}+{E}_{MT})\). Detailed results of the calculations are provided in Table S1. The combination type2 exhibits a binding energy of -18.788 kcal mol−1, which is higher than the − 17.650 kcal mol−1 of the combination type1. The results indicate that MT tends to bind with the end groups of BO-4Cl, which may have an impact on the crystallization of acceptors. To gain further insight into the interactions between BO-4Cl and MT, we performed non-covalent interaction (NCI) analysis using the Independent Gradient Model (IGM) method.22 The resulting IGM diagram was rendered using VMD.23 Figure S3 illustrates the NCI graphs of the two BO-4Cl:MT combinations from different perspectives. The blue color represents strong O-H…N and O-H…S hydrogen bonding interactions, while the green represents weak van der Waals attractions and the red represents strong repulsions between molecules. The presence of blue color indicates the formation of strong non-covalent interactions between BO-4Cl and MT.
Molecular properties. The thermogravimetric analysis (TGA) shown in Fig. 2a indicates a complete and instant volatilization of MT at 150℃. While at a typical TA temperature of 100℃, a complete volatilization of MT was observed after one-hour treatment, indicating its favorable volatility. The Fourier transform infrared (FT-IR) results in Fig. 2b show that the characteristic absorption of hydroxyl groups of MT, peated at 3269 cm− 1, disappears after TA treatment at 50℃ for 10 min. After depositing the PM6:BO-4Cl/MT film, we observed no detectable signals for the hydroxyl group, providing the conclusive evidence for the complete removal of MT during film formation. The UV-vis absorption spectra of pure films and blend films were investigated in Fig. 2c and S4. The absorption peaks of PM6/MT film show no significant changes, suggesting negligible intermolecular interactions between donors and additives. The principal absorption peak of BO-4Cl/MT film displays a redshift of 4 nm compared to that of BO-4Cl film, indicating that the additives primarily influence the intermolecular π-π interaction of acceptors, which is consistent with the computational results as shown in Fig. 1. Meanwhile, the primary peak of PM6:BO-4Cl/MT blend film exhibits a redshift from 818 to 820 nm, probably resulting from the optimization of the donor/acceptor molecular packing in the presence of additives. The analysis of spectral data indicates that MT molecules mainly interact with acceptors and can be readily removed, highlighting its potential as a promising candidate additive for photovoltaic devices.
Photovoltaic performance and charge dynamics. Inspired with the above analysis, OSCs with a configuration of ITO/2PACz/active layer/PFN-Br/Ag were fabricated (Fig. 3a). Figure 3b presents the current density-voltage (J-V) characteristics of PM6:BO-4Cl devices with and without MT under one sun illumination (AM 1.5 G, 100 mW cm− 2), and the corresponding photovoltaic parameters are summarized in Table 1. The PM6:BO-4Cl control device exhibits a moderate PCE of 18.09% with an open-circuit voltage (VOC) of 0.853 V, a short-circuit current density (JSC) of 27.55 mA cm− 2 and a fill factor (FF) of 76.95%. We optimized the weight ratio (33/67/100 wt%) of MT additive to BO-4Cl acceptor. The results, as illustrated in Figure S5 and Table S2, demonstrate that the device processed with 67 wt% MT additive yields a champion PCE of 18.81%, accompanied by a VOC of 0.853 V, a JSC of 28.35 mA cm− 2 and an FF of 77.76%. As far as we know, the reported efficiency is among the highest for PM6-based binary OSCs processed using non-halogenated solvents (Table S3).18,19,24–30 The external quantum efficiency (EQE) spectrum of PM6:BO-4Cl/MT device shows an increased response in the region of 400–800 nm (Fig. 3c), corresponding to the enhancement of JSC. The JSC value calculated from the integration of the EQE spectrum was 27.64 mA cm− 2, which agreed well with the value obtained from the J-V curve within a deviation of 2.6%.
Table 1
Detailed photovoltaic parameters of OSCs under illumination of AM 1.5G, 100 mW cm− 2.
Conditions
|
VOC (V)
|
JSC (mA cm− 2)
|
JSC EQE(mA cm− 2)
|
FF (%)
|
PCEmax (%)
|
PCEavga (%)
|
PM6:BO-4Cl
|
0.853
|
27.55
|
27.07
|
76.95
|
18.09
|
17.71 ± 0.28
|
PM6:BO-4Cl/MT
|
0.853
|
28.35
|
27.64
|
77.76
|
18.81
|
18.52 ± 0.19
|
a Average data were obtained from 16 independent devices.
We plotted the photocurrent density (Jph) versus the effective voltage (Veff) curves to gain a deeper understanding of the exciton dissociation and charge extraction processes. Jph can be defined as Jlight – Jdark, where Jlight and Jdark are the current densities under illumination and in the dark, respectively. Veff can be defined as V0 – Vbias, where V0 is the voltage at Jph = 0 and Vbias is the applied voltage. The exciton dissociation probability (Pdiss) was calculated as Jph/Jsat under short-circuit conditions, where Jsat is the saturation photocurrent density. Similarly, the charge collection probability (Pcoll) was obtained from Jph/Jsat under the maximum power output condition. As shown in Fig. 3d and Table S4, the PM6:BO-4Cl/MT device exhibits higher Pdiss (98.44%) and Pcoll (90.25%) compared to the PM6:BO-4Cl device (97.57% and 89.51%, respectively), indicating more efficient exciton dissociation and charge collection.
The dependence of JSC and VOC on the illumination intensity (Plight) was measured to investigate the charge recombination processes. The relationship between JSC and Plight can be expressed as JSC ∝ Plightɑ, where ɑ represents the exponential factor obtained from the slope of the linear fitting. In Fig. 3e, the PM6:BO-4Cl/MT device exhibits a higher ɑ value (0.998) compared to the PM6:BO-4Cl device (0.992), suggesting a weaker bimolecular recombination. As shown in Fig. 3f, the slopes of VOC versus Plight decrease from 1.19 kT/q for the PM6:BO-4Cl device to 1.12 kT/q for the PM6:BO-4Cl/MT device (where k is the Boltzmann’s constant, T is the temperature, and q is the elementary charge), indicating the dominance of biomolecular recombination in the PM6:BO-4Cl/MT device. The deviation from the expected kT/q value may be attributed to trap-assisted recombination derived from morphological traps or interfacial defects. These findings suggest that both bimolecular recombination and trap-assisted recombination are significantly suppressed in the MT-assisted devices, accounting for the high JSC and FF values obtained. Besides, we also performed transient photocurrent (TPC) measurements to gain insights into the carrier dynamics process. As depicted in Fig. 3g, the PM6:BO-4Cl/MT device exhibits a shorter carrier lifetime (0.15 µs) compared to the PM6:BO-4Cl device (0.19 µs), indicating more efficient carrier extraction, which agrees well with the observed charge collection probability in the PM6:BO-4Cl/MT device.
We subsequently investigated the charge transport properties by using the space charge limited current (SCLC) method. As shown in Fig. 3h and Table S5, the PM6:BO-4Cl device exhibits a hole mobility (µh) of 5.19 × 10− 4 cm2 V− 1 s− 1 and an electron mobility (µe) of 1.33 × 10− 3 cm2 V− 1 s− 1, resulting in a µh/µe value of 0.39. Upon the addition of MT to the binary system, the PM6:BO-4Cl/MT device demonstrates an increased µh of 5.88 × 10− 4 cm2 V− 1 s− 1 and µe of 1.42 × 10− 3 cm2 V− 1 s− 1. The µh/µe value slightly increases to 0.41, indicating a more balanced charge carrier mobility. To evaluate the trap state density (Nt) in different devices, the dark current-voltage (I-V) curves of the electron-only devices are plotted in Fig. 3i. At low bias voltage, the current exhibits an ohmic response within the linear region. As the voltage increases, the current enters the intermediate region, where traps become filled by the injected free carriers. The calculated values of Nt are found to be 8.09 × 1015 cm− 3 and 6.12 × 1015 cm− 3 for the PM6:BO-4Cl and PM6:BO-4Cl/MT devices, respectively. The balanced charge transport and reduced trap state density indicate that MT effectively suppresses charge accumulation and recombination in the active layers, contributing to the improved JSC and FF in the MT-assisted devices.
Morphology characterization and working mechanism. To explore the effect of MT on the morphology, we examined the miscibility of the blend films through contact angle measurements (Figure S6). According to the Owens-Wendt-Rabel-Kaelble (OWRK) model, the surface tensions (γ) are calculated to be 33.12 and 38.35 mN m− 1 for PM6 and BO-4Cl films, respectively. Upon the MT addition, the corresponding γ values change to 29.86 and 43.42 mN m− 1 for the respective films. The interaction parameters (χ) are further calculated according to the Flory-Huggins model, which can be expressed as \({{\chi }}_{\text{1,2}}\propto {(\sqrt{{\gamma }_{1}}-\sqrt{{\gamma }_{2}})}^{2}\). The calculated parameters are summarized in Table S6 and S7. The χ parameter of the BO-4Cl:MT blend (χ = 0.335) is significantly lower than that of the PM6:MT blend (χ = 1.033), indicating that MT exhibits better miscibility with the BO-4Cl molecule. Without the addition of MT, the PM6:BO-4Cl binary blend shows a χ parameter of 0.192, which is lower than that of PM6:(BO-4Cl/MT) blend (χ = 0.696) and BO-4Cl: (PM6/MT) blend (χ = 0.531). The higher χ parameters observed in the MT-assisted blends suggest that MT reduces the miscibility between the donor and acceptor, promoting appropriate phase separation in the binary blend film.
The film surface morphology properties were investigated by the atomic force microscopy (AFM) measurements (as shown in Fig. 4a, 4b and S7). There is no noticeable change in the images of PM6 and PM6/MT films, indicating that MT has a negligible effect on PM6. In contrast, the BO-4Cl/MT film displays uniform pores due to the volatilization of MT, resulting in an increased Rq value from 0.41 to 10.40 nm. These results confirm the uniform mixing of MT and BO-4Cl, aligning with the findings from contact angle measurements. Meanwhile, the PM6:BO-4Cl/MT film exhibits a finer nanofiber network with an Rq value of 1.20 nm comparing with 0.88 nm of the PM6:BO-4Cl film, facilitating more efficient exciton dissociation and charge transport.
In addition, we studied the crystallinity and crystal orientation of the blend films using grazing incident wide-angle X-ray scattering (GIWAXS). The 2D-GIWAXS patterns and the corresponding line cuts are shown in Fig. 4c-4e. The calculated values of d-spacing and crystal coherence length (CCL) for both films are summarized in Fig. 4f and Table S8. The PM6:BO-4Cl film exhibits a lamellar packing (100) diffraction peak at qxy = 0.272 Å−1 in the in-plane (IP) direction and a strong π-π stacking (010) diffraction peak at qz = 1.725 Å−1 along the out-of-plane (OOP) direction, manifesting a preferential face-on orientation. Similarly, the PM6:BO-4Cl/MT film shows the (100) and (010) diffraction peaks at qxy = 0.271 Å−1 and qz = 1.732 Å−1, respectively. The calculated d-spacing values31 estimated from the (100) peaks are 23.100 and 23.185 Å for PM6:BO-4Cl and PM6:BO-4Cl/MT films, respectively. The slightly increased d-spacing values of the lamellar packing may be attributed to the steric hindrance effect of MT. Moreover, the d-spacing values of π-π stacking at (010) peaks decrease from 3.642 to 3.628 Å, implying a more compact molecular packing. Upon the MT addition, the CCL values significantly increase from 38.209 to 71.581 Å for the (100) peaks and from 21.179 to 30.567 Å for the (010) peaks, respectively. These higher CCL values suggest fewer imperfections and dislocations, demonstrating the positive effect of MT in inducing higher crystallinity in the blend films.
Based on the analysis of the interaction between MT and PM6:BO-4Cl heterojunction, a working mechanism of morphological evolution influenced by the volatile additive is proposed (Fig. 4g). During the spin-coating process of the photovoltaic solution onto substrates, the films undergo an initial nucleation process. MT molecules predominantly localize around acceptors and form strong intermolecular interactions via hydrogen bonds. As the solvent evaporates, MT molecules become trapped between acceptors, effectively hindering the disordered aggregation of BO-4Cl molecules and reducing charge recombination. Due to the strong dipole force, MT can induce the ordered molecular arrangement and better crystallization of acceptors, facilitating the long-range charge transfer. Upon complete removal of MT, the blend films exhibit a refined nanofiber morphology with optimal phase separation.
Large-area modules. Furthermore, we manufactured large-area organic photovoltaic modules to evaluate the potential of MT-assisted OSCs for large-scale production. The P1/P2/P3 laser etching method32 was used to connect the anode of one cell to the cathode of the next cell (Fig. 5a). Seven sub-cells are connected in series to form the monolithic modules. Figure 5b presents a photo of the module with a side length of 6 cm. The total illumination area is 19.31 cm2, including the photoactive area and the inactive interconnecting area. Due to the laser patterning process, the geometric fill factor can be up to 98%. As depicted in Fig. 5c, the module achieves a remarkable champion PCE of 15.74%, accompanied by a VOC of 5.945 V, a JSC of 3.66 mA cm− 2 and an FF of 72.39%, which is one of the highest values reported to date for large-area organic photovoltaic modules processed with non-halogenated solvents (Fig. 5d and Table S9).13,25,33–46 The PCE was also certified with 14.95% efficiency by a third-party accrediting institute (Figure S8). Subsequently, the shelf stability of the champion module was tested in inert atmosphere, it maintained 86% efficiency after 3840 hours aging (Figure S9).