On the construction of diabatic and adiabatic potential energy surfaces based on ab initio valence bond theory

Hybrid Density Functional Valence Bond Method with Multistate Treatment

Recently, a hybrid density functional valence bond (VB) method, λ-DFVB(U), has been proposed and shown to give accuracy that is comparable to that of CASPT2 in calculations of atomization energies, atomic excitation energies, and reaction barriers, while its computational cost is approximately the same as the valence bond self-consistent-field (VBSCF) method. However, the interaction between electronic states is not included in λ-DFVB(U) since the last step of λ-DFVB(U) is not a diagonalization of the Hamiltonian matrix on the electronic state basis. Therefore, λ-DFVB(U) gives the wrong topology of the potential energy surfaces (PESs) near the conical intersection region. In the present paper, we propose a novel hybrid density functional VB method with multistate treatment, named λ-DFVB(MS), in which an effective Hamiltonian matrix is constructed on the basis of the diabatic states obtained by the valence-bond-based compression approach for the diabatization scheme, and the interaction between electronic states can be included through the diagonalization of the effective Hamiltonian matrix. Test calculations show that λ-DFVB(MS) gives the correct topology of the PESs near the conical intersection region. We also show that the VBSCF wave function with selected VB structures can be applied as a reference in λ-DFVB(MS).

Nonorthogonal Multireference Wave Function Description of Triplet-Triplet Energy Transfer Couplings

In this study, the use of self-consistent field quasi-diabats is investigated for calculation of triplet energy transfer diabatic coupling elements. It is proposed that self-consistent field quasi-diabats are particularly useful for studying energy transfer (EnT) processes because orbital relaxation in response to changes in electron configuration is implicitly built into the model. The conceptual model that is developed allows for the simultaneous evaluation of direct and charge-transfer mechanisms to establish the importance of the different possible EnT mechanisms. The method's performance is evaluated using two model systems: the ethylene dimer and ethylene with the methaniminium cation. While states that mediate the charge-transfer mechanism were found to be higher in energy than the states involved in the direct mechanism, the coupling elements that control the kinetics were found to be significantly larger in the charge-transfer mechanism. Subsequently, we discuss the advantage of the approach in the context of practical difficulties with the use of established approaches.

Leveling the Mountain Range of Excited-State Benchmarking through Multistate Density Functional Theory

The performance of multistate density functional theory (MSDFT) with nonorthogonal state interaction (NOSI) is assessed for 100 vertical excitation energies against the theoretical best estimates extracted to the full configuration interaction accuracy on the database developed by Loos et al. in 2018 (Loos2018). Two optimization techniques, namely, block-localized excitation and target state optimization, are examined along with two ways of estimating the transition density functional (TDF) for the correlation energy of the Hamiltonian matrix density functional. The results from the two optimization methods are similar. It was found that MSDFT-NOSI using the spin-multiplet degeneracy constraint for the TDF of spin-coupling interaction, along with the M06-2X functional, yields a root-mean-square error (RMSE) of 0.22 eV, which performs noticeably better than time-dependent density functional theory (DFT) at an RMSE of 0.43 eV using the same functional and basis set on the Loos2018 database. In comparison with wave function theory, NOSI has smaller errors than CIS(D∞), LR-CC2, and ADC(3) all of which have an RMSE of 0.28 eV, but somewhat greater than STEOM-CCSD (RMSE of 0.14 eV) and LR-CCSD (RMSE of 0.11 eV) wave function methods. In comparison with Kohn-Sham (KS) DFT calculations, the multistate DFT approach has little double counting of correlation. Importantly, there is no noticeable difference in the performance of MSDFT-NOSI on the valence, Rydberg, singlet, triplet, and double-excitation states. Although the use of another hybrid functional PBE0 leads to a greater RMSE of 0.36 eV, the deviation is systematic with a linear regression slope of 0.994 against the results with M06-2X. The present benchmark reveals that density functional approximations developed for KS-DFT for the ground state with a noninteracting reference may be adopted in MSDFT calculations in which the state interaction is key.

Multistate Energy Decomposition Analysis of Molecular Excited States

A multistate energy decomposition analysis (MS-EDA) method is described to dissect the energy components in molecular complexes in excited states. In MS-EDA, the total binding energy of an excimer or an exciplex is partitioned into a ground-state term, called local interaction energy, and excited-state contributions that include exciton excitation energy, superexchange stabilization, and orbital and configuration-state delocalization. An important feature of MS-EDA is that key intermediate states associated with different energy terms can be variationally optimized, providing quantitative insights into widely used physical concepts such as exciton delocalization and superexchange charge-transfer effects in excited states. By introducing structure-weighted adiabatic excitation energy as the minimum photoexcitation energy needed to produce an excited-state complex, the binding energy of an exciplex and excimer can be defined. On the basis of the nature of intermolecular forces through MS-EDA analysis, it was found that molecular complexes in the excited states can be classified into three main categories, including (1) encounter excited-state complex, (2) charge-transfer exciplex, and (3) intimate excimer or exciplex. The illustrative examples in this Perspective highlight the interplay of local excitation polarization, exciton resonance, and superexchange effects in molecular excited states. It is hoped that MS-EDA can be a useful tool for understanding photochemical and photobiological processes.

