Spin densities from subsystem density-functional theory: assessment and application to a photosynthetic reaction center complex model

Subsystem density-functional theory (DFT) is a powerful and efficient alternative to Kohn-Sham DFT for large systems composed of several weakly interacting subunits. Here, we provide a systematic investigation of the spin-density distributions obtained in subsystem DFT calculations for radicals in explicit environments. This includes a small radical in a solvent shell, a π-stacked guanine-thymine radical cation, and a benchmark application to a model for the special pair radical cation, which is a dimer of bacteriochlorophyll pigments, from the photosynthetic reaction center of purple bacteria. We investigate the differences in the spin densities resulting from subsystem DFT and Kohn-Sham DFT calculations. In these comparisons, we focus on the problem of overdelocalization of spin densities due to the self-interaction error in DFT. It is demonstrated that subsystem DFT can reduce this problem, while it still allows to describe spin-polarization effects crossing the boundaries of the subsystems. In practical calculations of spin densities for radicals in a given environment, it may thus be a pragmatic alternative to Kohn-Sham DFT calculations. In our calculation on the special pair radical cation, we show that the coordinating histidine residues reduce the spin-density asymmetry between the two halves of this system, while inclusion of a larger binding pocket model increases this asymmetry. The unidirectional energy transfer in photosynthetic reaction centers is related to the asymmetry introduced by the protein environment.     

A density functional theory for studying ionization processes in water clusters

A generalized Kohn-Sham (GKS) approach to density functional theory (DFT), based on the Baer-Neuhauser-Livshits range-separated hybrid, combined with ab initio motivated range-parameter tuning is used to study properties of water dimer and pentamer cations. The water dimer is first used as a benchmark system to check the approach. The present brand of DFT localizes the positive charge (hole), stabilizing the proton transferred geometry in agreement with recent coupled-cluster calculations. Relative energies of various conformers of the water dimer cation compare well with previously published coupled cluster results. The GKS orbital energies are good approximations to the experimental ionization potentials of the system. Low-lying excitation energies calculated from time-dependent DFT based on the present method compare well with recently published high-level "equation of motion-coupled-cluster" calculations. The harmonic frequencies of the water dimer cation are in good agreement with experimental and wave function calculations where available. The method is applied to study the water pentamer cation. Three conformers are identified: two are Eigen type and one is a Zundel type. The structure and harmonic vibrational structure are analyzed. The ionization dynamics of a pentamer water cluster at 0 K shows a fast <50 fs transient for transferring a proton from one of the water molecules, releasing a hydroxyl radical and creating a protonated tetramer carrying the excess hole.     

Legendre-transform functionals for spin-density-functional theory

We provide a rigorous proof that the Hohenberg-Kohn theorem holds for spin densities by extending Lieb's Legendre-transform formulation to spin densities. The resulting spin-density-functional theory resolves several troublesome issues. Most importantly, the present paper provides an explicit construction for the spin potentials at any point along the adiabatic connection curve, thus providing a formal basis for the use of exchange-correlation functionals of the spin density in the Kohn-Sham density-functional theory (DFT). The practical implications of this result for unrestricted Kohn-Sham DFT calculations is considered, and the existence of holes below the Fermi level is discussed. We argue that an orbital's energy tends to increase as its occupation number increases, which provides the basis for a computational algorithm for determining the occupation numbers in Kohn-Sham DFT and helps explain the origin of Hund's rules and holes below the Fermi level.     

Single molecule conductivity: the role of junction-orbital degeneracy in the artificially high currents predicted by ab initio approaches

