@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{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} } @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{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} } @phdthesis{Melani2019, author = {Melani, Giacomo}, title = {From structural fluctuations to vibrational spectroscopy of adsorbates on surfaces}, doi = {10.25932/publishup-44182}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-441826}, school = {Universit{\"a}t Potsdam}, pages = {119}, year = {2019}, abstract = {Aluminum oxide is an Earth-abundant geological material, and its interaction with water is of crucial importance for geochemical and environmental processes. Some aluminum oxide surfaces are also known to be useful in heterogeneous catalysis, while the surface chemistry of aqueous oxide interfaces determines the corrosion, growth and dissolution of such materials. In this doctoral work, we looked mainly at the (0001) surface of α-Al 2 O 3 and its reactivity towards water. In particular, a great focus of this work is dedicated to simulate and address the vibrational spectra of water adsorbed on the α-alumina(0001) surface in various conditions and at different coverages. In fact, the main source of comparison and inspiration for this work comes from the collaboration with the "Interfacial Molecular Spectroscopy" group led by Dr. R. Kramer Campen at the Fritz-Haber Institute of the MPG in Berlin. The expertise of our project partners in surface-sensitive Vibrational Sum Frequency (VSF) generation spectroscopy was crucial to develop and adapt specific simulation schemes used in this work. Methodologically, the main approach employed in this thesis is Ab Initio Molecular Dynamics (AIMD) based on periodic Density Functional Theory (DFT) using the PBE functional with D2 dispersion correction. The analysis of vibrational frequencies from both a static and a dynamic, finite-temperature perspective offers the ability to investigate the water / aluminum oxide interface in close connection to experiment. The first project presented in this work considers the characterization of dissociatively adsorbed deuterated water on the Al-terminated (0001) surface. This particular structure is known from both experiment and theory to be the thermodynamically most stable surface termination of α-alumina in Ultra-High Vacuum (UHV) conditions. Based on experiments performed by our colleagues at FHI, different adsorption sites and products have been proposed and identified for D 2 O. While previous theoretical investigations only looked at vibrational frequencies of dissociated OD groups by staticNormal Modes Analysis (NMA), we rather employed a more sophisticated approach to directly assess vibrational spectra (like IR and VSF) at finite temperature from AIMD. In this work, we have employed a recent implementation which makes use of velocity-velocity autocorrelation functions to simulate such spectral responses of O-H(D) bonds. This approach allows for an efficient and qualitatively accurate estimation of Vibrational Densities of States (VDOS) as well as IR and VSF spectra, which are then tested against experimental spectra from our collaborators. In order to extend previous work on unimolecularly dissociated water on α-Al 2 O 3 , we then considered a different system, namely, a fully hydroxylated (0001) surface, which results from the reconstruction of the UHV-stable Al-terminated surface at high water contents. This model is then further extended by considering a hydroxylated surface with additional water molecules, forming a two-dimensional layer which serves as a potential template to simulate an aqueous interface in environmental conditions. Again, employing finite-temperature AIMD trajectories at the PBE+D2 level, we investigated the behaviour of both hydroxylated surface (HS) and the water-covered structure derived from it (known as HS+2ML). A full range of spectra, from VDOS to IR and VSF, is then calculated using the same methodology, as described above. This is the main focus of the second project, reported in Chapter 5. In this case, comparison between theoretical spectra and experimental data is definitely good. In particular, we underline the nature of high-frequency resonances observed above 3700 cm -1 in VSF experiments to be associated with surface OH-groups, known as "aluminols" which are a key fingerprint of the fully hydroxylated surface. In the third and last project, which is presented in Chapter 6, the extension of VSF spectroscopy experiments to the time-resolved regime offered us the opportunity to investigate vibrational energy relaxation at the α-alumina / water interface. Specifically, using again DFT-based AIMD simulations, we simulated vibrational lifetimes for surface aluminols as experimentally detected via pump-probe VSF. We considered the water-covered HS model as a potential candidate to address this problem. The vibrational (IR) excitation and subsequent relaxation is performed by means of a non-equilibrium molecular dynamics scheme. In such a scheme, we specifically looked at the O-H stretching mode of surface aluminols. Afterwards, the analysis of non-equilibrium trajectories allows for an estimation of relaxation times in the order of 2-4 ps which are in overall agreement with measured ones. The aim of this work has been to provide, within a consistent theoretical framework, a better understanding of vibrational spectroscopy and dynamics for water on the α-alumina(0001) surface,ranging from very low water coverage (similar to the UHV case) up to medium-high coverages, resembling the hydroxylated oxide in environmental moist conditions.}, language = {en} }