“…For instance, classical Montecarlo (MC) and molecular dynamics (MD) are one of the best (non-atomistic and atomistic, respectively) approaches to model up to several nanometer large crystalline or amorphous organic and inorganic materials allowing obtaining rich information about molecular and nuclear position-related properties ( Zeng et al, 2003 ; Heinz and Ramezani-Dakhel, 2016 ; Ramakrishnan et al, 2017 ; Eckert et al, 2020 ). To compute electronic-related properties, we must perform first-principles calculations such as Density Functional Theory (DFT) and ab-initio MD (AIMD) to access static and dynamic information, respectively ( Zeng et al, 2003 ; Heinz and Ramezani-Dakhel, 2016 ; Ramakrishnan et al, 2017 ; Eckert et al, 2020 ). When performing first-principles calculations, due to the complexity of the h OI interface, one of the organic or inorganic phases is usually oversimplified to systems comprising organic monomer/inorganic surface ( Alexandre et al, 2010 ; Semoto et al, 2011 ; Hofmann et al, 2013 ; Motta et al, 2015 ; Wang et al, 2017a ; Pourrahimi et al, 2017 ; Liao et al, 2019 ) or organic oligomers/inorganic small clusters ( Mombrú et al, 2017a ; Ullah et al, 2017 ; Wang et al, 2020a ).…”