Computational investigation of α-SiO₂ surfaces as a support for Pd

The properties of a supported metal catalyst depend crucially on the interaction between the active metal and the support. A case in point is Pd supported on silica, Pd/SiO₂, which is widely used in oxidation catalysis. There is a need for a broad range of computational models that describe the interaction of Pd with silica surfaces so that active site models can be proposed and tested. In this work, we create well-defined, reproducible, periodic models of SiO₂ surfaces and investigate their interaction with Pd using dispersion-corrected DFT. We use crystalline α-SiO₂ as a useful starting point for creating and estimating the adsorption properties of metals on SiO₂ surfaces, which can represent the specific isolated functional groups present on more complex amorphous silica surfaces. We have modelled α-SiO₂ (001), (100) and (101) surfaces containing isolated siloxane and silanol functional groups and estimated their affinity towards the adsorption of Pd atoms regarding an isolated gaseous Pd atom and the fcc Pd solid. This provides additional information on the ease with which Pd can be dispersed on the surfaces in question. From our model, we characterise the surface energies of the α-SiO₂ (hkl) surfaces and calculate the geometries of the Pd₁/α-SiO₂ (hkl) adsorption site on each surface. We estimate that Pd₁(g) will prefer to adsorb close to strained four-membered siloxane rings or on a vicinal silanol group of α-SiO₂ (101).