Theoretical Investigation on Triplet Excitation Energy Transfer in Fluorene Dimer

Triplet-triplet energy transfer in fluorene dimer is investigated by combining rate theories with electronic structure calculations. The two key parameters for the control of energy transfer, electronic coupling and reorganization energy, are calculated based on the diabatic states constructed by the constrained density functional theory. The fluctuation of the electronic coupling is further revealed by molecular dynamics simulation. Succeedingly, the diagonal and off-diagonal fluctuations of the Hamiltonian are mapped from the correlation functions of those parameters, and the rate is then estimated both from the perturbationtheory and wavepacket diffusion method. The results manifest that both the static and dynamic fluctuations enhance the rate significantly, but the rate from the dynamic fluctuation is smaller than that from the static fluctuation.     

Low temperature electron transfer in strongly condensed phase

Electron transfer coupled to a collective vibronic degree of freedom is studied in strongly condensed phase and at lower temperatures where quantum fluctuations are essential. Based on an exact representation of the reduced density matrix of the electronic + reaction coordinate compound in terms of path integrals, recent findings on the overdamped limit in quantum dissipative systems are employed. This allows us to give a consistent generalization of the well-known Zusman equations to the quantum domain. Detailed conditions for the range of validity are specified. Using the Wigner transform these results are also extended to the quantum dynamics in full phase space. As an important application electronic transfer rates are derived that comprise adiabatic and nonadiabatic processes in the low temperature regime including nuclear tunneling. Accurate agreement with precise quantum Monte Carlo data is observed.     

Electronic coupling matrix elements from charge constrained density functional theory calculations using a plane wave basis set

We present a plane wave basis set implementation for the calculation of electronic coupling matrix elements of electron transfer reactions within the framework of constrained density functional theory (CDFT). Following the work of Wu and Van Voorhis [J. Chem. Phys. 125, 164105 (2006)], the diabatic wavefunctions are approximated by the Kohn-Sham determinants obtained from CDFT calculations, and the coupling matrix element calculated by an efficient integration scheme. Our results for intermolecular electron transfer in small systems agree very well with high-level ab initio calculations based on generalized Mulliken-Hush theory, and with previous local basis set CDFT calculations. The effect of thermal fluctuations on the coupling matrix element is demonstrated for intramolecular electron transfer in the tetrathiafulvalene-diquinone (Q-TTF-Q(-)) anion. Sampling the electronic coupling along density functional based molecular dynamics trajectories, we find that thermal fluctuations, in particular the slow bending motion of the molecule, can lead to changes in the instantaneous electron transfer rate by more than an order of magnitude. The thermal average, (<|H(ab)|(2)>)(1/2)=6.7 mH, is significantly higher than the value obtained for the minimum energy structure, |H(ab)|=3.8 mH. While CDFT in combination with generalized gradient approximation (GGA) functionals describes the intermolecular electron transfer in the studied systems well, exact exchange is required for Q-TTF-Q(-) in order to obtain coupling matrix elements in agreement with experiment (3.9 mH). The implementation presented opens up the possibility to compute electronic coupling matrix elements for extended systems where donor, acceptor, and the environment are treated at the quantum mechanical (QM) level.     

Influences of Quantum Mechanically Mixed Electronic and Vibrational Pigment States in 2D Electronic Spectra of Photosynthetic Systems: Strong Electronic Coupling Cases

