“…From the many possibilities available for the theoretical study of graphene, each of them with its own advantages, reliability, and range of application, we chose the procedure described in our previous work [ 70 ] and used density functional theory (DFT) [ 71 , 72 , 73 ] as the main tool for simulating the properties of graphene quantum dots through use of the the Gaussian 16 [ 74 ] package. For the exchange-correlation functional, local density approximation (LDA) [ 75 , 76 ] was used, in virtue of its better performance for graphitic systems and higher calculation speed compared with general gradient approximations (GGAs) or hybrid functionals such as B3LYP [ 77 , 78 , 79 , 80 , 81 ], and because it has been successfully used to analyze interactions in carbon nanostructures [ 82 , 83 , 84 , 85 , 86 , 87 ].…”