The diabatic picture of electron transfer, reaction barriers, and molecular dynamics

A systematic model study quantifying how conical intersection topography modulates photochemical reactions

Despite their important role in photochemistry and expected presence in most polyatomic molecules, conical intersections have been thoroughly characterized in a comparatively small number of systems. Conical intersections can confer molecular photoreactivity or photostability, often with remarkable efficacy, due to their unique structure: at a conical intersection, the adiabatic potential energy surfaces of two or more electronic states are degenerate, enabling ultrafast decay from an excited state without radiative emission, known as nonadiabatic transfer. Furthermore, the precise conical intersection topography determines fundamental properties of photochemical processes, including excited-state decay rate, efficacy, and molecular products that are formed. However, these relationships have yet to be defined comprehensively. In this article, we use an adaptable computational model to investigate a variety of conical intersection topographies, simulate resulting nonadiabatic dynamics, and calculate key photochemical observables. We varied the vibrational mode frequencies to modify conical intersection topography systematically in four primary classes of conical intersections and quantified the resulting rate, total yield, and product yield of nonadiabatic decay. The results reveal that higher vibrational mode frequencies reduce nonadiabatic transfer, but increase the transfer rate and resulting photoproduct formation. These trends can inform progress toward experimental control of photochemical reactions or tuning of molecules' photochemical properties based on conical intersections and their topography.

Uniform potential difference scheme to evaluate effective electronic couplings for superexchange electron transfer in donor-bridge-acceptor systems

This article proposes an ab initio quantum chemical method to evaluate the effective electronic coupling that determines the rate of superexchange electron transfer in donor-bridge-acceptor (D-B-A) systems. The method utilizes the fragment charge difference to define electronic diabatic states and to apply an electrostatic potential in a form of a uniform potential difference that mimics solvation effects on the relative energies of the electronic states. The two-state generalized Mulliken-Hush method is used to obtain the effective electronic coupling as the nondiagonal element of the effective Hamiltonian that is derived based on the Green's function approach and the quasi-degenerate perturbation theory. A theoretical basis is provided for the dependence of the calculated effective electronic coupling on the applied potential and for how to find the optimal potential to give the desired effective electronic coupling that coincides with the result of the minimum energy splitting method. The method is applied to typical D-B-A molecules and gives the effective electronic couplings in reasonable agreement with the experimental estimates.

Density-functional tight-binding: basic concepts and applications to molecules and clusters

The scope of this article is to present an overview of the Density Functional based Tight Binding (DFTB) method and its applications. The paper introduces the basics of DFTB and its standard formulation up to second order. It also addresses methodological developments such as third order expansion, inclusion of non-covalent interactions, schemes to solve the self-interaction error, implementation of long-range short-range separation, treatment of excited states via the time-dependent DFTB scheme, inclusion of DFTB in hybrid high-level/low level schemes (DFT/DFTB or DFTB/MM), fragment decomposition of large systems, large scale potential energy landscape exploration with molecular dynamics in ground or excited states, non-adiabatic dynamics. A number of applications are reviewed, focusing on -(i)- the variety of systems that have been studied such as small molecules, large molecules and biomolecules, bare orfunctionalized clusters, supported or embedded systems, and -(ii)- properties and processes, such as vibrational spectroscopy, collisions, fragmentation, thermodynamics or non-adiabatic dynamics. Finally outlines and perspectives are given.

Quasi-Diabatic Scheme for Nonadiabatic On-the-Fly Simulations

We use the quasi-diabatic (QD) propagation scheme to perform on-the-fly nonadiabatic simulations of the photodynamics of ethylene. The QD scheme enables a seamless interface between accurate diabatic-based quantum dynamics approaches and adiabatic electronic structure calculations, explicitly avoiding any efforts to construct global diabatic states or reformulate the diabatic dynamics approach to the adiabatic representation. Using the partial linearized path-integral approach and the symmetrical quasi-classical approach as the diabatic dynamics methods, the QD propagation scheme enables direct nonadiabatic simulation with complete active space self-consistent field on-the-fly electronic structure calculations. The population dynamics obtained from both approaches are in a close agreement with the quantum wavepacket-based method and outperform the widely used trajectory surface-hopping approach. Further analysis of the ethylene photodeactivation pathways demonstrates the correct predictions of competing processes of nonradiative relaxation mechanism through various conical intersections. This work provides the foundation of using accurate diabatic dynamics approaches and on-the-fly adiabatic electronic structure information to perform ab initio nonadiabatic simulation.