A priori evaluations, using Hartree-Fock self-consistent-field (SCF) theory or density-functional theory (DFT), of the current passing between two electrodes through a single bridging molecule result in predicted conductivities that may be up to one to two orders of magnitude larger than observed ones. We demonstrate that this is, in part, often due to the improper application of the computational methods. Conductivity is shown to arise from tunneling between junction states of the electrodes through the molecule; these states are inherently either quasi two-fold or four-fold degenerate and always comprise the (highest occupied molecular orbital) HOMO band at the Fermi energy of the system. Frequently, in previous cluster based molecular conduction calculations, closed-shell SCF or Kohn-Sham DFT methods have been applied to systems that we demonstrate to be intrinsically open shell in nature. Such calculations are shown to induce artificial HOMO-LUMO (LUMO-lowest unoccupied molecular orbital) band splittings that Landauer-based formalisms for steady-state conduction interpret as arising from extremely rapid through-molecule tunneling at the Fermi energy, hence, overestimating the low-voltage conductivity. It is demonstrated that these shortcomings can be eliminated, dramatically reducing calculated current magnitudes, through the alternate use of electronic-structure calculations based on the spin-restricted open-shell formalism and related multiconfigurational SCF of DFT approaches. Further, we demonstrate that most anomalies arising in DFT implementations arise through the use of hybrid density functionals such as B3LYP. While the enhanced band-gap properties of these functionals have made them the defacto standard in molecular conductivity calculations, we demonstrate that it also makes them particularly susceptible to open-shell anomalies.     

Accurate magnetic exchange couplings in transition-metal complexes from constrained density-functional theory

We demonstrate an accurate method for extracting Heisenberg exchange-coupling constants (J) from density-functional theory (DFT) calculations. We note that the true uncoupled low-spin state of a given molecule should be identified with the ground state of the system subject to a constraint on the spin density of the atoms. Using an efficient optimization strategy for constrained DFT we obtain these states directly, leading to a simple, physically motivated formula for J. Our method only depends on state energies and their associated electron densities and assigns no unphysical meaning to the Kohn-Sham determinant or individual orbitals. We study several bimetallic transition-metal complexes and find that the constrained DFT approach is competitive with, if not better than, the best broken symmetry DFT results. The success of constrained DFT in these cases appears to result from a balanced elimination of self-interaction error and static correlation from the simulation.     

Nonlinear optical property calculations by the long-range-corrected coupled-perturbed Kohn-Sham method

The long-range correction (LC) scheme for the exchange functional of density-functional theory (DFT) was combined with the coupled-perturbed Kohn-Sham (CPKS) method to calculate nonlinear optical response properties. By using this LC-CPKS method, we calculated the hyperpolarizabilities of typical molecules and the dipole moments, polarizabilities, and hyperpolarizabilities of push-pull pi-conjugated systems: p-nitroaniline, 4-amino-4'-nitrostilbene, and alpha,omega-nitroaminopolyenes. It was found that the LC scheme clearly improved the calculation of these optical properties for all of these systems, which have been significantly overestimated by conventional DFTs. We therefore concluded that the long-range exchange interaction played an important role in calculating the optical properties using the DFT formalism.     

DFT description of the electronic structure and spectromagnetic properties of strongly correlated electronic systems: NiII, CuII and ZnII o-dioxolene complexes

The spectroscopic and magnetic properties of dioxolene complexes of zinc, copper and nickel were studied by DFT calculations on model complexes of formulas [(NH(3))(4)M(II)(SQ)](+) (M=Zn, Ni; SQ=semiquinonato) and [(NH(3))(2)Cu(II)(SQ)](+). Standard approaches such as time-dependent DFT (TDDFT), the Slater transition state (STS), and broken symmetry (BS) were found to be unable to completely account for the physical properties of the systems, and complete active space-configuration interaction (CAS-CI) calculations based on the Kohn-Sham (KS) orbitals was applied. The CAS-CI energies, properly corrected with multireference perturbation theory (MR-PT), were found to be in good agreement with experimental data. We present here a calculation protocol that has a low CPU cost/accuracy ratio and seems to be very promising for interpreting the properties of strongly correlated electronic systems in complexes of real chemical size.     

Topological analysis of electron densities from Kohn-Sham and subsystem density functional theory

In this study, we compare the electron densities for a set of hydrogen-bonded complexes obtained with either conventional Kohn-Sham density functional theory (DFT) calculations or with the frozen-density embedding (FDE) method, which is a subsystem approach to DFT. For a detailed analysis of the differences between these two methods, we compare the topology of the electron densities obtained from Kohn-Sham DFT and FDE in terms of deformation densities, bond critical points, and the negative Laplacian of the electron density. Different kinetic-energy functionals as needed for the frozen-density embedding method are tested and compared to a purely electrostatic embedding. It is shown that FDE is able to reproduce the characteristics of the density in the bonding region even in systems such as the F-H-F(-) molecule, which contains one of the strongest hydrogen bonds. Basis functions on the frozen system are usually required to accurately reproduce the electron densities of supermolecular calculations. However, it is shown here that it is in general sufficient to provide just a few basis functions in the boundary region between the two subsystems so that the use of the full supermolecular basis set can be avoided. It also turns out that electron-density deformations upon bonding predicted by FDE lack directionality with currently available functionals for the nonadditive kinetic-energy contribution.     

