“…First, we develop a rigorous OET model based on the finiteelement method to address and quantify the multi-physics behaviors by coupling the optical, electrical, and thermal modules, [21,22] where the simulation details are provided in the Supporting Information and the related parameters used for this simulation are listed in Tables S1-S3 (Supporting Information). The microscopic energy conversion processes of photons, charge-carriers, and phonons of a PV cell are illustrated in Figure 1, which can be divided into six categories from the viewpoint of recombination sources: I) thermalization heat arising from the energy relaxation, i.e., photon-excited electrons (holes) with excess potential energy beyond bandgap (i.e., hν−E g , where ν is the frequency of incident light) return to the conduction (valence) band edge in picosecond timescales; [23] II) Joule heat caused by the motion of carriers under built-in electric-field within the depletion region; [24] III) bulk recombination heat contributing from Shockley-Read-Hall (SRH) and Auger recombinations; it is worth noting that although the radiation recombination is very important for V OC , it is not a source of heat generation, thus we do not discuss the radiation recombination when describing bulk recombination, and the Auger recombination is almost negligible compared with SRH recombination (Figure S1, Supporting Information). IV) surface recombination heat attributed to the carrier trapping effect by surface defects; V) Peltier heat at heterojunction interface due to the energy band offset; and VI) Peltier heat at the semiconductor/metal interface, where the transport carriers have to flow from conduction/valence band of semiconductor region to the quasi-Fermi level before being collected by the respective electrodes.…”