Two-Dimensional Analysis of the Diabatic Transition of a General Vectorial Physical Observable Based on Adiabatic-to-Diabatic Transformation

We present a full analysis of the magnitude and orientation of the diabatic transition matrix element of a general vectorial physical observable during the adiabatic-to-diabatic transformation. The diabatic transition is a function of the adiabatic-to-diabatic transformation angle and the two basic vectors of the adiabatic states, which are the off-diagonal matrix element and the difference between the two diagonal matrix elements. To the best of our knowledge, this is the first time that the transformation has been accomplished in a more general two-dimensional scale for a vectorial physical observable. All possible extreme values of a diabatic transition are deduced for systems with different features. By using an approximate diabatic transition dipole, the pilot implementation of the analysis produces an electronic coupling curve nearly identical to that obtained by the generalized Mulliken-Hush method for the testing molecule. Evidently, this complete analysis of a diabatic transition will be very useful in determining the adiabatic-to-diabatic transformation angle by using a physical observable and can also be used to evaluate the quality of various approximations for constructing the diabatic states.

Unified Model for Photophysical and Electro-Optical Properties of Green Fluorescent Proteins

Green fluorescent proteins (GFPs) have become indispensable imaging and optogenetic tools. Their absorption and emission properties can be optimized for specific applications. Currently, no unified framework exists to comprehensively describe these photophysical properties, namely the absorption maxima, emission maxima, Stokes shifts, vibronic progressions, extinction coefficients, Stark tuning rates, and spontaneous emission rates, especially one that includes the effects of the protein environment. In this work, we study the correlations among these properties from systematically tuned GFP environmental mutants and chromophore variants. Correlation plots reveal monotonic trends, suggesting that all these properties are governed by one underlying factor dependent on the chromophore's environment. By treating the anionic GFP chromophore as a mixed-valence compound existing as a superposition of two resonance forms, we argue that this underlying factor is defined as the difference in energy between the two forms, or the driving force, which is tuned by the environment. We then introduce a Marcus-Hush model with the bond length alternation vibrational mode, treating the GFP absorption band as an intervalence charge transfer band. This model explains all of the observed strong correlations among photophysical properties; related subtopics are extensively discussed in the Supporting Information. Finally, we demonstrate the model's predictive power by utilizing the additivity of the driving force. The model described here elucidates the role of the protein environment in modulating the photophysical properties of the chromophore, providing insights and limitations for designing new GFPs with desired phenotypes. We argue that this model should also be generally applicable to both biological and nonbiological polymethine dyes.

Lazy electrons in graphene

Within a tight-binding approximation, we numerically determine the time evolution of graphene electronic states in the presence of classically vibrating nuclei. There is no reliance on the Born-Oppenheimer approximation within the p-orbital tight-binding basis, although our approximation is "atomically adiabatic": the basis p-orbitals are taken to follow nuclear positions. Our calculations show that the strict adiabatic Born-Oppenheimer approximation fails badly. We find that a diabatic (lazy electrons responding weakly to nuclear distortions) Born-Oppenheimer model provides a much more accurate picture and suggests a generalized many-body Bloch orbital-nuclear basis set for describing electron-phonon interactions in graphene.

Toward the evaluation of intersystem crossing rates with variational relativistic methods

The change in electronic state from one spin multiplicity to another, known as intersystem crossing, occurs in molecules via the relativistic phenomenon of spin-orbit coupling. Current means of estimating intersystem crossing rates rely on the perturbative evaluation of spin-orbit coupling effects. This perturbative approach, valid in lighter atoms where spin-orbit coupling is weaker, is expected to break down for heavier elements where relativistic effects become dominant. Methods which incorporate spin-orbit effects variationally, such as the exact-two-component (X2C) method, will be necessary to treat this strong-coupling regime. We present a novel procedure which produces a diabatic basis of spin-pure electronic states coupled by spin-orbit terms, generated from fully variational relativistic calculations. This method is implemented within X2C using time-dependent density-functional theory and is compared to results from a perturbative relativistic study in the weak spin-orbit coupling regime. Additional calculations on a more strongly spin-orbit-coupled [UO₂Cl₄]^2- complex further illustrate the strengths of this method. This procedure will be valuable in the estimation of intersystem crossing rates within strongly spin-coupled species.

