Refine
Year of publication
Document Type
- Article (149)
- Postprint (9)
- Part of a Book (1)
- Conference Proceeding (1)
- Other (1)
Language
- English (161)
Is part of the Bibliography
- yes (161)
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)
- photooxidation (2)
- photoswitches (2)
- reduction (2)
- surface chemistry (2)
- thermochemistry (2)
- tight-binding (2)
- transition (2)
- visible-light (2)
- AM1 calculations (1)
- Ab initio quantum chemical methods and calculations (1)
- Adsorption (1)
- BLUF domains (1)
- Basis sets (1)
- CC2 calculations (1)
- Chemical calculations (1)
- Chemical dynamics (1)
- Classical MD (1)
- Cluster model (1)
- Coherent states (1)
- Computational chemistry (1)
- Correlation function (1)
- DNA origami (1)
- Dark reactions (1)
- Delta-Kohn Sham method (1)
- Density-matrix (1)
- Densityfunctional theory (1)
- Dimer (1)
- Dual-Responsiveness (1)
- Energy (1)
- Excimer (1)
- Generalized Langevin oscillator model (1)
- High-harmonic generation (1)
- Insulating films (1)
- Interchain interactions (1)
- Ionization (1)
- Isotope effect (1)
- Laser induced desorption (1)
- Lasers (1)
- Local density friction approximation (1)
- Manipulation of Emulsion Stability (1)
- Markov processes (1)
- Molecular dynamics with friction (1)
- Molecules (1)
- Monomers (1)
- Oligomers (1)
- Open quantum systems (1)
- PCM (1)
- Perylene (1)
- Photosynthesis (1)
- Poly vinyl alcohol (1)
- Quantum chemistry (1)
- Quantum mechanics (1)
- RuBisCO (1)
- Simulation of polymer XPS (1)
- Spiropyrane (1)
- Surface science (1)
- Surfactant (1)
- Switchable Surfactants (1)
- TD-DFT (1)
- Thermal broadening effects (1)
- Vibrational states (1)
- Vibronic spectrum (1)
- X-ray photoelectron spectroscopy (1)
- adamantane (1)
- azobenzene (1)
- azobenzenes (1)
- bismuth (1)
- blue-light sensor (1)
- configuration interactions (1)
- density functional theory (1)
- dipole approximation (1)
- dissipative dynamics (1)
- dynamics (1)
- electron dynamics (1)
- electronic wavepackets (1)
- excited states (1)
- fine-structure (1)
- flavin (1)
- gold (1)
- graphene (1)
- high harmonic generation (1)
- imines (1)
- ionization (1)
- isomerisation (1)
- isomerization (1)
- isotope effects (1)
- laser chemistry (1)
- micelles (1)
- models (1)
- molecular dynamics (1)
- multiphoton processes (1)
- nanoparticle dimers (1)
- optimal control theory (1)
- organic compounds (1)
- oxidation state (1)
- pH-Dependent Photoresponsivity (1)
- phonons (1)
- photochromism (1)
- photofragmentation (1)
- porphyrin (1)
- quantum control (1)
- quantum dynamics (1)
- rate constants (1)
- ray absorption-spectroscopy (1)
- resonance Raman (1)
- salts (1)
- scanning tunneling microscopy (1)
- scanning tunneling spectroscopy (1)
- scattering (1)
- spectra (1)
- spin state (1)
- states (1)
- strong field (1)
- substituent effects (1)
- surface (1)
- surface hopping (1)
- surface-enhanced Raman scattering (1)
- surfactants (1)
- thermal isomerization (1)
- time-dependent Schrödinger equation (1)
- time-dependent configuration interaction (1)
- ultrafast reactions (1)
- van der Waals forces (1)
- visible light (1)
- xanthenes (1)
Institute
CO oxidation on Ru(0001) is a long-standing example of a reaction that, being thermally forbidden in ultrahigh vacuum, can be activated by femtosecond laser pulses. In spite of its relevance, the precise dynamics of the photoinduced oxidation process as well as the reasons behind the dominant role of the competing CO photodesorption remain unclear. Here we use ab initio molecular dynamics with electronic friction that account for the highly excited and nonequilibrated system created by the laser to investigate both reactions. Our simulations successfully reproduce the main experimental findings: the existence of photoinduced oxidation and desorption, the large desorption to oxidation branching ratio, and the changes in the O K-edge X-ray absorption spectra attributed to the initial stage of the oxidation process. Now, we are able to monitor in detail the ultrafast CO desorption and CO oxidation occurring in the highly excited system and to disentangle what causes the unexpected inertness to the otherwise energetically favored oxidation.
Porous, layered materials containing sp(2)-hybridized carbon and nitrogen atoms, offer through their tunable properties, a versatile route towards tailormade catalysts for electrochemistry and photochemistry. A key molecule interacting with these quasi two-dimensional materials (2DM) is water, and a photo(electro)chemical key reaction catalyzed by them, is water splitting into H-2 and O-2, with the hydrogen evolution reaction (HER) and the oxygen evolution reaction (OER) as half reactions. The complexity of some C/N-based 2DM in contact with water raises special needs for their theoretical modelling, which in turn is needed for rational design of C/N-based catalysts. In this work, three classes of C/N-containing porous 2DM with varying pore sizes and C/N ratios, namely graphitic carbon nitride (g-C3N4), C2N, and poly(heptazine imides) (PHI), are studied with various computational methods. We elucidate the performance of different models and model chemistries (the combination of electronic structure method and basis set) for water and water fragment adsorption in the low-coverage regime. Further, properties related to the photo(electro)chemical activity like electrochemical overpotentials, band gaps, and optical excitation energies are in our focus. Specifically, periodic models will be tested vs. cluster models, and density functional theory (DFT) vs. wavefunction theory (WFT). This work serves as a basis for a systematic study of trends for the photo(electro)chemical activity of C/N-containing layered materials as a function of water content, pore size and density.
Recent molecular beam experiments have shown that water may adsorb molecularly or dissociatively on an α-Al2O3(0001) surface, with enhanced dissociation probability compared to “pinhole dosing”, i.e., adsorption under thermal equilibrium conditions. However, precise information on the ongoing reactions and their relative probabilities is missing. In order to shed light on molecular beam scattering for this system, we perform ab initio molecular dynamics calculations to simulate water colliding with α-Al2O3(0001). We find that single water molecules hitting a cold, clean surface from the gas phase are either reflected, molecularly adsorbed, or dissociated (so-called 1–2 dissociation only). A certain minimum translational energy (above 0.1 eV) seems to be required to enforce dissociation, which may explain the higher dissociation probability in molecular beam experiments. When the surface is heated and/or when refined surface and beam models are applied (preadsorption with water or water fragments, clustering and internal preexcitation in the beam), additional channels open, among them physisorption, water clustering on the surface, and so-called 1–4 and 1–4′ dissociation.
α-Al2O3 surfaces are common in a wide variety of applications and useful models of more complicated, environmentally abundant, alumino-silicate surfaces. While decades of work have clarified that all properties of these surfaces depend sensitively on the crystal face and the presence of even small amounts of water, quantitative insight into this dependence has proven challenging. Overcoming this challenge requires systematic study of the mechanism by which water interacts with various α-Al2O3 surfaces. Such insight is most easily gained for the interaction of small amounts of water with surfaces in ultra high vacuum. In this study, we continue our combined theoretical and experimental approach to this problem, previously applied to water interaction with the α-Al2O3 (0001) and (11̅02) surfaces, now to water interaction with the third most stable surface, that is, the (112̅0). Because we characterize all three surfaces using similar tools, it is straightforward to conclude that the (112̅0) is most reactive with water. The most important factor explaining its increased reactivity is that the high density of undercoordinated surface Al atoms on the (112̅0) surface allows the bidentate adsorption of OH fragments originating from dissociatively adsorbed water, while only monodentate adsorption is possible on the (0001) and (11̅02) surfaces: the reactivity of α-Al2O3 surfaces with water depends strongly, and nonlinearly, on the density of undercoordinated surface Al atoms.
In the light of recent intensity-voltage low energy electron diffraction (LEED-IV) experiments [Surf. Sci. 316, 92 (1994); Surf. Rev. Lett. 10, 487 (2003)], the electronic and geometric structure of a water bilayer adsorbed at the Ru(0001) surface are investigated through first-principles total energy calculations, using periodic slab geometries and gradient-corrected density functional theory (DFT). We consider five possible bilayer structures, all roughly consistent with the LEED-IV analysis (three intact structures and two half-dissociated), and a water single layer at Ru(0001). Adsorption energies and substrate-adsorbate geometry parameters are given and discussed in the light of the experiments. We also give a comparative analysis of the electron density redistribution (Delta rho) and of the dipole moment change (Delta mu) induced by water adsorption on the Ru(0001) surface. In agreement with Feibelman [Science 295, 99 (2002)], the half-dissociated structures are found to be more stable than the intact ones, and their adsorption geometries in better agreement with the LEED-IV data. However, the Delta rho analysis shows that a half-dissociated structure induces a Delta mu>0, which would be incompatible with the experimentally measured decrease of the work function following bilayer adsorption; the latter would be consistent, instead, with the Delta mu < 0 induced by the intact structures. It is the aim of this paper to compare various possible adsorption structures, most of them already considered previously, with one and the same method. For this purpose, thick slabs and restrictive computational parameters are chosen to generally address the accuracy and the limits of DFT in reproducing adsorption energies and bond lengths of water-metal interacting systems
Vibrationally resolved lowest-energy bands of the photoelectron spectra (PES) of adamantane, diamantane, and urotropine were simulated by a time-dependent correlation function approach within the harmonic approximation. Geometries and normal modes for neutral and cationic molecules were obtained from B3LYP hybrid density functional theory (DFT). It is shown that the simulated spectra reproduce the experimentally observed vibrational finestructure (or its absence) quite well. Origins of the finestructure are discussed and related to recurrences of autocorrelation functions and dominant vibrations. Remaining quantitative and qualitative errors of the DFT-derived PES spectra refer to (i) an overall redshift by ∼0.5 eV and (ii) the absence of satellites in the high-energy region of the spectra. The former error is shown to be due to the neglect of many-body corrections to ordinary Kohn-Sham methods, while the latter has been argued to be due to electron-nuclear couplings beyond the Born-Oppenheimer approximation [Gali et al., Nat. Commun. 7, 11327 (2016)].
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.
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 sp2/sp3 hybrid species with C[double bond, length as m-dash]C 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.
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 sp2/sp3 hybrid species with C[double bond, length as m-dash]C 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.
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.
Vibrationally resolved absorption and emission (fluorescence) spectra of perylene and its N-derivatives in gas phase and in solution (acetonitrile) were simulated using a time-dependent approach based on correlation functions determined by density functional theory. By systematically varying the number and position of N atoms, it is shown that the presence of nitrogen heteroatoms has a negligible effect on the molecular structure and geometric distortions upon electronic transitions, while spectral properties change: in particular the number of N atoms is important while their position is less decisive. Thus, the N-substitution can be used to fine-tune the optical properties of perylene-based molecules.
Diamondoids are hydrogen-saturated molecular motifs cut out of diamond, forming a class of materials with tunable optoelectronic properties. In this work, we extend previous work on neutral, closed-shell diamondoids by computing with hybrid density functional theory and time-dependent correlation functions vibrationally broadened absorption spectra of cations and radicals derived from the simplest diamondoid, adamantane, namely, the neutral 1- and 2-adamantyl radicals (C10H15), the 1- and 2-adamantyl cations (C10H15+), and the adamantane radical cation (C10H16+). For selected cases, we also report vibrationally broadened emission, photoelectron, and resonance Raman spectra. Furthermore, the effect of the damping factor on the vibrational fine-structure is studied. The following trends are found: (1) Low-energy absorptions of the adamantyl radicals and cations, and of the adamantane cation, are all strongly red-shifted with respect to adamantane; (2) also, emission spectra are strongly red-shifted, whereas photoelectron spectra are less affected for the cases studied; (3) vibrational fine-structures are reduced compared to those of adamantane; (4) the spectroscopic signals of 1- and 2-adamantyl species are significantly different from each other; and (5) reducing the damping factor has only a limited effect on the vibrational fine-structure in most cases. This suggests that removing hydrogen atoms and/or electrons from adamantane leads to new optoelectronic properties, which should be detectable by vibronic spectroscopy.
Using gradient- and dispersion-corrected density functional theory in connection with ab initio molecular dynamics and efficient, parametrized Velocity-Velocity Autocorrelation Function (VVAF) methodology, we study the vibrational spectra (Vibrational Sum Frequency, VSF, and infrared, IR) of hydroxylated alpha-Al2O3(0001) surfaces with and without additional water. Specifically, by considering a naked hydroxylated surface and the same surface with a particularly stable, "ice-like" hexagonal water later allows us to identify and disentangle main spectroscopic bands of OH bonds, their orientation and dynamics, and the role of water adsorption. In particular, we assign spectroscopic signals around 3700 cm(-1) as being dominated by perpendicularly oriented non-hydrogen bonded aluminol groups, with and without additional water. Furthermore, the thin water layer gives spectroscopic signals which are already comparable to previous theoretical and experimental findings for the solid/(bulk) liquid interface, showing that water molecules closest to the surface play a decisive role in the vibrational response of these systems. From a methodological point of view, the effects of temperature, anharmonicity, hydrogen-bonding, and structural dynamics are taken into account and analyzed, allowing us to compare the calculated IR and VSF spectra with the ones based on normal mode analysis and vibrational density of states. The VVAF approach employed in this work appears to be a computationally accurate yet feasible method to address the vibrational fingerprints and dynamical properties of water/metal oxide interfaces. Published by AIP Publishing.
Water can adsorb molecularly or dissociatively onto different sites of metal oxide surfaces. These adsorption sites can be disentangled using surface-sensitive vibrational spectroscopy. Here, we model Vibrational Sum Frequency (VSF) spectra for various forms of dissociated, deuterated water on a reconstructed, Al-terminated α-Al2O3(0001) surface at submonolayer coverages (the so-called 1-2, 1-4, and 1-4′ modes). Using an efficient scheme based on velocity-velocity autocorrelation functions, we go beyond previous normal mode analyses by including anharmonicity, mode coupling, and thermal surface motion in the framework of ab initio molecular dynamics. In this way, we calculate vibrational density of states curves, infrared, and VSF spectra. Comparing computed VSF spectra with measured ones, we find that relative frequencies of resonances are in quite good agreement and linewidths are reasonably well represented, while VSF intensities coincide not well. We argue that intensities are sensitively affected by local interactions and thermal fluctuations, even at such low coverage, while absolute peak positions strongly depend on the choice of the electronic structure method and on the appropriate inclusion of anharmonicity.
Carbon monoxide on copper surfaces continues to be a fascinating, rich microlab for many questions evolving in surface science. Recently, hot-electron mediated, femtosecond-laser pulse induced dynamics of CO molecules on Cu(100) were the focus of experiments [Inoue et al., Phys. Rev. Lett. 117, 186101 (2016)] and theory [Novko et al., Phys. Rev. Lett. 122, 016806 (2019)], unraveling details of the vibrational nonequilibrium dynamics on ultrashort (subpicoseconds) timescales. In the present work, full-dimensional time-resolved hot-electron driven dynamics are studied by molecular dynamics with electronic friction (MDEF). Dissipation is included by a friction term in a Langevin equation which describes the coupling of molecular degrees of freedom to electron-hole pairs in the copper surface, calculated from gradient-corrected density functional theory (DFT) via a local density friction approximation (LDFA). Relaxation due to surface phonons is included by a generalized Langevin oscillator model. The hot-electron induced excitation is described via a time-dependent electronic temperature, the latter derived from an improved two-temperature model. Our parameter-free simulations on a precomputed potential energy surface allow for excellent statistics, and the observed trends are confirmed by on-the-fly ab initio molecular dynamics with electronic friction (AIMDEF) calculations. By computing time-resolved frequency maps for selected molecular vibrations, instantaneous frequencies, probability distributions, and correlation functions, we gain microscopic insight into hot-electron driven dynamics and we can relate the time evolution of vibrational internal CO stretch-mode frequencies to measured data, notably an observed redshift. Quantitatively, the latter is found to be larger in MDEF than in experiment and possible reasons are discussed for this observation. In our model, in addition we observe the excitation and time evolution of large-amplitude low-frequency modes, lateral CO surface diffusion, and molecular desorption. Effects of surface atom motion and of the laser fluence are also discussed.
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.
Vibrational relaxation of adsorbates is a sensitive tool to probe energy transfer at gas/solid and liquid/solid interfaces. The most direct way to study relaxation dynamics uses time-resolved spectroscopy. Here we report on a non-equilibrium ab initio molecular dynamics (NE-AIMD) methodology to model vibrational relaxation of OH vibrations on a hydroxylated, water-covered alpha-Al2O3(0001) surface. In our NE-AIMD approach, after exciting selected O-H bonds their coupling to surface phonons and to the water adlayer is analyzed in detail, by following both the energy flow in time, as well as the time-evolution of Vibrational Density of States (VDOS) curves. The latter are obtained from Time-dependent Correlation Functions (TCFs) and serve as prototypical, generic representatives of time-resolved vibrational spectra. As most important results, (i) we find a few-picosecond lifetime of the excited modes and (ii) identify both hydrogen-bonded aluminols and water molecules in the adsorbed water layer as main dissipative channels, while the direct coupling to Al2O3 surface phonons is of minor importance on the timescales of interest. Our NE-AIMD/TCF methodology is powerful for complex adsorbate systems, in principle even reacting ones, and opens a way towards time-resolved vibrational spectroscopy.
Combining photochromism and nonlinear optical (NLO) properties of molecular switches-functionalized self-assembled monolayers (SAMs) represents a promising concept toward novel photonic and optoelectronic devices. Using second harmonic generation, density functional theory, and correlated wave function methods, we studied the switching abilities as well as the NLO contrasts between different molecular states of various fulgimide-containing SAMs on Si(111). Controlled variations of the linker systems as well as of the fulgimides enabled us to demonstrate very efficient reversible photoinduced ring-opening/closure reactions between the open and closed forms of the fulgimides. Thus, effective cross sections on the order of 10(-18) cm(-2) are observed. Moreover, the reversible switching is accompanied by pronounced NLO contrasts up to 32%. Further molecular engineering of the photochromic switches and the linker systems may even increase the NLO contrast upon switching.
Thermally stable photoswitches that are driven with low-energy light are rare, yet crucial for extending the applicability of photoresponsive molecules and materials towards, e.g., living systems. Combined ortho-fluorination and -amination couples high visible light absorptivity of o-aminoazobenzenes with the extraordinary bistability of o-fluoroazobenzenes. Herein, we report a library of easily accessible o-aminofluoroazobenzenes and establish structure-property relationships regarding spectral qualities, visible light isomerization efficiency and thermal stability of the cis-isomer with respect to the degree of o-substitution and choice of amino substituent. We rationalize the experimental results with quantum chemical calculations, revealing the nature of low-lying excited states and providing insight into thermal isomerization. The synthesized azobenzenes absorb at up to 600 nm and their thermal cis-lifetimes range from milliseconds to months. The most unique example can be driven from trans to cis with any wavelength from UV up to 595 nm, while still exhibiting a thermal cis-lifetime of 81 days. <br /> [GRAPHICS] <br /> .
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.
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.
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
The thermal Z -> E (back-) isomerization of azobenzenes is a prototypical reaction occurring in molecular switches. It has been studied for decades, yet its kinetics is not fully understood. In this paper, quantum chemical calculations are performed to model the kinetics of an experimental benchmark system, where a modified azobenzene (AzoBiPyB) is embedded in a metal-organic framework (MOF). The molecule can be switched thermally from cis to trans, under solvent-free conditions. We critically test the validity of Eyring transition state theory for this reaction. As previously found for other azobenzenes (albeit in solution), good agreement between theory and experiment emerges for activation energies and activation free energies, already at a comparatively simple level of theory, B3LYP/6-31G* including dispersion corrections. However, theoretical Arrhenius prefactors and activation entropies are in qualitiative disagreement with experiment. Several factors are discussed that may have an influence on activation entropies, among them dynamical and geometric constraints (imposed by the MOF). For a simpler model-Z -> E isomerization in azobenzene-a systematic test of quantum chemical methods from both density functional theory and wavefunction theory is carried out in the context of Eyring theory. Also, the effect of anharmonicities on activation entropies is discussed for this model system. Our work highlights capabilities and shortcomings of Eyring transition state theory and quantum chemical methods, when applied for the Z -> E (back-) isomerization of azobenzenes under solvent-free conditions.
We report on the experimental and theoretical investigation of a considerable increase in the rate for thermal cis -> trans isomerization of azobenzene-containing molecules in the presence of gold nanopartides. Experimentally, by means of UV vis spectroscopy, we studied a series of azobenzene-containing surfactants and 4-nitroazobenzene. We found that in the presence of gold,nanoparticles the thermal lifetime of the cis isomer of the azobenzenecontaining molecules was decreased by up to 3 orders of magnitude in comparison to the lifetime in solution without nanoparticles. The electron transfer between azobenzene-containing molecules and a surface of gold nanopartides is a possible reason to promote the thermal cis trans switching. To investigate the effect of electron attachment to, and withdrawal from, the azobenzene-containing molecules on the isomerization rate, we performed density functional theory calculations of activation energy barriers of the reaction together with Eyring's transition state theory calculations of the rates for azobenzene derivatives with donor and acceptor groups in para position of one of the phenyl rings, as well as for one of the azobenzene-containing surfactants. We found that activation barriers are greatly lowered for azobenzene-containing molecules, both upon electron attachment and withdrawal, which leads, in turn, to a dramatic increase in the thermal isomerization rate.
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
The quantum chemical description of the adsorption, vibrations, and reactions of molecules at periodic solid surfaces is frequently based on a methodological "standard model": density functional theory (DFT) in the generalized gradient approximation (GGA), using plane wave bases and three-dimensional supercells. Although the computationally efficient GGA functionals can be very successful, cases are known where they do not perform so well. Most importantly, activation energies for chemical reactions are typically underestimated, with the consequence of computed reaction rates being too large. In this work, we consider a well-studied model system: water or water fragments adsorbed on an Al-terminated alpha-Al2O3(0001) surface as a test bed for studying the performance of electronic structure methods, both from DFT and wave function theory. On the DFT side, we employ two GGA exchange correlation functionals: PW91 and PBE with and without dispersion corrections, whose results are then compared to those of hybrid functionals B3LYP and HSE06. Further, we follow a periodic wave function approach in the form of local second-order Moller-Plesset perturbation theory, LMP2, on a Hartree-Fock reference. En route, we address issues arising from the choice of the basis set. The key findings of our study are as follows: (i) DFT-GGA adsorption energies are in reasonable agreement with both hybrid-DFT and LMP2 values. In particular, the deviations between the relative energies, corresponding to different adsorption structures, are in the range of the error due to missing dispersion corrections or the basis set error. (ii) Harmonic DFT-GGA vibrational frequencies for oxygen hydrogen stretch modes are by several tens of wavenumbers red-shifted compared to corresponding hybrid-DFT values. The latter are in much better agreement with recent experimental data. (iii) The activation energy for a hydrogen diffusion reaction is grossly underestimated by GGA compared to hybrid-DFT or LMP2, which in turn are quite comparable.
In this paper a perturbation-theory study of vibrational lifetimes for the bending and stretching modes of hydrogen adsorbed on a Si(100) surface is presented. The hydrogen-silicon interaction is treated with a semiempirical bond-order potential. Calculations are performed for H-Si clusters of different sizes. The finite lifetime is due to vibration-phonon coupling, which is assumed to be linear or bilinear in the phonon and nonlinear in the H-Si stretching and bending modes. Lifetimes and vibrational transition rates are evaluated with one- and two-phonon processes taken into account. Temperature effects are also discussed. In agreement with the experiment and previous theoretical treatment it is found that the H-Si (upsilon(s)=1) stretching vibration decays on a nanosecond timescale, whereas for the H-Si (upsilon(b)=1) bending mode a picosecond decay is predicted. For higher-excited vibrations, simple scaling laws are found if the excitation energies are not too large. The relaxation mechanisms for the excited H-Si stretching and the H-Si bending modes are analyzed in detail.
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.
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
The electronic structure of the metal-organic interface of isolated ligand coated gold nanoparticles
(2022)
Light induced electron transfer reactions of molecules on the surface of noble metal nanoparticles (NPs) depend significantly on the electronic properties of the metal-organic interface. Hybridized metal-molecule states and dipoles at the interface alter the work function and facilitate or hinder electron transfer between the NPs and ligand. X-ray photoelectron spectroscopy (XPS) measurements of isolated AuNPs coated with thiolated ligands in a vacuum have been performed as a function of photon energy, and the depth dependent information of the metal-organic interface has been obtained. The role of surface dipoles in the XPS measurements of isolated ligand coated NPs is discussed and the binding energy of the Au 4f states is shifted by around 0.8 eV in the outer atomic layers of 4-nitrothiophenol coated AuNPs, facilitating electron transport towards the molecules. Moreover, the influence of the interface dipole depends significantly on the adsorbed ligand molecules. The present study paves the way towards the engineering of the electronic properties of the nanoparticle surface, which is of utmost importance for the application of plasmonic nanoparticles in the fields of heterogeneous catalysis and solar energy conversion.
In this work, we present theoretical simulations of laser-driven vibrational control of NO adsorbed on a gold surface. Our goal is to tailor laser pulses to selectively excite specific modes and vibrational eigenstates, as well as to favor photodesorption of the adsorbed molecule. To this end, various control schemes and algorithms are applied. For adsorbates at metallic surfaces, the creation of electron hole pairs in the substrate is known to play a dominant role in the transfer of energy from the system to the surroundings. These nonadiabatic couplings are included perturbatively in our reduced density matrix simulations using a generalization of the state-resolved position-dependent anharmonic rate model we recently introduced. An extension of the reduced density matrix is also proposed to provide a sound model for photodesorption in dissipative systems.
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.
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.
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.
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.
Azobenzenes easily photoswitch in solution, while their photoisomerization at surfaces is often hindered. In recent work, it was demonstrated by nonadiabatic molecular dynamics with trajectory surface hopping [Titov et al., J. Phys. Chem. Lett. 2016, 7, 3591-3596] that the experimentally observed suppression of trans -> cis isomerization yields in azobenzenes in a densely packed SAM (self-assembled monolayer) [Gahl et al., J. Am. Chem. Soc. 2010, 132, 1831-1838] is dominated by steric hindrance. In the present work, we systematically study by ground-state Langevin and nonadiabatic surface hopping dynamics, the effects of decreasing packing density on (i) UV/vis absorption spectra, (ii) trans -> cis isomerization yields, and (iii) excited-state lifetimes of photoexcited azobenzene. Within the quantum mechanics/ molecular mechanics models adopted here, we find that above a packing density of similar to 3 molecules/nm(2), switching yields are strongly reduced, while at smaller packing densities, the "monomer limit" is quickly approached. The UV/vis absorption spectra, on the other hand, depend on packing density over a larger range (down to at least similar to 1 molecule/nm(2)). Trends for excited-state lifetimes are less obvious, but it is found that lifetimes of pi pi* excited states decay monotonically with decreasing coverage. Effects of fluorination of the switches are also discussed for single, free molecules. Fluorination leads to comparatively large trans -> cis yields, in combination with long pi pi* lifetimes. Furthermore, for selected systems, also the effects of n pi* excitation at longer excitation wavelengths have been studied, which is found to enhance trans -> cis yields for free molecules but can lead to an opposite behavior in densely packed SAMs.
The alpha-Al2O3(0001) surface has been extensively studied because of its significance in both fundamental research and application. Prior work suggests that in ultra-high-vacuum (UHV), in the absence of water, the so-called Al-I termination is thermodynamically favored, while in ambient, in contact with liquid water, a Gibbsite-like layer is created. While the view of the alpha- Al2O3(0001)/H2O(l) interface appears relatively clear in theory, experimental characterization of this system has resulted in estimates of surface acidity, i.e., isoelectric points, that differ by 4 pH units and surface structure that in some reports has non-hydrogen-bonded surface aluminol (Al-OH) groups and in others does not. In this study, we employed vibrational sum frequency spectroscopy (VSFS) and density functional theory (DFT) simulation to study the surface phonon modes of the differently terminated alpha-Al2O3(0001) surfaces in both UHV and ambient. We find that, on either water dosing of the Al-I in UHV or heat-induced dehydroxylation of the Gibbsite-like in ambient, the surfaces do not interconvert. This observation offers a new explanation for disagreements in prior work on the alpha-Al2O3(0001)/liquid water interface -different preparation methods may create surfaces that do not interconvert-and shows that the surface phonon spectral response offers a novel probe of interfacial hydrogen bonding structure.
Complete sticking at low incidence energies and broad angular scattering distributions at higher energies are often observed in molecular beam experiments on gas-surface systems which feature a deep chemisorption well and lack early reaction barriers. Although CO binds strongly on Ru(0001), scattering is characterized by rather narrow angular distributions and sticking is incomplete even at low incidence energies. We perform molecular dynamics simulations, accounting for phononic (and electronic) energy loss channels, on a potential energy surface based on first-principles electronic structure calculations that reproduce the molecular beam experiments. We demonstrate that the mentioned unusual behavior is a consequence of a very strong rotational anisotropy in the molecule-surface interaction potential. Beyond the interpretation of scattering phenomena, we also discuss implications of our results for the recently proposed role of a precursor state for the desorption and scattering of CO from ruthenium.
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 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.
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.
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.
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.
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.
Recent experiments on laser-dissociation of aligned homonuclear diatomic molecules show an asymmetric forward-backward (spatial) electron-localization along the laser polarization axis. Most theoretical models attribute this asymmetry to interference effects between gerade and ungerade vibronic states. Presumably due to alignment, these models neglect molecular rotations and hence infer an asymmetric (post-dissociation) charge distribution over the two identical nuclei. In this paper, we question the equivalence that is made between spatial electron-localization, observed in experiments, and atomic electron-localization, alluded by these theoretical models. We show that (seeming) agreement between these models and experiments is due to an unfortunate omission of nuclear permutation symmetry, i.e., quantum statistics. Enforcement of the latter requires mandatory inclusion of the molecular rotational degree of freedom, even for perfectly aligned molecules. Unlike previous interpretations, we ascribe spatial electron-localization to the laser creation of a rovibronic wavepacket that involves field-free molecular eigenstates with opposite space-inversion symmetry i.e., even and odd parity. Space-inversion symmetry breaking would then lead to an asymmetric distribution of the (space-fixed) electronic density over the forward and backward hemisphere. However, owing to the simultaneous coexistence of two indistinguishable molecular orientational isomers, our analytical and computational results show that the post-dissociation electronic density along a specified space-fixed axis is equally shared between the two identical nuclei-a result that is in perfect accordance with the principle of the indistinguishability of identical particles. Published under an exclusive license by AIP Publishing.
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.
The time-dependent approach to electronic spectroscopy, as popularized by Heller and co-workers in the 1980s, is applied here in conjunction with linear-response, time-dependent density functional theory to study vibronic absorption and resonance Raman spectra of beta-carotene, with and without a solvent. Two-state models, the harmonic and the Condon approximations are used in order to do so. A new code has been developed which includes excited state displacements, vibrational frequency shifts, and Duschinsky rotation, i.e., mode mixing, for both non-adiabatic spectroscopies. It is shown that Duschinsky rotation has a pronounced effect on the resonance Raman spectra of beta-carotene. In particular, it can explain a recently found anomalous behaviour of the so-called nu(1) peak in resonance Raman spectra [N. Tschirner, M. Schenderlein, K. Brose, E. Schlodder, M. A. Mroginski, C. Thomsen, and P. Hildebrandt, Phys. Chem. Chem. Phys. 11, 11471 (2009)], which shifts with the change in excitation wavelength.
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.
The quest for "chemical accuracy" is becoming more and more demanded in the field of structure and kinetics of molecules at solid surfaces. In this paper, as an example, we focus on the barrier for hydrogen diffusion on a alpha-Al2O3 (0001) surface, aiming for a couple cluster singles, doubles, and perturbative triples [CCSD(T)]-level benchmark. We employ the density functional theory (DFT) optimized minimum and transition state structures reported by Heiden, Usvyat, and Saalfrank [J. Phys. Chem. C 123, 6675 (2019)]. The barrier is first evaluated at the periodic Hartree-Fock and local Moller-Plesset second-order perturbation (MP2) level of theory. The possible sources of errors are then analyzed, which includes basis set incompleteness error, frozen core, density fitting, local approximation errors, as well as the MP2 method error. Using periodic and embedded fragment models, corrections to these errors are evaluated. In particular, two corrections are found to be non-negligible (both from the chemical accuracy perspective and at the scale of the barrier value of 0.72 eV): the correction to the frozen core-approximation of 0.06 eV and the CCSD(T) correction of 0.07 eV. Our correlated wave function results are compared to barriers obtained from DFT. Among the tested DFT functionals, the best performing for this barrier is B3LYP-D3.
Surface-enhanced Raman scattering is a powerful approach to detect molecules at very low concentrations, even up to the single-molecule level. One important aspect of the materials used in such a technique is how much the signal is intensified, quantified by the enhancement factor (EF). Herein we obtained the EFs for gold nanoparticle dimers of 60 and 80 nm diameter, respectively, self-assembled using DNA origami nanotriangles. Cy5 and TAMRA were used as surface-enhanced Raman scattering (SERS) probes, which enable the observation of individual nanoparticles and dimers. EF distributions are determined at four distinct wavelengths based on the measurements of around 1000 individual dimer structures. The obtained results show that the EFs for the dimeric assemblies follow a log-normal distribution and are in the range of 10(6) at 633 nm and that the contribution of the molecular resonance effect to the EF is around 2, also showing that the plasmonic resonance is the main source of the observed signal. To support our studies, FDTD simulations of the nanoparticle's electromagnetic field enhancement has been carried out, as well as calculations of the resonance Raman spectra of the dyes using DFT. We observe a very close agreement between the experimental EF distribution and the simulated values.
Recent experiments using time- and angle-resolved two-photon photoemission (2PPE) spectroscopy at metal/polar adsorbate interfaces succeeded in time-dependent analysis of the process of electron solvation. A fully quantum mechanical, two-dimensional simulation of this process, which explicitly includes laser excitation, is presented here, confirming the origin of characteristic features, such as the experimental observation of an apparently negative dispersion. The inference of the spatial extent of the localized electron states from the angular dependence of the 2PPE spectra has been found to be non-trivial and system-dependent. (C) 2005 American Institute of Physics
Recent progress towards a quantum theory of laser-induced desorption and related phenomena is reviewed, for specific examples. These comprise the photodesorption of NO from Pt(111), the scanning tunnelling microscope and laser- induced desorption and switching of H at Si(100), and the electron stimulated desorption and dissociation of CO at Ru(0001). The theoretical methods used for nuclear dynamics range from open-system density matrix theory over nonadiabatically coupled multi-state models to electron-nuclear wavepackets. Also, aspects of time-dependent spectroscopy to probe ultrafast nonadiabatic processes at surfaces will be considered for the example of two-photon photoemission of solvated electrons in ice layers on Cu(111)
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.
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.
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.
We report on photoinduced remote control of work function and surface potential of a silicon surface modified with a photosensitive self-assembled monolayer consisting of chemisorbed azobenzene molecules (4-nitroazobenzene). Itwas found that the attachment of the organic monolayer increases the work function by hundreds of meV due to the increase in the electron affinity of silicon substrates. The change in the work function on UV light illumination is more pronounced for the azobenzene jacketed silicon substrate (ca. 250 meV) in comparison to 50 meV for the unmodified surface. Moreover, the photoisomerization of azobenzene results in complex kinetics of thework function change: immediate decrease due to light-driven processes in the silicon surface followed by slower recovery to the initial state due to azobenzene isomerization. This behavior could be of interest for electronic devices where the reaction on irradiation should be more pronounced at small time scales but the overall surface potential should stay constant over time independent of the irradiation conditions. Published by AIP Publishing.
Photosensitive azobenzene-containing surfactants have attracted great attention in past years because they offer a means to control soft-matter transformations with light. At concentrations higher than the critical micelle concentration (CMC), the surfactant molecules aggregate and form micelles, which leads to a slowdown of the photoinduced trans -> cis azobenzene isomerization. Here, we combine nonadiabatic dynamics simulations for the surfactant molecules embedded in the micelles with absorption spectroscopy measurements of micellar solutions to uncover the reasons responsible for the reaction slowdown. Our simulations reveal a decrease of isomerization quantum yields for molecules inside the micelles. We also observe a reduction of extinction coefficients upon micellization. These findings explain the deceleration of the trans -> cis switching in micelles of the azobenzene-containing surfactants.
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.
The question as to whether state-selective population of molecular vibrational levels by shaped infrared laser pulses is possible in a condensed phase environment is of central importance for such diverse fields as time-resolved spectroscopy, quantum computing, or "vibrationally mediated chemistry." This question is addressed here for a model system, representing carbon monoxide adsorbed on a Cu(100) surface. Three of the six vibrational modes are considered explicitly, namely, the CO stretch vibration, the CO-surface vibration, and a frustrated translation. Optimized infrared pulses for state-selective excitation of "bright" and "dark" vibrational levels are designed by optimal control theory in the framework of a Markovian open-system density matrix approach, with energy flow to substrate electrons and phonons, phase relaxation, and finite temperature accounted for. The pulses are analyzed by their Husimi "quasiprobability" distribution in time-energy space.
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.
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.
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 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.
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.
Recently, Nocera and co-workers (J. Am. Chem. Soc. 2018, 140, 13711) demonstrated that triaryl borate Lewis acids facilitate the direct electrochemical reduction of triphenylphosphine oxide (TPPO) to triphenylphosphine (TPP). In the present contribution, we report a quantum chemical study unravelling details of the reaction, which also supports the proposed ECrECi mechanism. Alternative electrochemical routes to TPPO reduction facilitated by other Lewis acids (CH3+), or by photocatalysis at semiconductor surfaces, are also briefly discussed.
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 cm2 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.
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.
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 cm2 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.
Vibrational dynamics of adsorbates near surfaces plays both an important role for applied surface science and as a model lab for studying fundamental problems of open quantum systems. We employ a previously developed model for the relaxation of a D-Si-Si bending mode at a D:Si(100)-(2 x 1) surface, induced by a "bath " of more than 2000 phonon modes [Lorenz and P. Saalfrank, Chem. Phys. 482, 69 (2017)], to extend previous work along various directions. First, we use a Hierarchical Effective Mode (HEM) model [Fischer et al., J. Chem. Phys. 153, 064704 (2020)] to study relaxation of higher excited vibrational states than hitherto done by solving a high-dimensional system-bath time-dependent Schrodinger equation (TDSE). In the HEM approach, (many) real bath modes are replaced by (much less) effective bath modes. Accordingly, we are able to examine scaling laws for vibrational relaxation lifetimes for a realistic surface science problem. Second, we compare the performance of the multilayer multiconfigurational time-dependent Hartree (ML-MCTDH) approach with that of the recently developed coherent-state-based multi-Davydov-D2 Ansatz [Zhou et al., J. Chem. Phys. 143, 014113 (2015)]. Both approaches work well, with some computational advantages for the latter in the presented context. Third, we apply open-system density matrix theory in comparison with basically "exact " solutions of the multi-mode TDSEs. Specifically, we use an open-system Liouville-von Neumann (LvN) equation treating vibration-phonon coupling as Markovian dissipation in Lindblad form to quantify effects beyond the Born-Markov approximation. Published under an exclusive license by AIP Publishing.
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 non-Born-Oppenheimer quantum dynamics of pp mu and pd mu molecular ions excited by ultrashort, superintense VUV laser pulses polarized along the molecular axis (z) is studied by the numerical solution of the time-dependent Schrodinger equation within a three-dimensional (3D) model, including the internuclear distance R and muon coordinates z and rho, a transversal degree of freedom. It is shown that in both pp mu and pd mu, muons approximately follow the applied laser field out of phase. After the end of the laser pulse, expectation values < z >, < p >, and < R > demonstrate "post-laser-pulse" oscillations in both pp mu and pd mu. In the case of pd mu, the post-laser-pulse oscillations of < z > and < R > appear as shaped "echo pulses." Power spectra, which are related to high-order harmonic generation (HHG), generated due to muonic and nuclear motion are calculated in the acceleration form. For pd mu it is found that there exists a unique characteristic frequency omega(pd mu)(osc) representing both frequencies of post-laser-pulse muonic oscillations and the frequency of nuclear vibrations, which manifest themselves by very sharp maxima in the corresponding power spectra of pd mu. The homonuclear pp p. ion does not possess such a unique characteristic frequency. The "exact" dynamics and power, and HHG spectra of the 3D model are compared with a Born-Oppenheimer, fixed-nuclei model featuring interesting differences: postpulse oscillations are absent and HHG spectra are affected indirectly or directly by nuclear motion.
The mechanism for signal transduction from the LOV-domains toward the kinase region of phototropin is still not well understood. We have performed molecular dynamics (MD) simulations and CONCOORD calculations on the LOV2 domain of Adiantum capillus-veneris, with the goal to detect possible differences between the two forms of the LOV domain which may not show up in the static crystal structures. Since no such clear differences are found in the MD simulations also, we suggest that the real, biologically active conformation of the LOV domain within the whole phototropin is different from the crystal structure of the isolated LOV domains. The MD simulations do offer, however, insight into details of the dynamics of the dark and illuminated LOV domains, which are discussed in the light of recent experiments
Modern laser technology and ultrafast spectroscopies have pushed the timescales for detecting and manipulating dynamical processes in molecules from the picosecond over femtosecond domains, to the attosecond regime (1 as = 10(-18) s). This way, real-time dynamics of electrons after their photoexcitation can be probed and manipulated. In particular, experiments are moving more and more from atomic and solid state systems to molecules, opening the fields of "molecular electron dynamics" and "attosecond chemistry." Also on the theory side, powerful quantum dynamical tools have been developed to rationalize experiments on ultrafast electron dynamics in molecular species. <br /> In this contribution, we concentrate on theoretical aspects of ultrafast electron dynamics in molecules, mostly driven by lasers. The dynamics will be described with the help of wavefunction-based ab initio methods such as time-dependent configuration interaction (TD-CI) or the multiconfigurational time-dependent Hartree-Fock (MCTDHF) methods. Besides a survey of the methods and their extensions toward, e.g., treatment of ionization, laser pulse optimization, and open quantum systems, two specific examples of applications will be considered: The creation and/or dynamical fate of electronic wavepackets, and the nonlinear optical response to laser pulse excitation in molecules by high harmonic generation (HHG).
The adsorption of molecules to the surface of carbon nanostructures opens a new field of hybrid systems with distinct and controllable properties. We present a microscopic study of the optical absorption in carbon nanotubes functionalized with molecular spiropyran photoswitches. The switching process induces a change in the dipole moment leading to a significant coupling to the charge carriers in the nanotube. As a result, the absorption spectra of functionalized tubes reveal a considerable redshift of transition energies depending on the switching state of the spiropyran molecule. Our results suggest that carbon nanotubes are excellent substrates for the optical readout of spiropyran-based molecular switches. The gained insights can be applied to other noncovalently functionalized one-dimensional nanostructures in an externally induced dipole field.
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.
With recent experimental advances in laser-driven electron dynamics in polyatomic molecules, the need arises for their reliable theoretical modelling. Among efficient, yet fairly accurate methods for many-electron dynamics are Time-Dependent Configuration Interaction Singles (TD-CIS) (a Wave Function Theory (WFT) method), and Real-Time Time-Dependent Density Functional Theory (RT-TD-DFT), respectively. Here we compare TD-CIS combined with extended Atomic Orbital (AO) bases, TD-CIS/AO, with RT-TD-DFT in a grid representation of the Kohn-Sham orbitals, RT-TD-DFT/Grid. Possible ionization losses are treated by complex absorbing potentials in energy space (for TD-CIS/AO) or real space (for RT-TD-DFT), respectively. The comparison is made for two test cases: (i) state-to-state transitions using resonant lasers (pi-pulses), i.e., bound electron motion, and (ii) large-amplitude electron motion leading to High Harmonic Generation (HHG). Test systems are a H-2 molecule and cis- and trans-1,2-dichlorethene, C2H2Cl2, (DCE). From time-dependent electronic energies, dipole moments and from HHG spectra, the following observations are made: first, for bound state-to-state transitions enforced by pi-pulses, TD-CIS nicely accounts for the expected population inversion in contrast to RT-TD-DFT, in agreement with earlier findings. Secondly, when using laser pulses under non-resonant conditions, dipole moments and lower harmonics in HHG spectra are obtained by TD-CIS/AO which are in good agreement with those obtained with RT-TD-DFT/Grid. Deviations become larger for higher harmonics and at low laser intensities, i.e., for low-intensity HHG signals. We also carefully test effects of basis sets for TD-CIS/AO and grid size for RT-TD-DFT/Grid, different exchange-correlation functionals in RT-TD-DFT, and absorbing boundaries. Finally, for the present examples, TD-CIS/AO is observed to be at least an order of magnitude more computationally efficient than RT-TD-DFT/Grid.
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.
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
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
The optical switching of molecular ensembles in a dissipative environment is a subject of various fields of chemical physics and physical chemistry. Here we try to switch arrays of molecules from a stable collective ground state to a state in which all molecules have been transferred to another stable higher-energy configuration. In our model switching proceeds through electronically excited intermediates which are coherently coupled to each other through dipolar interactions, and which decay incoherently within a finite lifetime by coupling to a dissipative environment. The model is quite general, but parameters are chosen to roughly resemble the all-trans -> all-cis isomerization of an array of azobenzene molecules on a surface. Using analytical and optimal control pulses and the concept of "laser distillation," we demonstrate that for various aggregates (dimers up to hexamers), controlled and complete switching should be possible.
In a recent paper [U. Lorenz and P. Saalfrank, Chem. Phys. 482, 69 (2017)], we proposed a robust scheme to set up a system-bath model Hamiltonian, describing the coupling of adsorbate vibrations (system) to surface phonons (bath), from first principles. The method is based on an embedded cluster approach, using orthogonal coordinates for system and bath modes, and an anharmonic phononic expansion of the system-bath interaction up to second order. In this contribution, we use this model Hamiltonian to calculate vibrational relaxation rates of H–Si and D–Si bending modes, coupled to a fully H(D)-covered Si(100)-(2×1) surface, at zero temperature. The D–Si bending mode has an anharmonic frequency lying inside the bath frequency spectrum, whereas the H–Si bending mode frequency is outside the bath Debye band. Therefore, in the present calculations, we only take into account one-phonon system-bath couplings for the D–Si system and both one- and two-phonon interaction terms in the case of H–Si. The computation of vibrational lifetimes is performed with two different approaches, namely, Fermi’s golden rule, and a generalized Bixon-Jortner model built in a restricted vibrational space of the adsorbate-surface zeroth-order Hamiltonian. For D–Si, the Bixon-Jortner Hamiltonian can be solved by exact diagonalization, serving as a benchmark, whereas for H–Si, an iterative scheme based on the recursive residue generation method is applied, with excellent convergence properties. We found that the lifetimes obtained with perturbation theory, albeit having almost the same order of magnitude—a few hundred fs for D–Si and a couple of ps for H–Si—, are strongly dependent on the discretized numerical representation of the bath spectral density. On the other hand, the Bixon-Jortner model is free of such numerical deficiencies, therefore providing better estimates of vibrational relaxation rates, at a very low computational cost. The results obtained with this model clearly show a net exponential decay of the time-dependent survival probability for the H–Si initial vibrational state, allowing an easy extraction of the bending mode “lifetime.” This is in contrast with the D–Si system, whose survival probability exhibits a non-monotonic decay, making it difficult to define such a lifetime. This different behavior of the vibrational decay is rationalized in terms of the power spectrum of the adsorbate-surface system. In the case of D–Si, it consists of several, non-uniformly distributed peaks around the bending mode frequency, whereas the H–Si spectrum exhibits a single Lorentzian lineshape, whose width corresponds to the calculated lifetime. The present work gives some insight into mechanisms of vibration-phonon coupling at surfaces. It also serves as a benchmark for multidimensional system-bath quantum dynamics, for comparison with approximate schemes such as reduced, open-system density matrix theory (where the bath is traced out and a Liouville-von Neumann equation is solved) or approximate wavefunction methods to solve the combined system-bath Schrödinger equation.
Whereas the reversible reduction of azobenzenes has been known for decades, their oxidation is destructive and as a result has been notoriously overlooked. Here, we show that a chain reaction leading to quantitative Z -> E isomerization can be initiated before reaching the destructive anodic peak potential. This hole-catalyzed pathway is accessible to all azobenzenes, without exception, and offers tremendous advantages over the recently reported reductive, radical-anionic pathway because it allows for convenient chemical initiation without the need for electrochemical setups and in the presence of air. In addition, catalytic amounts of metal-free sensitizers, such as methylene blue, can be used as excited-state electron acceptors, enabling a shift of the excitation wavelength to the far red of the azobenzene absorption (up to 660 nm) and providing quantum yields exceeding unity (up to 200%). Our approach will boost the efficiency and sensitivity of optically dense liquid-crystalline and solid photo-switchable materials.
Recent experiments and theory suggest that ground state properties and reactivity of molecules can be modified when placed inside a nanoscale cavity, giving rise to strong coupling between vibrational modes and the quantized cavity field. This is commonly thought to be caused either by a cavity-distorted Born-Oppenheimer ground state potential or by the formation of light-matter hybrid states, vibrational polaritons. Here, we systematically study the effect of a cavity on ground state properties and infrared spectra of single molecules, considering vibration-cavity coupling strengths from zero up to the vibrational ultrastrong coupling regime. Using single-mode models for Li-H and O-H stretch modes and for the NH3 inversion mode, respectively, a single cavity mode in resonance with vibrational transitions is coupled to position-dependent molecular dipole functions. We address the influence of the cavity mode on polariton ground state energies, equilibrium bond lengths, dissociation energies, activation energies for isomerization, and on vibro-polaritonic infrared spectra. In agreement with earlier work, we observe all mentioned properties being strongly affected by the cavity, but only if the dipole self-energy contribution in the interaction Hamiltonian is neglected. When this term is included, these properties do not depend significantly on the coupling anymore. Vibro-polaritonic infrared spectra, in contrast, are always affected by the cavity mode due to the formation of excited vibrational polaritons. It is argued that the quantized nature of vibrational polaritons is key to not only interpreting molecular spectra in cavities but also understanding the experimentally observed modification of molecular reactivity in cavities.
The response of the hydrogen molecular ion, H-2(+), to few-cycle laser pulses of different intensities is simulated. To treat the coupled electron-nuclear motion, we use adiabatic potentials computed with Gaussian-type basis sets together with a heuristic ionization model for the electron and a grid representation for the nuclei. Using this mixed-basis approach, the time-dependent Schrodinger equation is solved, either within the Born-Oppenheimer approximation or with nonadiabatic couplings included. The dipole response spectra are compared to all-grid-based solutions for the three-body problem, which we take as a reference to benchmark the Gaussian-type basis set approaches. Also, calculations employing the fixed-nuclei approximation are performed, to quantify effects due to nuclear motion. For low intensities and small ionization probabilities, we get excellent agreement of the dynamics using Gaussian-type basis sets with the all-grid solutions. Our investigations suggest that high harmonic generation (HHG) and high-frequency response, in general, can be reliably modeled using Gaussian-type basis sets for the electrons for not too high harmonics. Further, nuclear motion destroys electronic coherences in the response spectra even on the time scale of about 30 fs and affects HHG intensities, which reflect the electron dynamics occurring on the attosecond time scale. For the present system, non-Born-Oppenheimer effects are small. The Gaussian-based, nonadiabatically coupled, time-dependent multisurface approach to treat quantum electron-nuclear motion beyond the non-Born-Oppenheimer approximation can be easily extended to approximate wavefunction methods, such as time-dependent configuration interaction singles (TD-CIS), for systems where no benchmarks are available.
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.
Following excited-state chemical shifts in molecular ultrafast x-ray photoelectron spectroscopy
(2022)
The conversion of photon energy into other energetic forms in molecules is accompanied by charge moving on ultrafast timescales. We directly observe the charge motion at a specific site in an electronically excited molecule using time-resolved x-ray photoelectron spectroscopy (TR-XPS). We extend the concept of static chemical shift from conventional XPS by the excited-state chemical shift (ESCS), which is connected to the charge in the framework of a potential model. This allows us to invert TR-XPS spectra to the dynamic charge at a specific atom. We demonstrate the power of TR-XPS by using sulphur 2p-core-electron-emission probing to study the UV-excited dynamics of 2-thiouracil. The method allows us to discover that a major part of the population relaxes to the molecular ground state within 220–250 fs. In addition, a 250-fs oscillation, visible in the kinetic energy of the TR-XPS, reveals a coherent exchange of population among electronic states.
Following excited-state chemical shifts in molecular ultrafast x-ray photoelectron spectroscopy
(2022)
The conversion of photon energy into other energetic forms in molecules is accompanied by charge moving on ultrafast timescales. We directly observe the charge motion at a specific site in an electronically excited molecule using time-resolved x-ray photoelectron spectroscopy (TR-XPS). We extend the concept of static chemical shift from conventional XPS by the excited-state chemical shift (ESCS), which is connected to the charge in the framework of a potential model. This allows us to invert TR-XPS spectra to the dynamic charge at a specific atom. We demonstrate the power of TR-XPS by using sulphur 2p-core-electron-emission probing to study the UV-excited dynamics of 2-thiouracil. The method allows us to discover that a major part of the population relaxes to the molecular ground state within 220–250 fs. In addition, a 250-fs oscillation, visible in the kinetic energy of the TR-XPS, reveals a coherent exchange of population among electronic states.
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.
We simulate the femtosecond-laser-induced desorption dynamics of a diatomic molecule from a metal surface by including the effect of the electron and phonon excitations created by the laser pulse. Following previous models, the laser-induced surface excitation is treated through the two temperature model, while the multidimensional dynamics of the molecule is described by a classical Langevin equation, in which the friction and random forces account for the action of the heated electrons. In this work we propose the additional use of the generalized Langevin oscillator model to also include the effect of the energy exchange between the molecule and the heated surface lattice in the desorption dynamics. The model is applied to study the laser-induced desorption of O-2 from the Ag(110) surface, making use of a six-dimensional potential energy surface calculated within density functional theory. Our results reveal the importance of the phonon mediated process and show that, depending on the value of the electronic density in the surroundings of the molecule adsorption site, its inclusion can significantly enhance or reduce the desorption probabilities.
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 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.
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.
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.
We perform ab initio molecular dynamics simulations to study the femtosecond laser induced desorption of H-2, D-2, and HD from a H: D-saturated Ru(0001) surface. To this aim we have extended the ab initio molecular dynamics with electronic friction (AIMDEF) scheme to include a random force that is a function of a timedependent electronic temperature. The latter characterizes the action of the ultrashort laser pulse according to a two temperature model. This allows us to perform multidimensional, hot-electron driven reaction dynamics and investigate the dependence of the desorption yields on the relative H: D isotope concentration on the surface. Our AIMDEF simulations show that the desorption process takes place in the presence of a heated adsorbate system that clearly influences the desorption dynamics. The heating of the adsorbate system is more (less) pronounced the larger is the concentration of the lighter (heavier) isotope. As a result, we conclude that the presence of H on the surface favors the desorption of molecules, whereas the presence of D hampers it, in agreement with previous experimental observations in which the phenomenon of "dynamical promotion" of a surface reaction had been postulated.