Refine
Year of publication
Document Type
- Article (148)
- Postprint (9)
- Part of a Book (1)
- Conference Proceeding (1)
- Other (1)
Language
- English (160)
Is part of the Bibliography
- yes (160)
Keywords
- photochemistry (3)
- simulations (3)
- absorbtion fine-structure (2)
- augmented-wave method (2)
- density functional calculations (2)
- graphite (2)
- hydrogen (2)
- initio molecular-dynamics (2)
- molecules (2)
- oxidation (2)
Institute
Somewhat surprisingly, inverted ("O-down") CO adsorbates on NaCl(100) were recently observed experimentally after infrared vibrational excitation (Lau et al., Science, 2020, 367, 175-178). Here we characterize these species using periodic density functional theory and a quantum mechanical description of vibrations. We determine stationary points and minimum energy paths for CO inversion, for low (1/8 and 1/4 monolayers (ML)) and high (1 ML) coverages. Transition state theory is applied to estimate thermal rates for "C-down" to "O-down" isomerization and the reverse process. For the 1/4 ML p(1 x 1) structure, two-dimensional and three-dimensional potential energy surfaces and corresponding anharmonic vibrational eigenstates obtained from the time-independent nuclear Schrodinger equation are presented. We find (i) rather coverage-independent CO inversion energies (of about 0.08 eV or 8 kJ mol(-1) per CO) and corresponding classical activation energies for "C-down" to "O-down" isomerization (of about 0.15 eV or 14 kJ mol(-1) per CO); (ii) thermal isomerization rates at 22 K which are vanishingly small for the "C-down" to "O-down" isomerization but non-negligible for the back reaction; (iii) several "accidentally degenerate" pairs of eigenstates well below the barrier, each pair describing "C-down" to "O-down" localized states.
X-ray photoelectron spectroscopy (XPS) is a powerful tool for probing the local chemical environment of atoms near surfaces. When applied to soft matter, such as polymers, XPS spectra are frequently shifted and broadened due to thermal atom motion and by interchain interactions. We present a combined quantum mechanical QM/molecular dynamics (MD) simulation of X-ray photoelectron spectra of polyvinyl alcohol (PVA) using oligomer models in order to account for and quantify these effects on the XPS (C1s) signal. In our study, molecular dynamics at finite temperature were performed with a classical forcefield and by ab initio MD (AIMD) using the Car-Parrinello method. Snapshots along, the trajectories represent possible conformers and/or neighbouring environments, with different C1s ionization potentials for individual C atoms leading to broadened XPS peaks. The latter are determined by Delta-Kohn Sham calculations. We also examine the experimental practice of gauging XPS (C1s) signals of alkylic C-atoms in C-containing polymers to the C1s signal of polyethylene.
We find that (i) the experimental XPS (C1s) spectra of PVA (position and width) can be roughly represented by single-strand models, (ii) interchain interactions lead to red-shifts of the XPS peaks by about 0.6 eV, and (iii) AIMD simulations match the findings from classical MD semi-quantitatively. Further, (iv) the gauging procedure of XPS (C1s) signals to the values of PE, introduces errors of about 0.5 eV. (C) 2014 Elsevier B.V. All rights reserved.
In Near Edge X-Ray Absorption Fine Structure (NEXAFS) spectroscopy X-Ray photons are used to excite tightly bound core electrons to low-lying unoccupied orbitals of the system. This technique offers insight into the electronic structure of the system as well as useful structural information. In this work, we apply NEXAFS to two kinds of imidazolium based ionic liquids ([C(n)C(1)im](+)[NTf2](-) and [C(4)C(1)im](+)[I](-)). A combination of measurements and quantum chemical calculations of C K and N K NEXAFS resonances is presented. The simulations, based on the transition potential density functional theory method (TP-DFT), reproduce all characteristic features observed by the experiment. Furthermore, a detailed assignment of resonance features to excitation centers (carbon or nitrogen atoms) leads to a consistent interpretation of the spectra.
In Near Edge X-Ray Absorption Fine Structure (NEXAFS) spectroscopy X-Ray photons are used to excite tightly bound core electrons to low-lying unoccupied orbitals of the system. This technique offers insight into the electronic structure of the system as well as useful structural information. In this work, we apply NEXAFS to two kinds of imidazolium based ionic liquids ([CnC1im]+[NTf2]- and [C4C1im]+[I]-). A combination of measurements and quantum chemical calculations of C K and N K NEXAFS resonances is presented. The simulations, based on the transition potential density functional theory method (TP-DFT), reproduce all characteristic features observed by the experiment. Furthermore, a detailed assignment of resonance features to excitation centers (carbon or nitrogen atoms) leads to a consistent interpretation of the spectra.
A cationic surfactant containing a spiropyrane unit is prepared exhibiting a dual-responsive adjustability of its surface-active characteristics. The switching mechanism of the system relies on the reversible conversion of the non-ionic spiropyrane (SP) to a zwitterionic merocyanine (MC) and can be controlled by adjusting the pH value and via light, resulting in a pH-dependent photoactivity: While the compound possesses a pronounced difference in surface activity between both forms under acidic conditions, this behavior is suppressed at a neutral pH level. The underlying switching processes are investigated in detail, and a thermodynamic explanation based on a combination of theoretical and experimental results is provided. This complex stimuli-responsive behavior enables remote-control of colloidal systems. To demonstrate its applicability, the surfactant is utilized for the pH-dependent manipulation of oil-in-water emulsions.
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.
A novel quantum method to deal with typical system-bath dynamical problems is introduced. Subsystem discrete variable representation and bath coherent-state sets are used to write down a multiconfigurational expansion of the wave function of the whole system. With the help of the Dirac-Frenkel variational principle, simple equations of motion-a kind of Schrodinger-Langevin equation for the subsystem coupled to (pseudo) classical equations for the bath-are derived. True dissipative dynamics at all times is obtained by coupling the bath to a secondary, classical Ohmic bath, which is modeled by adding a friction coefficient in the derived pseudoclassical bath equations. The resulting equations are then solved for a number of model problems, ranging from tunneling to vibrational relaxation dynamics. Comparison of the results with those of exact, multiconfiguration time-dependent Hartree calculations in systems with up to 80 bath oscillators shows that the proposed method can be very accurate and might be of help in studying realistic problems with very large baths. To this end, its linear scaling behavior with respect to the number of bath degrees of freedom is shown in practice with model calculations using tens of thousands of bath oscillators.
A multi-reference study of the byproduct formation for a ring-closed dithienylethene photoswitch
(2015)
Photodriven molecular switches are sometimes hindered in their performance by forming byproducts which act as dead ends in sequences of switching cycles, leading to rapid fatigue effects. Understanding the reaction pathways to unwanted byproducts is a prerequisite for preventing them. This article presents a study of the photochemical reaction pathways for byproduct formation in the photochromic switch 1,2-bis-(3-thienyl)-ethene. Specifically, using single-and multi-reference methods the post-deexcitation reaction towards the byproduct in the electronic ground state S-0 when starting from the S-1-S-0 conical intersection (CoIn), is considered in detail. We find an unusual low-energy pathway, which offers the possibility for the formation of a dyotropic byproduct. Several high-energy pathways can be excluded with high probability.
A multi-reference study of the byproduct formation for a ring-closed dithienylethene photoswitch
(2015)
Photodriven molecular switches are sometimes hindered in their performance by forming byproducts which act as dead ends in sequences of switching cycles, leading to rapid fatigue effects. Understanding the reaction pathways to unwanted byproducts is a prerequisite for preventing them. This article presents a study of the photochemical reaction pathways for byproduct formation in the photochromic switch 1,2-bis-(3-thienyl)-ethene. Specifically, using single- and multi-reference methods the post-deexcitation reaction towards the byproduct in the electronic ground state S0 when starting from the S1–S0 conical intersection (CoIn), is considered in detail. We find an unusual low-energy pathway, which offers the possibility for the formation of a dyotropic byproduct. Several high-energy pathways can be excluded with high probability.
A multi-reference study of the byproduct formation for a ring-closed dithienylethene photoswitch
(2015)
Photodriven molecular switches are sometimes hindered in their performance by forming byproducts which act as dead ends in sequences of switching cycles, leading to rapid fatigue effects. Understanding the reaction pathways to unwanted byproducts is a prerequisite for preventing them. This article presents a study of the photochemical reaction pathways for byproduct formation in the photochromic switch 1,2-bis-(3-thienyl)-ethene. Specifically, using single- and multi-reference methods the post-deexcitation reaction towards the byproduct in the electronic ground state S0 when starting from the S1–S0 conical intersection (CoIn), is considered in detail. We find an unusual low-energy pathway, which offers the possibility for the formation of a dyotropic byproduct. Several high-energy pathways can be excluded with high probability.
We present a rigorous method to set up a system-bath Hamiltonian for the coupling of adsorbate vibrations (the system) to surface phonons (the bath). The Hamiltonian is straightforward to derive and exact up to second order in the environment coordinates, thus capable of treating one- and two-phonon contributions to vibration-phonon coupling. The construction of the Hamiltonian uses orthogonal coordinates for system and bath modes, is based on an embedded cluster approach, and generalizes previous Hamiltonians of a similar type, but avoids several (additional) approximations. While the parametrization of the full Hamiltonian is in principle feasible by a first principles quantum mechanical treatment, here we adopt in the spirit of a QM/MM model a combination of density functional theory (“QM”, for the system) and a semiempirical forcefield (“MM”, for the bath). We apply the Hamiltonian to a fully H-covered Si(100)-(2 × 1) surface, using Fermi’s Golden Rule to obtain vibrational relaxation rates of various H–Si bending modes of this system. As in earlier work it is found that the relaxation is dominated by two-phonon contributions because of an energy gap between the Si–H bending modes and the Si phonon bands. We obtain vibrational lifetimes (of the first excited state) on the order of 2 ps at K. The lifetimes depend only little on the type of bending mode (symmetric vs. antisymmetric, parallel vs. perpendicular to the Si2H2 dimers). They decrease by a factor of about two when heating the surface to 300 K. We also study isotope effects by replacing adsorbed H atoms by deuterium, D. The Si–D bending modes are shifted into the Si phonon band of the solid, opening up one-phonon decay channels and reducing the lifetimes to few hundred fs.
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.
We present a new global ground state potential energy surface (PES) for carbon monoxide at a coverage of 1/4, on a rigid Ru(0001) surface [Ru(0001)(2x2):CO]. All six adsorbate degrees of freedom are considered. For constructing the PES, we make use of more than 90 000 points calculated with periodic density functional theory using the RPBE exchange-correlation functional and an empirical van der Waals correction. These points are used for interpolation, utilizing a symmetry-adapted corrugation reducing procedure (CRP). Three different interpolation schemes with increasing accuracy have been realized, giving rise to three flavours of the CRP PES. The CRP PES yields in agreement with the DFT reference and experiments, the atop position of CO to be the most stable adsorption geometry, for the most accurate interpolation with an adsorption energy of 1.69 eV. The CRP PES shows that diffusion parallel to the surface is hindered by a barrier of 430 meV, and that dissociation is facilitated but still activated. As a first "real" application and further test of the new potential, the six-dimensional vibrational Schrodinger equation is solved variationally to arrive at fully coupled, anharmonic frequencies and vibrational wavefunctions for the vibrating, adsorbed CO molecule. Good agreement with experiment is found also here. Being analytical, the new PES opens an efficient way towards multidimensional dynamics. (C) 2014 AIP Publishing LLC.
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.
Ab initio calculations for XPS chemical shifts of poly(vinyl-trifluoroacetate) using trimer models
(2011)
X-ray photoelectron spectra (XPS) of the polymer poly(vinyl-trifluoroacetate) show C(1s) binding energy shifts which are unusual because they are influenced by atoms which are several bonds away from the probed atom. In this work, the influence of the trifluoroacetate substituent on the 1s ionization potential of the carbon atoms of the polyethylene chain is investigated theoretically using mono-substituted, diad and triad models of trimers representing the polymer. Carbon 1s ionization energies are calculated by the Hartree-Fock theory employing Koopmans' theorem. The influence of the configuration and conformation of the functional groups as well as the degree of substitution are found to be important determinants of XPS spectra. It is further found that the 1s binding energy correlates in a linear fashion, with the total electrostatic potential at the position of the probe atom, and depends not only on nearest neighbor effects. This may have implications for the interpretation of high-resolution XP spectra.
In this work, the adsorption and splitting of the water molecule by light and/or an external potential is investigated in the frame of (photo-) electrochemical cells using a rutile ruthenium dioxide anode. With the help of periodic density functional calculations, the adsorbed structures of H(2)O and some radicals involved in the splitting process (O, OH, OOH) are obtained and compared with the available experimental results. On the basis of these electronic-structure calculations, we use a method to calculate the stability of the reaction intermediates and conclude on the thermodynamical possibility of the water splitting reaction at the surface. We demonstrate that a moderate overpotential of 0.64 V is required for the reaction to take place at the RuO(2)(110) surface.
Graphitic carbon nitride, g-C₃N₄, is a promising organic photo-catalyst for a variety of redox reactions. In order to improve its efficiency in a systematic manner, however, a fundamental understanding of the microscopic interaction between catalyst, reactants and products is crucial. Here we present a systematic study of water adsorption on g-C₃N₄ by means of density functional theory and the density functional based tight-binding method as a prerequisite for understanding photocatalytic water splitting. We then analyze this prototypical redox reaction on the basis of a thermodynamic model providing an estimate of the overpotential for both water oxidation and H⁺ reduction. While the latter is found to occur readily upon irradiation with visible light, we derive a prohibitive overpotential of 1.56 eV for the water oxidation half reaction, comparing well with the experimental finding that in contrast to H₂ production O₂ evolution is only possible in the presence of oxidation cocatalysts.
Graphitic carbon nitride, g-C₃N₄, is a promising organic photo-catalyst for a variety of redox reactions. In order to improve its efficiency in a systematic manner, however, a fundamental understanding of the microscopic interaction between catalyst, reactants and products is crucial. Here we present a systematic study of water adsorption on g-C₃N₄ by means of density functional theory and the density functional based tight-binding method as a prerequisite for understanding photocatalytic water splitting. We then analyze this prototypical redox reaction on the basis of a thermodynamic model providing an estimate of the overpotential for both water oxidation and H⁺ reduction. While the latter is found to occur readily upon irradiation with visible light, we derive a prohibitive overpotential of 1.56 eV for the water oxidation half reaction, comparing well with the experimental finding that in contrast to H₂ production O₂ evolution is only possible in the presence of oxidation cocatalysts.
Graphitic carbon nitride, g-C3N4, is a promising organic photo-catalyst for a variety of redox reactions. In order to improve its efficiency in a systematic manner, however, a fundamental understanding of the microscopic interaction between catalyst, reactants and products is crucial. Here we present a systematic study of water adsorption on g-C3N4 by means of density functional theory and the density functional based tight-binding method as a prerequisite for understanding photocatalytic water splitting. We then analyze this prototypical redox reaction on the basis of a thermodynamic model providing an estimate of the overpotential for both water oxidation and H+ reduction. While the latter is found to occur readily upon irradiation with visible light, we derive a prohibitive overpotential of 1.56 eV for the water oxidation half reaction, comparing well with the experimental finding that in contrast to H-2 production O-2 evolution is only possible in the presence of oxidation cocatalysts.
We explore the photophysics of P(NDI2OD-T2), a high-mobility and air-stable n-type donor/acceptor polymer. Detailed steady-state UV-vis and photoluminescence (PL) measurements on solutions of P(NDI2OD-T2) reveal distinct signatures of aggregation. By performing quantum chemical calculations, we can assign these spectral features to unaggregated and stacked polymer chains. NMR measurements independently confirm the aggregation phenomena of P(NDI2OD-T2) in solution. The detailed analysis of the optical spectra shows that aggregation is a two-step process with different types of aggregates, which we confirm by time-dependent PL measurements. Analytical ultracentrifugation measurements suggest that aggregation takes place within the single polymer chain upon coiling. By transferring these results to thin P(NDI2OD-T2) films, we can conclude that film formation is mainly governed by the chain collapse, leading in general to a high aggregate content of similar to 45%. This process also inhibits the formation of amorphous and disordered P(NDI2OD-T2) films.
Pump-probe near edge X-ray absorption fine structure (PP-NEXAFS) spectra of molecules offer insight into valence-excited states, even if optically dark. In PP-NEXAFS spectroscopy, the molecule is "pumped" by UV or visible light enforcing a valence excitation, followed by an X-ray "probe" exciting core electrons into (now) partially empty valence orbitals. Calculations of PP-NEXAFS have so far been done by costly, correlated wavefunction methods which are not easily applicable to medium-sized or large molecules. Here we propose an efficient, first principles method based on density functional theory in combination with the transition potential and Delta SCF methodology (TP-DFT/Delta SCF) to compute molecular ground state and PP-NEXAFS spectra. We apply the method to n ->pi* pump/O-K-edge NEXAFS probe spectroscopy of thymine (for which both experimental and other theoretical data exist) and to n -> pi* or pi -> pi* pump/N-K-edge NEXAFS probe spectroscopies of trans-and cis-azobenzene. Published by AIP Publishing.
Blue light sensing using flavin (BLUF) protein photoreceptor domains change their hydrogen bond network after photoexcitation. To explore this phenomenon, BLUF domains from R. sphaeroides were simulated using Amber99 molecular dynamics (MD). Five starting configurations were considered, to study different BLUF proteins (AppA/BlrB), Trp conformations (Win/Wout), structure determination (X-ray/NMR), and finally, His protonation states. We found dependencies of the hydrogen bonds on almost all parameters. Our data show an especially strong correlation of the Trp position and hydrogen bonds involving Gln63. The latter is in some contradiction to earlier results (Obanayama et al., Photochem. Photobiol. 2008, 84 10031010). Possible origins and implications are discussed. Our calculations support conjectures that Gln63 is more flexible with Trp104 in Win position. Using snapshots from MD and time-dependent density functional theory, UV/vis spectra for the chromophore were determined, which account for molecular motion of the protein under ambient conditions. In accord with experiment, it is found that the UV/vis spectra of BLUF bound flavin are red-shifted and thermally broadened for all calculated p ? p* transitions, relative to gas phase flavin at T = 0 K. However, differences in the spectra between the various BLUF configurations cannot be resolved with the present approach. (c) 2012 Wiley Periodicals, Inc.
Recently, C K-edge Near Edge X-ray Absorption Fine Structure (NEXAFS) spectra of graphite (HOPG) surfaces have been measured for the pristine material, and for HOPG treated with either bromine or krypton plasmas (Lippitz et al., Surf. Sci., 2013, 611, L1). Changes of the NEXAFS spectra characteristic for physical (krypton) and/or chemical/physical modifications of the surface (bromine) upon plasma treatment were observed. Their molecular origin, however, remained elusive. In this work we study by density functional theory, the effects of selected point and line defects as well as chemical modifications on NEXAFS carbon K-edge spectra of single graphene layers. For Br-treated surfaces, also Br 3d X-ray Photoelectron Spectra (XPS) are simulated by a cluster approach, to identify possible chemical modifications. We observe that some of the defects related to plasma treatment lead to characteristic changes of NEXAFS spectra, similar to those in experiment. Theory provides possible microscopic origins for these changes.
Recently, C K-edge Near Edge X-ray Absorption Fine Structure (NEXAFS) spectra of graphite (HOPG) surfaces have been measured for the pristine material, and for HOPG treated with either bromine or krypton plasmas (Lippitz et al., Surf. Sci., 2013, 611, L1). Changes of the NEXAFS spectra characteristic for physical (krypton) and/or chemical/physical modifications of the surface (bromine) upon plasma treatment were observed. Their molecular origin, however, remained elusive. In this work we study by density functional theory, the effects of selected point and line defects as well as chemical modifications on NEXAFS carbon K-edge spectra of single graphene layers. For Br-treated surfaces, also Br 3d X-ray Photoelectron Spectra (XPS) are simulated by a cluster approach, to identify possible chemical modifications. We observe that some of the defects related to plasma treatment lead to characteristic changes of NEXAFS spectra, similar to those in experiment. Theory provides possible microscopic origins for these changes.
Recently, C K-edge Near Edge X-ray Absorption Fine Structure (NEXAFS) spectra of graphite (HOPG) surfaces have been measured for the pristine material, and for HOPG treated with either bromine or krypton plasmas (Lippitz et al., Surf. Sci., 2013, 611, L1). Changes of the NEXAFS spectra characteristic for physical (krypton) and/or chemical/physical modifications of the surface (bromine) upon plasma treatment were observed. Their molecular origin, however, remained elusive. In this work we study by density functional theory, the effects of selected point and line defects as well as chemical modifications on NEXAFS carbon K-edge spectra of single graphene layers. For Br-treated surfaces, also Br 3d X-ray Photoelectron Spectra (XPS) are simulated by a cluster approach, to identify possible chemical modifications. We observe that some of the defects related to plasma treatment lead to characteristic changes of NEXAFS spectra, similar to those in experiment. Theory provides possible microscopic origins for these changes.
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.
We study theoretically the quantum dynamics and spectroscopy of rovibrational polaritons formed in a model system composed of a single rovibrating diatomic molecule, which interacts with two degenerate, orthogonally polarized modes of an optical Fabry-Perot cavity. We employ an effective rovibrational Pauli-Fierz Hamiltonian in length gauge representation and identify three-state vibro-polaritonic conical intersections (VPCIs) between singly excited vibro-polaritonic states in a two-dimensional angular coordinate branching space. The lower and upper vibrational polaritons are of mixed light-matter hybrid character, whereas the intermediate state is purely photonic in nature. The VPCIs provide effective population transfer channels between singly excited vibrational polaritons, which manifest in rich interference patterns in rotational densities. Spectroscopically, three bright singly excited states are identified when an external infrared laser field couples to both a molecular and a cavity mode. The non-trivial VPCI topology manifests as pronounced multi-peak progression in the spectral region of the upper vibrational polariton, which is traced back to the emergence of rovibro-polaritonic light-matter hybrid states. Experimentally, ubiquitous spontaneous emission from cavity modes induces a dissipative reduction of intensity and peak broadening, which mainly influences the purely photonic intermediate state peak as well as the rovibro-polaritonic progression. Published under an exclusive license by AIP Publishing.
The interaction of water with a-alumina (i.e. alpha-Al2O3) surfaces is important in a variety of applications and a useful model for the interaction of water with environmentally abundant aluminosilicate phases. Despite its significance, studies of water interaction with alpha-Al2O3 surfaces other than the (0001) are extremely limited. Here we characterize the interaction of water (D2O) with a well defined alpha-Al2O3(1 (1) over bar 02) surface in UHV both experimentally, using temperature programmed desorption and surface-specific vibrational spectroscopy, and theoretically, using periodic-slab density functional theory calculations. This combined approach makes it possible to demonstrate that water adsorption occurs only at a single well defined surface site (the so-called 1-4 configuration) and that at this site the barrier between the molecularly and dissociatively adsorbed forms is very low: 0.06 eV. A subset of OD stretch vibrations are parallel to this dissociation coordinate, and thus would be expected to be shifted to low frequencies relative to an uncoupled harmonic oscillator. To quantify this effect we solve the vibrational Schrodinger equation along the dissociation coordinate and find fundamental frequencies red-shifted by more than 1500 cm(-1). Within the context of this model, at moderate temperatures, we further find that some fraction of surface deuterons are likely delocalized: dissociatively and molecularly absorbed states are no longer distinguishable.
The interaction of water with α-alumina (i.e. α-Al2O3) surfaces is important in a variety of applications and a useful model for the interaction of water with environmentally abundant aluminosilicate phases. Despite its significance, studies of water interaction with α-Al2O3 surfaces other than the (0001) are extremely limited. Here we characterize the interaction of water (D2O) with a well defined α-Al2O3(1[1 with combining macron]02) surface in UHV both experimentally, using temperature programmed desorption and surface-specific vibrational spectroscopy, and theoretically, using periodic-slab density functional theory calculations. This combined approach makes it possible to demonstrate that water adsorption occurs only at a single well defined surface site (the so-called 1–4 configuration) and that at this site the barrier between the molecularly and dissociatively adsorbed forms is very low: 0.06 eV. A subset of OD stretch vibrations are parallel to this dissociation coordinate, and thus would be expected to be shifted to low frequencies relative to an uncoupled harmonic oscillator. To quantify this effect we solve the vibrational Schrödinger equation along the dissociation coordinate and find fundamental frequencies red-shifted by more than 1500 cm−1. Within the context of this model, at moderate temperatures, we further find that some fraction of surface deuterons are likely delocalized: dissociatively and molecularly absorbed states are no longer distinguishable.
Azobenzene-based molecular photoswitches are becoming increasingly important for the development of photoresponsive, functional soft-matter material systems. Upon illumination with light, fast interconversion between a more stable trans and a metastable cis configuration can be established resulting in pronounced changes in conformation, dipole moment or hydrophobicity. A rational design of functional photosensitive molecules with embedded azo moieties requires a thorough understanding of isomerization mechanisms and rates, especially the thermally activated relaxation. For small azo derivatives considered in the gas phase or simple solvents, Eyring’s classical transition state theory (TST) approach yields useful predictions for trends in activation energies or corresponding half-life times of the cis isomer. However, TST or improved theories cannot easily be applied when the azo moiety is part of a larger molecular complex or embedded into a heterogeneous environment, where a multitude of possible reaction pathways may exist. In these cases, only the sampling of an ensemble of dynamic reactive trajectories (transition path sampling, TPS) with explicit models of the environment may reveal the nature of the processes involved. In the present work we show how a TPS approach can conveniently be implemented for the phenomenon of relaxation–isomerization of azobenzenes starting with the simple examples of pure azobenzene and a push–pull derivative immersed in a polar (DMSO) and apolar (toluene) solvent. The latter are represented explicitly at a molecular mechanical (MM) and the azo moiety at a quantum mechanical (QM) level. We demonstrate for the push–pull azobenzene that path sampling in combination with the chosen QM/MM scheme produces the expected change in isomerization pathway from inversion to rotation in going from a low to a high permittivity (explicit) solvent model. We discuss the potential of the simulation procedure presented for comparative calculation of reaction rates and an improved understanding of activated states.
The cis-trans isomerisation of N-benzylideneaniline (NBA) and derivatives containing a central C=N bond has been investigated experimentally and theoretically. Eight different NBA molecules in three different solvents were irradiated to enforce a photochemical trans (hv) -> cis isomerisation and the kinetics of the thermal backreaction cis (Delta)-> trans were determined by NMR spectroscopy measurements in the temperature range between 193 and 288 K. Theoretical calculations using density functional theory and Eyring transition-state theory were carried out for 12 different NBA species in the gas phase and three different solvents to compute thermal isomerisation rates of the thermal back reaction. While the computed absolute rates are too large, they reveal and explain experimental trends. Time-dependent density functional theory provides optical spectra for vertical transitions and excitation energy differences between trans and cis forms. Together with isomerisation rates, the latter can be used to identify "optimal switches" with good photochromicity and reasonable thermal stability.
In this work, we study the adsorption of CO from low to high coverage at a defect-free NaCl(100) surface by means of duster and periodic models, using highly accurate wave function-based QM:QM embedding as well as density functional theory. At low coverages, the most accurate methods predict a zero-point-corrected adsorption energy of around 13 kJ/mol, and the CO molecules are found to be oriented perpendicular to the surface. At higher coverages, lower-energy phases with nonparallel/upright, tilted orientations emerge. Besides the well-known p(2 x 1)/antiparallel phase (T/A), we find other tilted phases (tilted/irregular, T/I; tilted/spiral, T/S) as local minima. Vibrational frequencies for CO adsorbed on NaCl(100) and Davydov splittings of the C-O stretch vibration are also determined. The IR spectra are characteristic fingerprints for the relative orientation of CO molecules and may therefore be used as sensitive probes to distinguish parallel/upright from various tilted adsorption phases.
Comparing thermal wave function methods for multi-configuration time-dependent Hartree simulations
(2014)
We compare two methods for creating stochastic temperature wave functions that can be used for Multi-Configuration Time-Dependent Hartree (MCTDH) simulations. In the first method, the MCTDH coefficients are chosen randomly, while the other method uses a single Hartree product of random single-particle functions (SPFs). We find that using random SPFs dramatically improves convergence for a model system for surface sticking.
Here we apply and expand the knowledge developed in the case of the H atom to describe high-harmonic generation (HHG) for the H-2 molecule by using time-dependent configuration interaction with single excitations. The implications of using a finite atomic orbital basis set and the impact of a heuristic lifetime model which addresses ionisation losses are discussed. We also examine the influence of the angular momentum of the basis on the computed HHG spectra. Moreover, we discuss the impact of adding diffuse functions and ghost atoms in different geometrical configurations around the molecule. The effects of these additional centres on the HHG spectra are correlated with the physical interpretation of this nonlinear optical phenomenon as given by the three-step model, relating the maximal radial extent of the electron as predicted by the model to the radial extent of the Gaussian basis sets. [GRAPHICS] .
Laser-induced condensed phase reactions are often interpreted as nonequilibrium phenomena that go beyond conventional thermodynamics. Here, we show by Langevin dynamics and for the example of femtosecond-laser desorption of hydrogen from a ruthenium surface that light adsorbates thermalize rapidly due to ultrafast energy redistribution after laser excitation. Despite the complex reaction mechanism involving hot electrons in the surface region, all desorption product properties are characterized by equilibrium distributions associated with a single, unique temperature. This represents an example of ultrahot chemistry on the subpicosecond time scale.
Quantum mechanical tunnelling describes transmission of matter waves through a barrier with height larger than the energy of the wave(1). Tunnelling becomes important when the de Broglie wavelength of the particle exceeds the barrier thickness; because wavelength increases with decreasing mass, lighter particles tunnel more efficiently than heavier ones. However, there exist examples in condensed-phase chemistry where increasing mass leads to increased tunnelling rates(2). In contrast to the textbook approach, which considers transitions between continuum states, condensed-phase reactions involve transitions between bound states of reactants and products. Here this conceptual distinction is highlighted by experimental measurements of isotopologue-specific tunnelling rates for CO rotational isomerization at an NaCl surface(3,4), showing nonmonotonic mass dependence. A quantum rate theory of isomerization is developed wherein transitions between sub-barrier reactant and product states occur through interaction with the environment. Tunnelling is fastest for specific pairs of states (gateways), the quantum mechanical details of which lead to enhanced cross-barrier coupling; the energies of these gateways arise nonsystematically, giving an erratic mass dependence. Gateways also accelerate ground-state isomerization, acting as leaky holes through the reaction barrier. This simple model provides a way to account for tunnelling in condensed-phase chemistry, and indicates that heavy-atom tunnelling may be more important than typically assumed.
When new covalent organic frameworks (COFs) are designed, the main efforts are typically focused on selecting specific building blocks with certain geometries and properties to control the structure and function of the final COFs. The nature of the linkage (imine, boroxine, vinyl, etc.) between these building blocks naturally also defines their properties. However, besides the linkage type, the orientation, i.e., the constitutional isomerism of these linkages, has rarely been considered so far as an essential aspect. In this work, three pairs of constitutionally isomeric imine-linked donor-acceptor (D-A) COFs are synthesized, which are different in the orientation of the imine bonds (D-C=N-A (DCNA) and D-N=C-A (DNCA)). The constitutional isomers show substantial differences in their photophysical properties and consequently in their photocatalytic performance. Indeed, all DCNA COFs show enhanced photocatalytic H2 evolution performance than the corresponding DNCA COFs. Besides the imine COFs shown here, it can be concluded that the proposed concept of constitutional isomerism of linkages in COFs is quite universal and should be considered when designing and tuning the properties of COFs.
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 combine the stochastic pulse optimization (SPO) scheme with the time-dependent configuration interaction singles method in order to control the high frequency response of a simple molecular model system to a tailored femtosecond laser pulse. For this purpose, we use H-2 treated in the fixed nuclei approximation. The SPO scheme, as similar genetic algorithms, is especially suited to control highly non-linear processes, which we consider here in the context of high harmonic generation. Here, we will demonstrate that SPO can be used to realize a "non-harmonic" response of H2 to a laser pulse. Specifically, we will show how adding low intensity side frequencies to the dominant carrier frequency of the laser pulse and stochastically optimizing their contribution can create a high-frequency spectral signal of significant intensity, not harmonic to the carrier frequency. At the same time, it is possible to suppress the harmonic signals in the same spectral region, although the carrier frequency is kept dominant during the optimization. (C) 2016 AIP Publishing LLC.
The biconformational switching of single cyclooctadiene molecules chemisorbed on a Si(001) surface was explored by quantum chemical and quantum dynamical calculations and low-temperature scanning tunneling microscopy experiments. The calculations rationalize the experimentally observed switching driven by inelastic electron tunneling (IET) at 5 K. At higher temperatures, they predict a controllable crossover behavior between IET-driven and thermally activated switching, which is fully confirmed by experiment.
The switching of single cyclooctadiene molecules chemisorbed on a Si(100) surface between two stable conformations, can be achieved with a scanning tunneling microscope [Nacci , Phys. Rev. B 77, 121405(R) (2008)]. Recently, it was shown by quantum chemical and quantum dynamical simulations that major experimental facts can be explained by a single-mode model with switching enforced by inelastic electron tunneling (IET) excitations and perturbed by vibrational relaxation [Nacci , Nano Lett. 9, 2997 (2009)]. In the present paper, we extend the previous theoretical work in several respects: (1) The model is generalized to a two-mode description in which two C2H4 units of COD can move independently; (2) contributions of dipole and, in addition, (cation and anion) resonance-IET rates are considered; (3) the harmonic-linear vibrational relaxation model used previously is generalized to anharmonic vibrations. While the present models highlight generic aspects of IET-switching between two potential minima, they also rationalize specific experimental findings for COD/Si(100): (1) A single-electron excitation mechanism with a linear dependence of the switching rate on tunneling current I, (2) the capability to switch both at negative and positive sample biases, and (3) a crossover temperature around similar to 60 K from an IET-driven, T-independent atom tunneling regime, to classical over-the-barrier isomerization with exponential T-dependence at higher temperatures for a bias voltage of +1.5 V and an average tunneling current of 0.73 nA.
Diarylethene derivatives are photochromic molecular switches, undergoing a ring-opening/-closing reaction by illumination with light. The symmetry of the closed form is determined by the WoodWard Hoffinann rules according to which the reaction proceeds by corirotatory rotation -in that case. Here, we show by a cOrnbined approach of scanning tunneling microscopy (STM) and density functional theory (DFT) calculations that the Open isomer of 4,4'-(4,4'-(perfluorocydopent-1-ene-1,2-diyl)bis(5-methyl-thiophent-4,2,4-dipyridine) (PDTE) retains its open form upon adsorption on a Ag(111) surface. It caribe switched into a closed form, which we identify as the digrotatOly cydization product, by controlled manipulation 'With the STM tip, Evidence of an electric-field dependent switching-process 'is interpreted on the basis of a Simple electroStatic Model, which suggests that the reaction proceedS via an "upright" intermediate state. This pathway thus strongly differs from the switching reaction in solution.
High Harmonic Generation (HHG) is a nonlinear optical process that provides a tunable source for high-energy photons and ultrashort laser pulses. Recent experiments demonstrated that HHG spectroscopy may also be used as an analytical tool to discriminate between randomly oriented configurational isomers of polyatomic organic molecules, namely, between the cis- and trans-forms of 1,2-dichloroethene (DCE) [M. C. H. Wong et al., Phys. Rev. A 84, 051403 (2011)]. Here, we suggest as an economic and at the same time a reasonably accurate method to compute HHG spectra for polyatomic species, Time-Dependent Configuration Interaction Singles (TD-CIS) theory in combination with extended atomic orbital bases and different models to account for ionization losses. The HHG spectra are computed for aligned and unaligned cis- and trans-DCE. For the unaligned case, a coherent averaging over possible rotational orientations is introduced. Furthermore, using TD-CIS, possible differences between the HHG spectra of cis- and trans-DCE are studied. For aligned molecules, spectral differences between cis and trans emerge, which can be related to their different point group symmetries. For unaligned, randomly oriented molecules, we also find distinct HHG spectra in partial agreement with experiment. In addition to HHG response in the frequency space, we compute time-frequency HHG spectra to gain insight into which harmonics are emitted at which time. Further differences between the two isomers emerge, suggesting time-frequency HHG as another tool to discriminate configurational isomers.
An electronic friction approach based on Langevin dynamics is used to describe the multidimensional (six-dimensional) dynamics of femtosecond laser induced desorption of H-2 and D-2 from a H(D)-covered Ru(0001) surface. The paper extends previous reduced-dimensional models, using a similar approach. In the present treatment forces and frictional coefficients are calculated from periodic density functional theory (DFT) and essentially parameter-free, while the action of femtosecond laser pulses on the metal surface is treated by using the two-temperature model. Our calculations shed light on the performance and validity of various adiabatic, non-adiabatic, and Arrhenius/Kramers type kinetic models to describe hot-electron mediated photoreactions at metal surfaces. The multidimensional frictional dynamics are able to reproduce and explain known experimental facts, such as strong isotope effects, scaling of properties with laser fluence, and non-equipartitioning of vibrational, rotational, and translational energies of desorbing species. Further, detailed predictions regarding translations are made, and the question for the controllability of photoreactions at surfaces with the help of vibrational preexcitation is addressed.
In this paper, we perform many-electron dynamics using the time-dependent configuration-interaction method in its reduced density matrix formulation (rho-TDCI). Dissipation is treated implicitly using the Lindblad formalism. To include the effect of ionization on the state-resolved dynamics, we extend a recently introduced heuristic model for ionizing states to the rho-TDCI method, which leads to a reduced density matrix evolution that is not norm-preserving. We apply the new method to the laser-driven excitation of H(2) in a strongly dissipative environment, for which the state-resolve lifetimes are tuned to a few femtoseconds, typical for dynamics of adsorbate at metallic surfaces. Further testing is made on the laser-induced intramolecular charge transfer in a quinone derivative as a model for a molecular switch. A modified scheme to treat ionizing states is proposed to reduce the computational burden associated with the density matrix propagation, and it is thoroughly tested and compared to the results obtained with the former model. The new approach scales favorably (similar to N(2)) with the number of configurations N used to represent the reduced density matrix in the rho-TDCI method, as compared to a N(3) scaling for the model in its original form.
While azobenzenes readily photoswitch in solution, their photoisomerization in densely packed self-assembled monolayers (SAMs) can be suppressed. Reasons for this can be steric hindrance and/or electronic quenching, e.g., by exciton coupling. We address these possibilities by means of nonadiabatic molecular dynamics with trajectory surface hopping calculations, investigating the trans -> cis isomerization of azobenzene after excitation into the pi pi* absorption band. We consider a free monomer, an isolated dimer and a dimer embedded in a SAM-like environment of additional azobenzene molecules, imitating in this way the gradual transition from an unconstrained over an electronically coupled to an electronically coupled and sterically hindered, molecular switch. Our simulations reveal that in comparison to the single molecule the quantum yield of the trans -> cis photoisomerization is similar for the isolated dimer, but greatly reduced in the sterically constrained situation. Other implications of dimerization and steric constraints are also discussed.
A variety of azobenzenes were synthesized to study the behavior of their E and Z isomers upon electrochemical reduction. Our results show that the radical anion of the Z isomer is able to rapidly isomerize to the corresponding E configured counterpart with a dramatically enhanced rate as compared to the neutral species. Due to a subsequent electron transfer from the formed E radical anion to the neutral Z starting material the overall transformation is catalytic in electrons; i.e., a substoichiometric amount of reduced species can isomerize the entire mixture. This pathway greatly increases the efficiency of (photo)switching while also allowing one to reach photostationary state compositions that are not restricted to the spectral separation of the individual azobenzene isomers and their quantum yields. In addition, activating this radical isomerization pathway with photoelectron transfer agents allows us to override the intrinsic properties of an azobenzene species by triggering the reverse isomerization direction (Z -> E) by the same wavelength of light, which normally triggers E -> Z isomerization. The behavior we report appears to be general, implying that the metastable isomer of a photoswitch can be isomerized to the more stable one catalytically upon reduction, permitting the optimization of azobenzene switching in new as well as indirect ways.