Perspective: Fifty years of density-functional theory in chemical physics

H4: A challenging system for natural orbital functional approximations

The correct description of nondynamic correlation by electronic structure methods not belonging to the multireference family is a challenging issue. The transition of D(2h) to D(4h) symmetry in H4 molecule is among the most simple archetypal examples to illustrate the consequences of missing nondynamic correlation effects. The resurgence of interest in density matrix functional methods has brought several new methods including the family of Piris Natural Orbital Functionals (PNOF). In this work, we compare PNOF5 and PNOF6, which include nondynamic electron correlation effects to some extent, with other standard ab initio methods in the H4 D(4h)/D(2h) potential energy surface (PES). Thus far, the wrongful behavior of single-reference methods at the D(2h)-D(4h) transition of H4 has been attributed to wrong account of nondynamic correlation effects, whereas in geminal-based approaches, it has been assigned to a wrong coupling of spins and the localized nature of the orbitals. We will show that actually interpair nondynamic correlation is the key to a cusp-free qualitatively correct description of H4 PES. By introducing interpair nondynamic correlation, PNOF6 is shown to avoid cusps and provide the correct smooth PES features at distances close to the equilibrium, total and local spin properties along with the correct electron delocalization, as reflected by natural orbitals and multicenter delocalization indices.

How Much Can Density Functional Approximations (DFA) Fail? The Extreme Case of the FeO4 Species

A thorough theoretical study of the relative energies of various molecular Fe·4O isomers with different oxidation states of both Fe and O atoms is presented, comparing simple Hartree-Fock through many Kohn-Sham approximations up to extended coupled cluster and DMRG multiconfiguration benchmark methods. The ground state of Fe·4O is a singlet, hexavalent iron(VI) complex (1)C2v-[Fe(VI)O2](2+)(O2)(2-), with isomers of oxidation states Fe(II), Fe(III), Fe(IV), Fe(V), and Fe(VIII) all lying slightly higher within the range of 1 eV. The disputed existence of oxidation state Fe(VIII) is discussed for isolated FeO4 molecules. Density functional theory (DFT) at various DF approximation (DFA) levels of local and gradient approaches, Hartree-Fock exchange and meta hybrids, range dependent, DFT-D and DFT+U models do not perform better for the relative stabilities of the geometric and electronic Fe·4O isomers than within 1-5 eV. The Fe·4O isomeric species are an excellent testing and validation ground for the development of density functional and wave function methods for strongly correlated multireference states, which do not seem to always follow chemical intuition.

Analytical nuclear gradients for the range-separated many-body dispersion model of noncovalent interactions

An accurate treatment of the long-range electron correlation energy, including van der Waals (vdW) or dispersion interactions, is essential for describing the structure, dynamics, and function of a wide variety of systems. Among the most accurate models for including dispersion into density functional theory (DFT) is the range-separated many-body dispersion (MBD) method [A. Ambrosetti et al., J. Chem. Phys., 2014, 140, 18A508], in which the correlation energy is modeled at short-range by a semi-local density functional and at long-range by a model system of coupled quantum harmonic oscillators. In this work, we develop analytical gradients of the MBD energy with respect to nuclear coordinates, including all implicit coordinate dependencies arising from the partitioning of the charge density into Hirshfeld effective volumes. To demonstrate the efficiency and accuracy of these MBD gradients for geometry optimizations of systems with intermolecular and intramolecular interactions, we optimized conformers of the benzene dimer and isolated small peptides with aromatic side-chains. We find excellent agreement with the wavefunction theory reference geometries of these systems (at a fraction of the computational cost) and find that MBD consistently outperforms the popular TS and D3(BJ) dispersion corrections. To demonstrate the performance of the MBD model on a larger system with supramolecular interactions, we optimized the C60@C60H28 buckyball catcher host-guest complex. In our analysis, we also find that neglecting the implicit nuclear coordinate dependence arising from the charge density partitioning, as has been done in prior numerical treatments, leads to an unacceptable error in the MBD forces, with relative errors of ∼20% (on average) that can extend well beyond 100%.

Nanoparticle shape, thermodynamics and kinetics

Nanoparticles can be beautiful, as in stained glass windows, or they can be ugly as in wear and corrosion debris from implants. We estimate that there will be about 70,000 papers in 2015 with nanoparticles as a keyword, but only one in thirteen uses the nanoparticle shape as an additional keyword and research focus, and only one in two hundred has thermodynamics. Methods for synthesizing nanoparticles have exploded over the last decade, but our understanding of how and why they take their forms has not progressed as fast. This topical review attempts to take a critical snapshot of the current understanding, focusing more on methods to predict than a purely synthetic or descriptive approach. We look at models and themes which are largely independent of the exact synthetic method whether it is deposition, gas-phase condensation, solution based or hydrothermal synthesis. Elements are old dating back to the beginning of the 20th century-some of the pioneering models developed then are still relevant today. Others are newer, a merging of older concepts such as kinetic-Wulff constructions with methods to understand minimum energy shapes for particles with twins. Overall we find that while there are still many unknowns, the broad framework of understanding and predicting the structure of nanoparticles via diverse Wulff constructions, either thermodynamic, local minima or kinetic has been exceedingly successful. However, the field is still developing and there remain many unknowns and new avenues for research, a few of these being suggested towards the end of the review.