Quantum Trajectory Mean-Field Method for Nonadiabatic Dynamics in Photochemistry

The mixed quantum-classical dynamical approaches have been widely used to study nonadiabatic phenomena in photochemistry and photobiology, in which the time evolutions of the electronic and nuclear subsystems are treated based on quantum and classical mechanics, respectively. The key issue is how to deal with coherence and decoherence during the propagation of the two subsystems, which has been the subject of numerous investigations for a few decades. A brief description on Ehrenfest mean-field and surface-hopping (SH) methods is first provided, and then different algorithms for treatment of quantum decoherence are reviewed in the present paper. More attentions were paid to quantum trajectory mean-field (QTMF) method under the picture of quantum measurements, which is able to overcome the overcoherence problem. Furthermore, the combined QTMF and SH algorithm is proposed in the present work, which takes advantages of the QTMF and SH methods. The potential to extend the applicability of the QTMF method was briefly discussed, such as the generalization to other type of nonadiabatic transitions, the combination with multiscale computational models, and possible improvements on its accuracy and efficiency by using machine-learning techniques.

How To Extract Quantitative Information on Electronic Transitions from the Density Functional Theory "Black Box"

Electronic couplings and vertical excitation energies are crucial determinants of charge and excitation energy transfer rates in a broad variety of processes ranging from biological charge transfer to charge transport through inorganic materials, from molecular sensing to intracellular signaling. Density Functional Theory (DFT) is generally used to calculate these critical parameters, but the quality of the results is unpredictable because of the semiempirical nature of the available DFT approaches. This study identifies a small set of fundamental rules that enables accurate DFT computation of electronic couplings and vertical excitation energies in molecular complexes and materials. These rules are applied to predict efficient DFT approaches to coupling calculations. The result is an easy-to-use guide for reliable DFT descriptions of electronic transitions.

Solving the Trivial Crossing Problem While Preserving the Nodal Symmetry of the Wave Function

In an adiabatic mixed quantum-classical simulation, the avoided crossing of weakly coupled eigenstates can lead to unphysical discontinuities in wave function dynamics, otherwise known as the trivial crossing problem. A standard solution to the trivial crossing problem eliminates spatial discontinuities in wave function dynamics by imposing changes to the eigenstate of the wave function. In this paper, we show that this solution has the side effect of introducing transient discontinuities in the nodal symmetry of the wave function. We present an alternative solution to the trivial crossing problem that preserves both the spatial and nodal structure of the adiabatic wave function. By considering a model of exciton dynamics on conjugated polymer systems, we show that failure to preserve wave function symmetry yields exciton dynamics that depends unphysically on polymer system size. We demonstrate that our symmetry-preserving solution to the trivial crossing problem yields more realistic dynamics and can thus improve the accuracy of simulations of larger systems that are prone to the trivial crossing problem.

Nitration of Tyrosine Channels Photoenergy through a Conical Intersection in Water

Nitration of tyrosine occurs under oxidative stress in vivo. The product, 3-nitrotyrosine (3NY), has a dramatically decreased quantum yield and can be used as a molecular ruler. In this study, femtosecond transient absorption spectroscopy and quantum calculations were implemented to elucidate the photoinduced relaxation processes of anionic 3NY in water. Upon 400 nm excitation into an excited electronic state with notable charge-transfer (CT) character, a barrierless nitro-twisting motion rapidly (<100 fs) guides the chromophore into an adjacent twisted intramolecular CT state, therein reaching a sloped S₁/S₀ conical intersection on the ∼100 fs time scale. Once in the hot ground state, excess energy is further released through vibrational cooling with biexponential time constants of ∼140 and 680 fs in water. Nitro back-twisting occurs on longer time scales (∼1.1 and 9 ps in water), returning the system to original ground state. Systematic evaluations of excited-state potential energies of anionic 3NY were performed by density functional theory (DFT) and time-dependent DFT calculations, showing that intersystem crossing (ISC) from the first singlet state (S₁) to the first or second triplet state (T₁ or T₂) is unlikely. Inclusion of an explicit water molecule in calculations leads to improved mapping of the excited-state energy ordering of the second singlet state (S₂) and T₂, further diminishing ISC probability from S₁ and favoring an ultrafast internal conversion to S₀. These results provide deep insights into the highly efficient nonradiative decay of anionic 3NY in aqueous solution, with nitro-site-specific information that can help infer the characterization and potential optogenetic control of 3NY in protein environment.