Excimer Energies

A multistate energy decomposition analysis (MS-EDA) method is introduced for excimers using density functional theory. Although EDA has been widely applied to intermolecular interactions in the ground state, few methods are currently available for excited-state complexes. Here, the total energy of an excimer state is separated into exciton excitation energy ΔE_Ex(|Ψ_X·Ψ_Y⟩*), resulting from the state interaction between locally excited monomer states |Ψ_X^*·Ψ_Y⟩ and |Ψ_X·Ψ_Y^*⟩ , a superexchange stabilization energy ΔE_SE, originating from the mutual charge transfer between two monomers |Ψ_X⁺·Ψ_Y⟩ and |Ψ_X^-·Ψ_Y⁺⟩ , and an orbital-and-configuration delocalization term ΔE_OCD due to the expansion of configuration space and block-localized orbitals to the fully delocalized dimer system. Although there is no net charge transfer in symmetric excimer cases, the resonance of charge-transfer states is critical to stabilizing the excimer. The monomer localized excited and charge-transfer states are variationally optimized, forming a minimal active space for nonorthogonal state interaction (NOSI) calculations in multistate density functional theory to yield the intermediate states for energy analysis. The present MS-EDA method focuses on properties unique to excited states, providing insights into exciton coupling, superexchange and delocalization energies. MS-EDA is illustrated on the acetone and pentacene excimer systems; three configurations of the latter case are examined, including the optimized excimer, a stacked configuration of two pentacene molecules and the fishbone orientation. It is found that excited-state energy splitting is strongly dependent on the relative energies of the monomer excited states and the phase-matching of the monomer wave functions.

Target State Optimized Density Functional Theory for Electronic Excited and Diabatic States

A flexible self-consistent field method, called target state optimization (TSO), is presented for exploring electronic excited configurations and localized diabatic states. The key idea is to partition molecular orbitals into different subspaces according to the excitation or localization pattern for a target state. Because of the orbital-subspace constraint, orbitals belonging to different subspaces do not mix. Furthermore, the determinant wave function for such excited or diabatic configurations can be variationally optimized as a ground state procedure, unlike conventional ΔSCF methods, without the possibility of collapsing back to the ground state or other lower-energy configurations. The TSO method can be applied both in Hartree-Fock theory and in Kohn-Sham density functional theory (DFT). The density projection procedure and the working equations for implementing the TSO method are described along with several illustrative applications. For valence excited states of organic compounds, it was found that the computed excitation energies from TSO-DFT and time-dependent density functional theory (TD-DFT) are of similar quality with average errors of 0.5 and 0.4 eV, respectively. For core excitation, doubly excited states and charge-transfer states, the performance of TSO-DFT is clearly superior to that from conventional TD-DFT calculations. It is shown that variationally optimized charge-localized diabatic states can be defined using TSO-DFT in energy decomposition analysis to gain both qualitative and quantitative insights on intermolecular interactions. Alternatively, the variational diabatic states may be used in molecular dynamics simulation of charge transfer processes. The TSO method can also be used to define basis states in multistate density functional theory for excited states through nonorthogonal state interaction calculations. The software implementing TSO-DFT can be accessed from the authors.

Minimal Active Space for Diradicals Using Multistate Density Functional Theory

This work explores the electronic structure as well as the reactivity of singlet diradicals, making use of multistate density functional theory (MSDFT). In particular, we show that a minimal active space of two electrons in two orbitals is adequate to treat the relative energies of the singlet and triplet adiabatic ground state as well as the first singlet excited state in many cases. This is plausible because dynamic correlation is included in the first place in the optimization of orbitals in each determinant state via block-localized Kohn–Sham density functional theory. In addition, molecular fragment, i.e., block-localized Kohn–Sham orbitals, are optimized separately for each determinant, providing a variational diabatic representation of valence bond-like states, which are subsequently used in nonorthogonal state interactions (NOSIs). The computational procedure and its performance are illustrated on some prototypical diradical species. It is shown that NOSI calculations in MSDFT can be used to model bond dissociation and hydrogen-atom transfer reactions, employing a minimal number of configuration state functions as the basis states. For p- and s-types of diradicals, the closed-shell diradicals are found to be more reactive than the open-shell ones due to a larger diabatic coupling with the final product state. Such a diabatic representation may be useful to define reaction coordinates for electron transfer, proton transfer and coupled electron and proton transfer reactions in condensed-phase simulations.

Diabatic States of Molecules

