Advancing organic photovoltaic materials by machine learning-driven design with polymer-unit fingerprints

The feature-selection and ML-model enhancement in PCE Prediction for OPV materials

In this study, we began with 220 molecular descriptors calculated via RDKit from the SMILES strings of each macromolecule. Facilitating the development of machine learning models by analyzing feature-property relationships is an essential initial step. To identify the macromolecule properties most closely related to power conversion efficiency (PCE), we employed a feature selection method. This combined the 220 RDKit descriptors with 17 additional microscopic properties, including normalized HOMO level, LUMO level, bandgap (\({E}_{g}\)), molecular weight (\({M}_{w}\)), and number-average molecular weight (\({M}_{n}\)), totaling 237 features.

Model learning complexity is influenced by the correlation strength⁶⁰ between features and the target property, with molecular properties proving to enhance model accuracy⁶¹. Feature selection is one of the key steps in machine learning, selecting features from all features that positively impact the learning algorithm, which will reduce the difficulty of learning tasks and make the model more interpretable. As an integrated learning method based on decision tree, Random Forest (RF) regression has significant advantages in feature importance assessment. Additionally, Lasso regression, a linear model that incorporates an L1 regularization term (i.e., the sum of the absolute values of the variable coefficients) to mitigate overestimation of model performance, was utilized for efficient variable selection. In this section, we utilized the RF regression and Lasso regression to find the optimal feature subset. Figure 3a depicts the feature importance ranking of RF regression and Lasso regression. The weight ranking of each feature in the RF regression model is shown on the left side of Fig. 3a (blue bar chart), ignoring features with weights below 0.005. Top features include MaxPartialCharge, \({E}_{g}\), -HOMO_p, and -LUMO_n. The features with non-zero coefficients in Lasso regression are sorted, as shown on the right side (red bar chart) of Fig. 3a. Notably, the important features identified by Lasso regression are largely consistent with those selected for the trained model. Separate feature importance ranking plots for the two regression methods are also presented in Supporting Information. Reducing redundancy in the feature set while retaining informative elements can improve machine learning performance and mitigate overfitting.

a Feature importance ranking of RF regression and Lasso regression. b SHAP dependence plot of top 20 features of RF model. c SHAP interaction evaluation plot of top 5 features for P-type and N-type OPVs.

Through comprehensive consideration of the screening results of feature importance by RF and Lasso regression, we selected 17 features from total 237 features, including MaxPartialCharge, \({E}_{g}\)_n, -HOMO_p, -LUMO_n, \({M}_{w}\), M, \({M}_{n}\), PEOE_VSA9, \({E}_{{\rm{LL}}}^{{\rm{DA}}}\), -HOMO_n, SMR_VSA10, NumHeteroatoms, FpDensityMorgan1, -LUMO_p, \({E}_{g}\)_p, \({E}_{{\rm{HL}}}^{{\rm{DA}}}\) and PDI (see Table 1). Among them, the primary features for P-type materials are -HOMO_p, -LUMO_p, \({E}_{g}\)_p, while N-type materials include MaxPartialCharge, \({E}_{g}\)_n, -LUMO_n, \({M}_{w}\), M, \({M}_{n}\), PEOE_VSA9, -HOMO_n, SMR_VSA10, NumHeteroatoms, FpDensityMorgan1, PDI and other 220 features. Additionally, \({E}_{{\rm{HL}}}^{{\rm{DA}}}\) represents the energetic difference between HOMO of donor and LUMO of acceptor, while \({E}_{{\rm{LL}}}^{{\rm{DA}}}\) quantifies the energetic difference between the LUMO of the donor and the LUMO of the acceptor. The SHAP (Shapley Additive exPlanations) method was then used to analyze these features’ contributions to the PCE prediction model. We use this method as a feature selection criterion and extracted the top 17 features that contribute the most to the model according to the SHAP value. Each feature’s SHAP value is shown in Fig. 3b, delineating the marginal contribution to the model’s output. The complete SHAP evaluation diagram is shown in Figure S4. Additionally, to further analyze the correlation between features of the donor and acceptor and the PCE in OPVs, the Pearson correlation coefficient is calculated to measure the linear relationship between the input feature and output feature. As shown in Figure S5, it is not difficult to find that the characteristic values strongly correlated with PCE include MaxPartialCharge, Eg_n, -LUMO_n, -HOMO_p, etc. The correlations between MaxPartialCharge, Eg_n, -LUMO_n, -HOMO_p and PCE are -0.29, -0.46, 0.43, 0.46, respectively, as depicted in Figure S5. This further verifies the evaluation results of the random forest regression model.