In 2D electronic spectroscopy studies, long‐lived quantum beats have recently been observed in photosynthetic systems, and several theoretical studies have suggested that the beats are produced by quantum mechanically mixed electronic and vibrational states. Concerning the electronic‐vibrational quantum mixtures, the impact of protein‐induced fluctuations was examined by calculating the 2D electronic spectra of a weakly coupled dimer with the Franck‐Condon active vibrational modes in the resonant condition [Fujihashi et al., J. Chem. Phys.­ 2015 , 142, 212403.]. This analysis demonstrated that quantum mixtures of the vibronic resonance are rather robust under the influence of the fluctuations at cryogenic temperatures, whereas the mixtures are eradicated by the fluctuations at physiological temperatures. However, this conclusion cannot be generalized because the magnitude of the coupling inducing the quantum mixtures is proportional to the inter‐pigment electronic coupling. In this study, we explore the impact of the fluctuations on electronic‐vibrational quantum mixtures in a strongly coupled dimer with an off‐resonant vibrational mode. Toward this end, we calculate energy transfer dynamics and 2D electronic spectra of a model dimer that corresponds to the most strongly coupled bacteriochlorophyll molecules in the Fenna‐Matthews‐Olson complex in a numerically accurate manner. The quantum mixtures are found to be robust under the exposure of protein‐induced fluctuations at cryogenic temperatures, irrespective of the resonance. At 300 K, however, the quantum mixing is disturbed more strongly by the fluctuations, and therefore, the beats in the 2D spectra become obscure even in a strongly coupled dimer with a resonant vibrational mode. Further, the overall behaviors of the energy transfer dynamics are demonstrated to be dominated by the environment and coupling between the 0 0 vibronic transitions as long as the Huang‐Rhys factor of the vibrational mode is small. The electronic‐vibrational quantum mixtures do not necessarily play a significant role in electronic energy transfer dynamics despite contributing to the enhancement of long‐lived quantum beating in the 2D spectra.     

Non-exponentiality in electron transfer kinetics: Static versus dynamic disorder models

Non-exponential electron transfer kinetics in complex systems are often analyzed in terms of a quenched, static disorder model. In this work we present an alternative analysis in terms of a simple dynamic disorder model where the solvent is characterized by highly non-exponential dynamics. We consider both low and high barrier reactions. For the former, the main result is a simple analytical expression for the survival probability of the reactant. In this case, electron transfer, in the long time, is controlled by the solvent polarization relaxation—in agreement with the analyses of Rips and Jortner and of Nadler and Marcus. The short time dynamics is also non-exponential, but for different reasons. The high barrier reactions, on the other hand, show an interesting dynamic dependence on the electronic coupling element,V _el.     

Analysis of the vibronic fine structure in circularly polarized emission spectra from chiral molecular aggregates

Using a Frenkel-exciton model, the degree of circular polarization of the luminescence (g(lum)) from one-dimensional, helical aggregates of chromophoric molecules is investigated theoretically. The coupling between the electronic excitation and a local, intramolecular vibrational mode is taken into account. Analytical expressions for the fluorescence band shape and g(lum) are presented for the case of strong and weak electronic coupling between the chromophoric units. Results are compared to those from numerical calculations obtained using the three particle approximation. g(lum) for the 0-0 vibronic band is found to be independent of the relative strength of electronic coupling between chromophores and excitation-vibration coupling. It depends solely on the number of coherently coupled molecules. In contrast, for the higher vibronic transitions[g(lum)] decreases with decreasing strength of the electronic coupling. In the limit of strong electronic coupling, [g(lum)] is almost constant throughout the series of vibronic transitions but for weak coupling [g(lum)] becomes vanishingly small for all vibronic transitions except for the 0-0 transition. The results are interpreted in terms of dynamic localization of the excitation during the zero point vibrational motion in the excited state of the aggregate. It is concluded that circular polarization measurements provide an independent way to determine the coherence size and bandwidth of the lowest exciton state for chiral aggregates.     

Calculation of vibronic couplings for phenoxyl/phenol and benzyl/toluene self-exchange reactions: implications for proton-coupled electron transfer mechanisms

