Regional self-interaction correction of density functional theory

We propose a new simple scheme for self-interaction correction (SIC) of exchange functionals in the density functional theory. In the new scheme, exchange energies are corrected by substituting exchange self-interactions for exchange functionals in regions of self-interaction. To classify the regions of self-interaction, we take advantage of the property of the total kinetic energy density approaching the Weizs?cker density in the case of electrons in isolated orbitals. The scheme differs from conventional SIC methods in that it produces optimized molecular structures. Applying the scheme to the calculation of reaction energy barriers showed that it provides a clear improvement in cases where the barriers are underestimated by conventional "pure" functionals. In particular, we found that this scheme even reproduces a transition state that is not given by pure functionals.     

Theoretical calculation of bond dissociation energies and heats of formation for nitromethane and polynitromethanes with density functional theory

The C? NO₂ bond dissociation energies (BDEs) and the heats of formation (HOFs) of nitromethane and polynitromethanes (dinitromethane, trinitromethane, and tetranitromethane) system in gas phase at 298.15 K were calculated theoretically. Density functional theory (DFT) B3LYP, B3P86, B3PW91, and PBE0 methods in combination with different basis sets were employed. It was found that the C? NO₂ bond BDEs can be improved from B3LYP to B3PW91 to B3P86 or PBE0 functional. Levels of theory employing B3P86 and PBE0 functionals were found to be sufficiently reliable without the presence of diffusion functions. As the number of NO₂ groups on the same C atom increases, the PBE0 functional performs better than the B3P86 functional. Regarding the calculated HOFs, all four functionals can yield satisfactory results with deviations of <2 kcal mol^?1 from experimental ones for CH₂(NO₂)₂ and CH(NO₂)₃, when the diffusion functions are not augmented. For the C(NO₂)₄ molecule, the large basis sets augmented with polarization functions and diffusion functions are required to yield a good result. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Rapidly convergent procedure to solve the density profile equation in the classical density functional theory

An efficient recursive procedure to solve the density profile equation in the classical density functional theory (DFT) using an inverse Broyden method is described. The present iterative procedure is free of calculation of the Jacobian matrix, and its inversion unavoidable for the well-known Newton-Raphson (NR) method and its variants. Numerical calculation indicates that only the approximate solution and iterative matrix of the lower bulk density case are employed as the corresponding initial guesses of the higher bulk density case, the present recursive procedure can converge quickly to the physical solution with an accuracy of epsilon = 10(-14); therefore, the procedure provides an efficient numerical algorithm for the theory in which acquirement of a density profile of high accuracy is a key step. Extensive numerical calculation shows the advantage of the present inverse Broyden method over Broyles' mixing procedure and a modified Powell hybrid algorithm (a variation of the NR method).     

A density functional theory study on the electrocyclization of 1,2,4,6-heptatetraene analogues: converting a pericyclic to a pseudopericyclic reaction

A comprehensive B3LYP/6-31+G* study on the electrocyclization of 1,2,4,6-heptatetraene analogues was conducted. Starting from the cyclization of (2Z)-2,4,5-hexatrienal, a pericyclic disrotatory process favored by the assistance of a electron lone pair, we incorporated small modifications in its molecular structure to obtain a truly pseudopericyclic process. To this purpose electronegative atoms (fluorine and nitrogen) were added to give a more electrophilic character on the carbon atom which is attacked by the electron lone pair of the oxygen atom. The complete pathway for each reaction was determined, and changes in magnetic properties were monitored with a view to estimating the aromatization associated with each process. This information, together with the energetic and structural results, allowed us to classify the reactions as pseudopericyclic or pericyclic. Among all studied reactions only one was a truly pseudopericyclic process and another was a borderline case. The features of this unequivocally pseudopericyclic case were analyzed in depth.     

Application of the Sakurai-Sugiura projection method to core-excited-state calculation by time-dependent density functional theory

