530 Physik
Refine
Language
- English (13) (remove)
Is part of the Bibliography
- yes (13)
Keywords
- Computational chemistry (1)
- Energy (1)
- Molecules (1)
- Monomers (1)
- Oligomers (1)
- density functional theory (1)
- oxidation state (1)
- porphyrin (1)
- scanning tunneling microscopy (1)
- scanning tunneling spectroscopy (1)
Institute
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.
We introduce a thermofield-based formulation of the multilayer multiconfigurational time-dependent Hartree (MCTDH) method to study finite temperature effects on non-adiabatic quantum dynamics from a non-stochastic, wave function perspective. Our approach is based on the formal equivalence of bosonic many-body theory at zero temperature with a doubled number of degrees of freedom and the thermal quasi-particle representation of bosonic thermofield dynamics (TFD). This equivalence allows for a transfer of bosonic many-body MCTDH as introduced by Wang and Thoss to the finite temperature framework of thermal quasi-particle TFD. As an application, we study temperature effects on the ultrafast internal conversion dynamics in pyrazine. We show that finite temperature effects can be efficiently accounted for in the construction of multilayer expansions of thermofield states in the framework presented herein. Furthermore, we find our results to agree well with existing studies on the pyrazine model based on the pMCTDH method.
It has been experimentally demonstrated that reaction rates for molecules embedded in microfluidic optical cavities are altered when compared to rates observed under "ordinary" reaction conditions. However, precise mechanisms of how strong coupling of an optical cavity mode to molecular vibrations affects the reactivity and how resonance behavior emerges are still under dispute. In the present work, we approach these mechanistic issues from the perspective of a thermal model reaction, the inversion of ammonia along the umbrella mode, in the presence of a single-cavity mode of varying frequency and coupling strength. A topological analysis of the related cavity Born-Oppenheimer potential energy surface in combination with quantum mechanical and transition state theory rate calculations reveals two quantum effects, leading to decelerated reaction rates in qualitative agreement with experiments: the stiffening of quantized modes perpendicular to the reaction path at the transition state, which reduces the number of thermally accessible reaction channels, and the broadening of the barrier region, which attenuates tunneling. We find these two effects to be very robust in a fluctuating environment, causing statistical variations of potential parameters, such as the barrier height. Furthermore, by solving the time-dependent Schrodinger equation in the vibrational strong coupling regime, we identify a resonance behavior, in qualitative agreement with experimental and earlier theoretical work. The latter manifests as reduced reaction probability when the cavity frequency omega(c) is tuned resonant to a molecular reactant frequency. We find this effect to be based on the dynamical localization of the vibro-polaritonic wavepacket in the reactant well.
The oxidation and spin state of a metal-organic molecule determine its chemical reactivity and magnetic properties. Here, we demonstrate the reversible control of the oxidation and spin state in a single Fe porphyrin molecule in the force field of the tip of a scanning electron tunneling microscope. Within the regimes of half-integer and integer spin state, we can further track the evolution of the magnetocrystalline anisotropy. Our experimental results are corroborated by density functional theory and wave function theory. This combined analysis allows us to draw a complete picture of the molecular states over a large range of intramolecular deformations.
We perform ab initio molecular dynamics simulations to study the femtosecond laser induced desorption of H-2, D-2, and HD from a H: D-saturated Ru(0001) surface. To this aim we have extended the ab initio molecular dynamics with electronic friction (AIMDEF) scheme to include a random force that is a function of a timedependent electronic temperature. The latter characterizes the action of the ultrashort laser pulse according to a two temperature model. This allows us to perform multidimensional, hot-electron driven reaction dynamics and investigate the dependence of the desorption yields on the relative H: D isotope concentration on the surface. Our AIMDEF simulations show that the desorption process takes place in the presence of a heated adsorbate system that clearly influences the desorption dynamics. The heating of the adsorbate system is more (less) pronounced the larger is the concentration of the lighter (heavier) isotope. As a result, we conclude that the presence of H on the surface favors the desorption of molecules, whereas the presence of D hampers it, in agreement with previous experimental observations in which the phenomenon of "dynamical promotion" of a surface reaction had been postulated.
Following excited-state chemical shifts in molecular ultrafast x-ray photoelectron spectroscopy
(2022)
The conversion of photon energy into other energetic forms in molecules is accompanied by charge moving on ultrafast timescales. We directly observe the charge motion at a specific site in an electronically excited molecule using time-resolved x-ray photoelectron spectroscopy (TR-XPS). We extend the concept of static chemical shift from conventional XPS by the excited-state chemical shift (ESCS), which is connected to the charge in the framework of a potential model. This allows us to invert TR-XPS spectra to the dynamic charge at a specific atom. We demonstrate the power of TR-XPS by using sulphur 2p-core-electron-emission probing to study the UV-excited dynamics of 2-thiouracil. The method allows us to discover that a major part of the population relaxes to the molecular ground state within 220–250 fs. In addition, a 250-fs oscillation, visible in the kinetic energy of the TR-XPS, reveals a coherent exchange of population among electronic states.
Following excited-state chemical shifts in molecular ultrafast x-ray photoelectron spectroscopy
(2022)
The conversion of photon energy into other energetic forms in molecules is accompanied by charge moving on ultrafast timescales. We directly observe the charge motion at a specific site in an electronically excited molecule using time-resolved x-ray photoelectron spectroscopy (TR-XPS). We extend the concept of static chemical shift from conventional XPS by the excited-state chemical shift (ESCS), which is connected to the charge in the framework of a potential model. This allows us to invert TR-XPS spectra to the dynamic charge at a specific atom. We demonstrate the power of TR-XPS by using sulphur 2p-core-electron-emission probing to study the UV-excited dynamics of 2-thiouracil. The method allows us to discover that a major part of the population relaxes to the molecular ground state within 220–250 fs. In addition, a 250-fs oscillation, visible in the kinetic energy of the TR-XPS, reveals a coherent exchange of population among electronic states.
Recent experiments and theory suggest that ground state properties and reactivity of molecules can be modified when placed inside a nanoscale cavity, giving rise to strong coupling between vibrational modes and the quantized cavity field. This is commonly thought to be caused either by a cavity-distorted Born-Oppenheimer ground state potential or by the formation of light-matter hybrid states, vibrational polaritons. Here, we systematically study the effect of a cavity on ground state properties and infrared spectra of single molecules, considering vibration-cavity coupling strengths from zero up to the vibrational ultrastrong coupling regime. Using single-mode models for Li-H and O-H stretch modes and for the NH3 inversion mode, respectively, a single cavity mode in resonance with vibrational transitions is coupled to position-dependent molecular dipole functions. We address the influence of the cavity mode on polariton ground state energies, equilibrium bond lengths, dissociation energies, activation energies for isomerization, and on vibro-polaritonic infrared spectra. In agreement with earlier work, we observe all mentioned properties being strongly affected by the cavity, but only if the dipole self-energy contribution in the interaction Hamiltonian is neglected. When this term is included, these properties do not depend significantly on the coupling anymore. Vibro-polaritonic infrared spectra, in contrast, are always affected by the cavity mode due to the formation of excited vibrational polaritons. It is argued that the quantized nature of vibrational polaritons is key to not only interpreting molecular spectra in cavities but also understanding the experimentally observed modification of molecular reactivity in cavities.
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.
Recently, Nocera and co-workers (J. Am. Chem. Soc. 2018, 140, 13711) demonstrated that triaryl borate Lewis acids facilitate the direct electrochemical reduction of triphenylphosphine oxide (TPPO) to triphenylphosphine (TPP). In the present contribution, we report a quantum chemical study unravelling details of the reaction, which also supports the proposed ECrECi mechanism. Alternative electrochemical routes to TPPO reduction facilitated by other Lewis acids (CH3+), or by photocatalysis at semiconductor surfaces, are also briefly discussed.