The vibronic couplings for the phenoxyl/phenol and the benzyl/toluene self-exchange reactions are calculated with a semiclassical approach, in which all electrons and the transferring hydrogen nucleus are treated quantum mechanically. In this formulation, the vibronic coupling is the Hamiltonian matrix element between the reactant and product mixed electronic-proton vibrational wavefunctions. The magnitude of the vibronic coupling and its dependence on the proton donor-acceptor distance can significantly impact the rates and kinetic isotope effects, as well as the temperature dependences, of proton-coupled electron transfer reactions. Both of these self-exchange reactions are vibronically nonadiabatic with respect to a solvent environment at room temperature, but the proton tunneling is electronically nonadiabatic for the phenoxyl/phenol reaction and electronically adiabatic for the benzyl/toluene reaction. For the phenoxyl/phenol system, the electrons are unable to rearrange fast enough to follow the proton motion on the electronically adiabatic ground state, and the excited electronic state is involved in the reaction. For the benzyl/toluene system, the electrons can respond virtually instantaneously to the proton motion, and the proton moves on the electronically adiabatic ground state. For both systems, the vibronic coupling decreases exponentially with the proton donor-acceptor distance for the range of distances studied. When the transferring hydrogen is replaced with deuterium, the magnitude of the vibronic coupling decreases and the exponential decay with distance becomes faster. Previous studies designated the phenoxyl/phenol reaction as proton-coupled electron transfer and the benzyl/toluene reaction as hydrogen atom transfer. In addition to providing insights into the fundamental physical differences between these two types of reactions, the present analysis provides a new diagnostic for differentiating between the conventionally defined hydrogen atom transfer and proton-coupled electron transfer reactions.     

On quantifying the vibronic interaction between electronic states

An accurate estimation of the interstate vibronic coupling strength is of particular relevance for the treatment of nonadiabatic dynamics. This is not a trivial task because direct interactions between electronic states have to be separated from intrinsic frequency shifts. Surprisingly, this issue has not been discussed in detail in the literature so far. An analysis of the error dependence is given for two formulas derived from linear vibronic coupling theory. The difficulty in estimating the interstate coupling parameters is shown to originate from the initially unknown contribution of the diagonal quadratic coupling coefficients to the total vibronic coupling. An interpretation of the error analysis including a numerical case study is followed by a more general discussion of the different mechanisms that can shape adiabatic electronic potential energy functions. Qualitative criteria are formulated for the differentiation between interstate and intrastate vibronic coupling effects based on electronic structure information. These ideas are then applied to investigate vibronic coupling problems in pyrazine as well as trans- and cis-hexatriene.     

Coherence transfer processes in vibrational relaxation of polyatomic molecules in condensed media

The vibrational relaxation of a polyatomic molecule in a condensed host is studied by a consideration of two molecular vibrations. Relaxation processes, intermode coupling terms and vibrational frequency fluctuation contributions are retained. Population decay (T₁), dephasing (T₂), and coherence transfer rates are evaluated through second order in the limit where the host bath dynamics are rapid compared to these molecular timescales. The rates are expressed in terms of temperature and frequency dependent bath correlation functions. For the special case of a three level system (the ground state and ones where one of the two vibrational modes is excited) the important effects of anharmonicity are incorporated. It is shown that certain coherence transfer terms involve zero frequency bath correlation functions, so they should be larger than the high frequency ones which obey modified energy gap laws. A discussion is presented of the types of interactions which may contribute to these coherence transfer processes.     

Intermittent single-molecule interfacial electron transfer dynamics

We report on single-molecule studies of photosensitized interfacial electron transfer (ET) processes in Coumarin 343 (C343)-TiO(2) nanoparticles (NP) and Cresyl Violet (CV(+))-TiO(2) NP systems, using time-correlated single-photon counting coupled with scanning confocal fluorescence microscopy. Fluorescence intensity trajectories of individual dye molecules adsorbed on a semiconductor NP surface showed fluorescence fluctuations and blinking, with time constants distributed from milliseconds to seconds. The fluorescence fluctuation dynamics were found to be inhomogeneous from molecule to molecule and from time to time, showing significant static and dynamic disorders in the interfacial ET reaction dynamics. We attribute fluorescence fluctuations to the interfacial ET reaction rate fluctuations, associating redox reactivity intermittency with the fluctuations of molecule-TiO(2) electronic and Franck-Condon coupling. Intermittent interfacial ET dynamics of individual molecules could be characteristic of a surface chemical reaction strongly involved with and regulated by molecule-surface interactions. The intermittent interfacial reaction dynamics that likely occur among single molecules in other interfacial and surface chemical processes can typically be observed by single-molecule studies but not by conventional ensemble-averaged experiments.     