Quantitative simulations of electronically nonadiabatic molecular processes require both accurate dynamics algorithms and accurate electronic structure information. Direct semiclassical nonadiabatic dynamics is expensive due to the high cost of electronic structure calculations, and hence it is limited to small systems, limited ensemble averaging, ultrafast processes, and/or electronic structure methods that are only semiquantitatively accurate. The cost of dynamics calculations can be made manageable if analytic fits are made to the electronic structure data, and such fits are most conveniently carried out in a diabatic representation because the surfaces are smooth and the couplings between states are smooth scalar functions. Diabatic representations, unlike the adiabatic ones produced by most electronic structure methods, are not unique, and finding suitable diabatic representations often involves time-consuming nonsystematic diabatization steps. The biggest drawback of using diabatic bases is that it can require large amounts of effort to perform a globally consistent diabatization, and one of our goals has been to develop methods to do this efficiently and automatically. In this Feature Article, we introduce the mathematical framework of diabatic representations, and we discuss diabatization methods, including adiabatic-to-diabatic transformations and recent progress toward the goal of automatization.

Unified Formulation of Phase Space Mapping Approaches for Nonadiabatic Quantum Dynamics

Nonadiabatic dynamical processes are one of the most important quantum mechanical phenomena in chemical, materials, biological, and environmental molecular systems, where the coupling between different electronic states is either inherent in the molecular structure or induced by the (intense) external field. The curse of dimensionality indicates the intractable exponential scaling of calculation effort with system size and restricts the implementation of "numerically exact" approaches for realistic large systems. The phase space formulation of quantum mechanics offers an important theoretical framework for constructing practical approximate trajectory-based methods for quantum dynamics. This Account reviews our recent progress in phase space mapping theory: a unified framework for constructing the mapping Hamiltonian on phase space for coupled F-state systems where the renowned Meyer-Miller Hamiltonian model is a special case, a general phase space formulation of quantum mechanics for nonadiabatic systems where the electronic degrees of freedom are mapped onto constraint space and the nuclear degrees of freedom are mapped onto infinite space, and an isomorphism between the mapping phase space approach for nonadiabatic systems and that for nonequilibrium electron transport processes. While the zero-point-energy parameter is conventionally assumed to be positive, we show that the constraint implied in the conventional Meyer-Miller mapping Hamiltonian requires that such a parameter can be negative as well and lies in (-1/F, +∞) for each electronic degree of freedom. More importantly, the zero-point-energy parameter should be interpreted as a special case of a commutator matrix in the comprehensive phase space mapping Hamiltonian for nonadiabatic systems. From the rigorous formulation of mapping phase space, we propose approximate but practical trajectory-based nonadiabatic dynamics methods. The applications to both gas phase and condensed phase problems include the spin-boson model for condensed phase dissipative two-state systems, the three-state photodissociation models, the seven-site model of the Fenna-Matthews-Olson monomer in photosynthesis of green sulfur bacteria, the strongly coupled molecular/atomic matter-optical cavity systems designed for controlling and manipulating chemical dynamical processes, and the Landauer model for a quantum dot state coupled with two electrodes. In these applications the overall performance of our phase space mapping dynamics approach is superior to two prevailing trajectory-based methods, Ehrenfest dynamics and fewest switches surface hopping.

Graph-Theory-Based Molecular Fragmentation for Efficient and Accurate Potential Surface Calculations in Multiple Dimensions

We present a multitopology molecular fragmentation approach, based on graph theory, to calculate multidimensional potential energy surfaces in agreement with post-Hartree-Fock levels of theory but at the density functional theory cost. A molecular assembly is coarse-grained into a set of graph-theoretic nodes that are then connected with edges to represent a collection of locally interacting subsystems up to an arbitrary order. Each of the subsystems is treated at two levels of electronic structure theory, the result being used to construct many-body expansions that are embedded within an ONIOM scheme. These expansions converge rapidly with the many-body order (or graphical rank) of subsystems and capture many-body interactions accurately and efficiently. However, multiple graphs, and hence multiple fragmentation topologies, may be defined in molecular configuration space that may arise during conformational sampling or from reactive, bond breaking and bond formation, events. Obtaining the resultant potential surfaces is an exponential scaling proposition, given the number of electronic structure computations needed. We utilize a family of graph-theoretic representations within a variational scheme to obtain multidimensional potential surfaces at a reduced cost. The fast convergence of the graph-theoretic expansion with increasing order of many-body interactions alleviates the exponential scaling cost for computing potential surfaces, with the need to only use molecular fragments that contain a fewer number of quantum nuclear degrees of freedom compared to the full system. This is because the dimensionality of the conformational space sampled by the fragment subsystems is much smaller than the full molecular configurational space. Additionally, we also introduce a multidimensional clustering algorithm, based on physically defined criteria, to reduce the number of energy calculations by orders of magnitude. The molecular systems benchmarked include coupled proton motion in protonated water wires. The potential energy surfaces and multidimensional nuclear eigenstates obtained are shown to be in very good agreement with those from explicit post-Hartree-Fock calculations that become prohibitive as the number of quantum nuclear dimensions grows. The developments here provide a rigorous and efficient alternative to this important chemical physics problem.

