Reactivity indices and fluctuation formulas in density functional theory: Isomorphic ensembles and a new measure of local hardness

Expanding horizons in conceptual density functional theory: Novel ensembles and descriptors to decipher reactivity patterns

Conceptual density functional theory (CDFT) and the quantum reactivity descriptors stemming from it have proven to be valuable tools for understanding the chemical behavior of molecules. This article is presented as being intrinsically of dual character. In a first part, it briefly reviews, in a deliberately didactical way, the main ensembles in CDFT, while the second half presents two additional ensembles, where the chemical hardness acts as a natural variable, and their respective reactivity descriptors. The evaluation of these reactivity descriptors on common organic chemical reagents are presented and discussed.

Links among the Fukui potential, the alchemical hardness and the local hardness of an atom in a molecule

This paper presents a brief summary of the difficulty that resides in the definition of the elusive concept of local chemical hardness. We argue that a definition of local hardness should be useful to a reactivity principle and not just as a mere definition. We then continue with a formal discussion about the benefits and difficulties of using the Fukui potential, which is interpreted as an alchemical derivative (alchemical hardness), as descriptor of local hardness of molecules. Computational evidence shows that the alchemical hardness is at least as good a descriptor as the combination of other two well-stabilized descriptors of local hardness, such as the Fukui function and grand canonical local hardness. Although our results are auspicious for the alchemical hardness as descriptor of local hardness, we finish by calling the attention of the community on the importance of discussing the raison d'être of a local hardness function and its main characteristics. We suggest that an axiomatic construction of local hardness could be they way of constructing a local hardness which is both useful and free of arbitrariness.

N-Derivatives of Shannon entropy density as response functions

The exact first and second order partial derivatives of Shannon entropy density with respect to the number of electrons at constant external potential are introduced as new descriptors for prediction of the active sites of a molecule. The derivatives, which are a measure of the inhomogeneity of electron density, are calculated both exactly (from analytical forms) and approximately (using the finite difference method) for some molecular systems. According to the maximum entropy principle, the extreme value of the first order derivative on the surface of a given molecule should determine the active sites of the molecule in electrophilic and nucleophilic attack. The second order derivative indicates where the Shannon entropy is more concentrated or depleted during the electron exchange. Although these derivatives on the surfaces of helium and neon atoms are uniform, the corresponding values for argon, krypton and xenon atoms are not. This could explain the greater tendency of heavy noble gas atoms to form stable compounds. A dual descriptor is also defined as the difference between the left and right first order derivatives of Shannon entropy density, which allows one to simultaneously predict the preferable sites for electrophilic and nucleophilic attack over the system at point r. Therefore, the reactivity of an atom in a molecule requires the non-uniformity of the first and second order derivatives of Shannon entropy density on the surface of that atom.

Comment on “Revisiting the definition of local hardness and hardness kernel” by C. A. Polanco-Ramirez, M. Franco-Pérez, J. Carmona-Espíndola, J. L. Gázquez and P. W. Ayers, Phys. Chem. Chem. Phys., 2017, 19, 12355

In this comment we show that the derivation proposed Polanco-Ramirez et al. appears naturally in the Taylor expansion of the energy, showing that their whole construction is not artificially built.

Conceptual DFT analysis of the fragility spectra of atoms along the minimum energy reaction coordinate

Theoretical justification has been provided to the method for monitoring the sequence of chemical bonds' rearrangement along a reaction path, by tracing the evolution of the diagonal elements of the Hessian matrix. Relations between the divergences of Hellman-Feynman forces and the energy and electron density derivatives have been demonstrated. By the proof presented on the grounds of the conceptual density functional theory formalism, the spectral amplitude observed on the atomic fragility spectra [L. Komorowski et al., Phys. Chem. Chem. Phys. 18, 32658 (2016)] reflects selectively the electron density modifications in bonds of an atom. In fact the spectral peaks for an atom reveal changes of the electron density occurring with bonds creation, breaking, or varying with the reaction progress.

New Fukui, dual and hyper-dual kernels as bond reactivity descriptors