Constrained subsystem density functional theory

Constrained Subsystem Density Fucntional Theory (CSDFT) allows to compute diabatic states for charge transfer reactions using the machinery of the constrained DFT method, and at the same time is able to embed such diabatic states in a molecular environment via a subsystem DFT scheme.

Structural and relative energy assessments of DFT functionals and the MP2 method to describe the gas phase methylation of nitronates: [R1R2CNO2]− + CH3I

The structures and energetics of CH₃I + [R¹R²CNO²]⁻ gas phase C- and O-methylation reactions were computed with several functionals and ECP/basis sets and compared to CCSD(T)/CBS.

Quantum and classical dynamics of reactive scattering of H2 from metal surfaces

State-of-the-art theoretical models allow nowadays an accurate description of H₂/metal surface systems and phenomena relative to heterogeneous catalysis. Here we review the most relevant ones investigated during the last 10 years.

Recent application of calculations of metal complexes based on density functional theory

Recent application of density functional theory (DFT) for metal complexes is reviewed to show the achievements of DFT and the challenges for it, as well as the methods for selecting proper functionals.

The role of exchange–correlation functional on the description of multiferroic properties using density functional theory: the ATiO3 (A = Mn, Fe, Ni) case study

In this study, ab initio density functional theory calculations were performed on ATiO₃ (A = Mn, Fe, Ni) materials for multiferroic applications.

New approaches for the calibration of exchange-energy densities in local hybrid functionals

A new derivation of semi-local calibration functions is proposed, to ameliorate the gauge problem of exchange-energy densities in local hybrid functionals.

Positive and negative linear compressibility and electronic properties of energetic and porous hybrid crystals with nitrate anions

Within the framework of DFT-D calculations, the compressibility anisotropy of UN and DATN energetic nitrates and the Ag(en)N hybrid crystal is established.

Electronically non-adiabatic influences in surface chemistry and dynamics

Electronically nonadiabatic interactions between molecules and metal surfaces are now well known. But evidence that such interactions influence reaction rates is still scarce. This paper reviews research related to this topic and proposes pathways forward.

Gentlest ascent dynamics for calculating first excited state and exploring energy landscape of Kohn-Sham density functionals

We develop the gentlest ascent dynamics for Kohn-Sham density functional theory to search for the index-1 saddle points on the energy landscape of the Kohn-Sham density functionals. These stationary solutions correspond to excited states in the ground state functionals. As shown by various examples, the first excited states of many chemical systems are given by these index-1 saddle points. Our novel approach provides an alternative, more robust way to obtain these excited states, compared with the widely used ΔSCF approach. The method can be easily generalized to target higher index saddle points. Our results also reveal the physical interest and relevance of studying the Kohn-Sham energy landscape.

Density Functional Theory Study of Atomic Layer Deposition of Zinc Oxide on Graphene

The dissociation of zinc ions (Zn(2+)) from vapor-phase zinc acetylacetonate, Zn(C5H7O2)2, or Zn(acac)2 and its adsorption onto graphene oxide via atomic layer deposition (ALD) were studied using a quantum mechanics approach. Density functional theory (DFT) was used to obtain an approximate solution to the Schrödinger equation. The graphene oxide cluster model was used to represent the surface of the graphene film after pre-oxidation. In this study, the geometries of reactants, transition states, and products were optimized using the B3LYB/6-31G** level of theory or higher. Furthermore, the relative energies of the various intermediates and products in the gas-phase radical mechanism were calculated at the B3LYP/6-311++G** and MP2/6-311 + G(2df,2p) levels of theory. Additionally, a molecular orbital (MO) analysis was performed for the products of the decomposition of the Zn(acac)2 complex to investigate the dissociation of Zn(2+) and the subsequent adsorption of H atoms on the C5H7O2 cluster to form acetylacetonate enol. The reaction energies were calculated, and the reaction mechanism was accordingly proposed. A simulation of infrared (IR) properties was performed using the same approach to support the proposed mechanism via a complete explanation of bond forming and breaking during each reaction step.

Size-dependent properties of transition metal clusters: from molecules to crystals and surfaces--computational studies with the program ParaGauss