Extended Mulliken-Hush Method with Applications to the Theoretical Study of Electron Transfer

A novel adiabatic-to-diabatic (ATD) transformation strategy, namely, the extended Mulliken-Hush (XMH) method, is proposed to evaluate diabatic properties including electronic couplings, potential energy surfaces, and their crossings. The XMH method is developed by adopting our recently proposed ATD transformation formula of a general vectorial physical observable, in which a useful ATD transformation is further determined by using an auxiliary dipole between localized frontier orbitals as a simple approximation of the diabatic transition dipole. The XMH method is simple and practical that provides a flexible way to construct diabatic states. To some extent, it can be regarded as an extension of the generalized Mulliken-Hush (GMH) method since the latter takes a stronger approximation, in which the diabatic transition dipole is assumed to be vanishing. Test calculations on the HeH₂⁺ system show that the electronic couplings predicted by the XMH method are closer to the ones calculated by the valence bond block-diagonalization approach than the GMH ones since the XMH method takes into account both the magnitude and direction of the diabatic transition dipole, which is consistent with the properties of this molecule. In the study of electron transfer in the two kinds of donor-bridge-acceptor systems, the XMH method maintains the simplicity of the GMH method and gives reasonable results even when the latter fails, wherein the diabatic transition dipole is nearly perpendicular to the difference of the initial and final adiabatic dipoles. More importantly, the XMH method can be easily combined with high-level electronic structure methods, in which the properties of the ground and excited states may be more accurately calculated, and hence, one may expect that further development of the XMH method would result in a general computational model for studying electron transfer reactions.

Dynamic-then-Static Approach for Core Excitations of Open-Shell Molecules

Delta self-consistent-field methods are widely used in studies of electronically excited states. However, the nonaufbau determinants are generally spin-contaminated. Here, we describe a general approach for spin-coupling interactions of open-shell molecules, making use of multistate density functional theory (MSDFT). In particular, the effective exchange integrals that determine spin coupling are obtained by enforcing the multiplet degeneracy of the S+1 state in the M_S = S manifold. Consequently, they are consistent with the energy of the high-spin state that is adequately treated by Kohn-Sham density functional theory (DFT) and, thereby, free of double counting of correlation. The method was applied to core excitations of open-shell molecules and compared with those by spin-adapted time-dependent DFT. An excellent agreement with experiment was found employing the BLYP functional and aug-cc-pCVQZ basis set. Overall, MSDFT provides an effective combination of the strengths of DFT and wave function theory to achieve efficiency and accuracy.

Exact-two-component block-localized wave function: A simple scheme for the automatic computation of relativistic ΔSCF

Block-localized wave function is a useful method for optimizing constrained determinants. In this article, we extend the generalized block-localized wave function technique to a relativistic two-component framework. Optimization of excited state determinants for two-component wave functions presents a unique challenge because the excited state manifold is often quite dense with degenerate states. Furthermore, we test the degree to which certain symmetries result naturally from the ΔSCF optimization such as time-reversal symmetry and symmetry with respect to the total angular momentum operator on a series of atomic systems. Variational optimizations may often break the symmetry in order to lower the overall energy, just as unrestricted Hartree-Fock breaks spin symmetry. Overall, we demonstrate that time-reversal symmetry is roughly maintained when using Hartree-Fock, but less so when using Kohn-Sham density functional theory. Additionally, maintaining total angular momentum symmetry appears to be system dependent and not guaranteed. Finally, we were able to trace the breaking of total angular momentum symmetry to the relaxation of core electrons.

Electronic Couplings for Photoinduced Charge Transfer and Excitation Energy Transfer Based on Fragment Particle-Hole Densities

A new scheme is proposed to calculate the electronic couplings for photoinduced charge transfer and excitation energy transfer for both singlet and triplet states. In this scheme, the locally excited and charge-transfer states are constructed from the adiabatic ones by maximally localizing the particle (i.e., electron) and hole densities in terms of predefined molecular fragments. The construction process, after which the electronic couplings are directly obtained, is highly efficient and can be combined with various kinds of preliminary electronic structure calculations as long as the adiabatic excitation energies and transition densities are available. The method also applies to the systems with multiple charge or excitation centers. Its validity is demonstrated by the applications to the 6,13-dichloropentacene dimer and tetramer and the C₆₀-Zn porphyrin dyad. The results reveal that the environment has a strong impact on the electronic couplings and can even enlarge those for long-range charge transfer by several orders of magnitude.

Block-Localized Excitation for Excimer Complex and Diabatic Coupling

