TY - JOUR A1 - Fischer, Eric Wolfgang A1 - Werther, Michael A1 - Bouakline, Foudhil A1 - Grossmann, Frank A1 - Saalfrank, Peter T1 - Non-Markovian vibrational relaxation dynamics at surfaces JF - The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr N2 - 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. KW - phonons KW - Vibrational states KW - Chemical dynamics KW - Adsorption KW - Surface science KW - Open quantum systems KW - Density-matrix KW - Coherent states KW - Markov processes Y1 - 2022 U6 - https://doi.org/10.1063/5.0092836 SN - 0021-9606 SN - 1089-7690 SN - 1520-9032 VL - 156 IS - 21 PB - AIP Publishing CY - Melville ER - TY - JOUR A1 - Bouakline, Foudhil A1 - Saalfrank, Peter T1 - Seemingly asymmetric atom-localized electronic densities following laser-dissociation of homonuclear diatomics JF - The journal of chemical physics : bridges a gap between journals of physics and journals of chemistry N2 - Recent experiments on laser-dissociation of aligned homonuclear diatomic molecules show an asymmetric forward-backward (spatial) electron-localization along the laser polarization axis. Most theoretical models attribute this asymmetry to interference effects between gerade and ungerade vibronic states. Presumably due to alignment, these models neglect molecular rotations and hence infer an asymmetric (post-dissociation) charge distribution over the two identical nuclei. In this paper, we question the equivalence that is made between spatial electron-localization, observed in experiments, and atomic electron-localization, alluded by these theoretical models. We show that (seeming) agreement between these models and experiments is due to an unfortunate omission of nuclear permutation symmetry, i.e., quantum statistics. Enforcement of the latter requires mandatory inclusion of the molecular rotational degree of freedom, even for perfectly aligned molecules. Unlike previous interpretations, we ascribe spatial electron-localization to the laser creation of a rovibronic wavepacket that involves field-free molecular eigenstates with opposite space-inversion symmetry i.e., even and odd parity. Space-inversion symmetry breaking would then lead to an asymmetric distribution of the (space-fixed) electronic density over the forward and backward hemisphere. However, owing to the simultaneous coexistence of two indistinguishable molecular orientational isomers, our analytical and computational results show that the post-dissociation electronic density along a specified space-fixed axis is equally shared between the two identical nuclei-a result that is in perfect accordance with the principle of the indistinguishability of identical particles. Published under an exclusive license by AIP Publishing. Y1 - 2021 U6 - https://doi.org/10.1063/5.0049710 SN - 0021-9606 SN - 1089-7690 VL - 154 IS - 23 PB - American Institute of Physics CY - Melville ER - TY - JOUR A1 - Witzorky, Christoph A1 - Paramonov, Guennaddi A1 - Bouakline, Foudhil A1 - Jaquet, Ralph A1 - Saalfrank, Peter A1 - Klamroth, Tillmann T1 - Gaussian-type orbital calculations for high harmonic generation in vibrating molecules BT - Benchmarks for H-2(+) JF - Journal of chemical theory and computation N2 - The response of the hydrogen molecular ion, H-2(+), to few-cycle laser pulses of different intensities is simulated. To treat the coupled electron-nuclear motion, we use adiabatic potentials computed with Gaussian-type basis sets together with a heuristic ionization model for the electron and a grid representation for the nuclei. Using this mixed-basis approach, the time-dependent Schrodinger equation is solved, either within the Born-Oppenheimer approximation or with nonadiabatic couplings included. The dipole response spectra are compared to all-grid-based solutions for the three-body problem, which we take as a reference to benchmark the Gaussian-type basis set approaches. Also, calculations employing the fixed-nuclei approximation are performed, to quantify effects due to nuclear motion. For low intensities and small ionization probabilities, we get excellent agreement of the dynamics using Gaussian-type basis sets with the all-grid solutions. Our investigations suggest that high harmonic generation (HHG) and high-frequency response, in general, can be reliably modeled using Gaussian-type basis sets for the electrons for not too high harmonics. Further, nuclear motion destroys electronic coherences in the response spectra even on the time scale of about 30 fs and affects HHG intensities, which reflect the electron dynamics occurring on the attosecond time scale. For the present system, non-Born-Oppenheimer effects are small. The Gaussian-based, nonadiabatically coupled, time-dependent multisurface approach to treat quantum electron-nuclear motion beyond the non-Born-Oppenheimer approximation can be easily extended to approximate wavefunction methods, such as time-dependent configuration interaction singles (TD-CIS), for systems where no benchmarks are available. KW - Basis sets KW - Chemical calculations KW - Ionization KW - Lasers KW - Quantum mechanics Y1 - 2021 U6 - https://doi.org/10.1021/acs.jctc.1c00837 SN - 1549-9618 SN - 1549-9626 VL - 17 IS - 12 SP - 7353 EP - 7365 PB - American Chemical Society CY - Washington ER - TY - JOUR A1 - Fischer, Eric W. A1 - Werther, Michael A1 - Bouakline, Foudhil A1 - Saalfrank, Peter T1 - A hierarchical effective mode approach to phonon-driven multilevel vibrational relaxation dynamics at surfaces JF - The journal of chemical physics : bridges a gap between journals of physics and journals of chemistry N2 - 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. Y1 - 2020 U6 - https://doi.org/10.1063/5.0017716 SN - 0021-9606 SN - 1089-7690 VL - 153 IS - 6 PB - American Institute of Physics CY - Melville ER - TY - JOUR A1 - Bouakline, Foudhil A1 - Fischer, E. W. A1 - Saalfrank, Peter T1 - A quantum-mechanical tier model for phonon-driven vibrational relaxation dynamics of adsorbates at surfaces JF - The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr N2 - We present a quantum-mechanical tier model for vibrational relaxation of low-lying excited states of an adsorbate vibrational mode (system), coupled to surface phonons (bath), at zero temperature. The tier model, widely used in studies of intramolecular vibrational energy redistribution in polyatomics, is adapted here to adsorbate-surface systems with the help of an embedded cluster approach, using orthogonal coordinates for the system and bath modes, and a phononic expansion of their interaction. The key idea of the model is to organize the system-bath zeroth-order vibrational space into a hierarchical structure of vibrational tiers and keep therein only vibrational states that are sequentially generated from the system-bath initial vibrational state. Each tier is generated from the previous one by means of a successor operator, derived from the system-bath interaction Hamiltonian. This sequential procedure leads to a drastic reduction of the dimension of the system-bath vibrational space. We notably show that for harmonic vibrational motion of the system and linear system-bath couplings in the system coordinate, the dimension of the tier-model vibrational basis scales as similar to N-lxv. Here, N is the number of bath modes, l is the highest-order of the phononic expansion, and l is the size of the system vibrational basis. This polynomial scaling is computationally far superior to the exponential scaling of the original zeroth-order vibrational basis, similar to M-N, with M being the number of basis functions per bath mode. In addition, since each tier is coupled only to its adjacent neighbors, the matrix representation of the system-bath Hamiltonian in this new vibrational basis has a symmetric block-tridiagonal form, with each block being very sparse. This favors the combination of the tier-model with iterative Krylov techniques, such as the Lanczos algorithm, to solve the time-dependent Schrodinger equation for the full Hamiltonian. To illustrate the method, we study vibrational relaxation of a D-Si bending mode, coupled via two-and (mainly) one-phonon interactions to a fully D-covered Si(100)-(2 x 1) surface, using a recent first-principles system-bath Hamiltonian. The results of the tier model are compared with those obtained by the Lindblad formalism of the reduced density matrix. We find that the tier model provides much more information and insight into mechanisms of vibration-phonon couplings at surfaces, and gives more reliable estimates of the adsorbate vibrational lifetimes. Moreover, the tier model might also serve as a benchmark for other approximate quantum-dynamics methods, such as multiconfiguration wavefunction approaches. Published under license by AIP Publishing. Y1 - 2019 U6 - https://doi.org/10.1063/1.5099902 SN - 0021-9606 SN - 1089-7690 VL - 150 IS - 24 PB - American Institute of Physics CY - Melville ER - TY - JOUR A1 - Bouakline, Foudhil A1 - Lorenz, Ulrich J. A1 - Melani, Giacomo A1 - Paramonov, Guennaddi K. A1 - Saalfrank, Peter T1 - Isotopic effects in vibrational relaxation dynamics of H on a Si(100) surface JF - The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr N2 - 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ödinger equation. Y1 - 2017 U6 - https://doi.org/10.1063/1.4994635 SN - 0021-9606 SN - 1089-7690 VL - 147 IS - 14 PB - American Institute of Physics CY - Melville ER - TY - JOUR A1 - Bouakline, Foudhil A1 - Lüder, Franziska A1 - Martinazzo, Rocco A1 - Saalfrank, Peter T1 - Reduced and exact quantum dynamics of the vibrational relaxation of a molecular system interacting with a finite-dimensional bath JF - The journal of physical chemistry : A, Molecules, spectroscopy, kinetics, environment & general theory N2 - We investigate the vibrational relaxation of a Morse oscillator, nonlinearly coupled to a finite-dimensional bath of harmonic oscillators at zero temperature, using two different approaches: Reduced dynamics with the help of the Lindblad formalism of reduced density matrix theory in combination with Fermi's Golden Rule, and exact dynamics (within the chosen model). with the multiconfiguration time-dependent Hartree (MCTDH) method. Two different models have been constructed, the situation where the bath spectrum is exactly resonant with the anharmonic oscillator transition frequencies, and the case for which the subsystem is slightly off-resonant with the environment. At short times, reduced dynamics calculations describe the relaxation process qualitatively well but fail to reproduce recurrences observed with MCTDH for longer times. Lifetimes of all the vibrational levels of the Morse oscillator have been calculated, and both Lindblad and MCTDH. results show the same dependence of the lifetimes on the initial vibrational state quantum number. A prediction, which should be generic for adsorbate systems is a striking, sharp increase of lifetimes of the subsystem vibrational levels close to the dissociation This is contradictory with harmonic/linear extrapolation laws, which predict a monotonic decrease of the lifetime with initial vibrational quantum number. Y1 - 2012 U6 - https://doi.org/10.1021/jp304466u SN - 1089-5639 VL - 116 IS - 46 SP - 11118 EP - 11127 PB - American Chemical Society CY - Washington ER -