Molecular electronegativity in density functional theory (VI) --Atom-bond electronegativity equalization model

Ｒｅｃｅｎｔｌｙ,ｄｅｖｅｌｏｐｍｅｎｔａｎｄａｐｐｌｉｃａｔｉｏｎｏｆｄｅｎｓｉｔｙｆｕｎｃｔｉｏｎａｌｔｈｅｏｒｙ(ＤＦＴ)ａｒｅｑｕｉｔｅａｔｔｒａｃｔｉｖｅｔｏｓｃｉｅｎｔｉｓｔｓ’ａｔｔｅｎｔｉｏｎ.Ｔｈｅｒｅｌｅｖａｎｔｐａｐｅｒｓａｒｅｉｎｃｒｅａｓｉｎｇｖｅｒｙｑｕｉｃｋｌｙ[1—11].ＳｉｎｃｅｔｈｅｅｌｅｃｔｒｏｎｅｇａｔｉｖｉｔｙｃｏｎｃｅｐｔｗａｓｐｒｏｐｏｓｅｄｂｙＰａｕｌｉｎｇ[12]ｔｏｄｅｓｃｒｉｂｅｔｈｅｐｏｗｅｒｏｆａｎａｔｏｍｉｎａｍｏｌｅｃｕｌｅｔｏａｔｔｒａｃｔｅｌ…     

Molecular electronegativity in density functional theory (II) --Direct calculation of group electronegativity and the atomic charges in a group

On the basis of a more precise expression of the atomic effective electronegativity deduced from the density functional theory and electronegativity equalization principle, a new scheme for calculating the group electronegativity and the atomic charges in a group is proposed and programed, and various parameters of electronegativity and hardness are given for some common atoms. Through calculation, analysis and comparison of more than one hundred groups, it is shown that the results from this scheme are reasonable and may be extended.     

Studies of chemical hardness and Fukui function using the exact solution of the density functional theory

Neal's procedure has been applied to determine the electron density ρ(x) for the H₂ molecule. The chemical hardness has been calculated employing the ab initio and density functional theory methods and the values are found to be reasonably good. The principle of maximum hardness (PMH) was tested. Fukui functions and the distribution of electron density around the internuclear distance were studied employing the electron density of the H₂ molecule. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 4–10, 2001     

Electrophilicity kernel and its hierarchy through softness in conceptual density functional theory

In the context of long‐range density functional theory of softness kernel, the concepts of local and kernel electrophilicity are developed to obtain the local to global hierarchical criteria such as bilocal symmetry, asymptotic behavior, and integral local to global relationships. Further development into potential‐density dependence is also provided while considering the local plus nonlocal specification of softness kernel and the associated atomic scales for the derived kernel; local and global electrophilicity are provided whose reliability is judged through the periodical features along the periodic table especially for the alkali metal atoms and the halogen atoms. © 2013 Wiley Periodicals, Inc.     

Assessing the role of Hartree‐Fock exchange,correlation energy and long range corrections in evaluating ionization potential,and electron affinity in density functional theory

Accurate determination of ionization potentials (IPs), electron affinities (EAs), fundamental gaps (FGs), and HOMO, LUMO energy levels of organic molecules play an important role in modeling and predicting the efficiencies of organic photovoltaics, OLEDs etc. In this work, we investigate the effects of Hartree Fock (HF) Exchange, correlation energy, and long range corrections in predicting IP and EA in Hybrid Functionals. We observe increase in percentage of HF exchange results in increase of IPs and decrease in EAs. Contrary to the general expectations inclusion of both HF exchange and correlation energy (from the second order perturbation theory MP2) leads to poor prediction. Range separated Hybrid Functionals are found to be more reliable among various DFT Functionals investigated. DFT Functionals predict accurate IPs whereas post HF methods predict accurate EAs. © 2017 Wiley Periodicals, Inc.     

原子-键电负性均衡方法参数的确定与探讨