We describe a block-localized excitation (BLE) method to carry out constrained optimization of block-localized orbitals for constructing valence bond-like, diabatic excited configurations using multistate density functional theory (MSDFT). The method is an extension of the previous block-localized wave function method through a fragment-based ΔSCF approach to optimize excited determinants within a molecular complex. In BLE, both the number of electrons and the electronic spin of different fragments in a whole system can be constrained, whereas electrostatic, exchange, and polarization interactions among different blocks can be fully taken into account of. To avoid optimization collapse to unwanted states, a ΔSCF projection scheme and a maximum overlap of wave function approach have been presented. The method is illustrated by the excimer complex of two naphthalene molecules. With a minimum of eight spin-adapted configurational state functions, it was found that the inversion of L_a^- and L_b^- states near the optimal structure of the excimer complex is correctly produced, which is in quantitative agreement with DMRG-CASPT2 calculations and experiments. Trends in the computed transfer integrals associated with excited-state energy transfer both in the singlet and triplet states are discussed. The results suggest that MSDFT may be used as an efficient approach to treat intermolecular interactions in excited states with a minimal active space (MAS) for interpretation of the results and for dynamic simulations, although the selection of a small active space is often system dependent.

Generalization of Block-Localized Wave Function for Constrained Optimization of Excited Determinants

The block-localized wave function method is useful to provide insights on chemical bonding and intermolecular interactions through energy decomposition analysis. The method relies on block localization of molecular orbitals (MOs) by constraining the orbitals to basis functions within given blocks. Here, a generalized block-localized orbital (GBLO) method is described to allow both physically localized and delocalized MOs to be constrained in orbital-block definitions. Consequently, GBLO optimization can be conveniently tailored by imposing specific constraints. The GBLO method is illustrated by three examples: (1) constrained polarization response orbitals through dipole and quadrupole perturbation in a water dimer complex, (2) the ground and first excited-state potential energy curves of ethene about its C-C bond rotation, and (3) excitation energies of double electron excited states. Multistate density functional theory is used to determine the energies of the adiabatic ground and excited states using a minimal active space (MAS) comprising specifically charge-constrained and excited determinant configurations that are variationally optimized by the GBLO method. We find that the GBLO expansion that includes delocalized MOs in configurational blocks significantly reduces computational errors in comparison with physical block localization, and the computed ground- and excited-state energies are in good accordance with experiments and results obtained from multireference configuration interaction calculations.

A Novel Valence-Bond-Based Automatic Diabatization Method by Compression

A novel valence-bond-based automatic diabatization method by compression, called valence-bond-based compression approach for dibatization (VBCAD), is presented in this Letter. It is a "black-box" type method that provides an automatic diabatization from a classical valence bond (VB) perspective. In VBCAD, a model space projection is performed by an eigenvalue decomposition algorithm followed by dimensional reduction based on a sequence of Householder transformations. Our diabaticity criterion is implemented in a way that maximizes the diversity of VB structure weights between different diabatic states. Owing to the rigorous Householder transformations employed in this entire procedure, the invariance of the target eigensubspace is preserved. This is illustrated on two prototypical examples.

Coupled electron and proton transfer in the piperidine drug metabolism pathway by the active species of cytochromes P450

KS-DFT and MSDFT studies reveal a novel CEPT step that triggers ring contraction of piperidines by P450.

Diabatic States at Construction (DAC) through Generalized Singular Value Decomposition

A procedure, called generalized diabatic-at-construction (GDAC), is presented to transform adiabatic potential energy surfaces into a diabatic representation by generalized singular value decomposition. First, we use a set of localized, valence bond-like configuration state functions, called DAC, as the basis states. Then, the adiabatic ground and relevant excited states are determined using multistate density functional theory (MSDFT). GDAC differs in the opposite direction from traditional approaches based on adiabatic-to-diabatic transformation with certain property restraints. The method is illustrated with applications to a model first-order bond dissociation reaction of CH₃OCH₂Cl polarized by a solvent molecule, the ground- and first-excited-state potential energy surfaces near the minimum conical intersection for the ammonia dimer photodissociation, and the multiple avoided curve crossings in the dissociation of lithium hydride. The GDAC diabatization method may be useful for defining charge-localized states in studies of electron transfer and proton-coupled electron transfer reactions in proteins.

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.

Direct diabatization based on nonadiabatic couplings: the N/D method

We propose a new diabatization method that is direct, orbital-free, and adiabatic-equivalent based on directly calculated nonadiabatic couplings of states and the adiabatic energy gradients.

Spin-Multiplet Components and Energy Splittings by Multistate Density Functional Theory

Kohn-Sham density functional theory has been tremendously successful in chemistry and physics. Yet, it is unable to describe the energy degeneracy of spin-multiplet components with any approximate functional. This work features two contributions. (1) We present a multistate density functional theory (MSDFT) to represent spin-multiplet components and to determine multiplet energies. MSDFT is a hybrid approach, taking advantage of both wave function theory and density functional theory. Thus, the wave functions, electron densities and energy density-functionals for ground and excited states and for different components are treated on the same footing. The method is illustrated on valence excitations of atoms and molecules. (2) Importantly, a key result is that for cases in which the high-spin components can be determined separately by Kohn-Sham density functional theory, the transition density functional in MSDFT (which describes electronic coupling) can be defined rigorously. The numerical results may be explored to design and optimize transition density functionals for configuration coupling in multiconfigurational DFT.