The atom with the maximum partial charge (MaxPartialCharge) was found to contribute the most to the model and had an inhibitory effect on PCE. MaxPartialCharge refers to the local accumulation of positive or negative charge due to the uneven distribution of electron density between atoms in a molecule. The presence of a high MaxPartialCharge indicates poor electronic delocalization and low conjugation degree within the molecule, resulting in inefficient charge transport and thereby reducing the PCE. Notably, the molecular orbital energy levels significantly affect PCE. Higher HOMO_p values positively correlate with the PCE predicted by the model, while lower LUMO_n values have a negative correlation. The PCE is defined by \(\frac{{J}_{{sc}}{V}_{{oc}}{FF}}{{P}_{{in}}}\), where \({P}_{{in}}\) is the incident illumination power. Also, according to the empirical equation of \({V}_{{oc}}\) = \({e}^{-1}(\left|{{E}_{{HOMO}}}^{D}\right|-\left|{{E}_{{LUMO}}}^{A}\right|)\)-0.3 V (where e is the elementary charge), the alignment of donor’s HOMO and acceptor’s LUMO levels is crucial for estimating PCE. Figure 4a, b shows that \({V}_{{oc}}\) increases with the deepening of HOMO level of polymer donor (E_HOMO^D) and decreases with the deepening of LUMO level (E_LUMO^A), consistent with the above conclusion. A trade-off between achieving a small energy loss (\({E}_{{loss}}\)) (i.e., a high VOC) and a high charge generation efficiency (η_EQE) in OPV devices, which means they often suffer much larger energy losses (0.5–1.0 eV) than inorganic PV devices (0.3–0.4 eV). However, the emergence of NFAs has circumvented the V_OC-η_EQE trade-off, enabling the attainment of a higher PCE with a much smaller \({E}_{{loss}}\). In OPV systems, a LUMO offset of ∼0.3 eV between the donor and acceptor is required to ensure efficient electron transfer and subsequent dissociation into free charge carriers. This generates a charge transfer (CT) state that consists of a hole on the HOMO of the donor and an electron on the LUMO of the acceptor, with the energy of this CT state (E_CT) usually smaller than that of the narrowest band gap (\({E}_{g}\)). The V_OC of the resulting device was further improved due to the greater energy difference between the HOMO of the donor and the LUMO of the acceptor. Additionally, the LUMO energy level difference between the donor and the acceptor (ΔE₁) decreases, which are beneficial for reducing energy loss and improving PCE (as shown in Fig. 4f). To sum up, optimizing molecular orbital energy levels will become a key step to affect the performance of OPV devices.

a \({V}_{{oc}}\) versus −HOMO. b \({V}_{{oc}}\) versus −LUMO. c \({J}_{{sc}}\) versus \({{E}}_{g}\). d \({J}_{{sc}}\) versus \({M}_{w}\). e The data distribution of Eg for the N-type OPVs dataset. f Schematic illustration of band gap alignment between donor materials and NFAs. In state-of-the-art polymer: NFA systems, the \({J}_{{sc}}\) is jointly contributed by the large-band-gap polymer donor and narrow-band-gap NFA with complementary optical absorption profiles. The green arrows show the transfer of electrons upon photoexcitation. The red arrows show the transfer of holes upon photoexcitation.

