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.     

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.     

The spin-unrestricted molecular Kohn-Sham solution and the analogue of Koopmans's theorem for open-shell molecules

Spin-unrestricted Kohn-Sham (KS) solutions are constructed from accurate ab initio spin densities for the prototype doublet molecules NO(2), ClO(2), and NF(2) with the iterative local updating procedure of van Leeuwen and Baerends (LB). A qualitative justification of the LB procedure is given with a "strong" form of the Hohenberg-Kohn theorem. The calculated energies epsilon(isigma) of the occupied KS spin orbitals provide numerical support to the analogue of Koopmans' theorem in spin-density functional theory. In particular, the energies -epsilon(ibeta) of the minor spin (beta) valence orbitals of the considered doublet molecules correspond fairly well to the experimental vertical ionization potentials (VIPs) I(i) (1) to the triplet cationic states. The energy -epsilon(Halpha) of the highest occupied (spin-unpaired) alpha orbital is equal to the first VIP I(H) (0) to the singlet cationic state. In turn, the energies -epsilon(ialpha) of the major spin (alpha) valence orbitals of the closed subshells correspond to a fifty-fifty average of the experimental VIPs I(i) (1) and I(i) (0) to the triplet and singlet states. For the Li atom we find that the exact spin densities are represented by a spin-polarized Kohn-Sham system which is not in its ground state, i.e., the orbital energy of the lowest unoccupied beta spin orbital is lower than that of the highest occupied alpha spin orbital ("a hole below the Fermi level"). The addition of a magnetic field in the -z direction will shift the beta levels up so as to restore the Aufbau principle. This is an example of the nonuniqueness of the mapping of the spin density on the KS spin-dependent potentials discussed recently in the literature. The KS potentials may no longer go to zero at infinity, and it is in general the differences nu(ssigma)( infinity )-epsilon(isigma) that can be interpreted as (averages of) ionization energies. In total, the present results suggest the spin-unrestricted KS theory as a natural one-electron independent-particle model for interpretation and assignment of the experimental photoelectron spectra of open-shell molecules.     

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.     

Singlet-triplet energy gaps for diradicals from fractional-spin density-functional theory

Open-shell singlet diradicals are difficult to model accurately within conventional Kohn-Sham (KS) density-functional theory (DFT). These methods are hampered by spin contamination because the KS determinant wave function is neither a pure spin state nor an eigenfunction of the S(2) operator. Here we present a theoretical foray for using single-reference closed-shell ground states to describe diradicals by fractional-spin DFT (FS-DFT). This approach allows direct, self-consistent calculation of electronic properties using the electron density corresponding to the proper spin eigenfunction. The resulting FS-DFT approach is benchmarked against diradical singlet-triplet gaps for atoms and small molecules. We have also applied FS-DFT to the singlet-triplet gaps of hydrocarbon polyacenes.     

Non‐empirical calculations of NMR indirect carbon–carbon coupling constants: 1. Three‐membered rings

The carbon–carbon indirect nuclear spin–spin coupling constants in cyclopropane, aziridine and oxirane were investigated by means of ab initio calculations at the RPA, SOPPA and DFT/B3LYP levels. We found that the carbon–carbon couplings are by far dominated by the Fermi contact term. Our best SOPPA and DFT results are in a very good agreement with each other and with the experimental values, whereas calculations at the RPA level of theory strongly overestimate the carbon–carbon couplings. Significant differences in the basis set dependence of the calculated carbon–carbon coupling constants obtained with either wavefunction method, RPA or SOPPA, or the density functional method, DFT/B3LYP, are observed. The SOPPA results depend much more strongly on the quality of the basis set than the results of DFT/B3LYP calculations. The medium‐sized core‐valence basis sets cc‐pCVTZ and even cc‐pCVDZ were found to perform fairly well at the SOPPA level for the one‐bond carbon–carbon couplings investigated here. Copyright © 2002 John Wiley & Sons, Ltd.     

Extended treatment of charge response kernel comprising the density functional theory and charge regulation procedures

We propose an extended treatment of the charge response kernel (CRK), (partial differential Q(a)/partial differential V(b)), which describes the response of partial charges on atomic sites to external electrostatic potential, on the basis of the density functional theory (DFT) via the coupled perturbed Kohn-Sham equations. The present CRK theory incorporates regulation procedures in the definition of partial charges to avoid unphysical large fluctuation of the CRK on "buried" sites. The CRKs of some alcohol and organic molecules, methanol, ethanol, propanol, butanol, dimethylsulfoxide (DMSO), and tetrahydrofuran (THF) were calculated, demonstrating that the new CRK model at the DFT level has greatly improved the performance of accuracy in comparison with that at the Hartree-Fock level previously proposed. The CRK model was also applied to investigate spatial nonlocality of the charge response through alkyl chain sequences. The CRK model at the DFT level enables us to construct a nonempirical strategy for polarizable molecular modeling, with practical reliability and robustness.     