Diabatic-At-Construction Method for Diabatic and Adiabatic Ground and Excited States Based on Multistate Density Functional Theory

We describe a diabatic-at-construction (DAC) strategy for defining diabatic states to determine the adiabatic ground and excited electronic states and their potential energy surfaces using the multistate density functional theory (MSDFT). The DAC approach differs in two fundamental ways from the adiabatic-to-diabatic (ATD) procedures that transform a set of preselected adiabatic electronic states to a new representation. (1) The DAC states are defined in the first computation step to form an active space, whose configuration interaction produces the adiabatic ground and excited states in the second step of MSDFT. Thus, they do not result from a similarity transformation of the adiabatic states as in the ATD procedure; they are the basis for producing the adiabatic states. The appropriateness and completeness of the DAC active space can be validated by comparison with experimental observables of the ground and excited states. (2) The DAC diabatic states are defined using the valence bond characters of the asymptotic dissociation limits of the adiabatic states of interest, and they are strictly maintained at all molecular geometries. Consequently, DAC diabatic states have specific and well-defined physical and chemical meanings that can be used for understanding the nature of the adiabatic states and their energetic components. Here we present results for the four lowest singlet states of LiH and compare them to a well-tested ATD diabatization method, namely the 3-fold way; the comparison reveals both similarities and differences between the ATD diabatic states and the orthogonalized DAC diabatic states. Furthermore, MSDFT can provide a quantitative description of the ground and excited states for LiH with multiple strongly and weakly avoided curve crossings spanning over 10 Å of interatomic separation.

Potential energy construction in the diabatic picture for quantum mechanical rate constants of intermolecular proton transfer

We propose a simple method for potential construction in the diabatic picture and the estimation of thermal rate constants for intermolecular proton transfer reactions using quantum dynamics simulations carried out on the constructed potentials.

Beyond Kohn-Sham Approximation: Hybrid Multistate Wave Function and Density Functional Theory

A multistate density functional theory (MSDFT) is presented in which the energies and densities for the ground and excited states are treated on the same footing using multiconfigurational approaches. The method can be applied to systems with strong correlation and to correctly describe the dimensionality of the conical intersections between strongly coupled dissociative potential energy surfaces. A dynamic-then-static framework for treating electron correlation is developed to first incorporate dynamic correlation into contracted state functions through block-localized Kohn-Sham density functional theory (KSDFT), followed by diagonalization of the effective Hamiltonian to include static correlation. MSDFT can be regarded as a hybrid of wave function and density functional theory. The method is built on and makes use of the current approximate density functional developed in KSDFT, yet it retains its computational efficiency to treat strongly correlated systems that are problematic for KSDFT but too large for accurate WFT. The results presented in this work show that MSDFT can be applied to photochemical processes involving conical intersections.

Multistate Density Functional Theory for Effective Diabatic Electronic Coupling

Multistate density functional theory (MSDFT) is presented to estimate the effective transfer integral associated with electron and hole transfer reactions. In this approach, the charge-localized diabatic states are defined by block localization of Kohn-Sham orbitals, which constrain the electron density for each diabatic state in orbital space. This differs from the procedure used in constrained density functional theory that partitions the density within specific spatial regions. For a series of model systems, the computed transfer integrals are consistent with experimental data and show the expected exponential attenuation with the donor-acceptor separation. The present method can be used to model charge transfer reactions including processes involving coupled electron and proton transfer.

Electron transfer within a reaction path model calibrated by constrained DFT calculations: application to mixed-valence organic compounds

The quantum dynamics of electron transfer in mixed-valence organic compounds is investigated using a reaction path model calibrated by constrained density functional theory (cDFT). Constrained DFT is used to define diabatic states relevant for describing the electron transfer, to obtain equilibrium structures for each of these states and to estimate the electronic coupling between them. The harmonic analysis at the diabatic minima yields normal modes forming the dissipative bath coupled to the electronic states. In order to decrease the system-bath coupling, an effective one dimensional vibronic Hamiltonian is constructed by partitioning the modes into a linear reaction path which connects both equilibrium positions and a set of secondary vibrational modes, coupled to this reaction coordinate. Using this vibronic model Hamiltonian, dissipative quantum dynamics is carried out using Redfield theory, based on a spectral density which is determined from the cDFT results. In a first benchmark case, the model is applied to a series of mixed-valence organic compounds formed by two 1,4-dimethoxy-3-methylphenylene fragments linked by an increasing number of phenylene bridges. This allows us to examine the coherent electron transfer in extreme situations leading to a ground adiabatic state with or without a barrier and therefore to the trapping of the charge or to an easy delocalization.