Dynamic nature of the intramolecular electronic coupling mediated by a solvent molecule: a computational study

We present a combined Molecular Dynamics/Quantum Chemical study of the solvent-mediated electronic coupling between an electron donor and acceptor in a C-clamp molecule. We characterize the coupling fluctuations due to the solvent motion for different solvents (acetonitrile, benzene, 1,3-diisopropyl-benzene) for the charge separation and the charge recombination processes. The time scale for solvent-induced coupling fluctuation is approximately 0.1 ps. The effect of these fluctuations on the observed rate is discussed using a recently developed theoretical model. We show that, while the microscopic charge transfer process is very complicated and its computational modeling very subtle, the macroscopic phenomenology can be captured by the standard models. Analyzing the contribution to the coupling given by different solvent orbitals, we find that many solvent orbitals mediate the electron transfer and that paths through different solvent orbitals can interfere constructively or destructively. A relatively small subset of substrate-solvent configurations dominate contributions to solvent-mediated coupling. This subset of configurations is related to the electronic structure of the C-clamp molecule.     

Quantum versus classical electron transfer energy as reaction coordinate for the aqueous Ru2+/Ru3+ redox reaction

Applying density functional theory (DFT)-based molecular dynamics simulation methods we investigate the effect of explicit treatment of electronic structure on the solvation free energy of aqueous Ru²⁺ and Ru³⁺.Our approach is based on the Marcus theory of redox half reactions, focussing on the vertical energy gap for reduction or oxidation of a single aqua ion. We compare the fluctuations of the quantum and classical energy gap along the same equilibrium ab initio molecular dynamics trajectory for each oxidation state. The classical gap is evaluated using a standard point charge model for the charge distribution of the solvent molecules (water). The quantum gap is computed from the full DFT electronic ground state energies of reduced and oxidized species, thereby accounting for the delocalization of the electron in the donor orbital and reorganization of the electron cloud after electron transfer (ET). The fluctuations of the quantum ET energy are well approximated by gaussian statistics giving rise to parabolic free energy profiles. The curvature is found to be independent of the oxidation state in agreement with the linear response assumption underlying Marcus theory. By contrast, the diabatic free energy curves evaluated using the classical gap as order parameter, while also quadratic, are asymmetric reflecting the difference in oxidation state. The response of these two order parameters is further analysed by a comparison of the spectral density of the fluctuations and the corresponding reorganization free energies.     

Nonlinearity and irreversibility in the theory of the solvent effect of proton transfer dynamics. II. Dissipative system

The quantum-statistical theory of the dynamics of proton transfer in solvated H-bond complexes is formulated. The theory takes into account a H-bond complex interaction with an environmental fast electronic polarisation as well as the coupling to the environment slow degrees of freedom that are connected with vibrational-rotational motion. The equation of motion for a reduced density matrix is derived in the form of a nonlinear generalised master equation. For the dynamics of the proton transfer in a symmetric double-well potential, the kinetic equation for the reactant state probability density has also been derived and solved. Results for different environments and temperatures are presented.     

Nonadiabatic quantum dynamics of Br2 in solid Ar: A four-dimensional study of the B to C state predissociation

Quantum dynamics simulations are performed for a diatomics-in-molecules based model of Br₂ in solid Ar which incorporates four nuclear degrees of freedom and four electronic states. The nuclear motions comprise two large amplitude coordinates describing the Br₂ bond distance and an effective symmetry-preserving matrix mode. Two symmetry-lowering harmonic modes are added in the spirit of linear vibronic coupling theory. Initiating the dynamics on the B state by means of an ultrafast laser pulse, nonadiabatic transitions to the two degenerate C states are monitored and the effect of vibrational preexcitation in the electronic ground state is investigated.     

Charge transfer vibronic transitions in uranyl tetrachloride compounds