In the so-called scalable regime the size-dependent behavior of the physical and chemical properties of transition metal clusters is described by scaling relationships. For most quantities this scalable regime is reached for cluster sizes between a few tens and a few hundreds of atoms, hence for systems for which an accurate treatment by density functional theory is still feasible. Thus, by invoking scaling relations one is able to obtain properties of very large nanoparticles and even the bulk limit from the results of a series of smaller cluster models. In this invited review we illustrate this strategy by exploiting results from computational studies that mostly were carried out with the density functional theory software ParaGauss. We address mainly the size-dependent behavior of the properties of transition metal clusters. To this end, we first present benchmark studies probing various approximations that are used in such density functional calculations. Subsequently we show how physical insight may be gained by exploring less understood types of systems. These applications range from bare clusters to nanoislands and nanoalloys to adsorption complexes.

Water on BN doped benzene: a hard test for exchange-correlation functionals and the impact of exact exchange on weak binding

Density functional theory (DFT) studies of weakly interacting complexes have recently focused on the importance of van der Waals dispersion forces, whereas the role of exchange has received far less attention. Here, by exploiting the subtle binding between water and a boron and nitrogen doped benzene derivative (1,2-azaborine) we show how exact exchange can alter the binding conformation within a complex. Benchmark values have been calculated for three orientations of the water monomer on 1,2-azaborine from explicitly correlated quantum chemical methods, and we have also used diffusion quantum Monte Carlo. For a host of popular DFT exchange-correlation functionals we show that the lack of exact exchange leads to the wrong lowest energy orientation of water on 1,2-azaborine. As such, we suggest that a high proportion of exact exchange and the associated improvement in the electronic structure could be needed for the accurate prediction of physisorption sites on doped surfaces and in complex organic molecules. Meanwhile to predict correct absolute interaction energies an accurate description of exchange needs to be augmented by dispersion inclusive functionals, and certain non-local van der Waals functionals (optB88- and optB86b-vdW) perform very well for absolute interaction energies. Through a comparison with water on benzene and borazine (B3N3H6) we show that these results could have implications for the interaction of water with doped graphene surfaces, and suggest a possible way of tuning the interaction energy.

How close are the Slater and Becke-Roussel potentials in solids?

The Becke–Roussel (BR) potential [Phys. Rev. A 1989, 39, 3761] was proposed as an approximation to the Slater potential, which is the Coulomb potential generated by the exact exchange hole. In the present work, a detailed comparison between the Slater and BR potentials in solids is presented. It is shown that the two potentials usually lead to very similar results for the electronic structure; however, in a few cases, e.g., Si, Ge, or strongly correlated systems like NiO, the fundamental band gap or magnetic properties can differ markedly. Such differences should not be neglected when the computationally expensive Slater potential is replaced by the cheap semilocal BR potential in approximations to the exact-exchange Kohn–Sham potential, such as the one proposed by Becke and Johnson [J. Chem. Phys. 2006, 124, 221101].

Application of Koopmans' theorem for density functional theory to full valence-band photoemission spectroscopy modeling

In this work, Koopmans' theorem for Kohn-Sham density functional theory (KS-DFT) is applied to the photoemission spectra (PES) modeling over the entire valence-band. To examine the validity of this application, a PES modeling scheme is developed to facilitate a full valence-band comparison of theoretical PES spectra with experiments. The PES model incorporates the variations of electron ionization cross-sections over atomic orbitals and a linear dispersion of spectral broadening widths. KS-DFT simulations of pristine rubrene (5,6,11,12-tetraphenyltetracene) and potassium-rubrene complex are performed, and the simulation results are used as the input to the PES models. Two conclusions are reached. First, decompositions of the theoretical total spectra show that the dissociated electron of the potassium mainly remains on the backbone and has little effect on the electronic structures of phenyl side groups. This and other electronic-structure results deduced from the spectral decompositions have been qualitatively obtained with the anionic approximation to potassium-rubrene complexes. The qualitative validity of the anionic approximation is thus verified. Second, comparison of the theoretical PES with the experiments shows that the full-scale simulations combined with the PES modeling methods greatly enhance the agreement on spectral shapes over the anionic approximation. This agreement of the theoretical PES spectra with the experiments over the full valence-band can be regarded, to some extent, as a collective validation of the application of Koopmans' theorem for KS-DFT to valence-band PES, at least, for this hydrocarbon and its alkali-adsorbed complex.

Construction of exchange-correlation functionals through interpolation between the non-interacting and the strong-correlation limit