Mechanism of Spin-Exchange Internal Conversion: Practical Proxies for Diabatic and Nonadiabatic Couplings

Spin-exchange internal conversion (SEIC) is a general class of reactions having singlet fission and triplet fusion as particular cases. Based on a charge transfer (CT) mediated mechanism and analytical derivation with a model Hamiltonian, we propose proxies for estimating the coupling strength in both diabatic and adiabatic pictures for general SEIC reactions. In the diabatic picture, we demonstrated the existence of a bilinear relationship between the coupling strength and molecular orbital overlap, which provides a practical way to predict diabatic couplings. In the adiabatic picture, we showed that nonadiabatic couplings can be approximated by simple functions of the wave function CT coefficients. These approaches were verified through the investigation of singlet oxygen photosensitization, where both ¹Δ_g and ¹Σ_g oxygen states can be competitively generated by a triplet fusion reaction. The interplay between the CT-mediated mechanism, the spatial factors of the bimolecular complex, and the electronic structure of the oxygen molecule during the reaction explains the curiously small coupling to the ¹Σ_g state along specific incidence directions. The results from both the diabatic and adiabatic pictures provide a comprehensive understanding of the reaction mechanism, which applies to general SEIC problems.

Triplet Pair States in Singlet Fission

This account aims at providing an understanding of singlet fission, i.e., the photophysical process of a singlet state ( S₁) splitting into two triplet states (2 × T₁) in molecular chromophores. Since its discovery 50 years ago, the field of singlet fission has enjoyed rapid expansion in the past 8 years. However, there have been lingering confusion and debates on the nature of the all-important triplet pair intermediate states and the definition of singlet fission rates. Here we clarify the confusion from both theoretical and experimental perspectives. We distinguish the triplet pair state that maintains electronic coherence between the two constituent triplets, ¹(TT), from one which does not, ¹(T···T). Only the rate of formation of ¹(T···T) is defined as that of singlet fission. We present distinct experimental evidence for ¹(TT), whose formation may occur via incoherent and/or vibronic coherent mechanisms. We discuss the challenges in treating singlet fission beyond the dimer approximation, in understanding the often neglected roles of delocalization on singlet fission rates, and in realizing the much lauded goal of increasing solar energy conversion efficiencies with singlet fission chromophores.

Quantum interferences among Dexter energy transfer pathways

We explore Dexter coupling pathway interferences in non-covalent assemblies, employing a method that enables the assessment of Dexter coupling pathway strengths, interferences, and their physical origins in the context of one-particle and two-particle (i.e., coulombic) operators.

Corresponding Orbitals Derived from Periodic Bloch States for Electron Transfer Calculations of Transition Metal Oxides

An approach for modeling electron transfer in solids and at surfaces of iron-(oxyhydr)oxides and other redox active solids has been developed for electronic structure methods (i.e., plane-wave density functional theory) capable of performing calculations with periodic cells and large system sizes efficiently while at the same time being accurate enough to be used in the estimation of the electron-transfer coupling matrix element, V _AB, and the electron transfer transmission factor, κ_el. This method is an extension of the valence bond theory electron transfer method for molecules and clusters implemented by Dupuis and others and used extensively by Rosso and co-workers in which scaled corresponding orbitals derived from the Bloch states are used to calculate the off-diagonal matrix elements H _AB and S _AB. A key development of the present work is the formulation of algorithms to improve the accuracy of the integration of the exact exchange integral in periodic boundary conditions. This method is demonstrated on model systems for electron small polaron transfer in iron-(oxyhydr)oxides, including bare Fe²⁺-Fe³⁺ ions, and in [Fe³⁺(OH₂)₂ (OH^-)₂)] _nⁿ⁺ chains representing the common edge-sharing Fe octahedral motif in these materials.

Signatures of Herzberg-Teller coupling in three-dimensional electronic spectroscopy

