Constraints for hierarchies of many electron distribution functions

Some inequalities that constrain the reconstruction of k-electron distribution functions from lower-order distribution functions are presented. These inequalities are related to the N-representability conditions on electron distributions functions and they have two basic types: (1) general N-representability inequalities, which are very powerful but difficult to apply and (2) generalized “Davidson” inequalities, which are less powerful but which may be more facile in computational implementations. A constraint on the exchange-correlation hole is also presented.     

Universal mathematical identities in density functional theory: results from three different spin-resolved representations

This paper supersedes previous theoretical approaches to conceptual DFT because it provides a unified and systematic approach to all of the commonly considered formulations of conceptual DFT, and even provides the essential mathematical framework for new formulations. Global, local, and nonlocal chemical reactivity indicators associated with the "closed-system representation" ([N(alpha),N(beta),nu(alpha)(r),nu(beta)(r)]) of spin-polarized density functional theory (SP-DFT) are derived. The links between these indicators and the ones associated with the "open-system representation" ([mu(alpha),mu(beta),nu(alpha)(r),nu(beta)(r)]) are derived, including the spin-resolved Berkowitz-Parr identity. The Legendre transform to the "density representation" ([rho(alpha)(r),rho(beta)(r)]) is performed, and the spin-resolved Harbola-Chattaraj-Cedillo-Parr identities linking the density representation to the closed-system and open-system representations are derived. Taken together, these results provide the framework for understanding chemical reactions from both the electron-following perspective (using either the closed-system or the open-system representation) and electron-preceding perspective (density representation). A powerful matrix-vector notation is developed; with this notation, identities in conceptual DFT become universal. Specifically, this notation allows the fundamental identities in conventional (spin-free) conceptual DFT, the [N(alpha),N(beta)] representation, and the [N=N(alpha)+N(beta),N(S)=N(alpha)-N(beta)] representation to be written in exactly the same forms. In cases where spin transfer and electron transfer are coupled (e.g., radical+molecule reactions), we believe that the [N(alpha),N(beta)] representation may be more useful than the more common [N,N(S)] representation.     

Removing electrons can increase the electron density: a computational study of negative Fukui functions

Ab initio and density-functional theory calculations for a family of substituted acetylenes show that removing electrons from these molecules causes the electron density along the C-C bond to increase. This result contradicts the predictions of simple frontier molecular orbital theory, but it is easily explained using the nucleophilic Fukui function-provided that one is willing to allow for the Fukui function to be negative. Negative Fukui functions emerge as key indicators of redox-induced electron rearrangements, where oxidation of an entire molecule (acetylene) leads to reduction of a specific region of the molecule (along the bond axis, between the carbon atoms). Remarkably, further oxidization of these substituted acetylenes (one can remove as many as four electrons!) causes the electron density along the C-C bond to increase even more. This work provides substantial evidence that the molecular Fukui function is sometimes negative and reveals that this is due to orbital relaxation.     

Elucidating the hard/soft acid/base principle: a perspective based on half-reactions

A comprehensive analysis is presented for the acid-base double-exchange reaction as well as the associated acid-displacement and base-displacement "half-reactions" with the goal of elucidating the meaning of the hard/soft acid/base (HSAB) principle and the conditions for its validity. When electron-transfer effects are important and other effects are negligible, the HSAB principle is driven by the surpassing stability of the soft acid/soft base product. When electrostatic effects dominate the reactivity, the HSAB principle is driven by the surpassing stability of the hard acid/hard base product. Because electron-transfer effects favor soft/soft interactions, while electrostatic effects favor hard/hard interactions, acid-base exchange reactions may be used to determine whether a reagent's reactivity is dominated by electron-transfer or by electrostatic effects. Because electron-transfer and electrostatic considerations separately favor the HSAB principle whenever the electronic chemical potentials of the acids and bases involved in the reaction are similar, our analysis provides strong support for the HSAB principle. The electronic chemical potential measures the intrinsic strength of acids and bases.     

Woodward-Hoffmann rules in density functional theory: initial hardness response

The Woodward-Hoffmann rules for pericyclic reactions, a fundamental set of reactivity rules in organic chemistry, are formulated in the language of conceptual density functional theory (DFT). DFT provides an elegant framework to introduce chemical concepts and principles in a quantitative manner, partly because it is formulated without explicit reference to a wave function, on whose symmetry properties the Woodward-Hoffmann [J. Am. Chem. Soc. 87, 395 (1965)] rules are based. We have studied the initial chemical hardness response using a model reaction profile for two prototypical pericyclic reactions, the Diels-Alder cycloaddition of 1,3-butadiene to ethylene and the addition of ethylene to ethylene, both in the singlet ground state and in the first triplet excited state. For the reaction that is thermally allowed but photochemically forbidden, the initial hardness response is positive along the singlet reaction profile. (By contrast, for the triplet reaction profile, a negative hardness response is observed.) For the photochemically allowed, thermally forbidden reaction, the behavior of the chemical hardness along the initial stages of the singlet and triplet reaction profiles is reversed. This constitutes a first step in showing that chemical concepts from DFT can be invoked to explain results that would otherwise require invoking the phase of the wave function.     

Generalizations of the Hohenberg-Kohn theorem: I. Legendre transform constructions of variational principles for density matrices and electron distribution functions

Given a general, N-particle Hamiltonian operator, analogs of the Hohenberg-Kohn theorem are derived for functions that are more general than the particle density, including density matrices and the diagonal elements thereof. The generalization of Lieb's Legendre transform ansatz to the generalized Hohenberg-Kohn functional not only solves the upsilon-representability problem for these entities, but, more importantly, also solves the N-representability problem. Restricting the range of operators explored by the Legendre transform leads to a lower bound on the true functional. If all the operators of interest are incorporated in the restricted maximization, however, the variational principle dictates that exact results are obtained for the systems of interest. This might have important implications for practical work not only for density matrices but also for density functionals. A follow-up paper will present a useful alternative approach to the upsilon- and N-representability problems based on the constrained search formalism.