\({E}_{g}\) has a significant contribution to the in PCE and is negatively correlated with PCE⁶². Designing low-bandgap materials to match the solar spectrum is a common method to improve short-circuit current (\({J}_{{sc}}\)) and thus the PCE of OPV cells⁶³. Fig. 4c plots \({J}_{{sc}}\) as a function of the polymer Eg, showing that \({J}_{{sc}}\) tends to increase as \({E}_{g}\) decreases, since a narrower \({E}_{g}\) can harvest more energy from the sunlight. The use of narrow band gap NFAs can broaden the absorption spectrum of OPVs to the near-infrared region, reducing energy loss. This further validates that controlling the \({E}_{g}\) of the chemical structure within a relatively small range (~1.5-2.5 eV) can produce OPV materials with high PCE. Figure 4e, most \({E}_{g}\) values in our database fall within this range. Additionally, molecular weight (Mw) plays a critical role in enhancing PCE. Figure 4d indicates a non-uniform positive correlation between \({J}_{{sc}}\) and the logarithm of Mw. Increasing the Mw of an identical polymer backbone is a straightforward approach to improve the PCE. In fact, a high Mw is believed to enhance the PCE of polymer OPVs due to increased crystallinity and intercrystallite connectivity. Consequently, optimizing \({E}_{g}\) and HOMO-LUMO migration levels is a reasonable strategy for designing polymer molecules, benefiting the synthesis and application of OPV materials. The LUMO energy difference between the donor and acceptor (\({E}_{{\rm{LL}}}^{{\rm{DA}}}\)) is also crucial for PCE, a large difference can lead to significant energy loss (\({V}_{{oc}}\,\)Loss) at the D/A interface (Fig. 4f), whereas \({E}_{{\rm{HL}}}^{{\rm{DA}}}\) can roughly estimate the driving force to dissociate excitons in the D/A interface. An illustration of this is shown in Fig. 3b. Similarly, the polymer dispersion index (PDI) describes the molecular weight distribution of the polymer. Defined as \({PDI}=\frac{{M}_{w}}{{M}_{n}}\), where \({M}_{w}\) and \({M}_{n}\) represent the weight average and number average molecular weights respectively, it underscores the correlation between molecular weight and PCE. We also find that atomic contributions to the monomer surface area (VSA) or polarizability (or molecular refractivity) are crucial factors influencing a polymer’s PCE. PEOE_VSA9 descriptors, which combine partial charges and surface area, are significant for OPV’s PCE. A higher PEOE_VSA value indicates a greater positive impact on the predicted PCE value. Similar patterns are observed with other combined descriptors. For instance, SMR_VSA10, the total VSA of atoms within a specific range of molecular refractivity (MR), positively affects PCE. MR values are calculated for each atom type using Wildman and Crippen’s method. If the total VSA of atoms has an MR between 4 and ∞ (SMR_VSA10), a higher PCE is likely. Key contributing atom types include C doubly bonded to a heteroatom, aromatic C with a heteroatom neighbor, aromatic bridgehead C, and aromatic C = C. Additionally, NumHeteroatoms positively impacts PCE. These 2-D topological/topochemical properties provide insights into molecular surface interactions, while FpDensityMorgan1 generates similarity fingerprints based on atomic chemical and connectivity attributes, also positively affecting PCE. Overall, this work reveals the correlation between polymer PCE and its physicochemical descriptors, such as HOMO, LUMO, molecular weight, and molecular refractivity.

