Excited-State Energy Decomposition Analysis.

Multistate density functional theory for local and charge-transfer tripdoublet states from triplet-free radical interactions.

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.

Multistate density functional theory (MSDFT) employing a minimum active space (MAS) is presented to determine charge transfer (CT) and local excited states of bimolecular complexes. MSDFT is a hybrid wave function theory (WFT) and density functional theory, in which dynamic correlation is first incorporated in individual determinant configurations using a Kohn–Sham exchange-correlation functional. Then, nonorthogonal configuration-state interaction is performed to treat static correlation. Because molecular orbitals are optimized separately for each determinant by including Kohn–Sham dynamic correlation, a minimal number of configurations in the active space, essential to representing low-lying excited and CT states of interest, is sufficient to yield the adiabatic states. We found that the present MAS-MSDFT method provides a good description of covalent and CT excited states in comparison with experiments and high-level computational results. Because of the simplicity and interpretive capability through diabatic configuration weights, the method may be useful in dynamic simulations of CT and nonadiabatic processes.

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.

We report a rigorous formulation of multi-state density functional theory (MSDFT) that extends the Kohn-Sham (KS) energy functional for the ground state to a Hamiltonian matrix functional H[D] of the density matrix D in the space spanned by the lowest N adiabatic states. We establish a variational principle of MSDFT, which guarantees that the variational optimization results in a Hamiltonian matrix, whose eigenvalues are the lowest N eigen-energies of the system. We present an explicit expression of H[D] and introduce the correlation matrix functional. Akin to KS-DFT for the ground state, a universal multi-state correlation potential is derived for a two-state system as an illustrative example. This work shows that MSDFT is an exact density functional theory that treats the ground and excited states on an equal footing and provides a framework for practical applications and future developments of approximate functionals for excited states.

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.

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

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.

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.

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.

In this Perspective, we introduce a minimal active space (MAS) for the lowest N eigenstates of a molecular system in the framework of multistate density functional theory (MSDFT), consisting of no more than N2 nonorthgonal Slater determinants. In comparison with some methods in wave function theory in which one seeks to expand the ever increasing size of an active space to approximate the wave functions, it is possible to have an upper bound in MSDFT because the auxiliary states in a MAS are used to represent the exact N-dimensional matrix density function D(r). Here, we partition the total Hamiltonian matrix functional into an orbital-dependent part, including multistate kinetic energy Tms and Coulomb-exchange energy EHx plus an external potential energy ∫dr v(r)D(r), and a correlation matrix density functional . The latter accounts for the part of correlation energy not explicitly included in the minimal active space. A major difference from Kohn-Sham DFT is that state interactions are necessary to represent the N-matrix density D(r) in MSDFT, rather than a noninteracting reference state for the scalar ground-state density ρo(r). Two computational approaches are highlighted. We first derive a set of nonorthogonal multistate self-consistent-field (NOSCF) equations for the variational optimization of . We introduce the multistate correlation potential, as the functional derivative of , which includes both correlation effects within the MAS and that from the correlation matrix functional. Alternatively, we describe a nonorthogonal state interaction (NOSI) procedure, in which the determinant functions are optimized separately. Both computational methods are useful for determining the exact eigenstate energies and for constructing variational diabatic states, provided that the universal correlation matrix functional is known. It is hoped that this discussion would stimulate developments of approximate multistate density functionals both for the ground and excited 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.

We report a rigorous formulation of density functional theory for excited states, providing a theoretical foundation for a multistate density functional theory. We prove the existence of a Hamiltonian matrix functional of the multistate matrix density D(r) in the subspace spanned by the lowest N eigenstates. Here, D(r) is an N-dimensional matrix of state densities and transition densities. Then, a variational principle of the multistate subspace energy is established, whose minimization yields both the energies and densities of the individual N eigenstates. Furthermore, we prove that the N-dimensional matrix density D(r) can be sufficiently represented by N2 nonorthogonal Slater determinants, based on which an interacting active space is introduced for practical calculations. This work establishes that the ground and excited states can be treated on an equal footing in density functional theory.

An intramolecular excited charge transfer (CT) analysis of imidazole derivatives has been made. The determined electronic transition dipole moments has been used to estimate the electronic coupling interactions between the excited charge transfer singlet state ((1)CT) and the ground state (S0) or the locally excited state ((1)LE). The properties of excited (1)CT state imidazole derivatives have been exploited by the significant contribution of the electronic coupling interactions. The excited state intramolecular proton transfer (ESIPT) analysis has also been discussed.

Poly-CPDTBT, as typical low-band gap copolymers, have potential applications in organic bulk heterojunction solar cells. To have a clear picture of its excited-state processes, the first task is to understand their excited states, in particular, electronic character and relevant optical absorption. Herein, the low-lying singlet excited states of Poly-CPDTBT oligomers were investigated via Algebraic Diagrammatic Construction Second Order (ADC(2)) and time-dependent density functional theory (TDDFT) method with several functionals. Six CPDTBTN (N = 1–6) oligomers were taken as prototypes to study their excited states in detail. The results provide interesting clues to extrapolate the photophysical properties of such polymers with potential applications in photovoltaic materials. The result provided by ωB97XD functional gives good agreement with the experiment result. The vertical excitation energies of the four lowest excited states decrease almost linearly with increasing polymerization degree (N) for CPDTBTN (N = 1–6). The transition density analysis indicates that the local excitations (LE) and the short-distance charge transfer (CT) excitations between two adjacent CPDT and BT units are dominant for low-lying excited states for short oligomers. For the long-chain oligomers (trimer to hexamer), the transition density shows a ladder (or zigzag) pattern along the diagonal blocks at the planar geometry. For long oligomers, the whole chain is involved in the transitions, and the CT excitations only exist between two adjacent CPDT and BT units. The present work provides a valuable basis for understanding the excited-state processes of Poly-CPDTBT and other conjugated polymers that conduct solar energy conversions, which has great significance for the development of new solar cells.