Generalized X-Pol Theory and Charge Delocalization States

The mixed molecular orbital and valence bond (MOVB) method has been used to generalize the explicit polarization (X-Pol) potential to incorporate charge delocalization resonance effects in the framework of valence bond theory. In the original X-Pol method, a macromolecular system is partitioned into individual fragments or blocks, and the molecular orbitals of the system are strictly localized within each block. Consequently, these block-localized molecular orbitals (BLMOs) are nonorthogonal across different blocks. In the generalized X-Pol (GX-Pol) theory, we construct charge delocalization VB states by expanding the localization space from monomer blocks into pairwise delocalized blocks. Thus, the expansion of the basis space leads to charge delocalization between monomer pairs, and a series of pairwise delocalization states can be constructed. In general, L-body delocalized states can be analogously defined by grouping L monomer blocks into one. The Hartree product wave function for each state can be fully antisymmetrized, which introduces explicitly exchange repulsion among all blocks. The GX-Pol wave function is a linear combination of all L-body charge transfer (valence bond) states, which incorporates charge delocalization and their resonance as well as static correlation effects. The GX-Pol method provides a general and rigorous theory to incorporate charge delocalization explicitly into these fragment-based electronic structural methods for macromolecular systems.

Two Aromatic Rings Coupled a Sulfur-Containing Group to Favor Protein Electron Transfer by Instantaneous Formations of π∴S:π↔π:S∴π or π∴π:S↔π:π∴S Five-Electron Bindings

The cooperative interactions among two aromatic rings with a S-containing group are described, which may participate in electron hole transport in proteins. Ab initio calculations reveal the possibility for the formations of the π∴S:π↔π:S∴π and π∴π:S↔π:π∴S five-electron bindings in the corresponding microsurrounding structures in proteins, both facilitating electron hole transport as efficient relay stations. The relay functionality of these two special structures comes from their low local ionization energies and proper binding energies, which varies with the different aromatic amino acids, S-containing residues, and the arrangements of the same aromatic rings according to the local microsurroundings in proteins.

A New Type of Electron Relay Station in Proteins: Three-Piece S:Π∴S↔S∴Π:S Resonance Structure

A type of relay station for electron transfer in proteins, three-piece five-electron bonding, is introduced in this paper, which is also first proposed here. The ab initio calculations predict the formation of S:Π∴S↔S∴Π:S resonance binding with an aromatic ring located in the middle of two sulfur-containing groups, which may participate in electron-hole transport in proteins. These special structures can lower the local ionization energies to capture electron holes efficiently and may be easily formed and broken because of their proper binding energies. In addition, the UV-vis spectra provide evidence of the formations of the three-piece five-electron binding. The cooperation of three adjacent pieces may be advantage to promote electron transfer a longer distance.

A unified diabatic description for electron transfer reactions, isomerization reactions, proton transfer reactions, and aromaticity