The Sakurai-Sugiura projection (SS) method was implemented and numerically assessed for diagonalization of the Hamiltonian in time-dependent density functional theory (TDDFT). Since the SS method can be used to specify the range in which the eigenvalues are computed, it may be an efficient tool for use with eigenvalues in a particular range. In this article, the SS method is applied to core excited calculations for which the eigenvalues are located within a particular range, since the eigenvalues are unique to atomic species in molecules. The numerical assessment of formaldehyde molecule by TDDFT with core-valence Becke's three-parameter exchange (B3) plus Lee-Yang-Parr (LYP) correlation (CV-B3LYP) functional demonstrates that the SS method can be used to selectively obtain highly accurate eigenvalues and eigenvectors. Thus, the SS method is a new and powerful alternative for calculating core-excitation energies without high computation costs.     

Orbital functional theory of linear response and excitation

Orbital functional theory (OFT) is based on a rule that determines a single‐determinant reference state Φ for any exact N‐electron eigenstate Ψ. An OFT model postulates an explicit correlation energy functional E_c of occupied orbital functions {?_i} and occupation numbers {n_i}. The orbital Euler–Lagrange equations are analogous to Kohn–Sham equations, but do not in general contain local potential functions. Time‐dependent Hartree–Fock theory is generalized in OFT to a formally exact linear response theory that includes electronic correlation. In the exchange‐only limit, the theory reduces to the random‐phase approximation of many‐body theory. The formalism determines excitation energies. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

密度泛函理论研究香草醛与氨基硫脲的加成-缩合反应

针对氨基硫脲和香草醛亲核加成-缩合反应的反应物和生成物,用量子化学密度泛函方法在B3LYP/6-31G水平上进行几何优化,通过Mulliken电荷分析推测出香草醛具有较高的亲核加成反应活性.对优化后的构型进行振动分析,得到不同温度下反应物和生成物的热力学性质,据此计算出相应温度下的反应焓变、吉布斯自由能变和平衡常数.实验表明,在298K～1000K内,氨基硫脲与香草醛的亲核加成-缩合反应是一个放热反应,吉布斯自由能变成负值,平衡常数很大,说明反应在温和条件下即能自发进行.在此基础上采用含时密度泛函方法(TD-DFT)计算了分子激发态的电子跃迁能,得到对应激发态的吸收波长,所得结果与文献值基本吻合.     

Exploring odd-even effects of simple oligomer-like DRCNnT series: a study based on density functional theory/time-dependent density functional theory calculations

Odd-even effects of short-circuit current density and power conversion efficiency (PCE) are an interesting phenomenon in some organic solar cells. Although some explanations have been given, why they behave in such a way is still an open question. In the present work, we investigate a set of acceptor-donor-acceptor simple oligomer-like small molecules, named the DRCNnT (n = 5-9) series, to give an insight into this phenomenon because the solar cells based on them have high PCE (up to 10.08%) and show strong odd-even effects in experiments. By modeling the DRCNnT series and using density functional theory, we have studied the ground-state electronic structures of the DRCNnT (n = 5-9) series in condensed phase. The calculated results reproduce the experimental trends well. Furthermore, we find that the exciton-binding energies of the DRCNnT series may be one of the key parameters to explain this phenomenon because they also show odd-even effects. In addition, by studying the effects of alkyl branch and terminal group on odd-even effects of dipole moment, we find that eliminating one or two alkyl branches does not break the odd-even effects of dipole moments, but eliminating one or two terminal groups does. Finally, we conclude that removing one alkyl branch close to the terminal group of DRCN5T can decrease highest occupied molecular orbital (HOMO) energy (thus increasing open circuit voltage) and increase dipole moment (thus enhancing charge separation and short-circuit current). This could be a new and simple method to increase the PCE of DRCN5T-based solar cells.     

Exchange and correlation energy in density functional theory

Several different versions of density functional theory (DFT) that satisfy Hohenberg–Kohn theorems are characterized by different definitions of a reference or model state determined by an N‐electron ground state. A common formalism is developed in which exact Kohn–Sham equations are derived for standard Kohn–Sham theory, for reference‐state density functional theory, and for unrestricted Hartree–Fock (UHF) theory considered as an exactly soluble model Hohenberg–Kohn theory. A natural definition of exchange and correlation energy functionals is shown to be valid for all such theories. An easily computed necessary condition for the locality of exchange and correlation potentials is derived. While it is shown that in the UHF model of DFT the optimized effective potential (OEP) exchange satisfies this condition by construction, the derivation shows that this condition is not, in general, sufficient to define an exact local exchange potential. It serves as a test to eliminate proposed local potentials that are not exact for ground states. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 77: 521–525, 2000     

