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

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

Using the constrained DFT approach in generating diabatic surfaces and off diagonal empirical valence bond terms for modeling reactions in condensed phases

The empirical valence bond (EVB) model provides an extremely powerful way for modeling and analyzing chemical reactions in solutions and proteins. However, this model is based on the unverified assumption that the off diagonal elements of the EVB Hamiltonian do not change significantly upon transfer of the reacting system from one phase to another. This ad hoc assumption has been rationalized by its consistency with empirically observed linear free energy relationships, as well as by other qualitative considerations. Nevertheless, this assumption has not been rigorously established. The present work explores the validity of the above EVB key assumption by a rigorous numerical approach. This is done by exploiting the ability of the frozen density functional theory (FDFT) and the constrained density functional theory (CDFT) models to generate convenient diabatic states for QM/MM treatments, and thus to examine the relationship between the diabatic and adiabatic surfaces, as well as the corresponding effective off diagonal elements. It is found that, at least for the test case of S(N)()2 reactions, the off diagonal element does not change significantly upon moving from the gas phase to solutions and thus the EVB assumption is valid and extremely useful.     

Fast prediction of hydration free energies for SAMPL4 blind test from a classical density functional theory

We report the performance of a classical density functional theory (CDFT) in the competition for the solvation free-energy category of the SAMPL4 blind prediction event. The theoretical calculations were carried out with the TIP3P water model and different combinations of solute configurations and molecular force fields. In comparison with the experimental data, the blind test yields an average unsigned error of 2.38 kcal/mol and the root mean square deviation of 2.99 kcal/mol. Whereas these numbers are significantly larger than the best results from explicit-solvent MD simulations, we find that the theoretical performance is sensitive to both the molecular force fields and solute configurations and that a comparable level of accuracy can be achieved by a judicious selection of the solute configurations and the force-field parameters. Most importantly, CDFT reduces the computational cost of MD simulation by almost 3 orders of magnitude, making it very attractive for large-scale hydration free-energy calculations (e.g., screening the aqueous solubility of drug-like molecules).     

An Anderson impurity model for efficient sampling of adiabatic potential energy surfaces of transition metal complexes