Orbital-corrected orbital-free density functional theory

A new implementation of density functional theory (DFT), namely orbital-corrected orbital-free (OO) DFT, has been developed. With at most two non-self-consistent iterations, OO-DFT accomplishes the accuracy comparable to fully self-consistent Kohn-Sham DFT as demonstrated by its application on the cubic-diamond Si and the face-centered-cubic Ag systems. Our work provides a new impetus to further improve orbital-free DFT method and presents a robust means to significantly lower the cost associated with general applications of linear-scaling Kohn-Sham DFT methods on large systems of thousands of atoms within different chemical bonding environment.     

Multireference spin-adapted variant of density functional theory

A new Kohn-Sham formalism is developed for studying the lowest molecular electronic states of given space and spin symmetry whose densities are represented by weighted sums of several reference configurations. Unlike standard spin-density functional theory, the new formalism uses total spin conserving spin-density operators and spin-invariant density matrices so that the method is fully spin-adapted and solves the so-called spin-symmetry dilemma. The formalism permits the use of an arbitrary set of reference (noninteracting) configurations with any number of open shells. It is shown that the requirement of degeneracy of the total noninteracting energies of the reference configurations (or configuration state functions) is equivalent to the stationary condition of the exact energy relative to the weights of the configurations (or configuration state functions). Consequently, at any molecular geometry, the weights can be determined by minimization of the energy, and, for given reference weights, the Kohn-Sham orbitals can be determined. From this viewpoint, the developed theory can be interpreted as an analog of the multiconfiguration self-consistent field approach within density functional theory.     

Exact-exchange density functional theory for hyperpolarizabilities

Time-dependent density functional theory (TDDFT) employing the exact-exchange functional (TDDFTx) has been formulated using the optimized effective potential method for the beta static hyperpolarizabilities, where it reduces to coupled-perturbed Kohn-Sham theory. A diagrammatic technique is used to take the functional derivatives for the derivation of the adiabatic second kernel, which is required for the analytical calculation of the beta static hyperpolarizabilities with DFT. The derived formulas have been implemented using Gaussian basis sets. The structure of the adiabatic exact-exchange second kernel is described and numerical examples are presented. It is shown that no current DFT functional satisfies the correct properties of the second kernel. Not surprisingly, TDDFTx, which corrects the self-interaction error in standard DFT methods and has the correct long-range behavior, provides results close to those of time-dependent Hartree-Fock in the static limit.     

新一代密度泛函方法XYG3

大量的应用研究测试表明,以B3LYP为代表的传统密度泛函方法在反应能垒以及非键相互作用等重要性质的预测上存在困难,并且预测精度随研究体系的增大而不断变差。开发越来越精确的交换相关泛函是现代密度泛函理论发展的主线。近年来,引入未占轨道信息的新一代双杂化密度泛函方法的研究受到越来越多的关注。本篇综述回顾了这一领域已取得的一些进展;推导了双杂化泛函方法在Kohn-Sham密度泛函理论框架下的理论基础;依据各自特点,将现已提出的双杂化泛函划分为三种类别。本文着重测评了以XYG3为代表的一类双杂化泛函的表现。最后,就双杂化泛函的未来发展方向,提出了一些设想与建议。     

Calculation of nonadiabatic couplings with restricted open-shell Kohn-Sham density-functional theory

This paper generalizes the recently proposed approaches for calculating the derivative couplings between adiabatic states in density-functional theory (DFT) based on a Slater transition-state density to transitions such as singlet-singlet excitations, where a single-determinant ansatz is insufficient. The proposed approach is based on restricted open-shell Frank et al. [J. Chem. Phys. 108, 4060 (1998)] theory used to describe a spin-adapted Slater transition state. To treat the dependence of electron-electron interactions on the nuclear positions, variational linear-response density-functional perturbation theory is generalized to reference states with an orbital-dependent Kohn-Sham Hamiltonian and nontrivial occupation patterns. The methods proposed in this paper are not limited to the calculation of derivative coupling vectors, but can also be used for the calculation of other transition matrix elements. Moreover, they can be used to calculate the linear response of open-shell systems to arbitrary external perturbations in DFT.     