Electronic structures of 4d transition metal monoxides by density functional theory

Bond distances, dissociation energies, ionization potentials and electron affinities of 4d transition metal monoxides from YO to CdO and their positive and negative ions were studied by use of density functional methods B3LYP, BLYP, B3PW91, BPW91, B3P86, BP86, SVWN, MPW1PW91 and PBE1PBE. It was found that calculated properties are highly dependent on the functionals employed, especially for dissociation energy. For most neutral species, pure density functionals BLYP, BPW91 and BP86 have good performance in predicting dissociation energy than hybrid density functionals B3LYP, B3PW91 and B3P86. In addition, BLYP gives the largest bond distance compared with other density functional methods, while SVWN gives shortest bond distance, largest dissociation energy and electron affinity. For the ground state, the spin multiplicity of the charged species can be obtained by ± 1 of their corresponding neutral species.     

黄酮类化合物的密度泛函理论研究

在混合密度泛函B3LYP理论下,用6-31G*基函数研究了几种典型黄酮类化合物分子的几何结构、电子结构和分子的静电势,讨论了电子结构和分子活性部位的关系.     

Differentiability of Lieb functional in electronic density functional theory

A solid understanding of the Lieb functional F_L is important because of its centrality in the foundations of electronic density functional theory. A basic question is whether directional derivatives of F_L at an ensemble‐V‐representable density are given by (minus) the potential. A widely accepted purported proof that F_L is Gâteaux differentiable at EV‐representable densities would say, “yes.” But that proof is fallacious, as shown here. F_L is not Gâteaux differentiable in the normal sense, nor is it continuous. By means of a constructive approach, however, we are able to show that the derivative of F_L at an EV‐representable density ρ₀ in the direction of ρ₁ is given by the potential if ρ₀ and ρ₁ are everywhere strictly greater than zero, and they and the ground state wave function have square integrable derivatives through second order. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

The use of the hybrid density functional theory (DFT) approach for calculation of vibrational spectra: nitromethane

The molecular geometry of nitromethane was optimized and its force field and vibrational spectrum were calculated by the BECKE3LYP method. The accuracy of optimization of the geometry of MeNO₂ obtained by this method using the 6–311G(d,p) and 6–311++G(d,p) basis sets is not poorer than that obtained at the second-order Møller-Plesset level of perturbation theory (MP2). The vibrational frequencies of nitromethane and its d₁, d₂, and d₃ isotopomers obtained by the BECKE3LYP method are in much better agreement with the experimental data than those calculated at the MP2 level using the same basis set. The average absolute error of calculations performed without the use of any scaling factors is ～2% for frequencies; the maximum deviation is ～4%.     

Evaluation of density functional theory in the bond rupture of octane

Density functional methods at the 6-31G* level are applied to the rupture of n-octane into methyl–heptyl, ethyl–hexyl, propyl–pentyl, and butyl–butyl radical fragments. The energetics of the radicals at UMP3, UMP2/6-31G*//UHF/6-31G* (hereafter referred to as UMP), are compared to UB3LYP/6-31G* results (referred to as UB). Although the UMP approach matches additivity energies to within 5 kcal/mol, it fails to mimic the overall energetic trend. The UB energies agree with additivity estimates and trends to within 1–2 kcal/mol and radical entropies deviate by only 2 e.u. from available experimental data. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 154–167, 1998     

Implementation of empirical dispersion corrections to density functional theory for periodic systems