OPVs are composed of donor (electron-donating) and acceptor (electron-accepting) material, both of which are organic in nature. The performance of OPVs is largely determined by the properties of the donor/acceptor (D/A) materials. In other words, the design strategy of the D/A materials and the synergistic effect of their combination are crucial to the performance of OPV devices. Therefore, to reduce the need for trial and error experiments and achieve efficient device performance (including complementary absorption and highly balanced charge transport characteristics, among others), the search for and discovery of synergistic donor/acceptor (D/A) combinations is indispensable. Herein, when we compute SHAP interaction values for all features, the dimension of SHAP is 1343*5*5 (where 1343 is the sample size and 5 is the number of features), which is used to capture the interaction effect of pairs. Additionally, the selection of the five features of the interaction is based on the ranking of the marginal contribution rate in SHAP values. The color represents the characteristic value along the vertical axis (red for high values, blue for low values). The complete interaction evaluation diagram is presented in Figure S6. The feature selection criterion of the interaction is based on the SHAP interaction value distribution (intuitively, the prominence of the red and blue regions). Specifically, the SHAP interaction value is used to represent the influence of the interaction of the two features on the model prediction. In other words, the standout red and blue regions in the interaction diagram have large interaction values and are more suitable for feature combination, while those overlapping together have no obvious interaction effect. From Fig. 3c, the variable in the green rectangle is suitable for feature combination, as indicated by the standout red and blue regions, whereas the variable in the yellow rectangle is not suitable. It is evident that the interaction between MaxPartialCharge and \({{E}}_{g}\) is relatively obvious. The narrower \({{E}}_{g}\) is, the easier it is for electrons or holes to jump from the valence band to the conduction band, and the higher the intrinsic carrier concentration, which has a positive contribution to the current. Whereas the current actually characterizes the speed of charge flow, making charge transport more efficient, which in turn increases PCE. It follows that MaxPartialCharge and \({E}_{g}\) play a synergistic role, whether it is reducing MaxPartialCharge or reducing \({E}_{g}\), it can promote charge separation and transmission, and improve PCE. Additionally, it is more interesting that the HOMO of the P-type material and the LUMO of the N-type material interact, which is consistent with the OPV mechanism. Organic materials absorb light energy to generates tightly bounded electron-hole pairs, namely, excitons. Owing to the large binding energy of exciton, thermal separation of electron and hole is hardly possible at room temperature (around 20 °C). To separate the electron and hole, OPV utilizes the D/A interface to surpass such binding energy. The energy difference between the \({{E}}_{g}\) and the charge-transfer state energy (E_CT) provides the ΔE₁ for exciton dissociation, which is equal to the lowest unoccupied molecular orbital (LUMO) energy level difference between the donor and the acceptor (as shown in Fig. 4f). This optimal energy level difference helps to efficiently transfer excited electrons to the N-type material, minimizing charge recombination and boosting photovoltaic conversion efficiency. Additionally, aligning the energy levels of P-type and N-type materials enhances interface stability, facilitates efficient charge separation and transport, and minimizes energy loss, thereby improving device performance. The energy level difference between the HOMO and LUMO determines the generation of photocurrent-optimal discrepancies allow for maximum photon absorption and charge carrier production, thus enhancing current output and overall device efficiency.

In summary, effective interaction between P-type’s HOMO and N-type’s LUMO is a crucial factor in achieving high-performance OPV devices, and researching and optimizing this interaction can significantly enhance the performance and photovoltaic conversion efficiency of organic photovoltaic devices.

Explicable structure-property relationship analysis in OPVs

As outlined in METHODS section, the polymer units (PU) identified by PURS are collected into the polymer-unit library (as shown in Fig. 5a), which is organized by the number of rings and element types using PURS (Fig. 5b, c). Among them, by the number of rings can be divided into branch chain, mono ring, fused ring; and are then sorted by their element composition. This classification is essential to facilitate subsequent combinations of different PUs to develop new materials. More detailed information regarding the polymer units is available in the Supporting Information.

a polymer unit (PU) identified from the OPV database by PURS. b Classify by ring number. c Classify by element type.

To evaluate the marginal contribution of each PU to the PCE, we employed SHAP analysis on 260 donor materials and 1343 non-fullerene acceptor materials based on RF model. SHAP decomposes the prediction into the sum of contributions from each input feature, enabling the interpretability of the importance of each PU. A higher importance value indicates a greater reliance of the machine learning algorithm on a specific PU for determining the performance of an acceptor material.

