Schematic overview of the workflow. First the data are collected from the literature and divided into monolayer and multilayer coverage. Secondly, specific features are designated to the data, including atomic composition, functional groups, and valence electrons. Thirdly, the Gaussian process regression model is constructed and trained on the data to be able to predict BEs for new molecules.