Modelling of epitaxial graphene functionalization