The electronic and vibronic interactions of uranyl (UO(2))(2+) in three tetrachloride crystals have been investigated with spectroscopic experiments and theoretical modeling. Analysis and simulation of the absorption and photoluminescence spectra have resulted in a quantitative understanding of the charge transfer vibronic transitions of uranyl in the crystals. The spectra obtained at liquid helium temperature consist of extremely narrow zero-phonon lines (ZPL) and vibronic bands. The observed ZPLs are assigned to the first group of the excited states formed by electronic excitation from the 3σ ground state into the f(δ,?) orbitals of uranyl. The Huang-Rhys theory of vibronic coupling is modified successfully for simulating both the absorption and luminescence spectra. It is shown that only vibronic coupling to the axially symmetric stretching mode is Franck-Condon allowed, whereas other modes are involved through coupling with the symmetric stretching mode. The energies of electronic transitions, vibration frequencies of various local modes, and changes in the O═U═O bond length of uranyl in different electronic states and in different coordination geometries are evaluated in empirical simulations of the optical spectra. Multiple uranyl sites derived from the resolution of a superlattice at low temperature are resolved by crystallographic characterization and time- and energy-resolved spectroscopic studies. The present empirical simulation provides insights into fundamental understanding of uranyl electronic interactions and is useful for quantitative characterization of uranyl coordination.     

Nonadiabatic coupling and diabatic electronic population dynamics on 11A2 and 11B1 states of ozone molecule

In this work, we examine nonadiabatic population dynamics for 1¹B₁ and 1¹A₂ states of ozone molecule (O₃). In O₃, two lowest singlet excited states, ¹A₂ and ¹B₁, can be coupled. Thus, population transfer between them occurs through the seam involving these two states. At any point of the seam (conical intersection), the Born-Oppenheimer approximation breaks down, and it is necessary to investigate nonadiabatic dynamics. We consider a linear vibronic coupling Hamiltonian model and evaluate vibronic coupling constant, diabatic frequencies for three modes of O₃, bilinear and quadratic coupling constants for diabatic potentials, displacements, and Huang-Rhys coupling constants using ab initio calculations. The electronic structure calculations have been performed at the multireference configuration interaction and complete active space with second-order perturbation theory with a full-valence complete active space self-consistent field methods and augmented Dunning's standard correlation-consistent-polarized quadruple zeta basis set to determine ab initio potential energy surfaces for the ground state and first two excited states of O₃, respectively. We have chosen active space comprising 18 electrons distributed over 12 active orbitals. Our calculations predict the linear vibronic coupling constant 0.123 eV. We have obtained the population on the 1¹B₁ and 1¹A₂ excited electronic states for the first 500 fs after photoexcitation.     

Proton transfer mediated by the vibronic coupling in oxygen core ionized states of glyoxalmonoxime studied by infrared-X-ray pump-probe spectroscopy

The theory of IR-X-ray pump-probe spectroscopy beyond the Born-Oppenheimer approximation is developed and applied to the study of the dynamics of intramolecular proton transfer in glyoxalmonoxime leading to the formation of the tautomer 2-nitrosoethenol. Due to the IR pump pulses the molecule gains sufficient energy to promote a proton to a weakly bound well. A femtosecond X-ray pulse snapshots the wave packet route and, hence, the dynamics of the proton transfer. The glyoxalmonoxime molecule contains two chemically nonequivalent oxygen atoms that possess distinct roles in the hydrogen bond, a hydrogen donor and an acceptor. Core ionizations of these form two intersecting core-ionized states, the vibronic coupling between which along the OH stretching mode partially delocalizes the core hole, resulting in a hopping of the core hole from one site to another. This, in turn, affects the dynamics of the proton transfer in the core-ionized state. The quantum dynamical simulations of X-ray photoelectron spectra of glyoxalmonoxime driven by strong IR pulses demonstrate the general applicability of the technique for studies of intramolecular proton transfer in systems with vibronic coupling.     

Calculation of the vibronic structure of the X2Pi photoelectron spectra of XCN, X=F, Cl, and Br

