Filtern
Erscheinungsjahr
Dokumenttyp
- Wissenschaftlicher Artikel (126)
- Postprint (8)
- Teil eines Buches (Kapitel) (1)
- Sonstiges (1)
Sprache
- Englisch (136)
Gehört zur Bibliographie
- ja (136)
Schlagworte
- 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)
Institut
- Institut für Chemie (136) (entfernen)
Electron transport through molecules treated by LCAO-MO Green's functions with absorbing boundaries
(2004)
In this Letter, we present a method for calculating transport properties of molecular conductors using a time- independent scattering approach based on Green's functions with absorbing boundaries. The method, which has been used before for chemical reaction dynamics in a grid basis [Seideman, Miller, J. Chem. Phys. 96 (1992) 4412], is formulated here in an LCAO-MO form within simple Huckel theory and extended Huckel theory (EHT), respectively. Test calculations are for a quasi-one-dimensional atom chain. As a more realistic application, the organic molecules benzene- 1,4-dithiolate and biphenyl-4,4'-dithiolate between gold electrodes are studied. (C) 2004 Elsevier B.V. All rights reserved
In this contribution, recent advances in the theory of laser and, to a lesser extent, of scanning tunneling microscope (STM) induced cleavage of bonds between an adsorbate and a solid surface, will be reviewed. Special emphasis will be given to the quantum dynamics of electronically non-adiabatic reactions. (c) 2005 Elsevier Ltd. All rights reserved
We report simulations of laser-driven many-electron dynamics by means of the time-dependent configuration interaction singles (doubles) approach. The method accounts for the correlation of ground and excited states, is capable of describing explicitly time-dependent, nonlinear phenomena, and is systematically improvable. Lithium cyanide serves as a molecular test system in which the charge distribution and hence the dipole moment are shown to be switchable, in a controlled fashion, by (a series of) laser pulses which induce selective, state-to-state electronic transitions. One focus of our time-dependent calculations is the question of how fast the transition from the ionic ground state to a specific excited state that is embedded in a multitude of other states can be made, without creating an electronic wave packet. (c) 2005 American Institute of Physics
In this paper we present time-dependent, quantum-dynamical simulations of photoinduced processes at solid surfaces involving nonadiabatic transitions of electrons to and from short-lived intermediate excited states. In particular, two-photon photoemission (2PPE) spectra of naked metal surfaces and free-standing metal films are considered. One major problem in both cases is the presence of electron-electron scattering, which is treated here in various ways. The first way is to adopt an open-system density matrix approach, in which a single electron is weakly coupled to a "bath" of other electrons. The second approach is based on a many-electron Schrodinger equation, which is solved with the help of a time-dependent configuration interactions singles (TD-CIS) method
We apply the multiconfiguration time-dependent Hartree-Fock method to electronic structure calculations and show that quantum chemical information can be obtained with this explicitly time-dependent approach. Different equations of motion are discussed, as well as the numerical cost. The two-electron integrals are calculated using a natural potential expansion, of which we describe the convergence behavior in detail
In this paper we report dynamical simulations of laser-driven, coupled nuclear-electron dynamics for a molecule- surface system. Specifically, the laser desorption of a small molecule (NO) from a metal slab (Pt) in the so-called DIET limit (Desorption Induced by Electronic Transitions), is studied. The excitation of the metal electrons by a laser pulse followed by the formation of a negative ion resonance, its subsequent decay, and the simultaneous desorption of the molecule are all treated within a single quantum mechanical model. This model is based on an earlier theory of Harris and others [S. M. Harris, S. Holloway, and G. R. Darling, J. Chem. Phys. 102, 8235 (1995)], according to which a nuclear degree of freedom is coupled to an electronic one, both propagated on a single non-Born-Oppenheimer potential energy surface. The goals of the present contribution are (i) to make a conceptual connection of this model to the frequently adopted nonadiabatic "multi-state" models of photodesorption, (ii) to understand details of the desorption mechanism, (iii) to explicitly account for the laser pulse, and (iv) to study the photodesorption as a function of the thickness of the metal film, and the laser parameters. As an important methodological aspect we also present a highly efficient numerical scheme to propagate the wave packet in a problem-adapted diabatic basis
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.
The results of a quantum-mechanical study of vibrational relaxation of hydrogen adsorbed on a Si(100) surface with the multi-configurational time-dependent Hartree (MCTDH) method are presented. A two-dimensional subsystem is coupled non-linearly to a bath of harmonic oscillators (phonons of the Si bulk), and the relaxation of subsystem vibrations proceeds primarily via a two-phonon process. Characteristic times of the system evolution agree well with our previous perturbation theory study. The vibrational population decay is non-exponential, exhibiting pronounced recurrences due to finite bath size. The dependence of the lifetimes of the vibrational levels on the bath size and on the coupling details is investigated.
We report quantum chemical calculations, mostly based on density functional theory, on azobenzene and substituted azobenzenes as neutral molecules or ions, in ground and excited states. Both the cis and trans configurations are computed as well as the activation energies to transform one isomer into the other and the possible reaction paths and reaction surfaces along the torsion and inversion modes. All calculations are done for the isolated species, but results are discussed in light of recent experiments aiming at the switching of surface mounted azobenzenes by scanning tunneling microscopes.
An efficient method for the numerical solution of a non-Markovian, open-system density matrix equation of motion in coordinate representation is developed. We apply the scheme to model simulations of the laser-assisted O+H -> OH association reaction in an environment. The suggested approach is based on the application of the time-evolution operator to the "closed-system" part of the overall Hamiltonian and transformation of the open-system equation of motion to the Heisenberg picture suitable for numerical propagation. A dual role of the system-environment coupling with respect to the infrared (ir) laser-driven association of OH is demonstrated: the association probability is increased due to the coupling at relatively weak laser fields, but decreased at strong laser fields. Moreover, at a certain strength of the ir laser field, the association probability does not depend on the strength of the system-bath coupling at all.
A theoretical model for the selective subsurface absorption of atomic hydrogen in a Pd(111) surface by infrared (IR) laser pulses is presented. The dynamics of the adsorbate is studied within the reduced density matrix approach. Energy and phase relaxation of the hydrogen atom are treated using the semigroup formalism. The vibrational excitation leading to subsurface absorption is performed using rationally designed pulses as well as IR laser pulses optimized on- the-fly. It is shown that dissipation can be used as a tool to transfer population to an otherwise inaccessible state via a mechanism known as "laser distillation." We demonstrate that when the reaction path is generalized from a reduced one-dimensional to full three-dimensional treatment of the system, the laser control strategy can prove very different.
The BLUF (blue-light sensing using flavine) domain of the AppA photoreceptor protein from Rhodobacter sphaeroides was modelled by using quantum chemical chromophore plus amino acid models at the (TD-)B3LYP/6-31G* level of theory. The models were based on NMR structures, and further refined by CHARM force field molecular dynamics simulations. The goal is to explain the total redshift by about 10 nm in the UV/Vis spectra of BLUF domains after illumination, and to relate it to structural changes. For this purpose UV/Vis spectra of the available NMR structures were calculated and related to geometrical features. In particular, the hydrogen network embedding the central chromophore is discussed. Specifically, the position of a conserved glutamine, Q63, is found to be important in agreement with findings from previous works. Additionally, however, we find a systematic dependence also on the geometry of a conserved serine, S41. Based on a series of calculations with known structures and with artificial structural models, we argue that indeed the light-induced switching of both Q63 and S41 is necessary to explain the full similar to 10 nm redshift in the light (signalling) state of serine containing BLUF domains. Following or accompanying the double switching, two structurally highly important residues W104 and M106 exchange places, but do not affect the overall UV/ Vis properties of the chromophore.
Quantum chemical calculations of various azobenzene (AB) derivatives have been carried out with the goal to describe the energetics and kinetics of their thermal cis -> trans isomerization. The effects of substituents, in particular their type, number, and positioning, on activation energies have been systematically studied with the ultimate goal to tailor the switching process. Trends observed for mono- and disubstituted species are discussed. A polarizable continuum model is used to study, in an approximate fashion, the cis -> trans isomerization of azobenzenes in solution. The nature of the transition state(s) and its dependence on substituents and the environment is discussed. In particular for push-pull azobenzenes, the reaction mechanism is found to change from inversion in nonpolar solvents to rotation in polar solvents. Concerning kinetics, calculations based on the Eyring transition state theory give usually reliable activation energies and enthalpies when compared to experimentally determined values. Also, trends in the resulting rate constants are correct. Other computed properties such as activation entropies and thus preexponential rate factors are in only moderate agreement with experiment.
In this paper, we report simulations of laser-driven many-electron dynamics by means of the time-dependent configuration interaction singles (TD-CIS) approach. The method is capable of describing explicitly time-dependent phenomena beyond perturbation theory and is systematically improvable. In contrast to most time-dependent density functional methods it also allows us to treat long-range charge-transfer states properly. As an example, the laser-pulse induced charge transfer between a donor (ethylene) and an acceptor molecule (tetracyanoethylene, TCNE) is studied by means of TD-CIS. Also, larger aggregates consisting of several donors and/or acceptors are considered. It is shown that the charge distribution and hence the dipole moments of the systems under study are switchable by (a series of) laser pulses which induce selective, state-to-state electronic transitions.
We report simulations of laser-pulse driven many-electron dynamics by means of a simple, heuristic extension of the time-dependent configuration interaction singles (TD-CIS) approach. The extension allows for the treatment of ionizing states as nonstationary states with a finite, energy-dependent lifetime to account for above-threshold ionization losses in laser-driven many-electron dynamics. The extended TD-CIS method is applied to the following specific examples: (i) state-to-state transitions in the LiCN molecule which correspond to intramolecular charge transfer, (ii) creation of electronic wave packets in LiCN including wave packet analysis by pump-probe spectroscopy, and, finally, (iii) the effect of ionization on the dynamic polarizability of H-2 when calculated nonperturbatively by TD-CIS.
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 femtosecond-laser-induced, substrate-mediated associative desorption of molecular hydrogen and deuterium from a Ru(0001) surface in the so-called DIMET limit is studied theoretically. Two widely used models, a "quantum nonadiabatic" approach and a "classical adiabatic" one are employed and compared to each other. The quantum model is realized by the Monte Carlo wave packet (MCWP) method in the framework of open-system density matrix theory: The classical approach is realized with the help of (frictional) Langevin dynamics with stochastic forces. For both models the same ground-state potential energy surface is used and the same two-temperature model adopted to describe the hot- electron-driven desorption dynamics. Apart from these common features both models are different. Still, both account well for the main experimental findings (Wagner et al. Phys. Rev. B 2005, 72, 205404). In particular, an isotope effect in desorption probabilities, the energy content of the desorbing molecules, and the scaling of these observables with laser fluence are reproduced and explained. The similarity of the results obtained with both models is traced back to the fact that, in the present case, the photodynamics takes place dominantly in the ground electronic state because the electronically excited state is rapidly quenched. The short lifetime of the excited state has also the effect that photoreaction cross sections are typically very small. An IR+vis hybrid scheme, by which the adsorbate is vibrationally excited by IR photons prior to the heating of metal electrons with the vis pulse, is shown to successfully promote the reaction even for strongly coupled adsorbate-surface systems.
Near edge X-ray absorption. ne structure and X-ray photoelectron spectroscopy have been employed to follow the reversible trans to cis isomerization of tetra-tert-butyl-azobenzene (TBA) adsorbed on Au(111). For one monolayer the molecules adopt an adsorption geometry characteristic of the trans-TBA isomer. The azo-bridge (N = N) is aligned nearly parallel to the surface and the phenyl rings exhibit a planar orientation with a small tilt angle <= 4 degrees with respect to the surface normal. Illumination of the molecular layer at 455 nm triggers the trans to cis isomerization which is associated with a pronounced change of the geometrical and electronic structure. The N1s to pi* transition of the central azo-bridge shifts by 0.45 +/- 0.05 eV to higher photon energy and the transition dipole moment (TDM) is tilted by 59 +/- 5 degrees with respect to the surface normal. The pi-system of one phenyl ring is tilted by about 30 degrees with respect to the surface normal, while the second ring plane is oriented nearly perpendicular to the surface. This reorientation is supported by a shift and broadening of the C-H resonances associated with the tert-butyl legs of the molecule. These findings support a configuration of the photo-switched TBA molecule on Au(111) which is comparable to the cis-isomer of the free molecule. In the photo-stationary state 53 +/- 5% of the TBA molecules are switched to the cis configuration. Thermal activation induces the back reaction to trans-TBA.
Fluoroionophores of fluorophore-spacer-receptor format were prepared for detection of PdCl2 by fluorescence enhancement. The fluorophore probes 1-13 consist of a fluorophore group, in alkyl spacer and a dithiomaleonitrile PdCl2 receptor. First, varying the length of the alkylene spacer (compounds 1-3) revealed, dominant through-space pathway for oxidative photoinduced electron transfer (PET) in CH2-bridged dithiomaleonitrile fluoroionophores. Second. fluorescent probes 4-9 containing two anthracene or pyrene fragments connected through CH2 bridges to the dithiomaleonitrile unit were synthesized. Modulation of the oxidation potential (E-Ox) through electron-withdrawing or -donating groups on the anthracene moiety regulates file thermodynamic driving force for oxidative PET (Delta G(PET)) in bis(anthrylmethylthio)maleonitriles and therefore the fluorescence quantum yields (Phi(f)), too. The new concept was confirmed and transferred to pyrenyl ligands, and fluorescence enhancements (FE) greater than 3.2 in the presence of PdCl2 were achieved by 7 and 8 (FE=5.4 and 5.2). Finally, for comparison, monofluorophore ligands 10-13 were synthesized.
Fluoroionophores of fluorophore-spacer-receptor format were prepared for detection of PdCl2 by fluorescence enhancement. The fluorescent probes 1-13 consist of a fluorophore group, an alkyl spacer and a dithiomaleonitrile PdCl2 receptor. First, varying the length of the alkylene spacer (compounds 1-3) revealed a dominant through-space pathway for oxidative photoinduced electron transfer (PET) in CH2-bridged dithiomaleonitrile fluoroionophores. Second, fluorescent probes 4-9 containing two anthracene or pyrene fragments connected through CH2 bridges to the dithiomaleonitrile unit were synthesized. Modulation of the oxidation potential (EOx) through electron-withdrawing or -donating groups on the anthracene moiety regulates the thermodynamic driving force for oxidative PET (GPET) in bis(anthrylmethylthio)maleonitriles and therefore the fluorescence quantum yields (f), too. The new concept was confirmed and transferred to pyrenyl ligands, and fluorescence enhancements (FE) greater than 3.2 in the presence of PdCl2 were achieved by 7 and 8 (FE=5.4 and 5.2). Finally, for comparison, monofluorophore ligands 10-13 were synthesized.
We present a systematic study of the influence of energy and phase relaxation on dynamic polarizability simulations in the linear response regime. The nonperturbative approach is based on explicit electron dynamics using short laser pulses of low intensities. To include environmental effects on the property calculation, we use the time- dependent configuration-interaction method in its reduced density matrix formulation. Both energy dissipation and nonlocal pure dephasing are included. The explicit treatment of time-resolved electron dynamics gives access to the phase shift between the electric field and the induced dipole moment, which can be used to define a useful uncertainty measure for the dynamic polarizability. The nonperturbative treatment is compared to perturbation theory expressions, as applied to a simple model system, the rigid H-2 molecule. It is shown that both approaches are equivalent for low field intensities, but the time-dependent treatment provides complementary information on the phase of the induced dipole moment, which allows for the definition of an uncertainty associated with the computation of the dynamic polarizability in the linear response regime.
Stochastic approach to laser-induced ultrafast dynamics : the desorption of H-2/D-2 from Ru(0001)
(2010)
The desorption of molecular hydrogen and deuterium induced by femtosecond-laser pulses is studied theoretically for the so-called DIMET (Desorption Induced by Multiple Electronic Transitions) process. These investigations are based on nonadiabatic classical Monte Carlo trajectory (CMCT) simulations on a ground and an excited state potential energy surface, including up to all six adsorbate degrees of freedom. The focus is on the hot-electron mediated energy transfer from the surface to the molecule and back, and the energy partitioning between the different degrees of freedom of the desorbing molecules. We first validate for a two-mode model comprising the desorption mode and the internal vibrational coordinate, the classical Monte Carlo trajectory method by comparing with Monte Carlo wavepacket (MCWP) calculations arising from a fully quantum mechanical open-system density matrix treatment. We then proceed by extending the CMCT calculations to include all six nuclear degrees of freedom of the desorbing molecule. This allows for a detailed comparison between theory and experiment concerning isotope effects, energy partitioning (translational, vibrational, and rotational energies and their distributions), and the dependence of these properties on the laser fluence. The most important findings are as follows. (i) CMCT agrees qualitative with the MCWP scheme. (ii) The basic experimental features such as the large isotope effect, the non-linear increase of yield with laser fluence, translationally hot products (in the order of several 1000 K) and non-equipartitioning of translational and internal energies (E-trans > E- vib > E-rot) are well reproduced. (iii) Predictions concerning a strong angular dependence of translational energies at large observation angles are also made.
The nonadiabatic coupling of an adsorbate close to a metallic surface leads to electronic damping of adsorbate vibrations and line broadening in vibrational spectroscopy. Here, a perturbative treatment of the electronic contribution to the lifetime broadening serves as a building block for a new approach, in which anharmonic vibrational transition rates are calculated from a position-dependent coupling function. Different models for the coupling function will be tested, all related to embedding theory. The first two are models based on a scattering approach with (i) a jellium-type and (ii) a density functional theory based embedding density, respectively. In a third variant a further refined model is used for the embedding density, and a semiempirical approach is taken in which a scaling factor is chosen to match harmonic, single-site, first-principles transition rates, obtained from periodic density functional theory. For the example of hydrogen atoms on (adsorption) and below (subsurface absorption) a Pd(111) surface, lifetimes of and transition rates between vibrational levels are computed. The transition rates emerging from different models serve as input for the selective subsurface adsorption of hydrogen in palladium starting from an adsorption site, by using sequences of infrared laser pulses in a laser distillation scheme.
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.
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.
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.
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.
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.
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.
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.
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 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.
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.
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.
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.
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.
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.
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 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.
The role of electronic friction and, more generally, of nonadiabatic effects during dynamical processes at the gas/metal surface interface is still a matter of discussion. In particular, it is not clear if electronic nonadiabaticity has an effect under "mild" conditions, when molecules in low rovibrational states interact with a metal surface. In this paper, we investigate the role of electronic friction on the dissociative sticking and (inelastic) scattering of vibrationally and rotationally cold H-2 molecules at a Ru(0001) surface theoretically. For this purpose, classical molecular dynamics with electronic friction (MDEF) calculations are performed and compared to MD simulations without friction. The two H atoms move on a six-dimensional potential energy surface generated from gradient-corrected density functional theory (DFT), that is, all molecular degrees of freedom are accounted for. Electronic friction is included via atomic friction coefficients obtained from an embedded atom, free electron gas (FEG) model, with embedding densities taken from gradient-corrected DFT. We find that within this model, dissociative sticking probabilities as a function of impact kinetic energies and impact angles are hardly affected by nonadiabatic effects. If one accounts for a possibly enhanced electronic friction near the dissociation barrier, on the other hand, reduced sticking probabilities are observed, in particular, at high impact energies. Further, there is always an influence on inelastic scattering, in particular, as far as the translational and internal energy distribution of the reflected molecules is concerned. Additionally, our results shed light on the role played by the velocity distribution of the incident molecular beam for adsorption probabilities, where, in particular, at higher impact energies, large effects are found.
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.
The time-dependent approach to electronic spectroscopy, as popularized by Heller and coworkers in the 1980's, is applied here in conjunction with linear-response, time-dependent density functional theory to study vibronic absorption, emission and resonance Raman spectra of several diamondoids. Two-state models, the harmonic and the Condon approximations, are used for the calculations, making them easily applicable to larger molecules. The method is applied to nine pristine lower and higher diamondoids: adamantane, diamantane, triamantane, and three isomers each of tetramantane and pentamantane. We also consider a hybrid species “Dia = Dia” – a shorthand notation for a recently synthesized molecule comprising two diamantane units connected by a C[double bond, length as m-dash]C double bond. We resolve and interpret trends in optical and vibrational properties of these molecules as a function of their size, shape, and symmetry, as well as effects of “blending” with sp2-hybridized C-atoms. Time-dependent correlation functions facilitate the computations and shed light on the vibrational dynamics following electronic transitions.