Graphene presents different types of defects as impurities, vacancies and grain boundaries which concentration may depend on growing methods. In addition, experiments report a strong increase of defect concentration on BLG surfaces when increasing the pressure after PTM solidification, 10 as well as the formation of strong strain fields leading to distribution of local folds before PTM solidification. Motivated by such findings and the results of the previous section, we have investigated the effect of singleand double vacancies on BLG functionalization at 6.8 GPa by SPRINT metadynamics.