The vibronic structure of the photoelectron spectra of the X (2)Pi state of XCN(+) (X=F, Cl, and Br) has been calculated, assuming that the X (2)Pi state can be considered as an isolated electronic state. The Renner-Teller coupling of the two components of the (2)Pi state via the degenerate bending mode as well as spin-orbit coupling effects are taken into account. The two stretching modes are treated within the so-called linear vibronic-coupling model. The vibronic and spin-orbit parameters have been determined by accurate ab initio electronic-structure calculations. While spin-orbit effects are small in FCN(+), the large spin-orbit splitting of the X (2)Pi state of the BrCN(+) leads to a complete quenching of the Renner-Teller effect. The X (2)Pi state of the ClCN(+) is shown to be of particular interest: here the resonance condition for linear-relativistic Renner-Teller coupling is approximately fulfilled. This coupling mechanism leads to a significant intensity transfer to vibronic levels with odd quanta of the bending mode. The calculated spectrum indicates that this novel relativistic vibronic-coupling effect should be observable in high-resolution (electron energy resolution of the order of a few meV) photoelectron spectra of ClCN.     

Probing ground-state single-electron self-exchange across a molecule-metal interface

We have probed single-molecule redox reaction dynamics of hemin (chloride) adsorbed on Ag nanoparticle surfaces by single-molecule surface-enhanced Raman spectroscopy (SMSERS) combined with spectroelectrochemistry. Redox reaction at the molecule/Ag interface is identified and probed by the prominent fluctuations of the Raman frequency of a specific vibrational mode, ν(4), which is a typical marker of the redox state of the iron center in a hemin molecule. On the basis of the autocorrelation and cross-correlation analysis of the single-molecule Raman spectral trajectories and the control measurements of single-molecule spectroelectochemistry and electrochemical STM, we suggest that the single-molecule redox reaction dynamics at the hemin-Ag interface is primarily driven by thermal fluctuations. The spontaneous fluctuation dynamics of the single-molecule redox reaction is measured under no external electric potential across the molecule-metal interfaces, which provides a novel and unique approach to characterize the interfacial electron transfer at the molecule-metal interfaces. Our demonstrated approaches are powerful for obtaining molecular coupling and dynamics involved in interfacial electron transfer processes. The new information obtained is critical for a further understanding, design, and manipulation of the charge transfer processes at the molecule-metal interface or metal-molecule-metal junctions, which are fundamental elements in single-molecule electronics, catalysis, and solar energy conversion.     

Contributions of symmetric and asymmetric normal coordinates to the intervalence electronic absorption and resonance Raman spectra of a strongly coupled p-phenylenediamine radical cation

Resonance Raman spectroscopy, electronic absorption spectroscopy, and the time-dependent theory of spectroscopy are used to analyze the intervalence electron transfer properties of a strongly delocalized class III molecule, the tetraalkyl-p-phenylene diamine radical cation bis(3-oxo-9-azabicyclo[3.3.1]non-9-yl)benzene ((k33)(2)PD(+)). This molecule is a prototypical system for strongly coupled organic intervalence electron transfer spectroscopy. Resonance Raman excitation profiles in resonance with the lowest energy absorption band are measured. The normal modes of vibration that are most strongly coupled to the intervalence transition are identified and assigned by using UB3LYP/6-31G(d) calculations. Excited state distortions are obtained, and the resonance Raman intensities and excitation profiles are calculated by using the time-dependent theory of Raman spectroscopy. The most highly distorted normal modes are all totally symmetric, but intervalence electron transfer absorption spectra are usually interpreted in terms of a model based on coupling between potential surfaces that are displaced along an asymmetric normal coordinate. This model provides a convenient physical picture for the intervalence compound, but it is inadequate for explaining the spectra. The absorption spectrum arising from only the strongly coupled surfaces consists of a single narrow band, in contrast to the broad, vibronically structured experimental spectrum. The electronic absorption spectrum of (k33)(2)PD(+) is calculated by using exactly the same potential surfaces as those used for the Raman calculations. The importance of symmetric normal coordinates, in addition to the asymmetric coordinate, is discussed. The observed vibronic structure is an example of the missing mode effect; the spacing is interpreted in terms of the time-dependent overlaps in the time domain.