原子-键电负性均衡方法(ABEEM)是以密度泛函理论(DFT)和电负性均衡原理为基础发展而来，它明确地考虑了化学键是不引入任何实验数据的带纯理论性和计算的方法．使用统一标准并具有代表性和全面性地选择了200多个模型分子，利用可得到较准确结构的MP2／6-31G*优化结构，心／STO-3G单点计算得到Mulliken重叠布居，再用最小二乘法拟合得到许多主族元素在分子体系中的诸原子(包括单、双和叁键等不同成键状态)和化学键的ABEEM参数．所得到的原子的价态电负性可与已提出的其他电负性标度相比拟，计算CO得到的电荷负端为C(与从头计算的结果相反)，结果与实验相符，且原子电荷的正负不完全由原子电负性决定．     

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.     

Analytical gradients for subsystem density functional theory within the slater‐function‐based amsterdam density functional program

We present a new implementation of analytical gradients for subsystem density‐functional theory (sDFT) and frozen‐density embedding (FDE) into the Amsterdam Density Functional program (ADF). The underlying theory and necessary expressions for the implementation are derived and discussed in detail for various FDE and sDFT setups. The parallel implementation is numerically verified and geometry optimizations with different functional combinations (LDA/TF and PW91/PW91K) are conducted and compared to reference data. Our results confirm that sDFT‐LDA/TF yields good equilibrium distances for the systems studied here (mean absolute deviation: 0.09 Å) compared to reference wave‐function theory results. However, sDFT‐PW91/PW91k quite consistently yields smaller equilibrium distances (mean absolute deviation: 0.23 Å). The flexibility of our new implementation is demonstrated for an HCN‐trimer test system, for which several different setups are applied. © 2016 Wiley Periodicals, Inc.     

PyCDFT: A Python package for constrained density functional theory

We present PyCDFT, a Python package to compute diabatic states using constrained density functional theory (CDFT). PyCDFT provides an object-oriented, customizable implementation of CDFT, and allows for both single-point self-consistent-field calculations and geometry optimizations. PyCDFT is designed to interface with existing density functional theory (DFT) codes to perform CDFT calculations where constraint potentials are added to the Kohn–Sham Hamiltonian. Here, we demonstrate the use of PyCDFT by performing calculations with a massively parallel first-principles molecular dynamics code, Qbox, and we benchmark its accuracy by computing the electronic coupling between diabatic states for a set of organic molecules. We show that PyCDFT yields results in agreement with existing implementations and is a robust and flexible package for performing CDFT calculations. The program is available at https://dx.doi.org/10.5281/zenodo.3821097 .     

Correlation of bond metallicity measures to electronegativity for binary oxides

We have investigated alkali, alkaline‐earth, and rutile binary oxides within density functional theory (DFT) and Bader's atoms‐in‐molecules theory, focusing on properties of bond and ring critical points, and their relations to band gap and Pauling electronegativity. We find linear relations of kinetic energy density, electron density, and the gap divided by kinetic energy density at the bond critical points to the difference of Pauling electronegativities of the cation and oxygen anion. At the ring critical points of rutile compounds, we also find that some bond metallicity measures are linearly related to the difference of electronegativities. This study extends our knowledge about the relations between bond critical points, band gap, and electronegativity, but also shows for the first time a quantitative relation between quantities at the ring critical points and global properties of the compounds.     

Structure and stability of XeF6 isomers: density functional theory study and the maximum hardness principle

The geometrical configurations of the four possible isomers with C_3v, O_h, C_s and C_2v symmetry on the potential energy surface of the XeF₆ molecule are optimized by using DFT-LDA/NL. Their relative energies, vibration frequencies, electronic chemical potential and hardness have been calculated. It is found that the C_2v configuration has one imaginary frequency. The relative energies of the four isomers increase in order of C_3v, O_h, C_s and C_2v, and the hardness values in same order. The isomer stability obeys the maximum hardness principle (MHP), while their hardness values are very close to each other. It is quite evident that the very close hardness is the main reason for the structure fluxionality of XeF₆.     

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     

Acidity predictions in an electron propagator approach