We define three new linear response indices with promising applications for bond reactivity using the mathematical framework of τ-CRT (finite temperature chemical reactivity theory). The τ-Fukui kernel is defined as the ratio between the fluctuations of the average electron density at two different points in the space and the fluctuations in the average electron number and is designed to integrate to the finite-temperature definition of the electronic Fukui function. When this kernel is condensed, it can be interpreted as a site-reactivity descriptor of the boundary region between two atoms. The τ-dual kernel corresponds to the first order response of the Fukui kernel and is designed to integrate to the finite temperature definition of the dual descriptor; it indicates the ambiphilic reactivity of a specific bond and enriches the traditional dual descriptor by allowing one to distinguish between the electron-accepting and electron-donating processes. Finally, the τ-hyper dual kernel is defined as the second-order derivative of the Fukui kernel and is proposed as a measure of the strength of ambiphilic bonding interactions. Although these quantities have never been proposed, our results for the τ-Fukui kernel and for τ-dual kernel can be derived in zero-temperature formulation of the chemical reactivity theory with, among other things, the widely-used parabolic interpolation model.

Local chemical potential, local hardness, and dual descriptors in temperature dependent chemical reactivity theory

From the definition of a local chemical potential, well-behaved expressions for the local hardness and the dual descriptors are derived.

Revisiting the definition of local hardness and hardness kernel

Local hardness is redefined following similar rules to those of local softness. The new concept describes chemical trends correctly.

Fractional electron number, temperature, and perturbations in chemical reactions

The mathematical framework of conceptual density functional theory is extended to use the eigenstates and eigenvalues of perturbed subsystems. This unites, justifies, and extends, several previously proposed models.

New link between conceptual density functional theory and electron delocalization

In this paper we give a new definition of the softness kernel based on the exchange-correlation density. This new kernel is shown to correspond to the change of electron fluctuation upon external perturbation, thus helping to bridge the gap between conceptual density functional theory and some tools describing electron localization in molecules. With the aid of a few computational calculations on diatomics we illustrate the performance of this new computational tool.

Application of the electron density force to chemical reactivity

In this paper the concept of force experienced by the electron density is applied to chemical reactivity. The force is based upon the gradient of a local chemical potential. Closely related concepts such as force field lines and local electron flux are defined to provide insight in chemical reactivity. The time evolution of a molecular site density is also proposed. From the divergence of the force, the nucleophilic and electrophilic behaviour of atomic sites are characterized. Finally, the relations between the force and local conceptual DFT descriptors are also given.

Information carriers and (reading them through) information theory in quantum chemistry

This Perspective discusses the reduction of the electronic wave function via the second-order reduced density matrix to the electron density ρ(r), which is the key ingredient in density functional theory (DFT) as a basic carrier of information. Simplifying further, the 1-normalized density function turns out to contain essentially the same information as ρ(r) and is even of preferred use as an information carrier when discussing the periodic properties along Mendeleev's table where essentially the valence electrons are at stake. The Kullback-Leibler information deficiency turns out to be the most interesting choice to obtain information on the differences in ρ(r) or σ(r) between two systems. To put it otherwise: when looking for the construction of a functional F(AB) = F[ζ(A)(r),ζ(B)(r)] for extracting differences in information from an information carrier ζ(r) (i.e. ρ(r), σ(r)) for two systems A and B the Kullback-Leibler information measure ΔS is a particularly adequate choice. Examples are given, varying from atoms, to molecules and molecular interactions. Quantum similarity of atoms indicates that the shape function based KL information deficiency is the most appropriate tool to retrieve periodicity in the Periodic Table. The dissimilarity of enantiomers for which different information measures are presented at global and local (i.e. molecular and atomic) level leads to an extension of Mezey's holographic density theorem and shows numerical evidence that in a chiral molecule the whole molecule is pervaded by chirality. Finally Kullback-Leibler information profiles are discussed for intra- and intermolecular proton transfer reactions and a simple S(N)2 reaction indicating that the theoretical information profile can be used as a companion to the energy based Hammond postulate to discuss the early or late transition state character of a reaction. All in all this Perspective's answer is positive to the question of whether an even simpler carrier of information than the electron density function ρ(r) can be envisaged: the shape function, integrating to 1 by construction fulfils this role. On the other hand obtaining the information (or information difference) contained in one (or two) systems from ρ(r) or σ(r) can be most efficiently done by using information theory, the Kulback-Leibler information deficiency being at the moment (one of) the most advisable functionals.

