@phdthesis{Rietze2020,
  author    = {Rietze, Clemens},
  title     = {Optimierung und Analyse von molekularen Schaltern in komplexen Umgebungen: thermische Stabilit{\"a}t, Auslesbarkeit und Schaltbarkeit},
  doi       = {10.25932/publishup-45959},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-459594},
  school      = {Universit{\"a}t Potsdam},
  pages     = {131},
  year      = {2020},
  abstract  = {Seit Jahrzehnten stellen die molekularen Schalter ein wachsendes Forschungsgebiet dar. Im Rahmen dieser Dissertation stand die Verbesserung der thermischen Stabilit{\"a}t, der Auslesbarkeit und Schaltbarkeit dieser molekularen Schalter in komplexen Umgebungen mithilfe computergest{\"u}tzter Chemie im Vordergrund. Im ersten Projekt wurde die Kinetik der thermischen E → Z-Isomerisierung und die damit verbundene thermische Stabilit{\"a}t eines Azobenzol-Derivats untersucht. Daf{\"u}r wurde Dichtefunktionaltheorie (DFT) in Verbindung mit der Eyring-Theorie des {\"U}bergangszustandes (TST) angewendet. Das Azobenzol-Derivat diente als vereinfachtes Modell f{\"u}r das Schalten in einer komplexen Umgebung (hier in metallorganischen Ger{\"u}sten). Es wurden thermodynamische und kinetische Gr{\"o}ßen unter verschiedenen Einfl{\"u}ssen berechnet, wobei gute {\"U}bereinstimmungen mit dem Experiment gefunden wurden. Die hier verwendete Methode stellte einen geeigneten Ansatz dar, um diese Gr{\"o}ßen mit angemessener Genauigkeit vorherzusagen. Im zweiten Projekt wurde die Auslesbarkeit der Schaltzust{\"a}nde in Form des nichtlinearen optischen (NLO) Kontrastes f{\"u}r die Molek{\"u}lklasse der Fulgimide untersucht. Die daf{\"u}r ben{\"o}tigten dynamischen Hyperpolarisierbarkeiten unter Ber{\"u}cksichtigung der Elektronenkorrelation wurden mittels einer etablierten Skalierungsmethode berechnet. Es wurden verschiedene Fulgimide analysiert, wobei viele experimentelle Befunde best{\"a}tigt werden konnten. Dar{\"u}ber hinaus legte die theoretische Vorhersage f{\"u}r ein weiteres System nahe, dass insbesondere die Erweiterung des π-Elektronensystems ein vielversprechender Ansatz zur Verbesserung von NLO-Kontrasten darstellt. Die Fulgimide verf{\"u}gen somit {\"u}ber n{\"u}tzliche Eigenschaften, sodass diese in Zukunft als Bauelemente in photonischen und optoelektronischen Bereichen Anwendungen finden k{\"o}nnten. Im dritten Projekt wurde die E → Z-Isomerisierung auf ein quantenmechanisch (QM) behandeltes Dimer mit molekularmechanischer (MM) Umgebung und zwei Fluorazobenzol-Monomeren durch Molek{\"u}ldynamik simuliert. Dadurch wurde die Schaltbarkeit in komplexer Umgebung (hier selbstorgansierte Einzelschichten = SAMs) bzw. von Azobenzolderivaten analysiert. Mit dem QM/MM Modell wurden sowohl Van-der-Waals-Interaktionen mit der Umgebung als auch elektronische Kopplung (nur zwischen QM-Molek{\"u}len) ber{\"u}cksichtigt. Dabei wurden systematische Untersuchungen zur Packungsdichte durchgef{\"u}hrt. Es zeigte sich, dass bereits bei einem Molek{\"u}labstand von 4.5 {\AA} die Quantenausbeute (prozentuale Anzahl erfolgreicher Schaltprozesse) des Monomers erreicht wird. Die gr{\"o}ßten Quantenausbeuten wurden f{\"u}r die beiden untersuchten Fluorazobenzole erzielt. Es wurden die Effekte des Molek{\"u}labstandes und der Einfluss von Fluorsubstituenten auf die Dynamik eingehend untersucht, sodass der Weg f{\"u}r darauf aufbauende Studien geebnet ist.},
  language  = {de}
}
@phdthesis{Kar,
  author    = {Kar, Manaswita},
  title     = {Energy band gap tuning of halide perovskite materials from first principles},
  school      = {Universit{\"a}t Potsdam},
  abstract  = {Solar cells based on hybrid perovskites materials have become significantly important among the third generation photovoltaics over the last few years. The first solid state solar cell was reported in 2012. Over the years, the power conversion efficiencies of these devices have increased at a tremendous pace and this has made the perovskite solar cell devices a serious competitor in the well-established market of thin-film and wafer technologies. Over time, a large number of articles on this topic has been published in peer-reviewed journals. The presence of lead in the most efficient hybrid perovskite materials have raised questions about the possible toxicity of these devices and the extent of their environmental impact. Therefore, a lot of research has been devoted to finding alternative perovskite materials with similar or even better opto-electronic properties. An alternative strategy to improve the efficiency of thin film solar cells is to build efficient tandem cells by combining two or more perovskite materials with specifically tailored band gaps. The first step towards the development of perovskite-only tandem solar cells is to identify complementary hybrid perovskite materials with specific band gaps that maximize the efficiency of tandem solar cells. The optimal set of optical gaps for a tandem structure made of two materials is 1.9 eV and 1.0 eV. Since the electronic properties of hybrid perovskites are known to be strongly dependent on the composition and distortion of the crystal lattice, strong focus has been made towards the structure optimisation as well as the calculation of the energy band gaps of the materials using density functional theory (DFT). In an attempt to study the structure-property relationship of these perovskite materials and to find novel perovskite materials for future applications, researchers have employed computational screening procedures to study a large range of these materials by systematic replacement of the cations and anions from the prototypical perovskite. Density functional theory in particular is used as a theoretical tool, because of it's precision to determine the properties of materials and also it's computational viability in dealing with complex systems. In this thesis, the main focus is to do a systematic screening of the perovskite materials, of the composition ABX3 again by replacing the A-site, B-site and the X-site elements to find novel materials with band gaps suitable for application in tandem solar cells. As a first step towards contributing to this vibrant field of research, a high-throughput computational screening has been performed by replacing the metal and the halogen in the conventional CH3NH3PbI3 perovskites with homovalent metals and halogens to find materials in the desired range of band gaps that has already been mentioned earlier. This is achieved by performing a geometry optimisation on all the simulated structures followed by calculating their energy band gaps at the semilocal and the hybrid levels of theory. However, it is well known that the rotation of the organic cation CH3NH3 hinders the stability of these devices by the formation of hydrogen bonds between the hydrogen atoms of the cation and the halogens. This causes the materials to degrade under normal temperature and pressure conditions. As an attempt to prevent these devices from being unstable, a next step has been taken where the CH3NH3 cation has been replaced by inorganic cations of similar ionic radius. This is followed by another thorough screening, similar to the previous step. The stability of the materials has been determined by using the empirical Goldschmidt tolerance factor. As a last part of the thesis, a small proportion of the inorganic cation is mixed with CH3NH3 in order to form mixed-halide perovskites. These structures are optimised and their band gaps are calculated using density functional theory in order to predict materials suitable for single junction as well as tandem solar cell devices. It is expected that the contribution made through this thesis will be helpful for the progress of perovskite solar cells in terms of efficiencies and will also allow the community to explore the different properties these materials for further progress and development.},
  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}
}
@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}
}
@phdthesis{Lever2022,
  author    = {Lever, Fabiano},
  title     = {Probing the ultrafast dynamics of 2-Thiouracil with soft x-rays},
  doi       = {10.25932/publishup-55523},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-555230},
  school      = {Universit{\"a}t Potsdam},
  pages     = {129},
  year      = {2022},
  abstract  = {Understanding the changes that follow UV-excitation in thionucleobases is of great importance for the study of light-induced DNA lesions and, in a broader context, for their applications in medicine and biochemistry. Their ultrafast photophysical reactions can alter the chemical structure of DNA - leading to damages to the genetic code - as proven by the increased skin cancer risk observed for patients treated with thiouracil for its immunosuppressant properties. In this thesis, I present four research papers that result from an investigation of the ultrafast dynamics of 2-thiouracil by means of ultrafast x-ray probing combined with electron spectroscopy. A molecular jet in the gas phase is excited with a uv pulse and then ionized with x-ray radiation from a Free Electron Laser. The kinetic energy of the emitted electrons is measured in a magnetic bottle spectrometer. The spectra of the measured photo and Auger electrons are used to derive a picture of the changes in the geometrical and electronic configurations. The results allow us to look at the dynamical processes from a new perspective, thanks to the element- and site- sensitivity of x-rays. The custom-built URSA-PQ apparatus used in the experiment is described. It has been commissioned and used at the FL24 beamline of the FLASH2 FEL, showing an electron kinetic energy resolution of ∆E/E ~ 40 and a pump-probe timing resolution of 190 f s. X-ray only photoelectron and Auger spectra of 2-thiouracil are extracted from the data and used as reference. Photoelectrons following the formation a 2p core hole are identified, as well as resonant and non-resonant Auger electrons. At the L 1 edge, Coster-Kronig decay is observed from the 2s core hole. The UV-induced changes in the 2p photoline allow the study the electronic-state dynamics. With the use of an Excited-State Chemical Shift (ESCS) model, we observe a ultrafast ground-state relaxation within 250 f s. Furthermore, an oscillation with a 250 f s period is observed in the 2p binding energy, showing a coherent population exchange between electronic states. Auger electrons from the 2p core hole are analyzed and used to deduce a ultrafast C -S bond expansion on a sub 100 f s scale. A simple Coulomb-model, coupled to quantum chemical calculations, can be used to infer the geometrical changes in the molecular structure.},
  language  = {en}
}
@phdthesis{Fischer2022,
  author    = {Fischer, Eric Wolfgang},
  title     = {Quantum vibrational dynamics in complex environments: from vibrational strong coupling in molecular cavity QED to phonon-induced adsorbate relaxation},
  doi       = {10.25932/publishup-56721},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-567214},
  school      = {Universit{\"a}t Potsdam},
  pages     = {viii, 171},
  year      = {2022},
  abstract  = {Molecules are often naturally embedded in a complex environment. As a consequence, characteristic properties of a molecular subsystem can be substantially altered or new properties emerge due to interactions between molecular and environmental degrees of freedom. The present thesis is concerned with the numerical study of quantum dynamical and stationary properties of molecular vibrational systems embedded in selected complex environments. In the first part, we discuss "strong-coupling" model scenarios for molecular vibrations interacting with few quantized electromagnetic field modes of an optical Fabry-P{\´e}rot cavity. We thoroughly elaborate on properties of emerging "vibrational polariton" light-matter hybrid states and examine the relevance of the dipole self-energy. Further, we identify cavity-induced quantum effects and an emergent dynamical resonance in a cavity-altered thermal isomerization model, which lead to significant suppression of thermal reaction rates. Moreover, for a single rovibrating diatomic molecule in an optical cavity, we observe non-adiabatic signatures in dynamics due to "vibro-polaritonic conical intersections" and discuss spectroscopically accessible "rovibro-polaritonic" light-matter hybrid states. In the second part, we study a weakly coupled but numerically challenging quantum mechanical adsorbate-surface model system comprising a few thousand surface modes. We introduce an efficient construction scheme for a "hierarchical effective mode" approach to reduce the number of surface modes in a controlled manner. In combination with the multilayer multiconfigurational time-dependent Hartree (ML-MCTDH) method, we examine the vibrational adsorbate relaxation dynamics from different excited adsorbate states by solving the full non-Markovian system-bath dynamics for the characteristic relaxation time scale. We examine half-lifetime scaling laws from vibrational populations and identify prominent non-Markovian signatures as deviations from Markovian reduced system density matrix theory in vibrational coherences, system-bath entanglement and energy transfer dynamics. In the final part of this thesis, we approach the dynamics and spectroscopy of vibronic model systems at finite temperature by formulating the ML-MCTDH method in the non-stochastic framework of thermofield dynamics. We apply our method to thermally-altered ultrafast internal conversion in the well-known vibronic coupling model of pyrazine. Numerically beneficial representations of multilayer wave functions ("ML-trees") are identified for different temperature regimes, which allow us to access thermal effects on both electronic and vibrational dynamics as well as spectroscopic properties for several pyrazine models.},
  language  = {en}
}
@phdthesis{Bouakline2023,
  author    = {Bouakline, Foudhil},
  title     = {Manifestations of Quantum-Mechanical Effects in Molecular Reaction Dynamics},
  school      = {Universit{\"a}t Potsdam},
  pages     = {316},
  year      = {2023},
  abstract  = {This habilitation thesis summarises the research work performed by the author during the last quindecennial period. The dissertation reflects his main research interests, which revolve around quantum dynamics of small-sized molecular systems, including their interactions with electromagnetic radiation or dissipative environments. This covers various dynamical processes that involve bound-bound, bound-free, and free-free molecular transitions. The latter encompass light-triggered rovibrational or rovibronic dynamics in bound molecules, molecular photodissociation induced by weak or strong laser fields, state-to-state reactive and/or inelastic molecular collisions, and phonon-driven vibrational relaxation of adsorbates at solid surfaces. Although the dissertation covers different topics of molecular reaction dynamics, most of these studies focus on nuclear quantum effects and their manifestations in experimental measures. The latter are assessed through comparison between quantum and classical predictions, and/or direct confrontation of theory and experiment. Most well known quantum concepts and effects will be encountered in this work. Yet, almost all these quantum notions find their roots in the central pillar of quantum theory, namely, the quantum superposition principle. Indeed, quantum coherence is the main source of most quantum effects, including interference, entanglement, and even tunneling. Thus, the common and predominant theme of all the investigations of this thesis is quantum coherence, and the survival or quenching of subsequent interference effects in various molecular processes. The lion's share of the dissertation is devoted to two associated quantum concepts, which are usually overlooked in computational molecular dynamics, viz. the Berry phase and identical nuclei symmetry. The importance of the latter in dynamical molecular processes and their direct fingerprints in experimental observables also rely very much on quantum coherence and entanglement. All these quantum phenomena are thoroughly discussed within the four main topics that form the core of this thesis. Each topic is described in a separate chapter, where it is briefly summarised and then illustrated with three peer-reviewed publications. The first topic deals with the relevance of quantum coherence/interference in molecular collisions, with a focus on the hydrogen-exchange reaction, H+H2 --> H2+H, and its isotopologues. For these collision processes, the significance of interference of probability amplitudes arises because of the existence of two main scattering pathways. The latter could be inelastic and reactive scattering, direct and time-delayed scattering, or two encircling reaction paths that loop in opposite senses around a conical intersection (CI) of the H3 molecular system. Our joint theoretical-experimental investigations of these processes reveal strong interference and geometric phase (GP) effects in state-to-state reaction probabilities and differential cross sections. However, these coherent effects completely cancel in integral cross sections and reaction rate constants, due to efficient dephasing of interference between the different scattering amplitudes. As byproducts of these studies, we highlight the discovery of two novel scattering mechanisms, which contradict conventional textbook pictures of molecular reaction dynamics. The second topic concerns the effect of the Berry phase on molecular photodynamics at conical intersections. To understand this effect, we developed a topological approach that separates the total molecular wavefunction of an unbound molecular system into two components, which wind in opposite senses around the conical intersection. This separation reveals that the only effect of the geometric phase is to change the sign of the relative phase of these two components. This in turn leads to a shift in the interference pattern of the molecular system---a phase shift that is reminiscient of the celebrated Aharonov-Bohm effect. This procedure is numerically illustrated with photodynamics at model standard CIs, as well as strong-field dissociation of diatomics at light-induced conical intersections (LICIs). Besides the fundamental aspect of these studies, their findings allow to interpret and predict the effect of the GP on the state-resolved or angle-resolved spectra of pump-probe experimental schemes, particularly the distributions of photofragments in molecular photodissociation experiments. The third topic pertains to the role of the indistinguishability of identical nuclei in molecular reaction dynamics, with an emphasis on dynamical localization in highly symmetric molecules. The main object of these studies is whether nuclear-spin statistics allow dynamical localization of the electronic, vibrational, or even rotational density on a specific molecular substructure or configuration rather than on another one which is identical (indistinguishable). Group-theoretic analysis of the symmetrized molecular wavefunctions of these systems shows that nuclear permutation symmetry engenders quantum entanglement between the eigenstates of the different molecular degrees of freedom. This subsequently leads to complete quenching of dynamical localization over indistinguishable molecular substructures---an observation that is in sharp contradiction with well known textbook views of iconic molecular processes. This is illustrated with various examples of quantum dynamics in symmetric double-well achiral molecules, such as the prototypical umbrella inversion motion of ammonia, electronic Kekul{\´e} dynamics in the benzene molecule, and coupled electron-nuclear dynamics in laser-induced indirect photodissociation of the dihydrogen molecular cation. The last part of the thesis is devoted to the development of approximate wavefunction approaches for phonon-induced vibrational relaxation of adsorbates (system) at surfaces (bath). Due to the so-called 'curse of dimensionality', these system-bath complexes cannot be handled with standard wavefunction methods. To alleviate the exponential scaling of the latter, we developed approximate yet quite accurate numerical schemes that have a polynomial scaling with respect to the bath dimensionality. The corresponding algorithms combine symmetry-based reductions of the full vibrational Hilbert space and iterative Krylov techniques. These approximate wavefunction approaches resemble the 'Bixon-Jortner model' and the more general 'quantum tier model'. This is illustrated with the decay of H-Si (D-Si) vibrations on a fully H(D)-covered silicon surface, which is modelled with a phonon-bath of more than two thousand oscillators. These approximate methods allow reliable estimation of the adsorbate vibrational lifetimes, and provide some insight into vibration-phonon couplings at solid surfaces. Although this topic is mainly computational, the developed wavefunction approaches permit to describe quantum entanglement between the system and bath states, and to embody some coherent effects in the time-evolution of the (sub-)system, which cannot be accounted for with the widely used 'reduced density matrix formalism'.},
  language  = {en}
}