@article{FischerWertherBouaklineetal.2022,
  author    = {Fischer, Eric Wolfgang and Werther, Michael and Bouakline, Foudhil and Grossmann, Frank and Saalfrank, Peter},
  title     = {Non-Markovian vibrational relaxation dynamics at surfaces},
  series = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr},
  volume    = {156},
  journal   = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr},
  number    = {21},
  publisher = {AIP Publishing},
  address   = {Melville},
  issn      = {0021-9606},
  doi       = {10.1063/5.0092836},
  pages     = {16},
  year      = {2022},
  abstract  = {Vibrational dynamics of adsorbates near surfaces plays both an important role for applied surface science and as a model lab for studying fundamental problems of open quantum systems. We employ a previously developed model for the relaxation of a D-Si-Si bending mode at a D:Si(100)-(2 x 1) surface, induced by a "bath " of more than 2000 phonon modes [Lorenz and P. Saalfrank, Chem. Phys. 482, 69 (2017)], to extend previous work along various directions. First, we use a Hierarchical Effective Mode (HEM) model [Fischer et al., J. Chem. Phys. 153, 064704 (2020)] to study relaxation of higher excited vibrational states than hitherto done by solving a high-dimensional system-bath time-dependent Schrodinger equation (TDSE). In the HEM approach, (many) real bath modes are replaced by (much less) effective bath modes. Accordingly, we are able to examine scaling laws for vibrational relaxation lifetimes for a realistic surface science problem. Second, we compare the performance of the multilayer multiconfigurational time-dependent Hartree (ML-MCTDH) approach with that of the recently developed coherent-state-based multi-Davydov-D2 Ansatz [Zhou et al., J. Chem. Phys. 143, 014113 (2015)]. Both approaches work well, with some computational advantages for the latter in the presented context. Third, we apply open-system density matrix theory in comparison with basically "exact " solutions of the multi-mode TDSEs. Specifically, we use an open-system Liouville-von Neumann (LvN) equation treating vibration-phonon coupling as Markovian dissipation in Lindblad form to quantify effects beyond the Born-Markov approximation. Published under an exclusive license by AIP Publishing.},
  language  = {en}
}
@article{FischerWertherBouaklineetal.2020,
  author    = {Fischer, Eric W. and Werther, Michael and Bouakline, Foudhil and Saalfrank, Peter},
  title     = {A hierarchical effective mode approach to phonon-driven multilevel vibrational relaxation dynamics at surfaces},
  series = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistry},
  volume    = {153},
  journal   = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistry},
  number    = {6},
  publisher = {American Institute of Physics},
  address   = {Melville},
  issn      = {0021-9606},
  doi       = {10.1063/5.0017716},
  pages     = {15},
  year      = {2020},
  abstract  = {We discuss an efficient Hierarchical Effective Mode (HEM) representation of a high-dimensional harmonic oscillator bath, which describes phonon-driven vibrational relaxation of an adsorbate-surface system, namely, deuterium adsorbed on Si(100). Starting from the original Hamiltonian of the adsorbate-surface system, the HEM representation is constructed via iterative orthogonal transformations, which are efficiently implemented with Householder matrices. The detailed description of the HEM representation and its construction are given in the second quantization representation. The hierarchical nature of this representation allows access to the exact quantum dynamics of the adsorbate-surface system over finite time intervals, controllable via the truncation order of the hierarchy. To study the convergence properties of the effective mode representation, we solve the time-dependent Schrodinger equation of the truncated system-bath HEM Hamiltonian, with the help of the multilayer extension of the Multiconfigurational Time-Dependent Hartree (ML-MCTDH) method. The results of the HEM representation are compared with those obtained with a quantum-mechanical tier-model. The convergence of the HEM representation with respect to the truncation order of the hierarchy is discussed for different initial conditions of the adsorbate-surface system. The combination of the HEM representation with the ML-MCTDH method provides information on the time evolution of the system (adsorbate) and multiple effective modes of the bath (surface). This permits insight into mechanisms of vibration-phonon coupling of the adsorbate-surface system, as well as inter-mode couplings of the effective bath.},
  language  = {en}
}