Particle number and probability density functional theory and A-representability

In Hohenberg-Kohn density functional theory, the energy E is expressed as a unique functional of the ground state density rho(r): E = E[rho] with the internal energy component F(HK)[rho] being universal. Knowledge of the functional F(HK)[rho] by itself, however, is insufficient to obtain the energy: the particle number N is primary. By emphasizing this primacy, the energy E is written as a nonuniversal functional of N and probability density p(r): E = E[N,p]. The set of functions p(r) satisfies the constraints of normalization to unity and non-negativity, exists for each N; N = 1, ..., infinity, and defines the probability density or p-space. A particle number N and probability density p(r) functional theory is constructed. Two examples for which the exact energy functionals E[N,p] are known are provided. The concept of A-representability is introduced, by which it is meant the set of functions Psi(p) that leads to probability densities p(r) obtained as the quantum-mechanical expectation of the probability density operator, and which satisfies the above constraints. We show that the set of functions p(r) of p-space is equivalent to the A-representable probability density set. We also show via the Harriman and Gilbert constructions that the A-representable and N-representable probability density p(r) sets are equivalent.

Conceptual DFT: the chemical relevance of higher response functions

In recent years conceptual density functional theory offered a perspective for the interpretation/prediction of experimental/theoretical reactivity data on the basis of a series of response functions to perturbations in the number of electrons and/or external potential. This approach has enabled the sharp definition and computation, from first principles, of a series of well-known but sometimes vaguely defined chemical concepts such as electronegativity and hardness. In this contribution, a short overview of the shortcomings of the simplest, first order response functions is illustrated leading to a description of chemical bonding in a covalent interaction in terms of interacting atoms or groups, governed by electrostatics with the tendency to polarize bonds on the basis of electronegativity differences. The second order approach, well known until now, introduces the hardness/softness and Fukui function concepts related to polarizability and frontier MO theory, respectively. The introduction of polarizability/softness is also considered in a historical perspective in which polarizability was, with some exceptions, mainly put forward in non covalent interactions. A particular series of response functions, arising when the changes in the external potential are solely provoked by changes in nuclear configurations (the "R-analogues") are also systematically considered. The main part of the contribution is devoted to third order response functions which, at first sight, may be expected not to yield chemically significant information, as turns out to be for the hyperhardness. A counterexample is the dual descriptor and its R analogue, the initial hardness response, which turns out to yield a firm basis to regain the Woodward-Hoffmann rules for pericyclic reactions based on a density-only basis, i.e. without involving the phase, sign, symmetry of the wavefunction. Even the second order nonlinear response functions are shown possibly to bear interesting information, e.g. on the local and global polarizability. Its derivatives may govern the influence of charge on the polarizability, the R-analogues being the nuclear Fukui function and the quadratic and cubic force constants. Although some of the higher order derivatives may be difficult to evaluate a comparison with the energy expansion used in spectroscopy in terms of nuclear displacements, nuclear magnetic moments, electric and magnetic fields leads to the conjecture that, certainly cross terms may contain new, intricate information for understanding chemical reactivity.

Condensed-to-atoms electronic Fukui functions within the framework of spin-polarized density-functional theory

A simple formalism devised to calculate the condensed-to-atoms Fukui function [R. R. Contreras, P. Fuentealba, M. Galván, and P. Pérez, Chem. Phys. Lett. 304, 405 (1999)] has been further analyzed within a spin-polarized density-functional theory framework. The model is based on a frozen-core approximation to these local reactivity indices [M. Galván, A. Vela, and J. L. Gázquez, J. Phys. Chem. 92, 6470 (1988)], giving us an extended local reactivity description of systems based on the frontier spin-up and spin-down molecular orbitals. Degenerate molecular spin orbitals have been explicitly included in our model equations. Computational results for the nitric oxide (NO) and some simple carbene systems are presented in order to test the model. These quantities have been discussed in the context of changes both in charge density and spin density within the context of electron charge transfer or spin-polarization processes.