Using PUFp as input, we examined the characteristics of polymer units with substantial SHAP values across three RF models (P-type, N-type, and P/N interactions) (Fig. 6b–d) and labeled them as important PU. The chemical structures corresponding to the important PUs of P-type OPV materials are depicted in Fig. 6f, and serial No. refers to its index number in the PU library (Fig. 5). The No. 200 PU is benzo[1,2-b:4,5-b’] dithiophene-4,8-dione, which has a quinone resonance structure, giving the polymer a good plane, while further improving the electron absorption capacity. More importantly, the quinone resonance structure is beneficial to enhance the charge transfer within the D-A polymer molecule. PU No. 175 is a heterocyclic structure containing S atoms, which will increase the rigidity after entering the main chain of polymerization, so that the free spin of the molecular chain segment is limited, so that the polymer has excellent photoelectric properties. The No. 24 PU contain imide groups, the electron-withdrawing groups, which contribute to lowering the LUMO level of the polymer and facilitating electron injection into the conduction band.

a The generation strategy of PUFp. b–e The interpretations of the ML models for P-type, N-type, P-N classification and D-A polymer-unit interaction of N-type by the SHAP evaluation. The blue and red bars on the right denote the proportional relation between the units and the prediction values. f–i Chemical structures of the PUs and their roles are identified through the importance analysis.

The significant PUs for N-type OPVs is shown as insets in Fig. 6g. Key PUs include: PU No. 283 contains a thiazole structure, and the electrostatic attraction between the sulfur and nitrogen atoms in thiazoles is beneficial to forming a closer π-π packing structure, which is a common strategy in the D-A OPV design. The No. 305 PU is quinoxaline, as a well-known electron-deficient system, which can not only improve the coplanarity of the polymeric main chain, but also extend the length of the π-π conjugated system to a large extent and increase the intensity of π-π close packing, which is a promising acceptor unit at present. Halogenation of electron acceptor units, such as Nos.304 and 97 PUs, can enhance the intramolecular charge transfer (ICT) effect and reduce the band gap of small non-fullerene receptors, which is one of the more effective molecular design strategies.

In Fig. 6d, the characteristic interaction evaluation of polymer units of P-type and N-type OPV materials was carried out to construct characteristic combinations, and the important PU obtained was shown in Fig. 6h. The complete interactive evaluation diagram is shown in Figure S7. More information about the characteristic interaction evaluation of polymer units of P-type and N-type OPV materials can be found in the Supporting Information. Herein, when we compute SHAP interaction values for all features, the dimension of SHAP is 1343*7*7 (where 1343 is the sample size and 7 is the number of features), which is used to capture the interaction effect of pairs. From Fig. 6d, the variable in the green rectangle is suitable for the feature combination because of the red and blue parts that stand out, whereas the variable in the yellow rectangle is not suitable. As a result, we screened out five variables suitable for feature combination, whose sequence number combinations are Nos. 175 and 382, Nos. 175 and 304, Nos. 200 and 382, Nos. 24 and 382, and Nos. 135 and 304, respectively. The corresponding PU of each sequence number is given in Fig. 6h. For the No. 175, as a donor unit, 2-methylthiophene has a strong electron transfer effect, which enhances the conjugated plane gravity and reduces the π-π packing distance. Perylene diimide (PDI) plays a catalytic role. If the strong coplanar PDI unit is introduced into the main chain, the charge delocalization ability inside the molecule can be increased, the π-π packing distance can be reduced, and the PCE can be increased. For the Nos. 175 and 304 combinations, the large atomic radius and special atomic orbital arrangement of halogen atoms can disperse the electron cloud density. Conjugated polymers based on fluorine or chlorine substitution usually exhibit better FF and Voc. Introducing halogen atoms into the sealing of non-fullerene accepter materials can reduce the molecular energy level, enhance the intramolecular charge transfer, and enhance the molecular crystallization. Additionally, the introduction of two-dimensional conjugated side chains in PU No.304 can increase the molecular conjugated area, broaden the spectral absorption, promote the interaction between molecules, and facilitate the formation of nanoscale bicontinuous phase separation during the preparation of thin films to the donor-acceptor blend, thus showing good photovoltaic performance. Using PUFp as input, we analyze the Pearson coefficients for D/A materials, and the detailed thermal map is shown in Figures S8.

