@article{MelaniNagataWirthetal.2018, author = {Melani, Giacomo and Nagata, Yuki and Wirth, Jonas and Saalfrank, Peter}, title = {Vibrational spectroscopy of hydroxylated alpha-Al2O3(0001) surfaces with and without water}, series = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr}, volume = {149}, journal = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr}, number = {1}, publisher = {American Institute of Physics}, address = {Melville}, issn = {0021-9606}, doi = {10.1063/1.5023347}, pages = {10}, year = {2018}, abstract = {Using gradient- and dispersion-corrected density functional theory in connection with ab initio molecular dynamics and efficient, parametrized Velocity-Velocity Autocorrelation Function (VVAF) methodology, we study the vibrational spectra (Vibrational Sum Frequency, VSF, and infrared, IR) of hydroxylated alpha-Al2O3(0001) surfaces with and without additional water. Specifically, by considering a naked hydroxylated surface and the same surface with a particularly stable, "ice-like" hexagonal water later allows us to identify and disentangle main spectroscopic bands of OH bonds, their orientation and dynamics, and the role of water adsorption. In particular, we assign spectroscopic signals around 3700 cm(-1) as being dominated by perpendicularly oriented non-hydrogen bonded aluminol groups, with and without additional water. Furthermore, the thin water layer gives spectroscopic signals which are already comparable to previous theoretical and experimental findings for the solid/(bulk) liquid interface, showing that water molecules closest to the surface play a decisive role in the vibrational response of these systems. From a methodological point of view, the effects of temperature, anharmonicity, hydrogen-bonding, and structural dynamics are taken into account and analyzed, allowing us to compare the calculated IR and VSF spectra with the ones based on normal mode analysis and vibrational density of states. The VVAF approach employed in this work appears to be a computationally accurate yet feasible method to address the vibrational fingerprints and dynamical properties of water/metal oxide interfaces. Published by AIP Publishing.}, language = {en} } @article{BouaklineLorenzMelanietal.2017, author = {Bouakline, Foudhil and Lorenz, Ulrich J. and Melani, Giacomo and Paramonov, Guennaddi K. and Saalfrank, Peter}, title = {Isotopic effects in vibrational relaxation dynamics of H on a Si(100) surface}, series = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr}, volume = {147}, journal = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr}, number = {14}, publisher = {American Institute of Physics}, address = {Melville}, issn = {0021-9606}, doi = {10.1063/1.4994635}, pages = {11}, year = {2017}, abstract = {In a recent paper [U. Lorenz and P. Saalfrank, Chem. Phys. 482, 69 (2017)], we proposed a robust scheme to set up a system-bath model Hamiltonian, describing the coupling of adsorbate vibrations (system) to surface phonons (bath), from first principles. The method is based on an embedded cluster approach, using orthogonal coordinates for system and bath modes, and an anharmonic phononic expansion of the system-bath interaction up to second order. In this contribution, we use this model Hamiltonian to calculate vibrational relaxation rates of H-Si and D-Si bending modes, coupled to a fully H(D)-covered Si(100)-(2×1) surface, at zero temperature. The D-Si bending mode has an anharmonic frequency lying inside the bath frequency spectrum, whereas the H-Si bending mode frequency is outside the bath Debye band. Therefore, in the present calculations, we only take into account one-phonon system-bath couplings for the D-Si system and both one- and two-phonon interaction terms in the case of H-Si. The computation of vibrational lifetimes is performed with two different approaches, namely, Fermi's golden rule, and a generalized Bixon-Jortner model built in a restricted vibrational space of the adsorbate-surface zeroth-order Hamiltonian. For D-Si, the Bixon-Jortner Hamiltonian can be solved by exact diagonalization, serving as a benchmark, whereas for H-Si, an iterative scheme based on the recursive residue generation method is applied, with excellent convergence properties. We found that the lifetimes obtained with perturbation theory, albeit having almost the same order of magnitude—a few hundred fs for D-Si and a couple of ps for H-Si—, are strongly dependent on the discretized numerical representation of the bath spectral density. On the other hand, the Bixon-Jortner model is free of such numerical deficiencies, therefore providing better estimates of vibrational relaxation rates, at a very low computational cost. The results obtained with this model clearly show a net exponential decay of the time-dependent survival probability for the H-Si initial vibrational state, allowing an easy extraction of the bending mode "lifetime." This is in contrast with the D-Si system, whose survival probability exhibits a non-monotonic decay, making it difficult to define such a lifetime. This different behavior of the vibrational decay is rationalized in terms of the power spectrum of the adsorbate-surface system. In the case of D-Si, it consists of several, non-uniformly distributed peaks around the bending mode frequency, whereas the H-Si spectrum exhibits a single Lorentzian lineshape, whose width corresponds to the calculated lifetime. The present work gives some insight into mechanisms of vibration-phonon coupling at surfaces. It also serves as a benchmark for multidimensional system-bath quantum dynamics, for comparison with approximate schemes such as reduced, open-system density matrix theory (where the bath is traced out and a Liouville-von Neumann equation is solved) or approximate wavefunction methods to solve the combined system-bath Schr{\"o}dinger equation.}, language = {en} } @article{YueMelaniKirschetal.2022, author = {Yue, Yanhua and Melani, Giacomo and Kirsch, Harald and Paarmann, Alexander and Saalfrank, Peter and Campen, Richard Kramer and Tong, Yujin}, title = {Structure and Reactivity of a-Al2O3(0001) Surfaces: How Do Al-I and Gibbsite-like Terminations Interconvert?}, series = {The journal of physical chemistry / publ. weekly by the American Chemical Society. C, Energy, materials, and catalysis}, volume = {126}, journal = {The journal of physical chemistry / publ. weekly by the American Chemical Society. C, Energy, materials, and catalysis}, number = {31}, publisher = {American Chemical Society}, address = {Washington}, issn = {1932-7447}, doi = {10.1021/acs.jpcc.2c03743}, pages = {13467 -- 13476}, year = {2022}, abstract = {The alpha-Al2O3(0001) surface has been extensively studied because of its significance in both fundamental research and application. Prior work suggests that in ultra-high-vacuum (UHV), in the absence of water, the so-called Al-I termination is thermodynamically favored, while in ambient, in contact with liquid water, a Gibbsite-like layer is created. While the view of the alpha- Al2O3(0001)/H2O(l) interface appears relatively clear in theory, experimental characterization of this system has resulted in estimates of surface acidity, i.e., isoelectric points, that differ by 4 pH units and surface structure that in some reports has non-hydrogen-bonded surface aluminol (Al-OH) groups and in others does not. In this study, we employed vibrational sum frequency spectroscopy (VSFS) and density functional theory (DFT) simulation to study the surface phonon modes of the differently terminated alpha-Al2O3(0001) surfaces in both UHV and ambient. We find that, on either water dosing of the Al-I in UHV or heat-induced dehydroxylation of the Gibbsite-like in ambient, the surfaces do not interconvert. This observation offers a new explanation for disagreements in prior work on the alpha-Al2O3(0001)/liquid water interface -different preparation methods may create surfaces that do not interconvert-and shows that the surface phonon spectral response offers a novel probe of interfacial hydrogen bonding structure.}, language = {en} } @article{MelaniNagataCampenetal.2019, author = {Melani, Giacomo and Nagata, Yuki and Campen, Richard Kramer and Saalfrank, Peter}, title = {Vibrational spectra of dissociatively adsorbed D2O on Al-terminated alpha-Al2O3(0001) surfaces from ab initio molecular dynamics}, series = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr}, volume = {150}, journal = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr}, number = {24}, publisher = {American Institute of Physics}, address = {Melville}, issn = {0021-9606}, doi = {10.1063/1.5099895}, pages = {15}, year = {2019}, abstract = {Water can adsorb molecularly or dissociatively onto different sites of metal oxide surfaces. These adsorption sites can be disentangled using surface-sensitive vibrational spectroscopy. Here, we model Vibrational Sum Frequency (VSF) spectra for various forms of dissociated, deuterated water on a reconstructed, Al-terminated α-Al2O3(0001) surface at submonolayer coverages (the so-called 1-2, 1-4, and 1-4′ modes). Using an efficient scheme based on velocity-velocity autocorrelation functions, we go beyond previous normal mode analyses by including anharmonicity, mode coupling, and thermal surface motion in the framework of ab initio molecular dynamics. In this way, we calculate vibrational density of states curves, infrared, and VSF spectra. Comparing computed VSF spectra with measured ones, we find that relative frequencies of resonances are in quite good agreement and linewidths are reasonably well represented, while VSF intensities coincide not well. We argue that intensities are sensitively affected by local interactions and thermal fluctuations, even at such low coverage, while absolute peak positions strongly depend on the choice of the electronic structure method and on the appropriate inclusion of anharmonicity.}, language = {en} } @article{MelaniNagataSaalfrank2021, author = {Melani, Giacomo and Nagata, Yuki and Saalfrank, Peter}, title = {Vibrational energy relaxation of interfacial OH on a water-covered alpha-Al2O3(0001) surface}, series = {Physical chemistry, chemical physics : PCCP ; a journal of European chemical societies}, volume = {23}, journal = {Physical chemistry, chemical physics : PCCP ; a journal of European chemical societies}, number = {13}, publisher = {Royal Society of Chemistry}, address = {Cambridge}, issn = {1463-9076}, doi = {10.1039/d0cp03777j}, pages = {7714 -- 7723}, year = {2021}, abstract = {Vibrational relaxation of adsorbates is a sensitive tool to probe energy transfer at gas/solid and liquid/solid interfaces. The most direct way to study relaxation dynamics uses time-resolved spectroscopy. Here we report on a non-equilibrium ab initio molecular dynamics (NE-AIMD) methodology to model vibrational relaxation of OH vibrations on a hydroxylated, water-covered alpha-Al2O3(0001) surface. In our NE-AIMD approach, after exciting selected O-H bonds their coupling to surface phonons and to the water adlayer is analyzed in detail, by following both the energy flow in time, as well as the time-evolution of Vibrational Density of States (VDOS) curves. The latter are obtained from Time-dependent Correlation Functions (TCFs) and serve as prototypical, generic representatives of time-resolved vibrational spectra. As most important results, (i) we find a few-picosecond lifetime of the excited modes and (ii) identify both hydrogen-bonded aluminols and water molecules in the adsorbed water layer as main dissipative channels, while the direct coupling to Al2O3 surface phonons is of minor importance on the timescales of interest. Our NE-AIMD/TCF methodology is powerful for complex adsorbate systems, in principle even reacting ones, and opens a way towards time-resolved vibrational spectroscopy.}, language = {en} }