We present a model intended for rapid sampling of ground and excited state potential energy surfaces for first-row transition metal active sites. The method is computationally inexpensive and is suited for dynamics simulations where (1) adiabatic states are required "on-the-fly" and (2) the primary source of the electronic coupling between the diabatic states is the perturbative spin-orbit interaction among the 3d electrons. The model Hamiltonian we develop is a variant of the Anderson impurity model and achieves efficiency through a physically motivated basis set reduction based on the large value of the d-d Coulomb interaction U(d) and a Lanczos matrix diagonalization routine to solve for eigenvalues. The model parameters are constrained by fits to the partial density of states obtained from ab initio density functional theory calculations. For a particular application of our model we focus on electron transfer occurring between cobalt ions solvated by ammonium, incorporating configuration interaction between multiplet states for both metal ions. We demonstrate the capability of the method to efficiently calculate adiabatic potential energy surfaces and the electronic coupling factor we have calculated compares well to previous calculations and experiment. (     

Exact density functionals for two-electron systems in an external magnetic field

In principle, the extension of density functional theory (DFT) to Coulombic systems in a nonvanishing magnetic field is via current DFT (CDFT). Though CDFT is long established formally, relatively little is known with respect to any generally applicable, reliable approximate E(XC) and A(XC) functionals analogous with the workhorse approximate functionals (local density approximation and generalized gradient approximation) of ordinary DFT. Progress can be aided by having benchmark studies on a solvable correlated system. At zero field, the best-known finite system for such purposes is Hooke's atom. Recently we extended the exact ground state solutions for this two-electron system to certain combinations of nonzero external magnetic fields and confinement strengths. From those exact solutions, as well as high-accuracy numerical results for other field and confinement combinations, we construct the correlated electron density and paramagnetic current density, the exact Kohn-Sham orbitals, and the exact DFT and CDFT exchange-correlation energies and potentials. We compare with results from several widely used approximate functionals, all of which exhibit major qualitative failures, whether in CDFT or in naive application of ordinary DFT. We also illustrate how the CDFT vorticity variable nu is a computationally difficult quantity which may not be appropriate in practice to describe the external B field effects on E(XC) and A(XC).     

Ab initio calculations on the excited states of Na3 cluster to explore beyond Born-Oppenheimer theories: adiabatic to diabatic potential energy surfaces and nuclear dynamics

We perform ab initio calculation using quantum chemistry package (MOLPRO) on the excited states of Na(3) cluster and present the adiabatic PESs for the electronic states 2(2)E' and 1(2)A(1)', and the non-adiabatic coupling (NAC) terms among those states. Since the ab initio calculated NAC elements for the states 2(2)E' and 1(2)A(1)' demonstrate the numerical validity of so called "Curl Condition," such states closely form a sub-Hilbert space. For this subspace, we employ the NAC terms to solve the "adiabatic-diabatic transformation (ADT)" equations to obtain the functional form of the transformation angles and pave the way to construct the continuous and single valued diabatic potential energy surface matrix by exploiting the existing first principle based theoretical means on beyond Born-Oppenheimer treatment. Nuclear dynamics has been carried out on those diabatic surfaces to reproduce the experimental spectrum for system B of Na(3) cluster and thereby, to explore the numerical validity of the theoretical development on beyond Born-Oppenheimer approach for adiabatic to diabatic transformation.     

Direct calculation of coupled diabatic potential-energy surfaces for ammonia and mapping of a four-dimensional conical intersection seam

We used multiconfiguration quasidegenerate perturbation theory and the fourfold-way direct diabatization scheme to calculate ab initio potential-energy surfaces at 3600 nuclear geometries of NH3. The calculations yield the adiabatic and diabatic potential-energy surfaces for the ground and first electronically excited singlet states and also the diabatic coupling surfaces. The diabatic surfaces and coupling were fitted analytically to functional forms to obtain a permutationally invariant 2 x 2 diabatic potential-energy matrix. An analytic representation of the adiabatic potential-energy surfaces is then obtained by diagonalizing the diabatic potential-energy matrix. The analytic representation of the surfaces gives an analytic representation of the four-dimensional conical intersection seam which is discussed in detail.     

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

A theoretical model is presented for deriving effective diabatic states based on ab initio valence bond self-consistent field (VBSCF) theory by reducing the multiconfigurational VB Hamiltonian into an effective two-state model. We describe two computational approaches for the optimization of the effective diabatic configurations, resulting in two ways of interpreting such effective diabatic states. In the variational diabatic configuration (VDC) method, the energies of the diabatic states are variationally minimized. In the consistent diabatic configuration (CDC) method, both the configuration coefficients and orbital coefficients are simultaneously optimized to minimize the adiabatic ground-state energy in VBSCF calculations. In addition, we describe a mixed molecular orbital and valence bond (MOVB) approach to construct the CDC diabatic and adiabatic states for a chemical reaction. Note that the VDC-MOVB method has been described previously. Employing the symmetric S(N)2 reaction between NH(3) and CH(3)NH(3)(+) as a test system, we found that the results from ab initio VBSCF and from ab initio MOVB calculations using the same basis set are in good agreement, suggesting that the computationally efficient MOVB method is a reasonable model for VB simulations of condensed phase reactions. The results indicate that CDC and VDC diabatic states converge, respectively, to covalent and ionic states as the molecular geometries are distorted from the minimum of the respective diabatic state along the reaction coordinate. Furthermore, the resonance energy that stabilizes the energy of crossing between the two diabatic states, resulting in the transition state of the adiabatic ground-state reaction, has a strong dependence on the overlap integral between the two diabatic states and is a function of both the exchange integral and the total diabatic ground-state energy.     

Time-dependent density functional theory gradients in the Amsterdam density functional package: geometry optimizations of spin-flip excitations

An implementation of time-dependent density functional theory (TDDFT) energy gradients into the Amsterdam density functional theory program package (ADF) is described. The special challenges presented by Slater-type orbitals in quantum chemical calculation are outlined with particular emphasis on details that are important for TDDFT gradients. Equations for the gradients of spin-flip TDDFT excitation energies are derived. Example calculations utilizing the new implementation are presented. The results of standard calculations agree well with previous results. It is shown that starting from a triplet reference, spin-flip TDDFT can successfully optimize the geometry of the four lowest singlet states of CH₂ and three other isovalent species. Spin-flip TDDFT is used to calculate the potential energy curve of the breaking of the C?CC bond of ethane. The curve obtained is superior to that from a restricted density functional theory calculation, while at the same time the problems with spin contamination exhibited by unrestricted density functional theory calculations are avoided.     

Vibronic coupling and double excitations in linear response time-dependent density functional calculations: dipole-allowed states of N2

The present study serves two purposes. First, we evaluate the ability of present time-dependent density functional response theory (TDDFRT) methods to deal with avoided crossings, i.e., vibronic coupling effects. In the second place, taking the vibronic coupling effects into account enables us, by comparison to the configuration analysis in a recent ab initio study [J. Chem. Phys. 115, 6438 (2001)], to identify the neglect of double excitations as the prime cause of limited accuracy of these linear response based TDDFRT calculations for specific states. The "statistical averaging of (model) orbital potentials (SAOP)" Kohn-Sham potential is used together with the standard adiabatic local-density approximation (ALDA) for the exchange-correlation kernel. We use the N2 molecule as prototype, since the TDDFRT/SAOP calculations have already been shown to be accurate for the vertical excitations, while this molecule has a well-studied, intricate vibronic structure as well as significant double excitation nature in the lowest 1Pi(u) state at elongated bond lengths. A simple diabatizing scheme is employed to obtain a diabatic potential energy matrix, from which we obtain the absorption spectrum of N2 including vibronic coupling effects. Considering the six lowest dipole allowed transitions of 1Sigma(u)+ and 1Pi(u) symmetry, we observe a good general agreement and conclude that avoided crossings and vibronic coupling can indeed be treated satisfactorily on the basis of TDDFRT excitation energies. However, there is one state for which the accuracy of TDDFRT/ALDA clearly breaks down. This is the state for which the ab initio calculations find significant double excitation character. To deal with double excitation character is an important challenge for time-dependent density functional theory.     

Linking the historical and chemical definitions of diabatic states for charge and excitation energy transfer reactions in condensed phase

Marcus theory of electron transfer (ET) and Fo?rster theory of excitation energy transfer (EET) rely on the Condon approximation and the theoretical availability of initial and final states of ET and EET reactions, often called diabatic states. Recently [Subotnik et al., J. Chem. Phys. 130, 234102 (2009)], diabatic states for practical calculations of ET and EET reactions were defined in terms of their interactions with the surrounding environment. However, from a purely theoretical standpoint, the definition of diabatic states must arise from the minimization of the dynamic couplings between the trial diabatic states. In this work, we show that if the Condon approximation is valid, then a minimization of the derived dynamic couplings leads to corresponding diabatic states for ET reactions taking place in solution by diagonalization of the dipole moment matrix, which is equivalent to a Boys localization algorithm; while for EET reactions in solution, diabatic states are found through the Edmiston-Ruedenberg localization algorithm. In the derivation, we find interesting expressions for the environmental contribution to the dynamic coupling of the adiabatic states in condensed-phase processes. In one of the cases considered, we find that such a contribution is trivially evaluable as a scalar product of the transition dipole moment with a quantity directly derivable from the geometry arrangement of the nuclei in the molecular environment. Possibly, this has applications in the evaluation of dynamic couplings for large scale simulations.     

Optimal diabatic bases via thermodynamic bounds

Describing kinetic processes within a perturbation theory approach such as Fermi's golden rule requires an understanding of the initial and final states of the system. A number of different methods have been proposed for obtaining these diabatic-like states, but a robust criterion for evaluating their accuracy has not been established. Here, we approach the problem of determining the most appropriate set of diabatic states for use in incoherent rate expressions. We develop a method that rotates an initial set of diabats into an optimized set beginning with a zeroth-order diabatic Hamiltonian and choosing the rotation that minimizes the effect of non-diabatic terms on the thermodynamic free energy. The Gibbs-Bogoliubov (GB) bound on the Helmholtz free energy is thus used as the diabatic criterion. We first derive the GB free energy for a two site system and then find an expression general for any electronic system Hamiltonian. Efficient numerical methods are used to perform the minimization subject to orthogonality constraints, and we examine the resulting diabats for system Hamiltonians in various parameter regimes. The transition from localized to delocalized states is clearly seen in these calculations, and some interesting features are discussed.     

Two-state model based on the block-localized wave function method

The block-localized wave function (BLW) method is a variant of ab initio valence bond method but retains the efficiency of molecular orbital methods. It can derive the wave function for a diabatic (resonance) state self-consistently and is available at the Hartree-Fock (HF) and density functional theory (DFT) levels. In this work we present a two-state model based on the BLW method. Although numerous empirical and semiempirical two-state models, such as the Marcus-Hush two-state model, have been proposed to describe a chemical reaction process, the advantage of this BLW-based two-state model is that no empirical parameter is required. Important quantities such as the electronic coupling energy, structural weights of two diabatic states, and excitation energy can be uniquely derived from the energies of two diabatic states and the adiabatic state at the same HF or DFT level. Two simple examples of formamide and thioformamide in the gas phase and aqueous solution were presented and discussed. The solvation of formamide and thioformamide was studied with the combined ab initio quantum mechanical and molecular mechanical Monte Carlo simulations, together with the BLW-DFT calculations and analyses. Due to the favorable solute-solvent electrostatic interaction, the contribution of the ionic resonance structure to the ground state of formamide and thioformamide significantly increases, and for thioformamide the ionic form is even more stable than the covalent form. Thus, thioformamide in aqueous solution is essentially ionic rather than covalent. Although our two-state model in general underestimates the electronic excitation energies, it can predict relative solvatochromic shifts well. For instance, the intense pi-->pi* transition for formamide upon solvation undergoes a redshift of 0.3 eV, compared with the experimental data (0.40-0.5 eV).     

Current‐density functional study of the HeH+ molecular ion under a strong ultrashort magnetic field

The HeH⁺ molecular ion under an ultrashort magnetic field on the order of 10⁹ G is investigated through quantum fluid dynamics and a current‐density functional theory (CDFT) based approach, employing a vector exchange–correlation (XC) potential which depends on the electronic charge‐density as well as on the current‐density. The behavior of the exchange and correlation energies of the HeH⁺ ion is analyzed and compared with those obtained using an approach based on the time‐dependent density functional theory (TD‐DFT) under similar computational constraints but employing a scalar XC potential dependent only on the electronic charge‐density. The CDFT‐based approach yields exchange and correlation energies as well as TD electronic charge‐densities drastically different from those obtained using the TD‐DFT‐based approach particularly, at typical TD magnetic field strengths. This is attributed to the nonadiabatic effects induced by the vector XC potential of the CDFT in the oscillating charge‐density of the HeH⁺ ion, which are further explained in the terminology of quantum fluid dynamics. The vector XC potential of the CDFT‐based approach is observed to augment the magnetic interactions in the H₂ molecule and in the He ion, whereas it opposes the magnetic interactions in the HeH⁺ ion particularly, at the intermediate magnetic field strengths. © 2012 Wiley Periodicals, Inc.     

The electronic structure and energetics of V(+)-benzene half-sandwiches of different multiplicities: Comparative multireference and single-reference theoretical study

Multireference [complete active space self-consistent field (CASSCF) and multiconfigurational quasidegenerate perturbation theory (MCQDPT)] and single-reference ab initio (Moller-Plesset second order perturbation theory (MP2) and coupled clusters with singles, doubles and noniterative triples [CCSD(T)]) and density functional theory (PBE and B3LYP) electronic structure calculations of V(C(6)H(6))(+) half-sandwich in the states of different multiplicities are described and compared. Detailed analyses of the geometries and electronic structures of the all found states are given; adiabatic and diabatic dissociation energies are estimated. The lowest electronic state of V(C(6)H(6))(+) half-sandwich was found to be the quintet (5)B(2) state with a slightly deformed upside-down-boat-shaped benzene ring and d(4) configuration of V atom, followed by a triplet (3)A(2) state lying about 4 kcal/mol above. The lowest singlet state (1)A(1)(d(4)) lies much ( approximately 28 kcal/mol) higher. MCQDPT calculated adiabatic dissociation energy (53.6 kcal/mol) for the lowest (5)B(2)(d(4)) state agrees well with the current 56.4 (54.4) kcal/mol experimental estimate, giving a preference to the lower one. Compared to MCQDPT, B3LYP hybrid exchange-correlation functional provides the best results, while CCSD(T) performs usually worse. Gradient-corrected PBE calculations tend to systematically overestimate metal-benzene binding in the row quintet相似文献   

Theoretical Investigation on Triplet Excitation Energy Transfer in Fluorene Dimer

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