In Fig. 6e, the characteristic interaction evaluation of D-A polymer units of N-type OPV materials was carried out to construct characteristic combinations, and the important D-A PU obtained was shown in Fig. 6i. The complete interactive evaluation diagram is shown in Figure S9. From the interaction diagram, we can find that the units suitable for feature combination are: No.304 (A) and No.105 (D), No.151 (D) and No.382 (A), No.283 (A) and No.151 (D), No.304 (A) and No.197 (D), No.149 (A) and No.151 (D), No. 355(D) and No.305 (A), etc. It provides ideas for the next important PU combination and structure design. Using PUFp as input, we analyze the Pearson coefficients for D/A units (in type acceptor materials), and the detailed thermal map is shown in Figures S10. When there are multiple thieno[3,4-b] thiophene electron-absorbing units in the structure, the intramolecular charge transfer (ICT) can make the material better absorb sunlight and improve the photoelectric conversion efficiency. Then, the introduction of halogen atoms can enhance the ICT effect and reduce the band gap of non-fullerene acceptors, which is one of the most effective molecular design strategies. The combination of thieno[3,4-b] thiophene and thiazole forms a rigid conjugated plane with rich heteroatoms, which is conducive to electron delocalization, and is a promising PU. For D-A polymer unit, ICT is generated due to the push-pull electron interaction between D and A, which reduces the band gap and causes the absorption redshift. Meanwhile, π-bridge is often used between D and A to reduce steric hindrance and improve the molecular planarity of the polymer. More importantly, it can be found that these different types of polymer units are common building blocks in D/A polymer molecules for the synthesis of OPV materials. In brief, the optimization objectives are as follows:

• The copolymerization of donor unit and acceptor unit was used to reduce the energy level band gap and broaden the spectral absorption.

• The HOMO energy level is reduced by introducing electron pushing groups.

• Through the precise introduction of fluorine/chlorine atom substitution on the polymer skeleton, the regulation of molecular energy levels, absorption, film morphology and charge dynamics can be achieved, while improving the \({J}_{{sc}}\) and FF, thereby improving the PCE and reducing energy loss.

• By introducing conjugated side chains to construct two-dimensional molecules, the coplanarity of molecules are increased and the PCE is improved.

Design and Screening of High PCE OPV Acceptor Materials Based on Important Polymer Units

By combining important PU in N-type OPV materials, we designed new polymer molecules to test the accuracy and rapid screening capabilities of our framework. The top 20 important PUs in N-type materials were categorized into three groups: donor polymer units (D), acceptor polymer units (A), and branched chain (C), as shown in Fig. 7a. Among them, there are five types of donor polymer units, six types of acceptor polymer units, and nine types of branched chains. Without specific constraints, a vast space composed of numerous structures (~1,048,576) is generated. In Figure S12, distributions of the polymer-unit type are shown and the donor polymer units, acceptor polymer units, and branched chain categories were used as the axes for all OPVs in the studied database. Figure S12 shows many empty areas in both the N-type OPVs, and obviously, these unreported combinations generate a huge space composed of many structures. In other words, there are many new materials that have not been explored based on the combination of existing PU, leaving a lot of space to be explored. To reduce the number of unreported candidates OPV materials within this categorization, the range of D, A, and C is limited to macromolecule with a high PCE ( > 12) and the macromolecule composition of at least one macromolecule composition of type D-A or A-D. The machine learning-based scheme for high PCE prediction is shown in Fig. 7b. Using these qualifications, we generated 3336 acceptor material combinations that matched 260 donor materials and employed the trained RF model (shown in Figure S3) to predict their PCEs and identify the combination with the highest PCE. The example of screened high PCE OPV acceptor materials is shown in Fig. 7c, and PCE > 14 value about 2678 combinations are provided in Support information.

a The top 20 important PUs in N-type and classification, abbreviation definitions for the PUs: donor polymer unit (D), acceptor polymer unit (A) and branched chain group (C). b Scheme of high PCE combination prediction by ML. c New materials generated by the PU combination process.

