Multistate Density Functional Theory Research Articles

Combined Multistate and Kohn-Sham Density Functional Theory Studies of the Elusive Mechanism of N-Dealkylation of N,N-Dimethylanilines Mediated by the Biomimetic Nonheme Oxidant FeIV(O)(N4Py)(ClO4)2.

The oxidative C-H bond activation mediated by heme and nonheme enzymes and related biomimetics is one of the most interesting processes in bioinorganic and oxidative chemistry. However, the mechanisms of these reactions are still elusive and controversy due to the involvement of highly reactive metal-oxo intermediates with multiple spin states, despite extensive experimental efforts, especially for the N-dealkylation of N,N-dialkyalinines. In this work, we employed multistate density functional theory (MSDFT) and the Kohn-Sham DFT to investigate the mechanism of N-demethylation of N,N-dimethyalinines oxidized by the reaction intermediate FeIV(O)(N4Py)(ClO4)2. The Kohn-Sham DFT study demonstrated that the reaction proceeds via a rate-limiting hydrogen atom transfer (HAT) step and a subsequent barrier-free oxygen rebound step to form the carbinol product. The MSDFT investigation on the first C-H activation further showed that this step is an initial hydrogen atom abstraction that is highly correlated between CEPT and HAT, i.e., both CEPT and HAT processes make significant contributions to the mechanism before reaching the diabatic crossing point, then the valence bond character of the adiabatic ground state is switched to the CEPT product configuration. The findings from this work may be applicable to other hydrogen abstraction process.

The charger transfer electronic coupling in diabatic perspective: A multi-state density functional theory study

The multi-state density functional theory (MSDFT) provides a convenient way to estimate electronic coupling of charge transfer processes based on a diabatic representation. Its performance has been benchmarked against the HAB11 database with a mean unsigned error (MUE) of 17 meV between MSDFT and ab initio methods. The small difference may be attributed to different representations, diabatic from MSDFT and adiabatic from ab initio calculations. In this discussion, we conclude that MSDFT provides a general and efficient way to estimate the electronic coupling for charge-transfer rate calculations based on the Marcus-Hush model.

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.

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.

The Quantum Coherent Mechanism for Singlet Fission: Experiment and Theory

The absorption of one photon by a semiconductor material usually creates one electron-hole pair. However, this general rule breaks down in a few organic semiconductors, such as pentacene and tetracene, where one photon absorption may result in two electron-hole pairs. This process, where a singlet exciton transforms to two triplet excitons, can have quantum yields as high as 200%. Singlet fission may be useful to solar cell technologies to increase the power conversion efficiency beyond the so-called Shockley-Queisser limit. Through time-resolved two-photon photoemission (TR-2PPE) spectroscopy in crystalline pentacene and tetracene, our lab has recently provided the first spectroscopic signatures in singlet fission of a critical intermediate known as the multiexciton state (also called a correlated triplet pair). More importantly, we found that population of the multiexciton state rises at the same time as the singlet state on the ultrafast time scale upon photoexcitation. This observation does not fit with the traditional view of singlet fission involving the incoherent conversion of a singlet to a triplet pair. However, it provides an experimental foundation for a quantum coherent mechanism in which the electronic coupling creates a quantum superposition of the singlet and the multiexciton state immediately after optical excitation. In this Account, we review key experimental findings from TR-2PPE experiments and present a theoretical analysis of the quantum coherent mechanism based on electronic structural and density matrix calculations for crystalline tetracene lattices. Using multistate density functional theory, we find that the direct electronic coupling between singlet and multiexciton states is too weak to explain the experimental observation. Instead, indirect coupling via charge transfer intermediate states is two orders of magnitude stronger, and dominates the dynamics for ultrafast multiexciton formation. Density matrix calculation for the crystalline tetracene lattice satisfactorily accounts for the experimental observations. It also reveals the critical roles of the charge transfer states and the high dephasing rates in ensuring the ultrafast formation of multiexciton states. In addition, we address the origins of microscopic relaxation and dephasing rates, and adopt these rates in a quantum master equation description. We show the need to take the theoretical effort one step further in the near future by combining high-level electronic structure calculations with accurate quantum relaxation dynamics for large systems.

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.

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.

Block-Localized Density Functional Theory (BLDFT), Diabatic Coupling, and Their Use in Valence Bond Theory for Representing Reactive Potential Energy Surfaces.

A multistate density functional theory in the framework of the valence bond model is described. The method is based on a block-localized density functional theory (BLDFT) for the construction of valence-bond-like diabatic electronic states and is suitable for the study of electron transfer reactions and for the representation of reactive potential energy surfaces. The method is equivalent to a valence bond theory with the treatment of the localized configurations by using density functional theory (VBDFT). In VBDFT, the electron densities and energies of the valence bond states are determined by BLDFT. A functional estimate of the off-diagonal matrix elements of the VB Hamiltonian is proposed, making use of the overlap integral between Kohn-Sham determinants and the exchange-correlation functional for the ground state substituted with the transition (exchange) density. In addition, we describe an approximate approach, in which the off-diagonal matrix element is computed by wave function theory using block-localized Kohn-Sham orbitals. The key feature is that the electron density of the adiabatic ground state is not directly computed nor used to obtain the ground-state energy; the energy is determined by diagonalization of the multistate valence bond Hamiltonian. This represents a departure from the standard single-determinant Kohn-Sham density functional theory. The multistate VBDFT method is illustrated by the bond dissociation of H2+ and a set of three nucleophilic substitution reactions in the DBH24 database. In the dissociation of H2+, the VBDFT method yields the correct asymptotic behavior as the two protons stretch to infinity, whereas approximate functionals fail badly. For the S(N)2 nucleophilic substitution reactions, the hybrid functional B3LYP severely underestimates the barrier heights, while the approximate two-state VBDFT method overcomes the self-interaction error, and overestimates the barrier heights. Inclusion of the ionic state in a three-state model, VBDFT(3), significantly improves the computed barrier heights, which are found to be in accord with accurate results. The BLDFT method is a versatile theory that can be used to analyze conventional DFT results to gain insight into chemical bonding properties, and it is illustrated by examining the intricate energy contributions to the ion-dipole complex stabilization.