The coupling between electronic and nuclear variables is a key consideration in molecular dynamics and spectroscopy. However, simulations that include detailed vibronic coupling terms are challenging to perform, and thus a variety of approximations can be used to model and interpret experimental results. Recent work shows that these simplified models can be inadequate. It is therefore important to understand spectroscopic signals that can identify failures of those approximations. Here we use an extended response-function method to simulate coherent three-dimensional electronic spectroscopy (3D ES) and study the sensitivity of this method to the breakdown of the Franck-Condon approximation. The simulations include a coordinate-dependent transition dipole operator that produces nodes, phase shifts, and peak patterns in 3D ES that can be used to identify Herzberg-Teller coupling. Guided by the simulation results, we interpret measurements on a molecular aggregate.

On the Computational Characterization of Charge-Transfer Effects in Noncovalently Bound Molecular Complexes

Charge-transfer (CT) is an important binding force in the formation of intermolecular complexes, and there have been a variety of theoretical models proposed to quantify this effect. These approaches, which typically rely on a definition of a "CT-free" state based on a partition of the system, sometimes yield significantly different results for a given intermolecular complex. Two widely used definitions of the "CT-free" state, the absolutely localized molecular orbitals (ALMO) method (where only on-fragment orbital mixings are permitted) and the constrained density functional theory (CDFT) approach (where fragment electron populations are fixed), are carefully examined in this work. Natural bond orbital (NBO) and the regularized symmetry-adapted perturbation theory (SAPT) are also briefly considered. Results for the ALMO and CDFT definitions of CT are compared on a broad range of model systems, including hydrogen-bonding systems, borane complexes, metal-carbonyl complexes, and complexes formed by water and metal cations. For most of these systems, CDFT yields a much smaller equilibrium CT energy compared to that given by the ALMO-based definition. This is mainly because the CDFT population constraint does not fully inhibit CT, which means that the CDFT "CT-free" state is in fact CT-contaminated. Examples of this contamination include (i) matching forward and backward donation (e.g., formic acid dimer) and (ii) unidirectional CT without changing fragment populations. The magnitude of the latter effect is quantified in systems such as the water dimer by employing a 3-space density constraint in addition to the orbital constraint. Furthermore, by means of the adiabatic EDA, it is shown that several observable effects of CT, such as the "pyramidalization" of the planar BH₃ molecule upon the complexation with Lewis bases, already appear on the "CT-free" CDFT surface. These results reveal the essential distinctions between the ALMO and CDFT definitions of CT and suggest that the former is more consistent with accepted understanding of the role of CT in intermolecular binding.

Fragment-orbital tunneling currents and electronic couplings for analysis of molecular charge-transfer systems

In theoretical charge-transfer research, calculation of the electronic coupling element is crucial for examining the degree of the electronic donor-acceptor interaction. The tunneling current (TC), representing the magnitudes and directions of electron flow, provides a way of evaluating electronic couplings, along with the ability of visualizing how electrons flow in systems. Here, we applied the TC theory to π-conjugated organic dimer systems, in the form of our fragment-orbital tunneling current (FOTC) method, which uses the frontier molecular-orbitals of system fragments as diabatic states. For a comprehensive test of FOTC, we assessed how reasonable the computed electronic couplings and the corresponding TC densities are for the hole- and electron-transfer databases HAB11 and HAB7. FOTC gave 12.5% mean relative unsigned error with regard to the high-level ab initio reference. The shown performance is comparable with that of fragment-orbital density functional theory, which gave the same error by 20.6% or 13.9% depending on the formulation. In the test of a set of nucleobase π stacks, we showed that the original TC expression is also applicable to nondegenerate cases under the condition that the overlap between the charge distributions of diabatic states is small enough to offset the energy difference. Lastly, we carried out visual analysis on the FOTC densities of thiophene dimers with different intermolecular alignments. The result depicts an intimate topological connection between the system geometry and electron flow. Our work provides quantitative and qualitative grounds for FOTC, showing it to be a versatile tool in characterization of molecular charge-transfer systems.

Spatial and orientational dependence of electron transfer parameters in aggregates of iridium-containing host materials for OLEDs: coupling constrained density functional theory with molecular dynamics

The integration between molecular dynamics and constrained density functional theory allows to evaluate charge transport parameters in bulk organic semiconductors.

Electronic spectrum and characterization of diabatic potential energy surfaces for thiophenol

We present an accurate simulation of the UV spectrum and a diabatization of three singlet potential surfaces along four coordinates.