An alternative way to calculate vertical ionization potentials (VIP) and vertical electron affinity (VEA) is the application of Koopman's theorem, using the electron propagator theory. In the present work, the results of the application of this theorem using the electron propagator formalism have been compared with the experiment in order to validate different basis set. Using the basis set with the best performance, the acidity tendencies in some substituted acetic acid molecules have been analyzed by correlating the proton affinity (PA) with molecular electronegativity (χ) and hardness (η); these last indexes were obtained from the calculated VIP and VEA considering the finite difference approximation. The above correlations were compared with equivalent correlations using the energy of the frontier Hartree–Fock orbitals and the corresponding Kohn–Sham orbitals, which were calculated with the B3LYP‐DFT procedure. The results indicate that the electron propagator theory could be an interesting alternative to evaluate reactivity indexes, since this theory gives reliable values of VIP and VEA. It was also found that (i) the VIP values are very close to experiment, with only a 0.38% of error; (ii) acceptable results are inferred for VEA; (iii) a triple zeta quality function works quite well in these calculations, and particularly the 6‐311G(d,p) basis set is the best, as it had been reported; and (iv) using the depronation energy (DPE), good results were obtained in the correlations δDPE‐VEA and δDPE‐χ. The results tested that P3 approximation in the electron propagator approach can be a new and interesting alternative in predicting VIP, VEA, and some reactivity indexes, such as χ and η, at least for the compounds studied. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Electronic structure and spectroscopic properties of 6‐aminophenanthridine and its derivatives: Insights from density functional theory

6‐Aminophenanthridine (6AP) and its derivatives show important biological activities as antiprion compounds and inhibitors of the protein folding activity of the ribosome. Both of these activities depend on the RNA binding property of these compounds, which has been recently characterized by fluorescence spectroscopy. Hence, fundamental insights into the photophysical properties of 6AP compounds are highly important to understand their biological activities. In this work, we have calculated electronic structures and optical properties of 6AP and its three derivatives 6AP8CF₃, 6AP8Cl, and 6APi by density functional theory (DFT) and time‐dependent density functional theory (TDDFT). Our calculated spectra show a good agreement with the experimental absorption and fluorescence spectra, and thus, provide deep insights into the optical properties of the compounds. Furthermore, comparing the results obtained with four different hybrid functionals, we demonstrate that the accuracy of the functionals varies in the order B3LYP > PBE0 > M062X > M06HF. © 2015 Wiley Periodicals, Inc.     

The exact Fermi potential yielding the Hartree–Fock electron density from orbital‐free density functional theory

The exact expression for the Fermi potential yielding the Hartree–Fock electron density within an orbital‐free density functional formalism is derived. The Fermi potential, which is defined as that part of the potential that depends on the particles’ nature, is in this context given as the sum of the Pauli potential and the exchange potential. The exact exchange potential for an orbital‐free density functional formalism is shown to be the Slater potential.     

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.     

Fragmentation characteristics of hydroxycinnamic acids in ESI‐MSn by density functional theory

This work aims to analyze the electrospray ionization multistage mass spectrometry (ESI‐MSⁿ) fragmentation characteristics of hydroxycinnamic acids (HCAs) in negative ion mode. The geometric parameters, energies, natural bond orbitals and frontier orbitals of fragments were calculated by density functional theory (DFT) to investigate mass spectral fragmentation mechanisms. The results showed that proton transfer always occurred during fragmentation of HCAs; their quasi‐molecular ions ([M − H]⁻) existed in more than one form and were mainly with the lowest energy. The fragmentation characteristics included the followings: (1) according to the different substitution position of phenolic hydroxyl group, the ring contraction reaction by CO elimination from benzene was in an increasingly difficult order: m‐phenolic hydroxyl > p‐phenolic hydroxyl > o‐phenolic hydroxyl; and (2) ortho effect always occurred in o‐dihydroxycinnamic acids (o‐diHCAs), i.e. one phenolic hydroxyl group offered H⁺, which combined with the other one to lose H₂O. In addition, there was a nucleophilic reaction during ring contraction in diHCAs that oxygen atom attacked the carbon atom binding with the other phenolic hydroxyl to lose CO₂. The fragmentation characteristics and mechanism of HCAs could be used for analysis and identification of such compounds quickly and effectively, and as reference for structural analogues by ESI‐MS. Copyright © 2017 John Wiley & Sons, Ltd.     

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.