Away from generalized gradient approximation: orbital-dependent exchange-correlation functionals

The local-density approximation of density functional theory (DFT) is remarkably accurate, for instance, for geometries and frequencies, and the generalized gradient approximations have also made bond energies quite reliable. Sometimes, however, one meets with failure in individual cases. One of the possible routes towards better functionals would be the incorporation of orbital dependence (which is an implicit density dependency) in the functionals. We discuss this approach both for energies and for response properties. One possibility is the use of the Hartree-Fock-type exchange energy expression as orbital-dependent functional. We will argue that in spite of the increasing popularity of this approach, it does not offer any advantage over Hartree-Fock for energies. We will advocate not to apply the separation of exchange and correlation, which is so ingrained in quantum chemistry, but to model both simultaneously. For response properties the energies and shapes of the virtual orbitals are crucial. We will discuss the benefits that Kohn-Sham potentials can offer which are derived from either an orbital-dependent energy functional, including the exact-exchange functional, or which can be obtained directly as orbital-dependent functional. We highlight the similarity of the Hartree-Fock and Kohn-Sham occupied orbitals and orbital energies, and the essentially different meanings the virtual orbitals and orbital energies have in these two models. We will show that these differences are beneficial for DFT in the case of localized excitations (in a small molecule or in a fragment), but are detrimental for charge-transfer excitations. Again, orbital dependency, in this case in the exchange-correlation kernel, offers a solution.     

Orbital-dependent correlation energy in density-functional theory based on a second-order perturbation approach: success and failure

We have developed a second-order perturbation theory (PT) energy functional within density-functional theory (DFT). Based on PT with the Kohn-Sham (KS) determinant as a reference, this new ab initio exchange-correlation functional includes an exact exchange (EXX) energy in the first order and a correlation energy including all single and double excitations from the KS reference in the second order. The explicit dependence of the exchange and correlation energy on the KS orbitals in the functional fits well into our direct minimization approach for the optimized effective potential, which is a very efficient method to perform fully self-consistent calculations for any orbital-dependent functionals. To investigate the quality of the correlation functional, we have applied the method to selected atoms and molecules. For two-electron atoms and small molecules described with small basis sets, this new method provides excellent results, improving both second-order Moller-Plesset expression and any conventional DFT results significantly. For larger systems, however, it performs poorly, converging to very low unphysical total energies. The failure of PT based energy functionals is analyzed, and its origin is traced back to near degeneracy problems due to the orbital- and eigenvalue-dependent algebraic structure of the correlation functional. The failure emerges in the self-consistent approach but not in perturbative post-EXX calculations, emphasizing the crucial importance of self-consistency in testing new orbital-dependent energy functionals.     

Parameters for excess electron transfer in DNA. Estimation using unoccupied Kohn-Sham orbitals and TD DFT

We show that the energetics and electronic couplings for excess electron transfer (EET) can be accurately estimated by using unoccupied Kohn-Sham orbitals (UKSO) calculated for neutral pi stacks. To assess the performance of different DFT functionals, we use MS-PT2 results for seven pi stacks of nucleobases as reference data. The DFT calculations are carried out by using the local spin density approximation SVWN, two generalized gradient approximation functionals BP86 and BLYP, and two hybrid functionals B3LYP and BH&HLYP. Best estimations within the UKSO approach are obtained by the B3LYP and SVWN methods. TD DFT calculations provide less accurate values of the EET parameters as compared with the UKSO data. Also, the excess charge distribution in the radical anions is well described by the LUMOs of neutral systems. In contrast, spin-unrestricted DFT calculations of radical anions considerably overestimate delocalization of the excess electron. The excellent results obtained for the ground and excited states of the radical anions (excitation energy, transition dipole moment, electronic coupling, and excess electron distribution) by using UKSO of neutral dimers suggest an efficient strategy to calculate the EET parameters for DNA pi stacks.     