Drawing on the adiabatic connection of density functional theory, exchange-correlation functionals of Kohn-Sham density functional theory are constructed which interpolate between the extreme limits of the electron-electron interaction strength. The first limit is the non-interacting one, where there is only exchange. The second limit is the strong correlated one, characterized as the minimum of the electron-electron repulsion energy. The exchange-correlation energy in the strong-correlation limit is approximated through a model for the exchange-correlation hole that is referred to as nonlocal-radius model [L. O. Wagner and P. Gori-Giorgi, Phys. Rev. A 90, 052512 (2014)]. Using the non-interacting and strong-correlated extremes, various interpolation schemes are presented that yield new approximations to the adiabatic connection and thus to the exchange-correlation energy. Some of them rely on empiricism while others do not. Several of the proposed approximations yield the exact exchange-correlation energy for one-electron systems where local and semi-local approximations often fail badly. Other proposed approximations generalize existing global hybrids by using a fraction of the exchange-correlation energy in the strong-correlation limit to replace an equal fraction of the semi-local approximation to the exchange-correlation energy in the strong-correlation limit. The performance of the proposed approximations is evaluated for molecular atomization energies, total atomic energies, and ionization potentials.

DFTB3 Parametrization for Copper: The Importance of Orbital Angular Momentum Dependence of Hubbard Parameters

We report the parametrization of a density functional tight binding method (DFTB3) for copper in a spin-polarized formulation. The parametrization is consistent with the framework of 3OB for main group elements (ONCHPS) and can be readily used for biological applications that involve copper proteins/peptides. The key to our parametrization is to introduce orbital angular momentum dependence of the Hubbard parameter and its charge derivative, thus allowing the 3d and 4s orbitals to adopt different sizes and responses to the change of charge state. The parametrization has been tested by applying to a fairly broad set of molecules of biological relevance, and the properties of interest include optimized geometries, ligand binding energies, and ligand proton affinities. Compared to the reference QM level (B3LYP/aug-cc-pVTZ, which is shown here to be similar to the B97-1 and CCSD(T) results, in terms of many properties of interest for a set of small copper containing molecules), our parametrization generally gives reliable structural properties for both Cu(I) and Cu(II) compounds, although several exceptions are also noted. For energetics, the results are more accurate for neutral ligands than for charged ligands, likely reflecting the minimal basis limitation of DFTB3; the results generally outperform NDDO based methods such as PM6 and even PBE with the 6-31+G(d,p) basis. For all ligand types, single-point B3LYP calculations at DFTB3 geometries give results very close (∼1-2 kcal/mol) to the reference B3LYP values, highlighting the consistency between DFTB3 and B3LYP structures. Possible further developments of the DFTB3 model for a better treatment of transition-metal ions are also discussed. In the current form, our first generation of DFTB3 copper model is expected to be particularly valuable as a method that drives sampling in systems that feature a dynamical copper binding site.

London dispersion in molecular chemistry--reconsidering steric effects

London dispersion, which constitutes the attractive part of the famous van der Waals potential, has long been underappreciated in molecular chemistry as an important element of structural stability, and thus affects chemical reactivity and catalysis. This negligence is due to the common notion that dispersion is weak, which is only true for one pair of interacting atoms. For increasingly larger structures, the overall dispersion contribution grows rapidly and can amount to tens of kcal mol(-1) . This Review collects and emphasizes the importance of inter- and intramolecular dispersion for molecules consisting mostly of first row atoms. The synergy of experiment and theory has now reached a stage where dispersion effects can be examined in fine detail. This forces us to reconsider our perception of steric hindrance and stereoelectronic effects. The quantitation of dispersion energy donors will improve our ability to design sophisticated molecular structures and much better catalysts.

Density and Compressibility of Liquid Water and Ice from First-Principles Simulations with Hybrid Functionals

We determined the equilibrium density and compressibility of water and ice from first-principles molecular dynamics simulations using gradient-corrected (PBE) and hybrid (PBE0) functionals. Both functionals predicted the density of ice to be larger than that of water, by 15 (PBE) and 35% (PBE0). The PBE0 functional yielded a lower density of both ice and water with respect to PBE, leading to better agreement with experiment for ice but not for liquid water. Approximate inclusion of dispersion interactions on computed molecular-dynamics trajectories led to a substantial improvement of the PBE0 results for the density of liquid water, which, however, resulted to be slightly lower than that of ice.

The Representation and Parametrization of Orthogonal Matrices

Four representations and parametrizations of orthogonal matrices Q ∈ R(m×n) in terms of the minimal number of essential parameters {φ} are discussed: the exponential representation, the Householder reflector representation, the Givens rotation representation, and the rational Cayley transform representation. Both square n = m and rectangular n < m situations are considered. Two separate kinds of parametrizations are considered: one in which the individual columns of Q are distinct, the Stiefel manifold, and the other in which only span(Q) is significant, the Grassmann manifold. The practical issues of numerical stability, continuity, and uniqueness are discussed. The computation of Q in terms of the essential parameters {φ}, and also the extraction of {φ} for a given Q are considered for all of the parametrizations. The transformation of gradient arrays between the Q and {φ} variables is discussed for all representations. It is our hope that developers of new methods will benefit from this comparative presentation of an important but rarely analyzed subject.