The physical basis of the hard/soft acid/base principle

The dependence of molecular properties on the chemical hardness is explored. As the chemical hardness of a molecule increases, its size and polarizability typically decrease and its charge and electronegativity typically increase. On the basis of these properties, the interaction energy between hard and soft acids and bases is quantified, and the physical basis of the global and local hard/soft acid/base (HSAB) principles is elucidated.

Generalized nuclear Fukui functions in the framework of spin-polarized density-functional theory

An extension of Cohen's nuclear Fukui function is presented in the spin-polarized framework of density-functional theory (SP-DFT). The resulting new nuclear Fukui function indices PhiNalpha and PhiSalpha are intended to be the natural descriptors for the responses of the nuclei to changes involving charge transfer at constant multiplicity and also the spin polarization at constant number of electrons. These generalized quantities allow us to gain new insights within a perturbative scheme based on DFT. Calculations of the electronic and nuclear SP-DFT quantities are presented within a Kohn-Sham framework of chemical reactivity for a sample of molecules, including H2O, H2CO, and some simple nitrenes (NX) and phosphinidenes (PX), with X=H, Li, F, Cl, OH, SH, NH2, and PH2. Results have been interpreted in terms of chemical bonding in the context of Berlin's theorem, which provides a separation of the molecular space into binding and antibinding regions.

Nuclear reactivity indices within regional density functional theory

Regional chemical potential values-mu(R) have been obtained with the use of nuclear reactivity indices. Perturbational formulae use values of reactivity indices of isolated molecular fragments. The changes of the parameters (DeltaN(R),{ DeltaQ(i)}i epsilonR) within each fragment determine the value of the regional chemical potential after a chemical reaction. The computational scheme has been tested numerically along the chemical reaction path. We have studied a set of chemical reactions to obtain regional chemical potentials (mu(ts)R) and regional transfer potentials (T(ts)R) for transition states of the following chemical reactions: HF+CO=HFCO, HCl+CO=HClCO, HF+SiO=HFSiO and HF+GeO=HFGeO. The results are reasonable and encouraging. Values of these indices show the possible reactivity directions of the transition states examined.

On the importance of the “density per particle” (shape function) in the density functional theory

The central role of the shape function sigma(r) from the density functional theory (DFT), the ratio of the electron density rho(r) and the number of electrons N of the system (density per particle), is investigated. Moreover, its relationship with DFT based reactivity indices is established. In the first part, it is shown that an estimate for the chemical hardness can be obtained from the long range behavior of the shape function and its derivative with respect to the number of electrons at a fixed external potential. Next, the energy of the system is minimized with the constraint that the shape function should integrate to unity; the associated Lagrange multiplier is shown to be related to the electronic chemical potential micro of the system. Finally, the importance of the shape function for both molecular structure, reactivity, and similarity is outlined.

Variational principles for describing chemical reactions. Reactivity indices based on the external potential

In a recent paper [J. Am. Chem. Soc. 2000, 122, 2010], the authors explored variational principles that help one understand chemical reactivity on the basis of the changes in electron density associated with a chemical reaction. Here, similar methods are used to explore the effect changing the external potential has on chemical reactivity. Four new indices are defined: (1) a potential energy surface that results from the second-order truncation of the Taylor series in the external potential about some reference, Upsilon(R(1),R(2),.,R(M)()); (2) the stabilization energy for the equilibrium nuclear geometry (relative to some reference), Xi; (3) the flexibility, or "lability", of the molecule at equilibrium, Lambda; and (4) the proton hardness, Pi, which performs a role in the theory of Brönsted-Lowry acids and bases that is similar to the role of the chemical hardness in the theory of Lewis acids and bases. Applications considered include the orientation of a molecule in an external electric field, molecular association reactions, and reactions between Brönsted-Lowry acids and bases.