Organic photovoltaic (OPV) based on πconjugated molecules or polymers converts the solar energy to electricity utilizing the electron donor (D)/acceptor (A) bulk heterojunction structure to efficiently separate the bound excitons generated by light absorption into free charge carriers. [1] It has attracted immense research interest due to their numerous technology advantages such as inherent flexibility, portability, lightweight, low cost, and low toxicity, [2] as well as recently reported a high power conversion efficiency (PCE) of 19%. [3] However, to commercialize these technologies, further efficiency improvements are desired. This can be hardly accomplished by tedious and expensive trial-and-error methods based on chemical intuitions, labor-intensive synthesis, and characterization, especially when considering the large chemical space available for designing new organic donor or acceptor molecules and the numerous experimental parameters that can be tuned in the fabrication of the devices. [4] It is important to realize that the construction of a quantitative structure-property relationship (QSPR) model [5] would greatly benefit the search of new functional molecules and the optimization of experimental conditions, and, accordingly, it can become a key factor to improve the photovoltaic performance of OPVs. [5b] Furthermore, an interpretable QSPR model can be also used to unravel the fundamental working mechanism for a complex experimental behavior through its statistical correlation analysis with various microscopic physical parameters, which can hardly be achieved by blackbox predictions. Unfortunately, deriving such quantitative and interpretable QSPR models is highly nontrivial due to our current poor understanding of the complex photophysical phenomena in OPVs such as photoabsorption, exciton diffusion, exciton dissociation, charge generation and charge transport, as well as charge collection at electrodes. [1c,6] Earlier empirical QSPR models such as Scharber's model [7] for fullerene-based OPVs and the modified Scharber's model [8] for nonfullerene-based OPVs have been proposed for predicting OPV performance, yet, their accuracy for predicting PCE is unsatisfactory (with a linear Pearson correlation coefficient (r) between Scharber's PCE and experimental PCE at around only 0.1-0.2). [9] Driven by the development of artificial intelligence (AI) [10] algorithms, as well as the increased availability of a high-quality large dataset along with efficient data-mining techniques, [11]