Boosting renewable energies is the 7th goal in the 2030 Agenda for Sustainable development whose premise is to accelerate the transition towards a low-carbon competitive economy. Sunlight is the largest exploitable resource of energy and it has promoted an extraordinary growth of photovoltaic (PV) cell technology to supply the demand. Several PV manufacturers are now working on bringing solar cell devices based on organic-inorganic hybrid perovskite (PSCs) to market due to their easy fabrication procedures, low cost, and the fast improvement of their efficiency (PCE). In order to enhance the possibility for PSCs to be competitive throughout the market, a very recent technology introduces Fullerene (C60) and its derivatives -specially [6,6]-phenyl-C61-butyric acid methyl ester (PCBM)- in PSC devices to work as electron transport layer (ETL), additives, or to control the crystal growth, either in conventional or inverted devices. The use of fullerenes have eliminated/reduced the hysteresis,[1] improved the device stability in comparison to fullerene-free standard PSCs,[2] and favors the high-scale preparation because the high-temperature step.Nevertheless, important features in the state-of-the-art fullerene-perovskite devices remain unexplored. It is known that fullerenes passivate the perovskite defects, improving the optoelectronic properties, but it has not been systematically explored to date. It is not extensively investigated the effect of binding groups of fullerenes, how fullerenes passivate the perovskite surfaces and how passivation influence surface energetics and the corresponding carrier dynamics. In this work we use periodic density functional simulations to explore the interaction of fullerenes with different perovskite surfaces, considering as well possible defects (Iodide anti-sites) and vacancies (cation vacancies). Computing the binding energy, the adsorption energy site, bandgaps, and density of states we are able to compare these relevant descriptors with the experimental efficiency, in order to find the key descriptors that govern the fullerene-perovskite interactions. In this presentation, we will show our results for APbI3 (A= methylammonium, Cs) perovskite surfaces with the most common fullerenes like C60 and PC61BM. Our simulations show different binding energies depending on the perovskite and surface vacancy and the effect on the electronic structure will be discussed.1. J. Xu, et. al., Nat. Commun. 2015,6, 7081; Y. Zhong, et. al., Adv. Funct. Mater. 2020, 1908920.2. V. Ferguson, et. al., Energy Environ. Mater. 2019, 2, 107–118; T. Gatti et. al., Nano Energy, 2017, 41, 84−100. Figure 1