Book•
Density-functional theory of atoms and molecules
01 Jan 1989-
TL;DR: In this paper, a review of current studies in density functional theory and density matrix functional theory is presented, with special attention to the possible applications within chemistry, including the concept of an atom in a molecule, calculation of electronegativities from the Xα method, pressure, Gibbs-Duhem equation, Maxwell relations and stability conditions.
Abstract: Current studies in density functional theory and density matrix functional theory are reviewed, with special attention to the possible applications within chemistry. Topics discussed include the concept of electronegativity, the concept of an atom in a molecule, calculation of electronegativities from the Xα method, the concept of pressure, Gibbs-Duhem equation, Maxwell relations, stability conditions, and local density functional theory.
Citations
More filters
TL;DR: In this article, a semi-empirical exchange correlation functional with local spin density, gradient, and exact exchange terms was proposed. But this functional performed significantly better than previous functionals with gradient corrections only, and fits experimental atomization energies with an impressively small average absolute deviation of 2.4 kcal/mol.
Abstract: Despite the remarkable thermochemical accuracy of Kohn–Sham density‐functional theories with gradient corrections for exchange‐correlation [see, for example, A. D. Becke, J. Chem. Phys. 96, 2155 (1992)], we believe that further improvements are unlikely unless exact‐exchange information is considered. Arguments to support this view are presented, and a semiempirical exchange‐correlation functional containing local‐spin‐density, gradient, and exact‐exchange terms is tested on 56 atomization energies, 42 ionization potentials, 8 proton affinities, and 10 total atomic energies of first‐ and second‐row systems. This functional performs significantly better than previous functionals with gradient corrections only, and fits experimental atomization energies with an impressively small average absolute deviation of 2.4 kcal/mol.
87,732 citations
TL;DR: A new density functional of the generalized gradient approximation (GGA) type for general chemistry applications termed B97‐D is proposed, based on Becke's power‐series ansatz from 1997, and is explicitly parameterized by including damped atom‐pairwise dispersion corrections of the form C6 · R−6.
Abstract: A new density functional (DF) of the generalized gradient approximation (GGA) type for general chemistry applications termed B97-D is proposed. It is based on Becke's power-series ansatz from 1997 and is explicitly parameterized by including damped atom-pairwise dispersion corrections of the form C(6) x R(-6). A general computational scheme for the parameters used in this correction has been established and parameters for elements up to xenon and a scaling factor for the dispersion part for several common density functionals (BLYP, PBE, TPSS, B3LYP) are reported. The new functional is tested in comparison with other GGAs and the B3LYP hybrid functional on standard thermochemical benchmark sets, for 40 noncovalently bound complexes, including large stacked aromatic molecules and group II element clusters, and for the computation of molecular geometries. Further cross-validation tests were performed for organometallic reactions and other difficult problems for standard functionals. In summary, it is found that B97-D belongs to one of the most accurate general purpose GGAs, reaching, for example for the G97/2 set of heat of formations, a mean absolute deviation of only 3.8 kcal mol(-1). The performance for noncovalently bound systems including many pure van der Waals complexes is exceptionally good, reaching on the average CCSD(T) accuracy. The basic strategy in the development to restrict the density functional description to shorter electron correlation lengths scales and to describe situations with medium to large interatomic distances by damped C(6) x R(-6) terms seems to be very successful, as demonstrated for some notoriously difficult reactions. As an example, for the isomerization of larger branched to linear alkanes, B97-D is the only DF available that yields the right sign for the energy difference. From a practical point of view, the new functional seems to be quite robust and it is thus suggested as an efficient and accurate quantum chemical method for large systems where dispersion forces are of general importance.
23,058 citations
University of Udine1, International School for Advanced Studies2, National Research Council3, Massachusetts Institute of Technology4, University of Paris5, Princeton University6, University of Minnesota7, ParisTech8, University of Milan9, International Centre for Theoretical Physics10, University of Paderborn11, ETH Zurich12, École Polytechnique Fédérale de Lausanne13
TL;DR: QUANTUM ESPRESSO as discussed by the authors is an integrated suite of computer codes for electronic-structure calculations and materials modeling, based on density functional theory, plane waves, and pseudopotentials (norm-conserving, ultrasoft, and projector-augmented wave).
Abstract: QUANTUM ESPRESSO is an integrated suite of computer codes for electronic-structure calculations and materials modeling, based on density-functional theory, plane waves, and pseudopotentials (norm-conserving, ultrasoft, and projector-augmented wave). The acronym ESPRESSO stands for opEn Source Package for Research in Electronic Structure, Simulation, and Optimization. It is freely available to researchers around the world under the terms of the GNU General Public License. QUANTUM ESPRESSO builds upon newly-restructured electronic-structure codes that have been developed and tested by some of the original authors of novel electronic-structure algorithms and applied in the last twenty years by some of the leading materials modeling groups worldwide. Innovation and efficiency are still its main focus, with special attention paid to massively parallel architectures, and a great effort being devoted to user friendliness. QUANTUM ESPRESSO is evolving towards a distribution of independent and interoperable codes in the spirit of an open-source project, where researchers active in the field of electronic-structure calculations are encouraged to participate in the project by contributing their own codes or by implementing their own ideas into existing codes.
19,985 citations
TL;DR: Five practical examples involving a wide variety of systems and analysis methods are given to illustrate the usefulness of Multiwfn, a multifunctional program for wavefunction analysis.
Abstract: Multiwfn is a multifunctional program for wavefunction analysis. Its main functions are: (1) Calculating and visualizing real space function, such as electrostatic potential and electron localization function at point, in a line, in a plane or in a spatial scope. (2) Population analysis. (3) Bond order analysis. (4) Orbital composition analysis. (5) Plot density-of-states and spectrum. (6) Topology analysis for electron density. Some other useful utilities involved in quantum chemistry studies are also provided. The built-in graph module enables the results of wavefunction analysis to be plotted directly or exported to high-quality graphic file. The program interface is very user-friendly and suitable for both research and teaching purpose. The code of Multiwfn is substantially optimized and parallelized. Its efficiency is demonstrated to be significantly higher than related programs with the same functions. Five practical examples involving a wide variety of systems and analysis methods are given to illustrate the usefulness of Multiwfn. The program is free of charge and open-source. Its precompiled file and source codes are available from http://multiwfn.codeplex.com.
17,273 citations
Cites background from "Density-functional theory of atoms ..."
...have reinterpreted ELF in the view of kinetic energy,([25,26]) which made ELF also meaningful for Kohn–Sham DFT wavefunction.([13]) They indicated that D(r) reveals the excess kinetic energy density caused by Pauli repulsion, whereas D0(r) can be considered as Thomas–Fermi kinetic energy density....
[...]
...A similar function named spin polarization parameter([13]) is expressed as follows:...
[...]
TL;DR: In this article, a new coupling of Hartree-Fock theory with local density functional theory was proposed to improve the predictive power of the Hartree−Fock model for molecular bonding, and the results of tests on atomization energies, ionization potentials, and proton affinities were reported.
Abstract: Previous attempts to combine Hartree–Fock theory with local density‐functional theory have been unsuccessful in applications to molecular bonding. We derive a new coupling of these two theories that maintains their simplicity and computational efficiency, and yet greatly improves their predictive power. Very encouraging results of tests on atomization energies, ionization potentials, and proton affinities are reported, and the potential for future development is discussed.
13,853 citations
References
More filters
TL;DR: In this article, the ground state of an interacting electron gas in an external potential was investigated and it was proved that there exists a universal functional of the density, called F[n(mathrm{r})], independent of the potential of the electron gas.
Abstract: This paper deals with the ground state of an interacting electron gas in an external potential $v(\mathrm{r})$. It is proved that there exists a universal functional of the density, $F[n(\mathrm{r})]$, independent of $v(\mathrm{r})$, such that the expression $E\ensuremath{\equiv}\ensuremath{\int}v(\mathrm{r})n(\mathrm{r})d\mathrm{r}+F[n(\mathrm{r})]$ has as its minimum value the correct ground-state energy associated with $v(\mathrm{r})$. The functional $F[n(\mathrm{r})]$ is then discussed for two situations: (1) $n(\mathrm{r})={n}_{0}+\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{n}(\mathrm{r})$, $\frac{\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{n}}{{n}_{0}}\ensuremath{\ll}1$, and (2) $n(\mathrm{r})=\ensuremath{\phi}(\frac{\mathrm{r}}{{r}_{0}})$ with $\ensuremath{\phi}$ arbitrary and ${r}_{0}\ensuremath{\rightarrow}\ensuremath{\infty}$. In both cases $F$ can be expressed entirely in terms of the correlation energy and linear and higher order electronic polarizabilities of a uniform electron gas. This approach also sheds some light on generalized Thomas-Fermi methods and their limitations. Some new extensions of these methods are presented.
38,160 citations
TL;DR: The concept of electronegativity is defined in this paper as the negative of the chemical potential (the Lagrange multiplier for the normalization constraint) in the Hohenberg-Kohn density functional theory of the ground state: χ =−μ=−(∂E/∂N)v.
Abstract: Precision is given to the concept of electronegativity. It is the negative of the chemical potential (the Lagrange multiplier for the normalization constraint) in the Hohenberg–Kohn density functional theory of the ground state: χ=−μ=−(∂E/∂N)v. Electronegativity is constant throughout an atom or molecule, and constant from orbital to orbital within an atom or molecule. Definitions are given of the concepts of an atom in a molecule and of a valence state of an atom in a molecule, and it is shown how valence‐state electronegativity differences drive charge transfers on molecule formation. An equation of Gibbs–Duhem type is given for the change of electronegativity from one situation to another, and some discussion is given of certain relations among energy components discovered by Fraga.
2,632 citations
642 citations
TL;DR: In this article, the authors define the assumed Hermitian kernel FM[γ;x′,x]≡δEvM/δγ (x,x′) and define a set of coupled Euler equations for the representation coefficients and spin orbitals of a rank-M approximation to the exact ground-state density matrix.
Abstract: The Hohenberg–Kohn theorem implies the existence of an energy functional based solely on the first‐order reduced density matrix of the ground state of an atomic or molecular system. Application of the variational principle for the functional EvM[γ] generates a set of coupled Euler equations for the representation coefficients and spin orbitals of a rank‐M approximation to the exact ground‐state density matrix. Defining the (assumed Hermitian) kernel FM[γ;x′,x]≡δEvM/δγ (x,x′), the equations in an arbitrary representation for the approximate density matrix, γ (x,x′) =ΣijMψ1(x) γijψj* (x′) are the following: FdξFM[γ;x′,ξ]ψi(ξ) =μψi(x′); i=1, 2,...,M; FdξFM[γ;x′,ξ] ΣMiψi(ξ) γ ij=ΣMiψi(x′) λij; j=1, 2, ...,M. The quantity μ is the chemical potential of the system of interest and the λij are a set of M2 Lagrange multipliers constraining the orthonormality of the spin orbital basis {ψ}. The coefficients γij must be chosen such that FM has the degenerate eigenvalue spectrum, FiiM=μ, i=1, 2,...,M, for all partiall...
245 citations