Time-dependent density-functional theory calculations of triplet-triplet absorption

We present density-functional theory calculations of triplet-triplet absorption by three different approaches based on time-dependent density-functional theory (DFT): unrestricted DFT linear response, open-shell restricted DFT linear response applied to the triplet state, and quadratic response with triplet excitations applied to the ground state. Comparison is also made with corresponding results obtained by Hartree-Fock and multiconfiguration self-consistent-field response theory. Two main conclusions concerning triplet-triplet transitions are drawn in this study: First, the very good agreement between unrestricted and restricted DFT results indicates that spin contamination of the triplet state is not a serious problem when computing triplet-triplet spectra of common organic molecules. Second, DFT response calculations of triplet-triplet transitions can be affected by triplet instability problems, especially for the combination of DFT quadratic response with functionals containing fractional exact Hartree-Fock exchange.     

Spin-spin contributions to the zero-field splitting tensor in organic triplets, carbenes and biradicals-a density functional and ab initio study

An evaluation study for the direct dipolar electron spin-spin (SS) contribution to the zero-field splitting (ZFS) tensor in electron paramagnetic resonance (EPR) spectroscopy is presented. Calculations were performed on a wide variety of organic systems where the SS contribution to the ZFS dominates over the second-order spin-orbit coupling (SOC) contribution. Calculations were performed using (hybrid) density functional theory (DFT), as well as complete active space self-consistent field (CASSCF) wave functions. In the former case, our implementation is an approximation, because we use the two-particle reduced spin-density matrix of the noninteracting reference system. In the latter case, the SS contribution is approximated by a mean-field method which, nevertheless, gives accurate results, compared to the approximation free computation of the SS part in a CASSCF framework. For the case of the triplet dioxygen molecule, it was shown that restricted open-shell density functional theory (RODFT), as well as CASSCF, can provide accurate spin-spin couplings while spin-unrestricted DFT leads to much larger errors. Furthermore, 15 organic radicals, including several 1,3 and 1,5 diradicals, dinitroxide biradicals, and even a chlorophyll a model system, were examined as test cases to demonstrate the accuracy and efficiency of our approach within a DFT framework. Accurate D values with root-mean-square deviations of 0.0035 cm(-1) were obtained. Furthermore, all trends, including those due to substituent effects, were correctly reproduced. In a different set of calculations, the polyacenes benzene, naphthalene, anthracene, and tetracene were studied. Applying DFT, the absolute D values were noticeably underestimated, but it was possible to correctly reproduce the trend to smaller D values with larger size of the systems. Finally, it was demonstrated that our approach is also well-suited for the study of carbenes. The smaller organic radicals of this work were also studied, through the use of CASSCF wave functions. This was a special advantage in the case of the triplet polyacenes, where the CASSCF approach gave better results than the DFT method. In comparing spin-restricted and spin-unrestricted results, it was shown through a natural orbital analysis and comparison to high-level ab initio calculations that even small amounts of spin polarization introduced by the unrestricted calculations lead to large deviations between the unrestricted Kohn-Sham (UKS) and restricted open-shell Kohn-Sham (ROKS) approaches. It is challenging to understand why the ROKS results show much better correlation with the experimental data.     

Natural transition orbitals for complex two-component excited state calculations

While the natural transition orbital (NTO) method has allowed electronic excitations from time-dependent Hartree-Fock and density functional theory to be viewed in a traditional orbital picture, the extension to multicomponent molecular orbitals such as those used in relativistic two-component methods or generalized Hartree-Fock (GHF) or generalized Kohn-Sham (GKS) is less straightforward due to mixing of spin-components and the inherent inclusion of spin-flip transitions in time-dependent GHF/GKS. An extension of single-component NTOs to the two-component framework is presented, in addition to a brief discussion of the practical aspects of visualizing two-component complex orbitals. Unlike the single-component analog, the method explicitly describes the spin and frequently obtains solutions with several significant orbital pairs. The method is presented using calculations on a mercury atom and a CrO₂Cl₂ complex.     