A way is found for describing general chemical reactions using diabatic multi-state and “twin-state” models. (Image adapted with permission from https://www.flickr.com/photos/cybaea/64638988/).

Bond angle variations in XH3 [X = N, P, As, Sb, Bi]: the critical role of Rydberg orbitals exposed using a diabatic state model

The origins of the observed bond angles in XH₃ and XH₃⁺ are explained using high-level calculations and a simple diabatic model.

The Third Dimension of a More O'Ferrall-Jencks Diagram for Hydrogen Atom Transfer in the Isoelectronic Hydrogen Exchange Reactions of (PhX)(2)H(•) with X = O, NH, and CH(2)

A critical element in theoretical characterization of the mechanism of proton-coupled electron transfer (PCET) reactions, including hydrogen atom transfer (HAT), is the formulation of the electron and proton localized diabatic states, based on which a More O'Ferrall-Jencks diagram can be represented to determine the step-wise and concerted nature of the reaction. Although the More O'Ferrall-Jencks diabatic states have often been used empirically to develop theoretical models for PCET reactions, the potential energy surfaces for these states have never been determined directly based on first principles calculations using electronic structure theory. The difficulty is due to a lack of practical method to constrain electron and proton localized diabatic states in wave function or density functional theory calculations. Employing a multistate density functional theory (MSDFT), in which the electron and proton localized diabatic configurations are constructed through block-localization of Kohn-Sham orbitals, we show that distinction between concerted proton-electron transfer (CPET) and HAT, which are not distinguishable experimentally from phenomenological kinetic data, can be made by examining the third dimension of a More O'Ferrall-Jencks diagram that includes both the ground and excited state potential surfaces. In addition, we formulate a pair of effective two-state valence bond models to represent the CPET and HAT mechanisms. We found that the lower energy of the CPET and HAT effective diabatic states at the intersection point can be used as an energetic criterion to distinguish the two mechanisms. In the isoelectronic series of hydrogen exchange reaction in (PhX)(2)H(•), where X = O, NH, and CH(2), there is a continuous transition from a CPET mechanism for the phenoxy radical-phenol pair to a HAT process for benzyl radical and toluene, while the reaction between PhNH(2) and PhNH(•) has a mechanism intermediate of CPET and HAT. The electronically nonadiabatic nature of the CPET mechanism in the phenol system can be attributed to the overlap interactions between the ground and excited state surfaces, resulting in roughly orthogonal minimum energy paths on the adiabatic ground and excited state potential energy surfaces. On the other hand, the minimum energy path on the adiabatic ground state for the HAT mechanism coincides with that on the excited state, producing a large electronic coupling that separates the two surfaces by more than 120 kcal/mol.

Nonorthogonality Problem and Effective Electronic Coupling Calculation: Application to Charge Transfer in π-Stacks Relevant to Biochemistry and Molecular Electronics

A recently proposed method for the calculation of the effective electronic coupling (or charge-transfer integral) in a two-state system is discussed and related to other methods in the literature. The theoretical expression of the coupling is exact within the two-state model and applies to the general case where the charge transfer (CT) process involves nonorthogonal initial and final diabatic (localized) states. In this work, it is shown how this effective electronic coupling is also the one to be used in a suitable extension of Rabi's formula to the nonorthogonal representation of two-state dynamical problems. The formula for the transfer integral is inspected in the regime of long-range CT and applied to CT reactions in redox molecular systems of interest to biochemistry and/or to molecular electronics: the guanine-thymine stack from regular B-DNA, the polyaromatic perylenediimide stack, and the quinol-semiquinone couple. The calculations are performed within the framework of the Density Functional Theory (DFT), using hybrid exchange-correlation (XC) density functionals, which also allowed investigation of the appropriateness of such hybrid-DFT methods for computing electronic couplings. The use of the recently developed M06-2X and M06-HF density functionals in appropriate ways is supported by the results of this work.

Energy decomposition analysis based on a block-localized wavefunction and multistate density functional theory

An interaction energy decomposition analysis method based on the block-localized wavefunction (BLW-ED) approach is described. The first main feature of the BLW-ED method is that it combines concepts of valence bond and molecular orbital theories such that the intermediate and physically intuitive electron-localized states are variationally optimized by self-consistent field calculations. Furthermore, the block-localization scheme can be used both in wave function theory and in density functional theory, providing a useful tool to gain insights on intermolecular interactions that would otherwise be difficult to obtain using the delocalized Kohn-Sham DFT. These features allow broad applications of the BLW method to energy decomposition (BLW-ED) analysis for intermolecular interactions. In this perspective, we outline theoretical aspects of the BLW-ED method, and illustrate its applications in hydrogen-bonding and π-cation intermolecular interactions as well as metal-carbonyl complexes. Future prospects on the development of a multistate density functional theory (MSDFT) are presented, making use of block-localized electronic states as the basis configurations.

A Non-Orthogonal Block-Localized Effective Hamiltonian Approach for Chemical and Enzymatic Reactions

The effective Hamiltonian-molecular orbital and valence bond (EH-MOVB) method based on non-orthogonal block-localized fragment orbitals has been implemented into the program CHARMM for molecular dynamics simulations of chemical and enzymatic reactions, making use of semiempirical quantum mechanical models. Building upon ab initio MOVB theory, we make use of two parameters in the EH-MOVB method to fit the barrier height and the relative energy between the reactant and product state for a given chemical reaction to be in agreement with experiment or high-level ab initio or density functional results. Consequently, the EH-MOVB method provides a highly accurate and computationally efficient QM/MM model for dynamics simulation of chemical reactions in solution. The EH-MOVB method is illustrated by examination of the potential energy surface of the hydride transfer reaction from trimethylamine to a flavin cofactor model in the gas phase. In the present study, we employed the semiempirical AM1 model, which yields a reaction barrier that is more than 5 kcal/mol too high. We use a parameter calibration procedure for the EH-MOVB method similar to that employed to adjust the results of semiempirical and empirical models. Thus, the relative energy of these two diabatic states can be shifted to reproduce the experimental energy of reaction, and the barrier height is optimized to reproduce the desired (accurate) value by adding a constant to the off-diagonal matrix element. The present EH-MOVB method offers a viable approach to characterizing solvent and protein-reorganization effects in the realm of combined QM/MM simulations.

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

Diabatic states have a long history in chemistry, beginning with early valence bond pictures of molecular bonding and extending through the construction of model potential energy surfaces to the modern proliferation of methods for computing these elusive states. In this review, we summarize the basic principles that define the diabatic basis and demonstrate how they can be applied in the specific context of constrained density functional theory. Using illustrative examples from electron transfer and chemical reactions, we show how the diabatic picture can be used to extract qualitative insight and quantitative predictions about energy landscapes. The review closes with a brief summary of the challenges and prospects for the further application of diabatic states in chemistry.