“…Regarding the computational predictions of GDDS, molecular dynamics (MD) simulation and quantum mechanics (QM) programs are powerful to accurately investigate the non-bonded interaction between the graphene/GO and ligand for a given GDDS (Gráfová et al, 2010;Ramraj and Hillier, 2010;Calero et al, 2013;Cho et al, 2013;Guo et al, 2013;Vovusha et al, 2013Vovusha et al, , 2018Mudedla et al, 2014;Vincent and Hillier, 2014;Wang et al, 2015;Mahdavi et al, 2016Mahdavi et al, , 2020Krepel and Hod, 2017;Safdari et al, 2017;Ajala et al, 2019;Alkathiri et al, 2019;Azhagiya Singam et al, 2019;Mason et al, 2019). However, both MD and QM methods are usually time-consuming, need to know a reasonable initial binding mode between the graphene/GO and ligand before calculations and tend to entrap the graphene/GO and ligand in a local minimum.…”