Linear response and quasiparticle calculations as probes of the Kohn-Sham eigenvalues in metals

The Kohn-Sham eigenvalues were formally introduced into density functional theory as Lagrange multipliers in the implementation of the minimum principle for the total energy of a many-electron system. No general results are available concerning the physical significance of these one-electron eigenvalues (with the exception of the highest occupied level, which equals the Fermi energy). Recent ab initio calculations of dynamical response in metals make explicit use of the Kohn-Sham band structure, and associated wave functions, through the use of spectral representations. This opens up the possibility of examining the significance of the eigenvalues at an “empirical” level, i.e., through direct comparison with the results of spectroscopic measurements. A particularly interesting example is afforded by new inelastic x-ray scattering experiments on A1. For a special wave vector transfer, q_o ≈︂ 1.5k_F, the measured spectrum provides a direct mapping of the Kohn-Sham noninteracting spectrum. For a range of wave vectors about q_o, the bare Kohn-Sham spectrum still reproduces all the main features of the measurements; this suggests that, in this metal, the Kohn-Sham eigenvalues are good approximations to the quasiparticle energies. We also discuss the interplay between Kohn-Sham bands and the energy of the “anomalous” plasmon in Cs, whose dispersion bears a signature of the excited-state band structure. Finally, and in a more formal framework, we outline the results of a first-principles comparison between quasiparticle amplitudes and Kohn-Sham wave functions at a jellium surface; the latter turn out to be excellent approximations to the former. © 1996 John Wiley & Sons, Inc.     

Evaluation of <S2> in density functional theory

The evaluation of in density functional theory (DFT) is considered. Wang et al. [J. Chem. Phys. 102, 3477 (1995)] have derived an approximate, local density expression for and in the present study their formula is evaluated using densities from unrestricted Hartree-Fock (UHF) and a range of DFT exchange-correlation functionals. The results are compared with those obtained by evaluating the conventional UHF expression using the Kohn-Sham orbitals, which is appropriate for the noninteracting system. A generalized gradient approximation for is then proposed and investigated.     

Intracule densities in the strong-interaction limit of density functional theory

The correlation energy in density functional theory can be expressed exactly in terms of the change in the probability of finding two electrons at a given distance r(12) (intracule density) when the electron-electron interaction is multiplied by a real parameter lambda varying between 0 (Kohn-Sham system) and 1 (physical system). In this process, usually called adiabatic connection, the one-electron density is (ideally) kept fixed by a suitable local one-body potential. While an accurate intracule density of the physical system can only be obtained from expensive wavefunction-based calculations, being able to construct good models starting from Kohn-Sham ingredients would highly improve the accuracy of density functional calculations. To this purpose, we investigate the intracule density in the lambda --> infinity limit of the adiabatic connection. This strong-interaction limit of density functional theory turns out to be, like the opposite non-interacting Kohn-Sham limit, mathematically simple and can be entirely constructed from the knowledge of the one-electron density. We develop here the theoretical framework and, using accurate correlated one-electron densities, we calculate the intracule densities in the strong interaction limit for few atoms. Comparison of our results with the corresponding Kohn-Sham and physical quantities provides useful hints for building approximate intracule densities along the adiabatic connection of density functional theory.     

Density functional restricted-unrestricted approach for nonlinear properties: application to electron paramagnetic resonance parameters of square planar copper complexes

The density functional restricted-unrestricted approach for treatments of spin polarization effects in molecular properties using spin restricted Kohn-Sham theory has been extended from linear to nonlinear properties. It is shown that the spin polarization contribution to a nonlinear property has the form of a quadratic response function that includes the zero-order Kohn-Sham operator, in analogy to the lower order case where the spin polarization correction to an expectation value has the form of a linear response function. The developed approach is used to formulate new schemes for computation of electronic g-tensors and hyperfine coupling constants, which include spin polarization effects within the framework of spin restricted Kohn-Sham theory. The proposed computational schemes are in the present work employed to study the spin polarization effects on electron paramagnetic resonance spin Hamiltonian parameters of square planar copper complexes. The obtained results indicate that spin polarization gives rise to sizable contributions to the hyperfine coupling tensor of copper in all investigated complexes, while the electronic g-tensors of these complexes are only marginally affected by spin polarization and other factors, such as choice of exchange-correlation functional or molecular structures, will have more pronounced impact on the accuracy of the results.