Refine
Year of publication
Document Type
- Article (126)
- Postprint (8)
- Part of a Book (1)
- Other (1)
Language
- English (136)
Is part of the Bibliography
- yes (136) (remove)
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)
- 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)
- Thermal broadening effects (1)
- Vibrational states (1)
- Vibronic spectrum (1)
- X-ray photoelectron spectroscopy (1)
- adamantane (1)
- azobenzene (1)
- azobenzenes (1)
- bismuth (1)
- configuration interactions (1)
- dipole approximation (1)
- dissipative dynamics (1)
- dynamics (1)
- electron dynamics (1)
- electronic wavepackets (1)
- excited states (1)
- fine-structure (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)
- multiphoton processes (1)
- nanoparticle dimers (1)
- optimal control theory (1)
- organic compounds (1)
- pH-Dependent Photoresponsivity (1)
- phonons (1)
- photochromism (1)
- photofragmentation (1)
- quantum control (1)
- quantum dynamics (1)
- rate constants (1)
- ray absorption-spectroscopy (1)
- resonance Raman (1)
- salts (1)
- scattering (1)
- spectra (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
- Institut für Chemie (136) (remove)
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 switching of single cyclooctadiene molecules chemisorbed on a Si(100) surface between two stable conformations, can be achieved with a scanning tunneling microscope [Nacci , Phys. Rev. B 77, 121405(R) (2008)]. Recently, it was shown by quantum chemical and quantum dynamical simulations that major experimental facts can be explained by a single-mode model with switching enforced by inelastic electron tunneling (IET) excitations and perturbed by vibrational relaxation [Nacci , Nano Lett. 9, 2997 (2009)]. In the present paper, we extend the previous theoretical work in several respects: (1) The model is generalized to a two-mode description in which two C2H4 units of COD can move independently; (2) contributions of dipole and, in addition, (cation and anion) resonance-IET rates are considered; (3) the harmonic-linear vibrational relaxation model used previously is generalized to anharmonic vibrations. While the present models highlight generic aspects of IET-switching between two potential minima, they also rationalize specific experimental findings for COD/Si(100): (1) A single-electron excitation mechanism with a linear dependence of the switching rate on tunneling current I, (2) the capability to switch both at negative and positive sample biases, and (3) a crossover temperature around similar to 60 K from an IET-driven, T-independent atom tunneling regime, to classical over-the-barrier isomerization with exponential T-dependence at higher temperatures for a bias voltage of +1.5 V and an average tunneling current of 0.73 nA.
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.
When new covalent organic frameworks (COFs) are designed, the main efforts are typically focused on selecting specific building blocks with certain geometries and properties to control the structure and function of the final COFs. The nature of the linkage (imine, boroxine, vinyl, etc.) between these building blocks naturally also defines their properties. However, besides the linkage type, the orientation, i.e., the constitutional isomerism of these linkages, has rarely been considered so far as an essential aspect. In this work, three pairs of constitutionally isomeric imine-linked donor-acceptor (D-A) COFs are synthesized, which are different in the orientation of the imine bonds (D-C=N-A (DCNA) and D-N=C-A (DNCA)). The constitutional isomers show substantial differences in their photophysical properties and consequently in their photocatalytic performance. Indeed, all DCNA COFs show enhanced photocatalytic H2 evolution performance than the corresponding DNCA COFs. Besides the imine COFs shown here, it can be concluded that the proposed concept of constitutional isomerism of linkages in COFs is quite universal and should be considered when designing and tuning the properties of COFs.
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)].
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.
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.
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.
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.
Graphitic carbon nitride, g-C3N4, is a promising organic photo-catalyst for a variety of redox reactions. In order to improve its efficiency in a systematic manner, however, a fundamental understanding of the microscopic interaction between catalyst, reactants and products is crucial. Here we present a systematic study of water adsorption on g-C3N4 by means of density functional theory and the density functional based tight-binding method as a prerequisite for understanding photocatalytic water splitting. We then analyze this prototypical redox reaction on the basis of a thermodynamic model providing an estimate of the overpotential for both water oxidation and H+ reduction. While the latter is found to occur readily upon irradiation with visible light, we derive a prohibitive overpotential of 1.56 eV for the water oxidation half reaction, comparing well with the experimental finding that in contrast to H-2 production O-2 evolution is only possible in the presence of oxidation cocatalysts.
Graphitic carbon nitride, g-C₃N₄, is a promising organic photo-catalyst for a variety of redox reactions. In order to improve its efficiency in a systematic manner, however, a fundamental understanding of the microscopic interaction between catalyst, reactants and products is crucial. Here we present a systematic study of water adsorption on g-C₃N₄ by means of density functional theory and the density functional based tight-binding method as a prerequisite for understanding photocatalytic water splitting. We then analyze this prototypical redox reaction on the basis of a thermodynamic model providing an estimate of the overpotential for both water oxidation and H⁺ reduction. While the latter is found to occur readily upon irradiation with visible light, we derive a prohibitive overpotential of 1.56 eV for the water oxidation half reaction, comparing well with the experimental finding that in contrast to H₂ production O₂ evolution is only possible in the presence of oxidation cocatalysts.
Graphitic carbon nitride, g-C₃N₄, is a promising organic photo-catalyst for a variety of redox reactions. In order to improve its efficiency in a systematic manner, however, a fundamental understanding of the microscopic interaction between catalyst, reactants and products is crucial. Here we present a systematic study of water adsorption on g-C₃N₄ by means of density functional theory and the density functional based tight-binding method as a prerequisite for understanding photocatalytic water splitting. We then analyze this prototypical redox reaction on the basis of a thermodynamic model providing an estimate of the overpotential for both water oxidation and H⁺ reduction. While the latter is found to occur readily upon irradiation with visible light, we derive a prohibitive overpotential of 1.56 eV for the water oxidation half reaction, comparing well with the experimental finding that in contrast to H₂ production O₂ evolution is only possible in the presence of oxidation cocatalysts.
In this work, the adsorption and splitting of the water molecule by light and/or an external potential is investigated in the frame of (photo-) electrochemical cells using a rutile ruthenium dioxide anode. With the help of periodic density functional calculations, the adsorbed structures of H(2)O and some radicals involved in the splitting process (O, OH, OOH) are obtained and compared with the available experimental results. On the basis of these electronic-structure calculations, we use a method to calculate the stability of the reaction intermediates and conclude on the thermodynamical possibility of the water splitting reaction at the surface. We demonstrate that a moderate overpotential of 0.64 V is required for the reaction to take place at the RuO(2)(110) surface.
The interaction of water with α-alumina (i.e. α-Al2O3) surfaces is important in a variety of applications and a useful model for the interaction of water with environmentally abundant aluminosilicate phases. Despite its significance, studies of water interaction with α-Al2O3 surfaces other than the (0001) are extremely limited. Here we characterize the interaction of water (D2O) with a well defined α-Al2O3(1[1 with combining macron]02) surface in UHV both experimentally, using temperature programmed desorption and surface-specific vibrational spectroscopy, and theoretically, using periodic-slab density functional theory calculations. This combined approach makes it possible to demonstrate that water adsorption occurs only at a single well defined surface site (the so-called 1–4 configuration) and that at this site the barrier between the molecularly and dissociatively adsorbed forms is very low: 0.06 eV. A subset of OD stretch vibrations are parallel to this dissociation coordinate, and thus would be expected to be shifted to low frequencies relative to an uncoupled harmonic oscillator. To quantify this effect we solve the vibrational Schrödinger equation along the dissociation coordinate and find fundamental frequencies red-shifted by more than 1500 cm−1. Within the context of this model, at moderate temperatures, we further find that some fraction of surface deuterons are likely delocalized: dissociatively and molecularly absorbed states are no longer distinguishable.
The interaction of water with a-alumina (i.e. alpha-Al2O3) surfaces is important in a variety of applications and a useful model for the interaction of water with environmentally abundant aluminosilicate phases. Despite its significance, studies of water interaction with alpha-Al2O3 surfaces other than the (0001) are extremely limited. Here we characterize the interaction of water (D2O) with a well defined alpha-Al2O3(1 (1) over bar 02) surface in UHV both experimentally, using temperature programmed desorption and surface-specific vibrational spectroscopy, and theoretically, using periodic-slab density functional theory calculations. This combined approach makes it possible to demonstrate that water adsorption occurs only at a single well defined surface site (the so-called 1-4 configuration) and that at this site the barrier between the molecularly and dissociatively adsorbed forms is very low: 0.06 eV. A subset of OD stretch vibrations are parallel to this dissociation coordinate, and thus would be expected to be shifted to low frequencies relative to an uncoupled harmonic oscillator. To quantify this effect we solve the vibrational Schrodinger equation along the dissociation coordinate and find fundamental frequencies red-shifted by more than 1500 cm(-1). Within the context of this model, at moderate temperatures, we further find that some fraction of surface deuterons are likely delocalized: dissociatively and molecularly absorbed states are no longer distinguishable.
Diarylethene derivatives are photochromic molecular switches, undergoing a ring-opening/-closing reaction by illumination with light. The symmetry of the closed form is determined by the WoodWard Hoffinann rules according to which the reaction proceeds by corirotatory rotation -in that case. Here, we show by a cOrnbined approach of scanning tunneling microscopy (STM) and density functional theory (DFT) calculations that the Open isomer of 4,4'-(4,4'-(perfluorocydopent-1-ene-1,2-diyl)bis(5-methyl-thiophent-4,2,4-dipyridine) (PDTE) retains its open form upon adsorption on a Ag(111) surface. It caribe switched into a closed form, which we identify as the digrotatOly cydization product, by controlled manipulation 'With the STM tip, Evidence of an electric-field dependent switching-process 'is interpreted on the basis of a Simple electroStatic Model, which suggests that the reaction proceedS via an "upright" intermediate state. This pathway thus strongly differs from the switching reaction in solution.
Here we apply and expand the knowledge developed in the case of the H atom to describe high-harmonic generation (HHG) for the H-2 molecule by using time-dependent configuration interaction with single excitations. The implications of using a finite atomic orbital basis set and the impact of a heuristic lifetime model which addresses ionisation losses are discussed. We also examine the influence of the angular momentum of the basis on the computed HHG spectra. Moreover, we discuss the impact of adding diffuse functions and ghost atoms in different geometrical configurations around the molecule. The effects of these additional centres on the HHG spectra are correlated with the physical interpretation of this nonlinear optical phenomenon as given by the three-step model, relating the maximal radial extent of the electron as predicted by the model to the radial extent of the Gaussian basis sets. [GRAPHICS] .
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.
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.
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.
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.
A theoretical model for the selective subsurface absorption of atomic hydrogen in a Pd(111) surface by infrared (IR) laser pulses is presented. The dynamics of the adsorbate is studied within the reduced density matrix approach. Energy and phase relaxation of the hydrogen atom are treated using the semigroup formalism. The vibrational excitation leading to subsurface absorption is performed using rationally designed pulses as well as IR laser pulses optimized on- the-fly. It is shown that dissipation can be used as a tool to transfer population to an otherwise inaccessible state via a mechanism known as "laser distillation." We demonstrate that when the reaction path is generalized from a reduced one-dimensional to full three-dimensional treatment of the system, the laser control strategy can prove very different.
The nonadiabatic coupling of an adsorbate close to a metallic surface leads to electronic damping of adsorbate vibrations and line broadening in vibrational spectroscopy. Here, a perturbative treatment of the electronic contribution to the lifetime broadening serves as a building block for a new approach, in which anharmonic vibrational transition rates are calculated from a position-dependent coupling function. Different models for the coupling function will be tested, all related to embedding theory. The first two are models based on a scattering approach with (i) a jellium-type and (ii) a density functional theory based embedding density, respectively. In a third variant a further refined model is used for the embedding density, and a semiempirical approach is taken in which a scaling factor is chosen to match harmonic, single-site, first-principles transition rates, obtained from periodic density functional theory. For the example of hydrogen atoms on (adsorption) and below (subsurface absorption) a Pd(111) surface, lifetimes of and transition rates between vibrational levels are computed. The transition rates emerging from different models serve as input for the selective subsurface adsorption of hydrogen in palladium starting from an adsorption site, by using sequences of infrared laser pulses in a laser distillation scheme.
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.
In this paper, we perform many-electron dynamics using the time-dependent configuration-interaction method in its reduced density matrix formulation (rho-TDCI). Dissipation is treated implicitly using the Lindblad formalism. To include the effect of ionization on the state-resolved dynamics, we extend a recently introduced heuristic model for ionizing states to the rho-TDCI method, which leads to a reduced density matrix evolution that is not norm-preserving. We apply the new method to the laser-driven excitation of H(2) in a strongly dissipative environment, for which the state-resolve lifetimes are tuned to a few femtoseconds, typical for dynamics of adsorbate at metallic surfaces. Further testing is made on the laser-induced intramolecular charge transfer in a quinone derivative as a model for a molecular switch. A modified scheme to treat ionizing states is proposed to reduce the computational burden associated with the density matrix propagation, and it is thoroughly tested and compared to the results obtained with the former model. The new approach scales favorably (similar to N(2)) with the number of configurations N used to represent the reduced density matrix in the rho-TDCI method, as compared to a N(3) scaling for the model in its original form.
We investigate the effect of intermode coupling and anharmonicity on the excitation and relaxation dynamics of CO on Cu(100). The nonadiabatic coupling of the adsorbate to the surface is treated perturbatively using a position-dependent state-resolved transition rate model. Using the potential energy surface of Marquardt et al. [J. Chem. Phys. 132, 074108 (2010)], which provides an accurate description of intermode interactions, we propose a four-dimensional model that represents simultaneously the diffusion and the desorption of the adsorbate. The system is driven by both rational and optimized infrared laser pulses to favor either selective mode and state excitations or lateral displacement along the diffusion coordinate. The dissipative dynamics is simulated using the reduced density matrix in its Lindblad form. We show that coupling between the degrees of freedom, mediated by the creation and annihilation of electron-hole pairs in the metal substrate, significantly affects the system excitation and relaxation dynamics. In particular, the angular degrees of freedom appear to play an important role in the energy redistribution among the molecule-surface vibrations. We also show that coherent excitation using simple IR pulses can achieve population transfer to a specific target to some extent but does not allow enforcement of the directionality to the diffusion motion.
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.
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.
Molecular photoswitches such as azobenzenes, which undergo photochemical trans <-> cis isomerizations, are often mounted for possible applications on a surface and/or surrounded by other switches, for example, in self-assembled monolayers. This may suppress the isomerization cross section due to possible steric reasons, or, as recently speculated, by exciton coupling to. neighboring switches, leading to ultrafast electronic quenching (Gahl et al., J. Am. Chem. Soc. 2010, 132, 1831). The presence of exciton coupling has been anticipated from a blue shift of the optical absorption band, compared to molecules in solution. From the theory side the need arises to properly analyze and quantify the change of absorption spectra of interacting and adsorbed switches. In particular, suitable methods should be identified, and effects of intermolecule and molecule surface interactions on spectra should be disentangled. In this paper by means of time-dependent Hartree-Fock. (TD-HF), various flavors of time-dependent density functional theory (TD-DFT), and the correlated wave function based, coupled cluster (CC2) method we investigated the 4-nitroazobenzene molecule as an:example: The low-lying singlet excited states in the isolated trans monomer and dieter as well as their composites with a silicon pentamantane nanocluster, which serves also as a crude model for a silicon surface, were determined. As most important results we found that (i) HF, CC2, range-separated density functionals, or global hybrids with large amount of exact exchange are able to describe exciton (Davydov) splitting properly, while hybrids with small amount of exact exchange fail producing spurious charge transfer. (ii) The exciton splitting in a free dimer would lead to a blue shift of the absorption signal; however, this effect is almost nullified or even overcompensated by the shift arising from van der Waals interactions between the two molecules. (iii) Adsorption on the Si "surface" leads to a further, strong red shift for the present system. (iv) At a next-nearest neighbor distance (of similar to 3.6 angstrom), the exciton splitting is similar to 0.3 eV, with or without "surface", suggesting a rapid quenching of the molecular pi ->pi* excitation. At larger distances, exciton splitting decreases rapidly.
While azobenzenes readily photoswitch in solution, their photoisomerization in densely packed self-assembled monolayers (SAMs) can be suppressed. Reasons for this can be steric hindrance and/or electronic quenching, e.g., by exciton coupling. We address these possibilities by means of nonadiabatic molecular dynamics with trajectory surface hopping calculations, investigating the trans -> cis isomerization of azobenzene after excitation into the pi pi* absorption band. We consider a free monomer, an isolated dimer and a dimer embedded in a SAM-like environment of additional azobenzene molecules, imitating in this way the gradual transition from an unconstrained over an electronically coupled to an electronically coupled and sterically hindered, molecular switch. Our simulations reveal that in comparison to the single molecule the quantum yield of the trans -> cis photoisomerization is similar for the isolated dimer, but greatly reduced in the sterically constrained situation. Other implications of dimerization and steric constraints are also discussed.
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.
Somewhat surprisingly, inverted ("O-down") CO adsorbates on NaCl(100) were recently observed experimentally after infrared vibrational excitation (Lau et al., Science, 2020, 367, 175-178). Here we characterize these species using periodic density functional theory and a quantum mechanical description of vibrations. We determine stationary points and minimum energy paths for CO inversion, for low (1/8 and 1/4 monolayers (ML)) and high (1 ML) coverages. Transition state theory is applied to estimate thermal rates for "C-down" to "O-down" isomerization and the reverse process. For the 1/4 ML p(1 x 1) structure, two-dimensional and three-dimensional potential energy surfaces and corresponding anharmonic vibrational eigenstates obtained from the time-independent nuclear Schrodinger equation are presented. We find (i) rather coverage-independent CO inversion energies (of about 0.08 eV or 8 kJ mol(-1) per CO) and corresponding classical activation energies for "C-down" to "O-down" isomerization (of about 0.15 eV or 14 kJ mol(-1) per CO); (ii) thermal isomerization rates at 22 K which are vanishingly small for the "C-down" to "O-down" isomerization but non-negligible for the back reaction; (iii) several "accidentally degenerate" pairs of eigenstates well below the barrier, each pair describing "C-down" to "O-down" localized states.
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.
We combine the stochastic pulse optimization (SPO) scheme with the time-dependent configuration interaction singles method in order to control the high frequency response of a simple molecular model system to a tailored femtosecond laser pulse. For this purpose, we use H-2 treated in the fixed nuclei approximation. The SPO scheme, as similar genetic algorithms, is especially suited to control highly non-linear processes, which we consider here in the context of high harmonic generation. Here, we will demonstrate that SPO can be used to realize a "non-harmonic" response of H2 to a laser pulse. Specifically, we will show how adding low intensity side frequencies to the dominant carrier frequency of the laser pulse and stochastically optimizing their contribution can create a high-frequency spectral signal of significant intensity, not harmonic to the carrier frequency. At the same time, it is possible to suppress the harmonic signals in the same spectral region, although the carrier frequency is kept dominant during the optimization. (C) 2016 AIP Publishing LLC.
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.
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.
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.
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.
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.
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.
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
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.
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).
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 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.
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.
A cationic surfactant containing a spiropyrane unit is prepared exhibiting a dual-responsive adjustability of its surface-active characteristics. The switching mechanism of the system relies on the reversible conversion of the non-ionic spiropyrane (SP) to a zwitterionic merocyanine (MC) and can be controlled by adjusting the pH value and via light, resulting in a pH-dependent photoactivity: While the compound possesses a pronounced difference in surface activity between both forms under acidic conditions, this behavior is suppressed at a neutral pH level. The underlying switching processes are investigated in detail, and a thermodynamic explanation based on a combination of theoretical and experimental results is provided. This complex stimuli-responsive behavior enables remote-control of colloidal systems. To demonstrate its applicability, the surfactant is utilized for the pH-dependent manipulation of oil-in-water emulsions.
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.
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.
In this paper, we present quantum dynamical calculations on electron correlation dynamics in atoms and molecules using explicitly time-dependent ab initio configuration interaction theory. The goals are (i) to show that in which cases it is possible to switch off the electronic correlation by ultrashort laser pulses, and (ii) to understand the temporal evolution and the time scale on which it reappears. We characterize the appearance of correlation through electron-electron scattering when starting from an uncorrelated state, and we identify pathways for the preparation of a Hartree-Fock state from the correlated, true ground state. Exemplary results for noble gases, alkaline earth elements, and selected molecules are provided. For Mg we show that the uncorrelated state can be prepared using a shaped ultrashort laser pulse.
We 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 biconformational switching of single cyclooctadiene molecules chemisorbed on a Si(001) surface was explored by quantum chemical and quantum dynamical calculations and low-temperature scanning tunneling microscopy experiments. The calculations rationalize the experimentally observed switching driven by inelastic electron tunneling (IET) at 5 K. At higher temperatures, they predict a controllable crossover behavior between IET-driven and thermally activated switching, which is fully confirmed by experiment.
Azobenzene-based molecular photoswitches are becoming increasingly important for the development of photoresponsive, functional soft-matter material systems. Upon illumination with light, fast interconversion between a more stable trans and a metastable cis configuration can be established resulting in pronounced changes in conformation, dipole moment or hydrophobicity. A rational design of functional photosensitive molecules with embedded azo moieties requires a thorough understanding of isomerization mechanisms and rates, especially the thermally activated relaxation. For small azo derivatives considered in the gas phase or simple solvents, Eyring’s classical transition state theory (TST) approach yields useful predictions for trends in activation energies or corresponding half-life times of the cis isomer. However, TST or improved theories cannot easily be applied when the azo moiety is part of a larger molecular complex or embedded into a heterogeneous environment, where a multitude of possible reaction pathways may exist. In these cases, only the sampling of an ensemble of dynamic reactive trajectories (transition path sampling, TPS) with explicit models of the environment may reveal the nature of the processes involved. In the present work we show how a TPS approach can conveniently be implemented for the phenomenon of relaxation–isomerization of azobenzenes starting with the simple examples of pure azobenzene and a push–pull derivative immersed in a polar (DMSO) and apolar (toluene) solvent. The latter are represented explicitly at a molecular mechanical (MM) and the azo moiety at a quantum mechanical (QM) level. We demonstrate for the push–pull azobenzene that path sampling in combination with the chosen QM/MM scheme produces the expected change in isomerization pathway from inversion to rotation in going from a low to a high permittivity (explicit) solvent model. We discuss the potential of the simulation procedure presented for comparative calculation of reaction rates and an improved understanding of activated states.
The 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.
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.
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.
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.
A multi-reference study of the byproduct formation for a ring-closed dithienylethene photoswitch
(2015)
Photodriven molecular switches are sometimes hindered in their performance by forming byproducts which act as dead ends in sequences of switching cycles, leading to rapid fatigue effects. Understanding the reaction pathways to unwanted byproducts is a prerequisite for preventing them. This article presents a study of the photochemical reaction pathways for byproduct formation in the photochromic switch 1,2-bis-(3-thienyl)-ethene. Specifically, using single- and multi-reference methods the post-deexcitation reaction towards the byproduct in the electronic ground state S0 when starting from the S1–S0 conical intersection (CoIn), is considered in detail. We find an unusual low-energy pathway, which offers the possibility for the formation of a dyotropic byproduct. Several high-energy pathways can be excluded with high probability.
A multi-reference study of the byproduct formation for a ring-closed dithienylethene photoswitch
(2015)
Photodriven molecular switches are sometimes hindered in their performance by forming byproducts which act as dead ends in sequences of switching cycles, leading to rapid fatigue effects. Understanding the reaction pathways to unwanted byproducts is a prerequisite for preventing them. This article presents a study of the photochemical reaction pathways for byproduct formation in the photochromic switch 1,2-bis-(3-thienyl)-ethene. Specifically, using single- and multi-reference methods the post-deexcitation reaction towards the byproduct in the electronic ground state S0 when starting from the S1–S0 conical intersection (CoIn), is considered in detail. We find an unusual low-energy pathway, which offers the possibility for the formation of a dyotropic byproduct. Several high-energy pathways can be excluded with high probability.
A multi-reference study of the byproduct formation for a ring-closed dithienylethene photoswitch
(2015)
Photodriven molecular switches are sometimes hindered in their performance by forming byproducts which act as dead ends in sequences of switching cycles, leading to rapid fatigue effects. Understanding the reaction pathways to unwanted byproducts is a prerequisite for preventing them. This article presents a study of the photochemical reaction pathways for byproduct formation in the photochromic switch 1,2-bis-(3-thienyl)-ethene. Specifically, using single-and multi-reference methods the post-deexcitation reaction towards the byproduct in the electronic ground state S-0 when starting from the S-1-S-0 conical intersection (CoIn), is considered in detail. We find an unusual low-energy pathway, which offers the possibility for the formation of a dyotropic byproduct. Several high-energy pathways can be excluded with high probability.
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 novel quantum method to deal with typical system-bath dynamical problems is introduced. Subsystem discrete variable representation and bath coherent-state sets are used to write down a multiconfigurational expansion of the wave function of the whole system. With the help of the Dirac-Frenkel variational principle, simple equations of motion-a kind of Schrodinger-Langevin equation for the subsystem coupled to (pseudo) classical equations for the bath-are derived. True dissipative dynamics at all times is obtained by coupling the bath to a secondary, classical Ohmic bath, which is modeled by adding a friction coefficient in the derived pseudoclassical bath equations. The resulting equations are then solved for a number of model problems, ranging from tunneling to vibrational relaxation dynamics. Comparison of the results with those of exact, multiconfiguration time-dependent Hartree calculations in systems with up to 80 bath oscillators shows that the proposed method can be very accurate and might be of help in studying realistic problems with very large baths. To this end, its linear scaling behavior with respect to the number of bath degrees of freedom is shown in practice with model calculations using tens of thousands of bath oscillators.
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.
The cis-trans isomerisation of N-benzylideneaniline (NBA) and derivatives containing a central C=N bond has been investigated experimentally and theoretically. Eight different NBA molecules in three different solvents were irradiated to enforce a photochemical trans (hv) -> cis isomerisation and the kinetics of the thermal backreaction cis (Delta)-> trans were determined by NMR spectroscopy measurements in the temperature range between 193 and 288 K. Theoretical calculations using density functional theory and Eyring transition-state theory were carried out for 12 different NBA species in the gas phase and three different solvents to compute thermal isomerisation rates of the thermal back reaction. While the computed absolute rates are too large, they reveal and explain experimental trends. Time-dependent density functional theory provides optical spectra for vertical transitions and excitation energy differences between trans and cis forms. Together with isomerisation rates, the latter can be used to identify "optimal switches" with good photochromicity and reasonable thermal stability.
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.
We present a rigorous method to set up a system-bath Hamiltonian for the coupling of adsorbate vibrations (the system) to surface phonons (the bath). The Hamiltonian is straightforward to derive and exact up to second order in the environment coordinates, thus capable of treating one- and two-phonon contributions to vibration-phonon coupling. The construction of the Hamiltonian uses orthogonal coordinates for system and bath modes, is based on an embedded cluster approach, and generalizes previous Hamiltonians of a similar type, but avoids several (additional) approximations. While the parametrization of the full Hamiltonian is in principle feasible by a first principles quantum mechanical treatment, here we adopt in the spirit of a QM/MM model a combination of density functional theory (“QM”, for the system) and a semiempirical forcefield (“MM”, for the bath). We apply the Hamiltonian to a fully H-covered Si(100)-(2 × 1) surface, using Fermi’s Golden Rule to obtain vibrational relaxation rates of various H–Si bending modes of this system. As in earlier work it is found that the relaxation is dominated by two-phonon contributions because of an energy gap between the Si–H bending modes and the Si phonon bands. We obtain vibrational lifetimes (of the first excited state) on the order of 2 ps at K. The lifetimes depend only little on the type of bending mode (symmetric vs. antisymmetric, parallel vs. perpendicular to the Si2H2 dimers). They decrease by a factor of about two when heating the surface to 300 K. We also study isotope effects by replacing adsorbed H atoms by deuterium, D. The Si–D bending modes are shifted into the Si phonon band of the solid, opening up one-phonon decay channels and reducing the lifetimes to few hundred fs.
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.
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.
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.
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 /> .
Action spectroscopy has emerged as an analytical tool to probe excited states in the gas phase. Although comparison of gas-phase absorption properties with quantum-chemical calculations is, in principle, straightforward, popular methods often fail to describe many molecules of interest-such as xanthene analogues. We, therefore, face their nano-and picosecond laser-induced photofragmentation with excited-state computations by using the CC2 method and time-dependent density functional theory (TDDFT). Whereas the extracted absorption maxima agree with CC2 predictions, the TDDFT excitation energies are blueshifted. Lowering the amount of Hartree-Fock exchange in the DFT functional can reduce this shift but at the cost of changing the nature of the excited state. Additional bandwidth observed in the photofragmentation spectra is rationalized in terms of multiphoton processes. Observed fragmentation from higher-lying excited states conforms to intense excited-to-excited state transitions calculated with CC2. The CC2 method is thus suitable for the comparison with photofragmentation in xanthene analogues.
Ab initio calculations for XPS chemical shifts of poly(vinyl-trifluoroacetate) using trimer models
(2011)
X-ray photoelectron spectra (XPS) of the polymer poly(vinyl-trifluoroacetate) show C(1s) binding energy shifts which are unusual because they are influenced by atoms which are several bonds away from the probed atom. In this work, the influence of the trifluoroacetate substituent on the 1s ionization potential of the carbon atoms of the polyethylene chain is investigated theoretically using mono-substituted, diad and triad models of trimers representing the polymer. Carbon 1s ionization energies are calculated by the Hartree-Fock theory employing Koopmans' theorem. The influence of the configuration and conformation of the functional groups as well as the degree of substitution are found to be important determinants of XPS spectra. It is further found that the 1s binding energy correlates in a linear fashion, with the total electrostatic potential at the position of the probe atom, and depends not only on nearest neighbor effects. This may have implications for the interpretation of high-resolution XP spectra.
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
Electron transport through molecules treated by LCAO-MO Green's functions with absorbing boundaries
(2004)
In this Letter, we present a method for calculating transport properties of molecular conductors using a time- independent scattering approach based on Green's functions with absorbing boundaries. The method, which has been used before for chemical reaction dynamics in a grid basis [Seideman, Miller, J. Chem. Phys. 96 (1992) 4412], is formulated here in an LCAO-MO form within simple Huckel theory and extended Huckel theory (EHT), respectively. Test calculations are for a quasi-one-dimensional atom chain. As a more realistic application, the organic molecules benzene- 1,4-dithiolate and biphenyl-4,4'-dithiolate between gold electrodes are studied. (C) 2004 Elsevier B.V. All rights reserved
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.
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.
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 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.
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.
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
alpha-Al2O3 surfaces are common in both engineered applications and the environment. Much prior work indicates that their properties, e.g., reactivity, polarity, and charge, change dramatically on interaction with water. Perhaps the simplest question that can be asked of alpha-Al2O3/water interaction is how a single water molecule interacts with the most stable alpha-Al2O3 surface: the alpha-Al2O3(0001). Over the last 15 years, a series of theoretical studies have found that water dissociatively adsorbs on alpha-Al2O3(0001) through two channels. However, to our knowledge no experimental evidence of these dissociation pathways has appeared. By combining sample preparation via supersonic molecular beam dosing, sample characterization via coherent, surface specific vibrational spectroscopy and electronic structure theory, we report the first experimental observation of reaction products of each, theoretically predicted, dissociation channel. These results thus overcome a 15 year old experiment/theory disconnect and make possible a variety of intriguing experiments that promise to provide significant new insights into water/Al2O3 and water/oxide interaction more generally.
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.
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.
α-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.
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.
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.
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.
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.
A variety of azobenzenes were synthesized to study the behavior of their E and Z isomers upon electrochemical reduction. Our results show that the radical anion of the Z isomer is able to rapidly isomerize to the corresponding E configured counterpart with a dramatically enhanced rate as compared to the neutral species. Due to a subsequent electron transfer from the formed E radical anion to the neutral Z starting material the overall transformation is catalytic in electrons; i.e., a substoichiometric amount of reduced species can isomerize the entire mixture. This pathway greatly increases the efficiency of (photo)switching while also allowing one to reach photostationary state compositions that are not restricted to the spectral separation of the individual azobenzene isomers and their quantum yields. In addition, activating this radical isomerization pathway with photoelectron transfer agents allows us to override the intrinsic properties of an azobenzene species by triggering the reverse isomerization direction (Z -> E) by the same wavelength of light, which normally triggers E -> Z isomerization. The behavior we report appears to be general, implying that the metastable isomer of a photoswitch can be isomerized to the more stable one catalytically upon reduction, permitting the optimization of azobenzene switching in new as well as indirect ways.
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.
We present a new global ground state potential energy surface (PES) for carbon monoxide at a coverage of 1/4, on a rigid Ru(0001) surface [Ru(0001)(2x2):CO]. All six adsorbate degrees of freedom are considered. For constructing the PES, we make use of more than 90 000 points calculated with periodic density functional theory using the RPBE exchange-correlation functional and an empirical van der Waals correction. These points are used for interpolation, utilizing a symmetry-adapted corrugation reducing procedure (CRP). Three different interpolation schemes with increasing accuracy have been realized, giving rise to three flavours of the CRP PES. The CRP PES yields in agreement with the DFT reference and experiments, the atop position of CO to be the most stable adsorption geometry, for the most accurate interpolation with an adsorption energy of 1.69 eV. The CRP PES shows that diffusion parallel to the surface is hindered by a barrier of 430 meV, and that dissociation is facilitated but still activated. As a first "real" application and further test of the new potential, the six-dimensional vibrational Schrodinger equation is solved variationally to arrive at fully coupled, anharmonic frequencies and vibrational wavefunctions for the vibrating, adsorbed CO molecule. Good agreement with experiment is found also here. Being analytical, the new PES opens an efficient way towards multidimensional dynamics. (C) 2014 AIP Publishing LLC.
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.