Additionally, we mapped the key building blocks from these highlighted polymer units and compared them to the structures of high-performance OPV acceptor materials. Our chemical structure analysis revealed that PUs like No.283 and No.149 were prevalent in over 14% of high PCE polymer acceptor materials. Firstly, the electrostatic attraction between sulfur and nitrogen atoms in the thiazoles (Fig. 7a) promotes tighter π-π stacking, which is a common strategy for designing D-A type acceptor materials. Meanwhile, quinoline enhances the coplanarity of the main chain, and its nitrogen atoms usually share electrons through covalent bonds with empty orbital electrons in other elements. Furthermore, halogenation of electron-accepting units can enhance intramolecular charge transfer (ICT) effects and reduce the bandgaps of non-fullerene acceptors.

As shown in Fig. 7c, the structure contains chlorine-containing fused rings, which broaden light absorption and contribute to higher short-circuit currents (\({J}_{{sc}}\)). The inclusion of strong electron donor groups like alkoxy chloride in the polymer backbone improves both the processability and photoelectric properties of the conjugated polymer. Additionally, chlorination is easier to synthesize compared to fluorination. Studies have shown that molecular design involving chlorination can expand light absorption and improve output voltage. This enables modification of non-radiative energy losses in OPV cells through chemical modification of the photoactive material, providing an opportunity to design efficient OPV materials with low bandgap-voltage offsets.

More importantly, we visualized the top 1000 of the 3336 combinations (targeting the acceptor material) for which PCE > 12 had been predicted. The violin plot, a data visualization that combines features of a boxplot and a kernel density map, shows how the data is distributed. Here, red represents all OPV acceptor materials, green represents A-D-A type OPV materials, and blue represents A-DA’-D-A type OPV materials. Figure 8b shows the distribution of predicted PCE values for all top 1000 designed and screened OPV acceptor materials, A-D-A type OPV materials, and A-DA’-D-A type OPV materials, respectively. The density curve illustrates the distribution of PCE under three categories of classification. The wider parts indicate more concentrated data, while the narrower parts indicate relatively fewer data points. Notably, the overall distribution of PCE values for A-DA’-D-A type OPV materials is higher than for A-D-A type OPV materials, indicating that A-DA’-D-A type OPV materials exhibit better structural properties and are more suitable as candidate structures for OPV materials. This finding provides a reliable strategy and guideline for the design of OPV materials.

a The representative OPV material acceptor structure. b Predicted PCE value distribution for the top 1000 combinations.

In summary, by leveraging advanced machine learning (ML) technology, we studied polymers to model highly optimized, efficient, and stable polymer structures for organic photovoltaic (OPV) cells. A significant amount of photovoltaic property data was collected from reported experimental studies and used to train ML models. We developed five models using RF, MLP, KNN, KRR, and SVM algorithms, with the RF regression model demonstrating the best predictive ability. Various representations of acceptor molecules, including descriptors, MACCS, and polymer unit fingerprint (PUFp), were employed to build ML models for predicting the corresponding OPV PCE class. The results indicate that PUFp with a length greater than 600 bits provides the best representation of acceptor molecules. In feature-property analysis, the polymers’ highest occupied molecular orbital (HOMO), lowest unoccupied molecular orbital (LUMO), molecular weight (\({M}_{w}\)), and band gap (\({E}_{g}\)) emerged as the most decisive descriptors. A library of 413 polymer units was constructed, and key polymer units affecting NFA (non-fullerene acceptor) materials were identified. More importantly, by combining these key polymer units in N-type OPV materials, new polymer molecules were designed to test the accuracy and rapid screening capabilities of our framework. Our research for the relationship between feature/structure and PCE can accelerate the design of new acceptor materials, thus advancing the development of high-PCE OPVs. Our methodology offers a promising approach for screening and designing new polymer acceptors for OPVs and can be applied to a wide range of donor materials, thereby accelerating the development of high-performance OPVs.