Optimized effective potential method for individual low-lying excited states

This paper presents an optimized effective potential (OEP) approach based on density functional theory (DFT) for individual excited states that implements a simple method of taking the necessary orthogonality constraints into account. The amended Kohn-Sham (KS) equations for orbitals of excited states having the same symmetry as the ground one are proposed. Using a variational principle with some orthogonality constraints, the OEP equations determining a local exchange potential for excited states are derived. Specifically, local potentials are derived whose KS determinants minimize the total energies and are simultaneously orthogonal to the determinants for states of lower energies. The parametrized form of an effective DFT potential expressed as a direct mapping of the external potential is used to simplify the OEP integral equations. A performance of the presented method is examined by exchange-only calculations of excited state energies for simple atoms and molecules.     

Cluster Amplitudes and Their Interplay with Self-Consistency in Density Functional Methods

Density functional theory (DFT) provides convenient electronic structure methods for the study of molecular systems and materials. Regular Kohn-Sham DFT calculations rely on unitary transformations to determine the ground-state electronic density, ground state energy, and related properties. However, for dissociation of molecular systems into open-shell fragments, due to the self-interaction error present in a large number of density functional approximations, the self-consistent procedure based on the this type of transformation gives rise to the well-known charge delocalization problem. To avoid this issue, we showed previously that the cluster operator of coupled-cluster theory can be utilized within the context of DFT to solve in an alternative and approximate fashion the ground-state self-consistent problem. This work further examines the application of the singles cluster operator to molecular ground state calculations. Two approximations are derived and explored: i) A linearized scheme of the quadratic equation used to determine the cluster amplitudes. ii) The effect of carrying the calculations in a non-self-consistent field fashion. These approaches are found to be capable of improving the energy and density of the system and are quite stable in either case. The theoretical framework discussed in this work could be used to describe, with an added flexibility, quantum systems that display challenging features and require expanded theoretical methods.     

Increasing the applicability of density functional theory. II. Correlation potentials from the random phase approximation and beyond

Density functional theory (DFT) results are mistrusted at times due to the presence of an unknown exchange correlation functional, with no practical way to guarantee convergence to the right answer. The use of a known exchange correlation functional based on wave-function theory helps to alleviate such mistrust. The exchange correlation functionals can be written exactly in terms of the density-density response function using the adiabatic-connection and fluctuation-dissipation framework. The random phase approximation (RPA) is the simplest approximation for the density-density response function. Since the correlation functional obtained from RPA is equivalent to the direct ring coupled cluster doubles (ring-CCD) correlation functional, meaning only Coulomb interactions are included, one can bracket RPA between many body perturbation theory (MBPT)-2 and CCD with the latter having all ring, ladder, and exchange contributions. Using an optimized effective potential strategy, we obtain correlation potentials corresponding to MBPT-2, RPA (ring-CCD), linear-CCD, and CCD. Using the suitable choice of the unperturbed Hamiltonian, Kohn-Sham self-consistent calculations are performed. The spatial behavior of the resulting potentials, total energies, and the HOMO eigenvalues are compared with the exact values for spherical atoms. Further, we demonstrate that the self-consistent eigenvalues obtained from these consistent potentials used in ab initio dft approximate all principal ionization potentials as demanded by ionization potential theorem.     

Performance of multi-configurational calculations for a 1,4-bis(phenylethynyl)benzene derivative conjugated molecule.

The theoretical challenge of finding a single method that quantitatively reproduces both the experimental low-lying excitation energies and the torsional barrier of a prototypical conjugated molecule, which could act as a molecular wire, has been addressed here. The results indicate that this goal can be reasonably achieved when multi-reference perturbation theory up to second order (MRMP2) based on a complete active space self-consistent field (CASSCF) wave function using large active spaces is used. The results obtained were also used to compare with less expensive Kohn-Sham (KS) density functional theory (DFT) calculations when applied to these properties. The results obtained with BLYP and B3LYP exchange-correlation functionals indicate that quantitative agreement with all the experimental data cannot be obtained with this methodology, with a clear dependence on the exchange-correlation form selected. We thus encourage a careful testing of pure and hybrid density functionals whenever KS DFT is used for the rational design of conjugated materials for charge conduits.