include the strong electron-electron interactions and strong electron-phonon couplings as well as the complicated donor/ acceptor (D/A) interface morphology, which are fundamentally different from inorganic semiconductors. [23,24] Therefore, the accurate simulating of OPVs requires high-level theoretical methods in quantum chemistry, quantum dynamics and statistical mechanics, and in recent years there have been substantial progress for the theoretical understanding of many microscopic processes such as charge transport, [25] exciton dissociation, [26,27] singlet fission, [28][29][30] etc. These methods are useful for few benchmark systems but cannot be used to explore the chemical space, i.e., screen a large number of candidate materials.The predictive power of a theoretical model by high-throughput virtual screening has been explored in the recent years. [31][32][33][34][35][36][37][38][39][40] For example, in the Harvard Clean Energy Project (CEP), [41] Aspuru-Guzik and co-workers have screened ≈2.3 million compounds by employing the Scharber model [42,43] to find out efficient new donor molecules for OPVs. Here, the computed energies of the highest occupied molecular orbital (HOMO)/the lowest unoccupied molecular orbital (LUMO) of organic molecules are calibrated to experimentally determined values, and then employing an averaging scheme, energies of HOMO/LUMO are estimated to use them as input in the Scharber model. However, the prediction of PCE of an OPV is much more challenging compared to the energies of orbitals. Very recently, the same research team [44] noticed that the Scharber model, widely used in these virtual screening works, is not accurate enough for the prediction of PCE, and calibration of PCE by considering molecular similarity, molecular weight and band gap energy leads to improvement in correlation between experimental (49 OPVs) and calculated PCE (Pearson's coefficient (r) is increased from 0.30 to 0.43). Further, the Scharber model is only optimized for the PC 61 BM acceptor [42] and a more general model is desired to take account of nonfullerene molecules.In an OPV, energy conversion is accomplished by four consecutive steps: i) absorption of photons and exciton formation, ii) exciton diffusion to D/A interface, iii) exciton dissociation, and iv) transport of holes/electrons to the respective electrode. Taking account of all these processes, the efficiency of a device depends on many microscopic properties of the organic material such as optical gap, charge-carrier mobility, ionization To design efficient materials for organic photovoltaics (OPVs), it is essential to identify the largest number of parameters that control their properties and build a model using these parameters (known as descriptors) for the prediction of the power conversion efficiency (PCE). By constructing a dataset for 280 small molecule OPV systems, it is found that for all high-performing devices, frontier molecular orbitals of donor molecules are nearly degenerated and in such cases, orbitals other than just high...