A recently developed empirical dispersion correction (Grimme et al., J. Chem. Phys. 2010, 132, 154104) to standard density functional theory (DFT‐D3) is implemented in the plane‐wave program package VASP. The DFT‐D3 implementation is compared with an implementation of the earlier DFT‐D2 version (Grimme, J. Comput. Chem. 2004, 25, 1463; Grimme, J. Comput. Chem. 2006, 27, 1787). Summation of empirical pair potential terms is performed over all atom pairs in the reference cell and over atoms in shells of neighboring cells until convergence of the dispersion energy is obtained. For DFT‐D3, the definition of coordination numbers has to be modified with respect to the molecular version to ensure convergence. The effect of three‐center terms as implemented in the original molecular DFT‐D3 version is investigated. The empirical parameters are taken from the original DFT‐D3 version where they had been optimized for a reference set of small molecules. As the coordination numbers of atoms in bulk and surfaces are much larger than in the reference compounds, this effect has to be discussed. The results of test calculations for bulk properties of metals, metal oxides, benzene, and graphite indicate that the original parameters are also suitable for solid‐state systems. In particular, the interlayer distance in bulk graphite and lattice constants of molecular crystals is considerably improved over standard functionals. With the molecular standard parameters (Grimme et al., J. Chem. Phys. 2010, 132, 154104; Grimme, J. Comput. Chem. 2006, 27, 1787) a slight overbinding is observed for ionic oxides where dispersion should not contribute to the bond. For simple adsorbate systems, such as Xe atoms and benzene on Ag(111), the DFT‐D implementations reproduce experimental results with a similar accuracy as more sophisticated approaches based on perturbation theory (Rohlfing and Bredow, Phys. Rev. Lett. 2008, 101, 266106). © 2012 Wiley Periodicals, Inc.     

左旋葡聚糖热解机理的密度泛函理论研究

采用密度泛函理论B3LYP/6-31++G(d,p)方法,对纤维素热解的主要产物左旋葡聚糖的热解反应机理进行了理论计算分析,设计了四种可能的热解反应途径, 对各种反应的反应物、产物和过渡态的结构进行了能量梯度全优化。计算结果表明,左旋葡聚糖开环成链状中间体时,首先,左旋葡聚糖中的两个半缩醛键C(1)-O(7)和C(6)-O(8)断裂,经过渡态TS₁形成中间体IM₁,同时,C(6)-O(7)结合成键使C(5)-C(6)-O(7)形成环状结构,该反应的能垒较高,为296.53 kJ/mol,然后IM₁经过渡态TS₂转变为中间体IM₂,该反应的能垒为234.09 kJ/mol;对IM₂设计了四条可能的反应路径,反应路径2和3能垒较低,是IM₂最可能的热解反应途径;在反应路径1和4中都包含了脱羰基反应,其反应能垒较高,不易发生。     

The formulation of a self-consistent constricted variational density functional theory for the description of excited states

We outline here a self-consistent approach to the calculation of transition energies within density functional theory. The method is based on constricted variational theory (CV-DFT). It constitutes in the first place an improvement over a previous scheme [T. Ziegler, M. Seth, M. Krykunov, J. Autschbach, F. Wang, Chem. Phys. 130 (2009) 154102] in that it includes terms in the variational parameters to any desired order n including n = ∞. For n = 2, CV(n)-DFT is similar to TD-DFT. Adiabatic TD-DFT becomes identical to CV(2)-DFT after the Tamm-Dancoff approximation is applied to both theories. We have termed the new scheme CV(n)-DFT. In the second place, the scheme can be implemented self-consistently, SCF-CV(n)-DFT. The procedure outlined here could also be used to formulate a SCF-CV(n) Hartree-Fock theory. The approach is further kindred to the ΔSCF-DFT procedures predating TD-DFT and we describe how adiabatic TD-DFT and ΔSCF-DFT are related through different approximations to SCF-CV(n)-DFT.     

氢化可的松和表氢化可的松的密度泛函理论研究

采用密度泛函理论B3LYP/6-31G和B3LYP/6-311G*方法优化了氢化可的松和表氢化可的松的几何结构,利用优化的结构得到了氢化可的松和表氢化可的松的原子净电荷、总能量及前沿分子轨道组成.基于简谐振动分析求得了氢化可的松和表氢化可的松的红外光谱频率和强度,由统计热力学分析得到了热力学函数;进而确定了氢化可的松和...