@article{HeidenUsvyatSaalfrank2019,
  author    = {Heiden, Sophia and Usvyat, Denis and Saalfrank, Peter},
  title     = {Theoretical Surface Science Beyond Gradient-Corrected Density Functional Theory},
  series = {The journal of physical chemistry : C, Nanomaterials and interfaces},
  volume    = {123},
  journal   = {The journal of physical chemistry : C, Nanomaterials and interfaces},
  number    = {11},
  publisher = {American Chemical Society},
  address   = {Washington},
  issn      = {1932-7447},
  doi       = {10.1021/acs.jpcc.9b00407},
  pages     = {6675 -- 6684},
  year      = {2019},
  abstract  = {The quantum chemical description of the adsorption, vibrations, and reactions of molecules at periodic solid surfaces is frequently based on a methodological "standard model": density functional theory (DFT) in the generalized gradient approximation (GGA), using plane wave bases and three-dimensional supercells. Although the computationally efficient GGA functionals can be very successful, cases are known where they do not perform so well. Most importantly, activation energies for chemical reactions are typically underestimated, with the consequence of computed reaction rates being too large. In this work, we consider a well-studied model system: water or water fragments adsorbed on an Al-terminated alpha-Al2O3(0001) surface as a test bed for studying the performance of electronic structure methods, both from DFT and wave function theory. On the DFT side, we employ two GGA exchange correlation functionals: PW91 and PBE with and without dispersion corrections, whose results are then compared to those of hybrid functionals B3LYP and HSE06. Further, we follow a periodic wave function approach in the form of local second-order Moller-Plesset perturbation theory, LMP2, on a Hartree-Fock reference. En route, we address issues arising from the choice of the basis set. The key findings of our study are as follows: (i) DFT-GGA adsorption energies are in reasonable agreement with both hybrid-DFT and LMP2 values. In particular, the deviations between the relative energies, corresponding to different adsorption structures, are in the range of the error due to missing dispersion corrections or the basis set error. (ii) Harmonic DFT-GGA vibrational frequencies for oxygen hydrogen stretch modes are by several tens of wavenumbers red-shifted compared to corresponding hybrid-DFT values. The latter are in much better agreement with recent experimental data. (iii) The activation energy for a hydrogen diffusion reaction is grossly underestimated by GGA compared to hybrid-DFT or LMP2, which in turn are quite comparable.},
  language  = {en}
}
@phdthesis{Heiden2018,
  author    = {Heiden, Sophia L.},
  title     = {Water at α-alumina surfaces},
  doi       = {10.25932/publishup-42636},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-426366},
  school      = {Universit{\"a}t Potsdam},
  pages     = {124},
  year      = {2018},
  abstract  = {The (0001) surface of α-Al₂O₃ is the most stable surface cut under UHV conditions and was studied by many groups both theoretically and experimentally. Reaction barriers computed with GGA functionals are known to be underestimated. Based on an example reaction at the (0001) surface, this work seeks to improve this rate by applying a hybrid functional method and perturbation theory (LMP2) with an atomic orbital basis, rather than a plane wave basis. In addition to activation barriers, we calculate the stability and vibrational frequencies of water on the surface. Adsorption energies were compared to PW calculations and confirmed PBE+D2/PW stability results. Especially the vibrational frequencies with the B3LYP hybrid functional that have been calculated for the (0001) surface are in good agreement with experimental findings. Concerning the barriers and the reaction rate constant, the expectations are fully met. It could be shown that recalculation of the transition state leads to an increased barrier, and a decreased rate constant when hybrid functionals or LMP2 are applied. Furthermore, the molecular beam scattering of water on (0001) surface was studied. In a previous work by Hass the dissociation was studied by AIMD of molecularly adsorbed water, referring to an equilibrium situation. The experimental method to obtaining this is pinhole dosing. In contrast to this earlier work, the dissociation process of heavy water that is brought onto the surface from a molecular beam source was modeled in this work by periodic ab initio molecular dynamics simulations. This experimental method results in a non-equilibrium situation. The calculations with different surface and beam models allow us to understand the results of the non-equilibrium situation better. In contrast to a more equilibrium situation with pinhole dosing, this gives an increase in the dissociation probability, which could be explained and also understood mechanistically by those calculations. In this work good progress was made in understanding the (1120) surface of α-Al₂O₃ in contact with water in the low-coverage regime. This surface cut is the third most stable one under UHV conditions and has not been studied to a great extent yet. After optimization of the clean, defect free surface, the stability of different adsorbed species could be classified. One molecular minimum and several dissociated species could be detected. Starting from these, reaction rates for various surface reactions were evaluated. A dissociation reaction was shown to be very fast because the molecular minimum is relatively unstable, whereas diffusion reactions cover a wider range from fast to slow. In general, the (112‾0) surface appears to be much more reactive against water than the (0001) surface. In addition to reactivity, harmonic vibrational frequencies were determined for comparison with the findings of the experimental "Interfacial Molecular Spectroscopy" group from Fritz-Haber institute in Berlin. Especially the vibrational frequencies of OD species could be assigned to vibrations from experimental SFG spectra with very good agreement. Also, lattice vibrations were studied in close collaboration with the experimental partners. They perform SFG spectra at very low frequencies to get deep into the lattice vibration region. Correspondingly, a bigger slab model with greater expansion perpendicular to the surface was applied, considering more layers in the bulk. Also with the lattice vibrations we could obtain reasonably good agreement in terms of energy differences between the peaks.},
  language  = {en}
}
@article{HeidenWirthCampenetal.2018,
  author    = {Heiden, Sophia and Wirth, Jonas and Campen, Richard Kramer and Saalfrank, Peter},
  title     = {Water molecular beam scattering at alpha-Al2O3(0001)},
  series = {The journal of physical chemistry : C, Nanomaterials and interfaces},
  volume    = {122},
  journal   = {The journal of physical chemistry : C, Nanomaterials and interfaces},
  number    = {27},
  publisher = {American Chemical Society},
  address   = {Washington},
  issn      = {1932-7447},
  doi       = {10.1021/acs.jpcc.8b04179},
  pages     = {15494 -- 15504},
  year      = {2018},
  abstract  = {Recent molecular beam experiments have shown that water may adsorb molecularly or dissociatively on an α-Al2O3(0001) surface, with enhanced dissociation probability compared to "pinhole dosing", i.e., adsorption under thermal equilibrium conditions. However, precise information on the ongoing reactions and their relative probabilities is missing. In order to shed light on molecular beam scattering for this system, we perform ab initio molecular dynamics calculations to simulate water colliding with α-Al2O3(0001). We find that single water molecules hitting a cold, clean surface from the gas phase are either reflected, molecularly adsorbed, or dissociated (so-called 1-2 dissociation only). A certain minimum translational energy (above 0.1 eV) seems to be required to enforce dissociation, which may explain the higher dissociation probability in molecular beam experiments. When the surface is heated and/or when refined surface and beam models are applied (preadsorption with water or water fragments, clustering and internal preexcitation in the beam), additional channels open, among them physisorption, water clustering on the surface, and so-called 1-4 and 1-4′ dissociation.},
  language  = {en}
}
@article{HeidenYueKirschetal.2018,
  author    = {Heiden, Sophia and Yue, Yanhua and Kirsch, Harald and Wirth, Jonas A. and Saalfrank, Peter and Campen, Richard Kramer},
  title     = {Water dissociative adsorption on α-Al2O3(112̅0) is controlled by surface site undercoordination, density, and topology},
  series = {The journal of physical chemistry / publ. weekly by the American Chemical Society : C, Nanomaterials and interfaces},
  volume    = {122},
  journal   = {The journal of physical chemistry / publ. weekly by the American Chemical Society : C, Nanomaterials and interfaces},
  number    = {12},
  publisher = {American Chemical Society},
  address   = {Washington},
  issn      = {1932-7447},
  doi       = {10.1021/acs.jpcc.7b10410},
  pages     = {6573 -- 6584},
  year      = {2018},
  abstract  = {α-Al2O3 surfaces are common in a wide variety of applications and useful models of more complicated, environmentally abundant, alumino-silicate surfaces. While decades of work have clarified that all properties of these surfaces depend sensitively on the crystal face and the presence of even small amounts of water, quantitative insight into this dependence has proven challenging. Overcoming this challenge requires systematic study of the mechanism by which water interacts with various α-Al2O3 surfaces. Such insight is most easily gained for the interaction of small amounts of water with surfaces in ultra high vacuum. In this study, we continue our combined theoretical and experimental approach to this problem, previously applied to water interaction with the α-Al2O3 (0001) and (11̅02) surfaces, now to water interaction with the third most stable surface, that is, the (112̅0). Because we characterize all three surfaces using similar tools, it is straightforward to conclude that the (112̅0) is most reactive with water. The most important factor explaining its increased reactivity is that the high density of undercoordinated surface Al atoms on the (112̅0) surface allows the bidentate adsorption of OH fragments originating from dissociatively adsorbed water, while only monodentate adsorption is possible on the (0001) and (11̅02) surfaces: the reactivity of α-Al2O3 surfaces with water depends strongly, and nonlinearly, on the density of undercoordinated surface Al atoms.},
  language  = {en}
}