Refine
Year of publication
Document Type
- Article (125)
- Postprint (8)
- Part of a Book (1)
- Other (1)
Language
- English (135) (remove)
Is part of the Bibliography
- yes (135)
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
- Institut für Chemie (135) (remove)
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] .
Fluorination of the hydroxylated alpha-Al2O3 (0001) surface is studied using periodic density functional theory calculations. On the basis of a hypothetical reaction substituting surface hydroxyl groups with fluorine atoms, we find surface fluorination to be strongly exergonic but kinetically hindered. Fluorinated surface areas turn out to be rather hydrophobic as compared to hydroxylated areas, suggesting fluorination as a potential route for tuning oxide surface properties such as hydrophilicity.
Molecular photoswitches such as azobenzenes, which undergo photochemical trans <-> cis isomerizations, are often mounted for possible applications on a surface and/or surrounded by other switches, for example, in self-assembled monolayers. This may suppress the isomerization cross section due to possible steric reasons, or, as recently speculated, by exciton coupling to. neighboring switches, leading to ultrafast electronic quenching (Gahl et al., J. Am. Chem. Soc. 2010, 132, 1831). The presence of exciton coupling has been anticipated from a blue shift of the optical absorption band, compared to molecules in solution. From the theory side the need arises to properly analyze and quantify the change of absorption spectra of interacting and adsorbed switches. In particular, suitable methods should be identified, and effects of intermolecule and molecule surface interactions on spectra should be disentangled. In this paper by means of time-dependent Hartree-Fock. (TD-HF), various flavors of time-dependent density functional theory (TD-DFT), and the correlated wave function based, coupled cluster (CC2) method we investigated the 4-nitroazobenzene molecule as an:example: The low-lying singlet excited states in the isolated trans monomer and dieter as well as their composites with a silicon pentamantane nanocluster, which serves also as a crude model for a silicon surface, were determined. As most important results we found that (i) HF, CC2, range-separated density functionals, or global hybrids with large amount of exact exchange are able to describe exciton (Davydov) splitting properly, while hybrids with small amount of exact exchange fail producing spurious charge transfer. (ii) The exciton splitting in a free dimer would lead to a blue shift of the absorption signal; however, this effect is almost nullified or even overcompensated by the shift arising from van der Waals interactions between the two molecules. (iii) Adsorption on the Si "surface" leads to a further, strong red shift for the present system. (iv) At a next-nearest neighbor distance (of similar to 3.6 angstrom), the exciton splitting is similar to 0.3 eV, with or without "surface", suggesting a rapid quenching of the molecular pi ->pi* excitation. At larger distances, exciton splitting decreases rapidly.
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.
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.
In recent simulations of femtosecond laser induced desorption of molecular oxygen from the Ag(110) surface, it has been shown that depending on the properties (depth and electronic environment) of the well in which 02 is adsorbed, the desorption can be either induced dominantly by hot electrons or via excitations of phonons. In this work we explore whether the ratios between the desorption yields from different adsorption wells can be tuned by changing initial surface temperature and laser pulse properties. We show that the initial surface temperature is an important parameter, and that by using low initial surface temperatures the electronically mediated process can be favored. In contrast, laser properties seem to have only a modest influence on the results. (C) 2016 Elsevier B.V. All rights reserved.
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.
Action spectroscopy has emerged as an analytical tool to probe excited states in the gas phase. Although comparison of gas-phase absorption properties with quantum-chemical calculations is, in principle, straightforward, popular methods often fail to describe many molecules of interest-such as xanthene analogues. We, therefore, face their nano-and picosecond laser-induced photofragmentation with excited-state computations by using the CC2 method and time-dependent density functional theory (TDDFT). Whereas the extracted absorption maxima agree with CC2 predictions, the TDDFT excitation energies are blueshifted. Lowering the amount of Hartree-Fock exchange in the DFT functional can reduce this shift but at the cost of changing the nature of the excited state. Additional bandwidth observed in the photofragmentation spectra is rationalized in terms of multiphoton processes. Observed fragmentation from higher-lying excited states conforms to intense excited-to-excited state transitions calculated with CC2. The CC2 method is thus suitable for the comparison with photofragmentation in xanthene analogues.
A Langevin model accounting for all six molecular degrees of freedom is applied to femtosecond-laser induced, hot-electron driven dynamics of Ru(0001)(2 x 2): CO. In our molecular dynamics with electronic friction approach, a recently developed potential energy surface based on gradient-corrected density functional theory accounting for van der Waals interactions is adopted. Electronic friction due to the coupling of molecular degrees of freedom to electron-hole pairs in the metal are included via a local density friction approximation, and surface phonons by a generalized Langevin oscillator model. The action of ultrashort laser pulses enters through a substrate-mediated, hot-electron mechanism via a time-dependent electronic temperature (derived from a two-temperature model), causing random forces acting on the molecule. The model is applied to laser induced lateral diffusion of CO on the surface, "hot adsorbate" formation, and laser induced desorption. Reaction probabilities are strongly enhanced compared to purely thermal processes, both for diffusion and desorption. Reaction yields depend in a characteristic (nonlinear) fashion on the applied laser fluence, as well as branching ratios for various reaction channels. Computed two-pulse correlation traces for desorption and other indicators suggest that aside from electron-hole pairs, phonons play a non-negligible role for laser induced dynamics in this system, acting on a surprisingly short time scale. Our simulations on precomputed potentials allow for good statistics and the treatment of long-time dynamics (300 ps), giving insight into this system which hitherto has not been reached. We find generally good agreement with experimental data where available and make predictions in addition. A recently proposed laser induced population of physisorbed precursor states could not be observed with the present low-coverage model.
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.
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 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.
Nonlinear optical response of photochromic azobenzene-functionalized self-assembled monolayers
(2015)
The combination of photochromic and nonlinear optical (NLO) properties of azobenzene-functionalized self-assembled monolayers (SAMs) constitutes an intriguing step towards novel photonic and optoelectronic devices. By utilizing the second-order NLO process of second harmonic generation (SHG), supported by density-functional theory and correlated wave function method calculations, we demonstrate that the photochromic interface provides the necessary prerequisites en route towards possible future technical applications: we find a high NLO contrast on the order of 16% between the switching states. These are furthermore accessible reversibly and with high efficiencies in terms of cross sections on the order of 10(-18) cm(2) for both photoisomerization reactions, i.e., drivable by means of low-power LED light sources. Finally, both photostationary states (PSSs) are thermally stable at ambient conditions.
Optical properties of modified diamondoids have been studied theoretically using vibrationally resolved electronic absorption, emission and resonance Raman spectra. A time-dependent correlation function approach has been used for electronic two-state models, comprising a ground state (g) and a bright, excited state (e), the latter determined from linear-response, time-dependent density functional theory (TD-DFT). The harmonic and Condon approximations were adopted. In most cases origin shifts, frequency alteration and Duschinsky rotation in excited states were considered. For other cases where no excited state geometry optimization and normal mode analysis were possible or desired, a short-time approximation was used. The optical properties and spectra have been computed for (i) a set of recently synthesized sp(2)/sp(3) hybrid species with CQC double-bond connected saturated diamondoid subunits, (ii) functionalized (mostly by thiol or thione groups) diamondoids and (iii) urotropine and other C-substituted diamondoids. The ultimate goal is to tailor optical and electronic features of diamondoids by electronic blending, functionalization and substitution, based on a molecular-level understanding of the ongoing photophysics.
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.
Incorporating photochromic molecules into organic/inorganic hybrid materials may lead to photoresponsive systems. In such systems, the second-order nonlinear properties can be controlled via external stimulation with light at an appropriate wavelength. By creating photochromic molecular switches containing self-assembled monolayers on Si(111), we can demonstrate efficient reversible switching, which is accompanied by a pronounced modulation of the nonlinear optical (NLO) response of the system. The concept of utilizing functionalized photoswitchable Si surfaces could be a way for the generation of two-dimensional NLO switching materials, which are promising for applications in photonic and optoelectronic devices.
Oxide/water interfaces are ubiquitous in a wide variety of applications and the environment. Despite this ubiquity, and attendant decades of study, gaining molecular level insight into water/oxide interaction has proven challenging. In part, this challenge springs from a lack of tools to concurrently characterize changes in surface structure (i.e., water/oxide interaction from the perspective of the solid) and O-H population and local environment (i.e., water/oxide interaction from the water perspective). Here, we demonstrate the application of surface specific vibrational spectroscopy to the characterization of the interaction of the paradigmatic alpha-Al2O3(0001) surface and water. By probing both the interfacial Al-O (surface phonon) and O-H spectral response, we characterize this interaction from both perspectives. Through electronic structure calculation, we assign the interfacial Al-O response and rationalize its changes on surface dehydroxylation and reconstruction. Because our technique is all-optical and interface specific, it is equally applicable to oxide surfaces in vacuum, ambient atmospheres and at the solid/liquid interface. Application of this approach to additional alumina surfaces and other oxides thus seems likely to significantly expand our understanding of how water meets oxide surfaces and thus the wide variety of phenomena this interaction controls. (C) 2015 AIP Publishing LLC.
System-bath problems in physics and chemistry are often described by Markovian master equations. However, the Markov approximation, i.e., neglect of bath memory effects is not always justified, and different measures of non-Markovianity have been suggested in the literature to judge the validity of this approximation. Here we calculate several computable measures of non-Markovianity for the non-trivial problem of a harmonic oscillator coupled to a large number of bath oscillators. The Multi Configurational Time Dependent Hart ree nietliod is used to provide a numerically converged solution of the system-bath Schrodinger equation, from which the appropriate quantities can be calculated. In particular, we consider measures based on trace-distances and quantum discord for a variety of initial states. These quantities have proven useful in the case of two-level and other small model systems Tpically encountered in quantum optics; but are less straightforward to interpret for the more complex model systems that are relevant for chemical physics.
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.
Femtosecond-laser pulse driven non-adiabatic spectroscopy and dynamics in molecular and condensed phase systems continue to be a challenge for theoretical modelling. One of the main obstacles is the "curse of dimensionality" encountered in non-adiabatic, exact wavepacket propagation. A possible route towards treating complex molecular systems is via semiclassical surface-hopping schemes, in particular if they account not only for non-adiabatic post-excitation dynamics but also for the initial optical excitation. One such approach, based on initial condition filtering, will be put forward in what follows. As a simple test case which can be compared with exact wavepacket dynamics, we investigate the influence of the different parameters determining the shape of a laser pulse (e.g., its finite width and a possible chirp) on the predissociation dynamics of a NaI molecule, upon photoexcitation of the A(0(+)) state. The finite-pulse effects are mapped into the initial conditions for semiclassical surface-hopping simulations. The simulated surface-hopping diabatic populations are in qualitative agreement with the quantum mechanical results, especially concerning the subpicosend photoinduced dynamics, the main deviations being the relative delay of the non-adiabatic transitions in the semiclassical picture. Likewise, these differences in the time-dependent electronic populations calculated via the semiclassical and the quantum methods are found to have a mild influence on the overall probability density distribution. As a result, the branching ratios between the bound and the dissociative reaction channels and the time-evolution of the molecular wavepacket predicted by the semiclassical method agree with those computed using quantum wavepacket propagation. Implications for more challenging molecular systems are given. (C) 2015 AIP Publishing LLC.
The Photoinduced E -> Z Isomerization of Bisazobenzenes: A Surface Hopping Molecular Dynamics Study
(2015)
The photoinduced E -> Z isomerization of azobenzene is a prototypical example of molecular switching. On the way toward rigid molecular rods such as those for opto-mechanical applications, multiazobenzene structures have been suggested in which several switching units are linked together within the same molecule (Bleger et al., J. Phys. Chem. B 2011, 115, 9930-9940). Large differences in the switching efficiency of multiazobenzenes have been observed, depending on whether the switching units are electronically decoupled or not. In this paper we study, on a time-resolved molecular level, the E -> Z isomerization of the simplest multiazobenzene, bisazobenzene (BAB). Two isomers (ortho- and para-BAB), differing only in the connectivity of two azo groups on a shared phenyl ring will be considered.To do so, nonadiabatic semiclassical dynamics after photo-excitation of the isomers are studied by employing an "on-the-fly", fewest switches surface hopping approach. States and couplings are calculated by Configuration Interaction (CI) based on a semiempirical (AM1) Hamiltonian (Persico and co-workers, Chem. Eur. J. 2004, 10, 2327-2341). In the case of para-BAB, computed quantum yields for photoswitching are drastically reduced compared to pristine azobenzene, due to electronic coupling of both switching units. A reason for this (apart from altered absorption spectra and reduced photochromicity) is the drastically reduced lifetimes of electronically excited states which are transiently populated. In contrast for meta-connected species, electronic subsystems are largely decoupled, and computed quantum yields are slightly higher than that for pristine azobenzene because of new isomerization channels. In this case we can also distinguish between single- and double-switch events and we find a cooperative effect: The isomerization of a single azo group is facilitated if the other azo group is already in the Z-configuration.
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.
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.
alpha-Al2O3 surfaces are common in both engineered applications and the environment. Much prior work indicates that their properties, e.g., reactivity, polarity, and charge, change dramatically on interaction with water. Perhaps the simplest question that can be asked of alpha-Al2O3/water interaction is how a single water molecule interacts with the most stable alpha-Al2O3 surface: the alpha-Al2O3(0001). Over the last 15 years, a series of theoretical studies have found that water dissociatively adsorbs on alpha-Al2O3(0001) through two channels. However, to our knowledge no experimental evidence of these dissociation pathways has appeared. By combining sample preparation via supersonic molecular beam dosing, sample characterization via coherent, surface specific vibrational spectroscopy and electronic structure theory, we report the first experimental observation of reaction products of each, theoretically predicted, dissociation channel. These results thus overcome a 15 year old experiment/theory disconnect and make possible a variety of intriguing experiments that promise to provide significant new insights into water/Al2O3 and water/oxide interaction more generally.
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.
Improving the photochemical properties of molecular photoswitches is crucial for the development of light-responsive systems in materials and life sciences. ortho-Fluoroazobenzenes are a new class of rationally designed photochromic azo compounds with optimized properties, such as the ability to isomerize with visible light only, high photoconversions, and unprecedented robust bistable character. Introducing sigma-electron-withdrawing F atoms ortho to the N=N unit leads to both an effective separation of the n -> pi* bands of the E and Z isomers, thus offering the possibility of using these two transitions for selectively inducing E/Z iso-merizations, and greatly enhanced thermal stability of the Z isomers. Additional para-electron-withdrawing groups (EWGs) work in concert with ortho-F atoms, giving rise to enhanced separation of the n -> pi* transitions. A comprehensive study of the effect of substitution on the key photochemical properties of ortho-fluoroazobenzenes is reported herein. In particular, the position, number, and nature of the EWGs have been varied, and the visible light photoconversions, quantum yields of isomerization, and thermal stabilities have been measured and rationalized by DFT calculations.
Using density functional theory and Ab Initio Molecular Dynamics with Electronic Friction (AIMDEF), we study the adsorption and dissipative vibrational dynamics of hydrogen atoms chemisorbed on free-standing lead films of increasing thickness. Lead films are known for their oscillatory behaviour of certain properties with increasing thickness, e.g., energy and electron spill-out change in discontinuous manner, due to quantum size effects [G. Materzanini, P. Saalfrank, and P.J.D. Lindan, Phys. Rev. B 63, 235405 (2001)]. Here, we demonstrate that oscillatory features arise also for hydrogen when chemisorbed on lead films. Besides stationary properties of the adsorbate, we concentrate on finite vibrational lifetimes of H-surface vibrations. As shown by AIMDEF, the damping via vibration-electron hole pair coupling dominates clearly over the vibration-phonon channel, in particular for high-frequency modes. Vibrational relaxation times are a characteristic function of layer thickness due to the oscillating behaviour of the embedding surface electronic density. Implications derived from AIMDEF for frictional many-atom dynamics, and physisorbed species will also be given. (C) 2014 AIP Publishing LLC.
Based on the analysis of optical absorption spectra, it has recently been speculated that the excitonic coupling between individual azobenzene-functionalized alkanethiols arranged in a self-assembled monolayer (SAM) on a gold surface could be strong enough to hinder collective trans-cis isomerization-on top of steric hindrance [Gahl et al., J. Am. Chem. Soc., 2010, 132, 1831]. Using models of SAMs of increasing complexity (dimer, linear N-mers, and two-dimensionally arranged N-mers) and density functional theory on the (TD-) B3LYP/6-31G* level, we determine optical absorption spectra, the nature and magnitude of excitonic couplings, and the corresponding spectral shifts. It is found that at inter-monomer distances of about 20 angstrom and above, TD-B3LYP excitation frequencies (and signal intensities) can be well described by the frequently used point-dipole approximation. Further, calculated blue shifts in optical absorption spectra account for the experimental observations made for azobenzene/gold SAMs, and hint to the fact that they can indeed be responsible for reduced switching probability in densely packed self-assembled structures.
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.
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.
In blue-light photoreceptors using flavin (BLUF), the signaling state is formed already within several 100 ps after illumination, with only small changes of the absorption spectrum. The accompanying structural evolution can, in principle, be monitored by femtosecond stimulated Raman spectroscopy (FSRS). The method is used here to characterize the excited-state properties of riboflavin and flavin adenine dinucleotide in polar solvents. Raman modes are observed in the range 90-1800 cm(-1) for the electronic ground state S-0 and upon excitation to the S-1 state, and modes >1000 cm(-1) of both states are assigned with the help of quantum-chemical calculations. Line shapes are shown to depend sensitively on resonance conditions. They are affected by wavepacket motion in any of the participating electronic states, resulting in complex amplitude modulation of the stimulated Raman spectra. Wavepackets in S-1 can be marked, and thus isolated, by stimulated-emission pumping with the picosecond Raman pulses. Excited-state absorption spectra are obtained from a quantitative comparison of broadband transient fluorescence and absorption. In this way, the resonance conditions for FSRS are determined. Early differences of the emission spectrum depend on excess vibrational energy, and solvation is seen as dynamic Stokes shift of the emission band. The ne state is evidenced only through changes of emission oscillator strength during solvation. S-1 quenching by adenine is seen with all methods in terms of dynamics, not by spectral intermediates.
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.
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.
A strategy to optimize the photoswitching efficiency of rigid, linear multiazobenzene constructs is presented. It consists of introducing large dihedral angles between azobenzene moieties linked via aryl-aryl connections in their para positions. Four bisazobenzenes exhibiting different dihedral angles as well as three single azobenzene reference compounds have been synthesized, and their switching behavior has been studied as well as experimentally and theoretically analyzed. As the dihedral angle between the two azobenzene units increases and consequently the electronic conjugation decreases, the photochromic characteristics improve, finally leading to individual azobenzene switches operating independently in the case of the perpendicular ortho,ortho,ortho',ortho'-tetramethyl biphenyl linker. The electronic decoupling leads to efficient separation of the absorption spectra of the involved switching states and hence by choosing the appropriate irradiation wavelength, an almost quantitative E -> Z photoisomerization up to 97% overall Z-content can be achieved. In addition, thermal Z -> E isomerization processes become independent of each other with increasing decoupling. The electronic decoupling could furthermore be proven by electrochemistry. The experimental data are supported by theory, and calculations additionally provide mechanistic insight into the preferred pathway for the thermal Z,Z -> Z,E -> E,E isomerization via inversion on the inner N-atoms. Our decoupling approach outlined herein provides the basis for constructing rigid rod architectures composed of multiple azobenzene photochromes, which display practically quantitative photoswitching properties, a necessary prerequisite to achieve highly efficient transduction of light energy directly into motion.
Flavins are chromophores in light-gated enzymes and therefore central in many photobiological processes. To unravel the optical excitation process as the initial, elementary step towards signal transduction, detailed ultrafast (femtosecond) experiments probing the photo-activation of flavins have been carried out recently [Weigel et al., J. Phys. Chem. B, 2011, 115, 3656-3680.]. The present paper contributes to a further understanding and interpretation of these experiments by studying the post-excitation vibrational dynamics of riboflavin (RF) and microsolvated riboflavin, RF center dot 4H(2)O, using first principles non-adiabatic molecular dynamics. By analyzing the characteristic atom motions and calculating time-resolved stimulated emission spectra following pi pi* excitation, it is found that after optical excitation C-N and C-C vibrations in the isoalloxazine rings of riboflavin set in. The Franck-Condon (vertically excited) state decays within about 10 fs, in agreement with experiment. Anharmonic coupling leads to Intramolecular Vibrational energy Redistribution (IVR) on the timescale of about 80-100 fs, first to (other) C-C stretching modes of the isoalloxazine rings, then by energy spread over the whole molecule, including low-frequency in-plane modes. The IVR is accompanied by a red-shift and broadening of the emission spectrum. When RF is microsolvated with four water molecules, an overall redshift of optical spectra by about 20 nm is observed but the relaxation dynamics is only slightly affected. For several trajectories, a tendency for hydrogen transfer from water to flavin-nitrogen (N-5) was found.
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.
Ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO), the most important enzyme for the assimilation of carbon into biomass, features a well-known isotope effect with regards to the CO2 carbon atom. This kinetic isotope effect alpha = k (12)/k (13) for the carboxylation step of the RuBisCO reaction sequence, and its microscopic origin, was investigated with the help of cluster models and quantum chemical methods [B3LYP/6-31G(d,p)]. We use a recently proposed model for the RuBisCO active site, in which a water molecule remains close to the reaction center during carboxylation of ribulose-1,5-bisphosphate [B. Kannappan, J.E. Gready, J. Am. Chem. Soc. 130 (2008), 15063]. Alternative active-site models and/or computational approaches were also tested. An isotope effect alpha for carboxylation is found, which is reasonably close to the one measured for the overall reaction, and which originates from a simple frequency shift of the bending vibration of (CO2)-C-12 compared to (CO2)-C-13. The latter is the dominant mode for the product formation at the transition state.
Selective excitation of molecule-surface vibrations in H2 and D2 dissociatively adsorbed on Ru(0001)
(2012)
In this contribution we report about the selective vibrational excitation of H2 and D2 on Ru(0001) as an example for nonadiabatic coupling of an open quantum system to a dissipative environment. We investigate the possibility of achieving state-selective vibrational excitations of H2 and D2 adsorbed on a Ru(0001) surface using picosecond infrared laser pulses. The systems behavior is explored using pulses that are rationally designed and others that are optimized using a time-local variant of Optimal Control Theory. The effects of dissipation on the laser-driven dynamics are studied using the reduced-density matrix formalism. The non-adiabatic couplings between adsorbate and surface are computed perturbatively, for which our recently introduced state-resolved anharmonic rate model is used. It is shown that mode- and state-selective excitation can be achieved in the absence of dissipation when using optimized laser pulses. The inclusion of dissipation in the model reduces the state selectivity and the population transfer yield to highly excited states. In this case, mode activation is most effectively realized by a rational pulse of carefully chosen duration rather than by a locally optimized pulse.
We investigate the effect of intermode coupling and anharmonicity on the excitation and relaxation dynamics of CO on Cu(100). The nonadiabatic coupling of the adsorbate to the surface is treated perturbatively using a position-dependent state-resolved transition rate model. Using the potential energy surface of Marquardt et al. [J. Chem. Phys. 132, 074108 (2010)], which provides an accurate description of intermode interactions, we propose a four-dimensional model that represents simultaneously the diffusion and the desorption of the adsorbate. The system is driven by both rational and optimized infrared laser pulses to favor either selective mode and state excitations or lateral displacement along the diffusion coordinate. The dissipative dynamics is simulated using the reduced density matrix in its Lindblad form. We show that coupling between the degrees of freedom, mediated by the creation and annihilation of electron-hole pairs in the metal substrate, significantly affects the system excitation and relaxation dynamics. In particular, the angular degrees of freedom appear to play an important role in the energy redistribution among the molecule-surface vibrations. We also show that coherent excitation using simple IR pulses can achieve population transfer to a specific target to some extent but does not allow enforcement of the directionality to the diffusion motion.
Angle-resolved two-photon photoemission and high-resolution electron energy loss spectroscopy are employed to derive the electronic structure of a subnanometer atomically precise quasi-one-dimensional graphene nanoribbon (GNR) on Au(111). We resolved occupied and unoccupied electronic bands including their dispersion and determined the band gap, which possesses an unexpectedly large value of 5.1 eV. Supported by density functional theory calculations for the idealized infinite polymer and finite size oligomers, an unoccupied nondispersive electronic state with an energetic position in the middle of the band gap of the GNR could be identified. This state resides at both ends of the ribbon (end state) and is only found in the finite sized systems, i.e., the oligomers.
The scanning tunnelling microscope (STM)-induced switching of a single cyclooctadiene molecule between two stable conformations chemisorbed on a Si(100) surface is investigated using an above threshold model including a neutral ground state and an ionic excited state potential. Switching was recently achieved experimentally with an STM operated at cryogenic temperatures (Nacci et al 2008 Phys. Rev. B 77 121405(R)) and rationalized by a below threshold model using just a single potential energy surface (Nacci et al 2009 Nano Lett. 9 2997).
In the present paper, we show that experimental key findings on the inelastic electron tunnelling (IET) switching can also be rationalized using an above threshold density matrix model, which includes, in addition to the neutral ground state potential, an anionic or cationic excited potential. We use one and two-dimensional potential energy surfaces. Furthermore, the influence of two key parameters of the density matrix description, namely the electronic lifetime of the ionic resonance and the vibrational lifetimes, on the ground state potential are discussed.
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.
With ongoing miniaturization of electronic devices, the need for individually addressable, switchable molecules arises. An example are azobenzenes on surfaces which have been shown to be switchable between trans and cis forms. Here, we examine the "direct" (rather than substrate-mediated) channel of the trans -> cis photoisomerization after pi pi* excitation of tetra-tert-butyl-azobenzene physisorbed on surfaces mimicking Au(111) and Bi(111), respectively. In spirit of the direct channel, the electronic structure of the surface is neglected, the latter merely acting as a rigid platform which weakly interacts with the molecule via Van-der-Waals forces. Starting from thermal ensembles which represent the trans-form, sudden excitations promote the molecules to pi pi*-excited states which are non-adiabatically coupled among themselves and to a n pi*-excited and the ground state, respectively. After excitation, relaxation to the ground state by internal conversion takes place, possibly accompanied by isomerization. The process is described here by "on the fly" semiclassical surface hopping dynamics in conjunction with a semiempirical Hamiltonian (AM1) and configuration-interaction type methods. It is found that steric constraints imposed by the substrate lead to reduced but non-vanishing, trans -> cis reaction yields and longer internal conversion times than for the isolated molecule. Implications for recent experiments for azobenzenes on surfaces are discussed.
The chemistry of water on alpha-alumina kinetics and nuclear quantum effects from first principles
(2012)
Water adsorption on an alumina (alpha-Al2O3) surface is studied here from first principles using periodic density functional theory in the generalized gradient approximation. Two different coverage regimes, low and high, are considered. For the low-coverage regime (with a coverage of 1/4 with respect to the number of coordinatively unsaturated Al sites), possible reactions at the surface such as dissociation, rotation, and diffusion of water and its fragments are investigated, using first principles thermodynamics and kinetics. A microkinetic model is set up with rates calculated from Eyring's transition state theory in order to cover a wide range of time scales. Special emphasis of this study is on the magnitude of quantum effects and on anharmonic corrections, particularly for reactions and dynamics. These have often been neglected in the past for water/alumina systems but can influence the system. This is particularly true for processes involving hydrogen atoms, where, for example, tunneling corrections to reaction rates are found to be important even at room temperature. For a higher-coverage regime (with a coverage of 2 ML), hydrogen dynamics becomes even more complex and is characterized, e.g., by concerted atom motion, strong anharmonicity, and delocalization. In this regime, classical molecular dynamics becomes questionable as well as quantum mechanical treatments based on the harmonic approximation.
The photochemistry as well as electrochemistry of novel donor-acceptor bis(morpholinothiazolyl)maleimides has been investigated. Proper substitution of these diarylethene-type molecular switches leads to the unique situation in which their ring-closure can only be accomplished electrochemically, while ring-opening can only be achieved photochemically. Hence, these switches operate with orthogonal stimuli, i.e. redox potential and light, respectively. The switch system could be optimized by introducing trifluoromethyl groups at the reactive carbon atoms in order to avoid by-product formation during oxidative ring closure. Both photochemical and electrochemical pathways were investigated for methylated, trifluoromethylated, and nonsymmetrical bis(morpholinothiazolyl) maleimides as well as the bis(morpholinothiazolyl) cyclopentene reference compound. With the aid of the nonsymmetrical "mixed" derivative, the mechanism of electrochemically driven ring closure could be elucidated and seems to proceed via a dicationic intermediate generated by two-fold oxidation. All experimental work has been complemented by density functional theory that provides detailed insights into the thermodynamics of the ring-open and closed forms, the nature of their excited states, and the reactivity of their neutral as well as ionized species in different electronic configurations. The particular diarylethene systems described herein could serve in multifunctional (logic) devices operated by different stimuli (inputs) and may pave the way to converting light into electrical energy via photoinduced "pumping" of redox-active meta-stable states.
In this paper, we present quantum dynamical calculations on electron correlation dynamics in atoms and molecules using explicitly time-dependent ab initio configuration interaction theory. The goals are (i) to show that in which cases it is possible to switch off the electronic correlation by ultrashort laser pulses, and (ii) to understand the temporal evolution and the time scale on which it reappears. We characterize the appearance of correlation through electron-electron scattering when starting from an uncorrelated state, and we identify pathways for the preparation of a Hartree-Fock state from the correlated, true ground state. Exemplary results for noble gases, alkaline earth elements, and selected molecules are provided. For Mg we show that the uncorrelated state can be prepared using a shaped ultrashort laser pulse.
We investigate the recombinative desorption of hydrogen and deuterium from a Ru(0001) surface initiated by femtosecond laser pulses. We adopt a quantum mechanical two-state model including three molecular degrees of freedom to describe the dynamics within the desorption induced by electronic transition (DIET) limit. The energy distributions as well as the state-resolved and ensemble properties of the desorbed molecules are analyzed in detail by using the time-energy method. Our results shed light on the experimentally observed 1) large isotopic effects regarding desorption yields and translational energies and 2) the nonequal energy partitioning into internal and translational modes. In particular, it is shown that a single temperature is sufficient to characterize the energy distributions for all degrees of freedom. Further, we confirm that quantization effects play an important role in the determination of the energy partitioning.