Local and linear chemical reactivity response functions at finite temperature in density functional theory

Spin-Polarized Conceptual Density Functional Theory from the Convex Hull

We present a new, nonarbitrary, internally consistent, and unambiguous framework for spin-polarized conceptual density-functional theory (SP-DFT). We explicitly characterize the convex hull of energy, as a function of the number of electrons and their spin, as the only accessible ground states in spin-polarized density functional theory. Then, we construct continuous linear and quadratic models for the energy. The nondifferentiable linear model exactly captures the simplicial geometry of the complex hull about the point of interest and gives exact representations for the conceptual DFT reactivity indicators. The continuous quadratic energy model is the paraboloid of maximum curvature, which most tightly encloses the point of interest and neighboring vertices. The quadratic model is invariant to the choice of coordinate system (i.e., {N, S} vs {Nα, Nβ}) and reduces to a sensible formulation of spin-free conceptual DFT in the appropriate limit. Using the quadratic model, we generalize the Parr functions {P+(r), P-(r)} (and their derivatives with respect to number of electrons) to this new spin-polarized framework, integrating the Parr function concept into the context of (spin-polarized) conceptual DFT, and extending it to include higher-order effects.

Exploring the Reactivity of Donor-Acceptor Systems through a Combined Conceptual and Constrained DFT Approach

In the context of the conceptual density functional theory (cDFT) and based on the computational efficiency of the constrained DFT (CDFT), we demonstrate that chemical reactivity can be governed by the difference between the local interacting chemical potentials of the reactants (referred as Edual), in agreement with Sanderson's equalization principle. In a proof-of-concept study, we investigated illustrative examples involving typical non-covalent donor-acceptor systems and reactive systems are provided. For the selected systems, our approach reveals significant mimicking between Edual and the DFT-computed intermolecular interaction energy profiles. We further evaluate the influence of the Coulomb and exchange-correlation contributions in Edual. These latter results suggest that numerous potential energy surfaces of clusters can be explored using a Sanderson-like model only based on classical interactions between molecular orbitals domains. To conclude, this study achieved a deeper understanding of the principles of cDFT and assessed, in a wider context, its efficiency in predicting the chemical reactivity.

Properties of the density functional response kernels and its implications on chemistry

An overview of mathematical properties of the non-local second order derivatives of the canonical, grand canonical, isomorphic, and grand isomorphic ensembles is given. The significance of their positive or negative semidefiniteness and the implications of these properties for atoms and molecules are discussed. Based on this property, many other interesting properties can be derived, such as the expansion in eigenfunctions, bounds on the diagonal and off-diagonal elements, and the eigenvalues of these kernels. We also prove Kato's theorem for the softness kernel and linear response and the dissociation limit of the linear responses as the sum of the linear responses of the individual fragments when dissociating a system into two non-interacting molecular fragments. Finally, strategies for the practical calculation of these kernels, their eigenfunctions, and their eigenvalues are discussed.

From Density Functional Theory to Conceptual Density Functional Theory and Biosystems

The position of conceptual density functional theory (CDFT) in the history of density functional theory (DFT) is sketched followed by a chronological report on the introduction of the various DFT descriptors such as the electronegativity, hardness, softness, Fukui function, local version of softness and hardness, dual descriptor, linear response function, and softness kernel. Through a perturbational approach they can all be characterized as response functions, reflecting the intrinsic reactivity of an atom or molecule upon perturbation by a different system, including recent extensions by external fields. Derived descriptors such as the electrophilicity or generalized philicity, derived from the nature of the energy vs. N behavior, complete this picture. These descriptors can be used as such or in the context of principles such as Sanderson’s electronegativity equalization principle, Pearson’s hard and soft acids and bases principle, the maximum hardness, and more recently, the minimum electrophilicity principle. CDFT has known an ever-growing use in various subdisciplines of chemistry: from organic to inorganic chemistry, from polymer to materials chemistry, and from catalysis to nanotechnology. The increasing size of the systems under study has been coped with thanks to methodological evolutions but also through the impressive evolution in software and hardware. In this flow, biosystems entered the application portfolio in the past twenty years with studies varying (among others) from enzymatic catalysis to biological activity and/or the toxicity of organic molecules and to computational peptidology. On the basis of this evolution, one can expect that “the best is yet to come”.

Molecular interactions from the density functional theory for chemical reactivity: Interaction chemical potential, hardness, and reactivity principles

In the first paper of this series, the authors derived an expression for the interaction energy between two reagents in terms of the chemical reactivity indicators that can be derived from density functional perturbation theory. While negative interaction energies can explain reactivity, reactivity is often more simply explained using the “|dμ| big is good” rule or the maximum hardness principle. Expressions for the change in chemical potential (μ) and hardness when two reagents interact are derived. A partial justification for the maximum hardness principle is that the terms that appear in the interaction energy expression often reappear in the expression for the interaction hardness, but with opposite sign.

Molecular Interactions From the Density Functional Theory for Chemical Reactivity: The Interaction Energy Between Two-Reagents

Reactivity descriptors indicate where a reagent is most reactive and how it is most likely to react. However, a reaction will only occur when the reagent encounters a suitable reaction partner. Determining whether a pair of reagents is well-matched requires developing reactivity rules that depend on both reagents. This can be achieved using the expression for the minimum-interaction-energy obtained from the density functional reactivity theory. Different terms in this expression will be dominant in different circumstances; depending on which terms control the reactivity, different reactivity indicators will be preferred.

From the Electron Density Gradient to the Quantitative Reactivity Indicators: Local Softness and the Fukui Function

Important reactivity measures such as the local softness, the Fukui function, and the global hardness have been calculated directly from first principles with the use of the electron density function, beyond the finite difference approximation. Our recently derived density gradient theorem and the principle of nearsightedness of the electronic matter have been instrumental in obtaining the original, albeit approximate, result on the local softness of an atom. By integration of the local softness s(r), we obtain the global softness S and the Fukui function f(r) = s(r)/S. Local and global softness values have also been calculated analytically for the basic hydrogenic orbitals; the general relation to the atomic number S = σZ ^-2 has been demonstrated, with constants σ characteristic for each orbital type. Global hardness η = 1/S calculated for atoms and ions has been favorably tested against its conventional measure given by the finite difference approximation: (I - A). Calculated test results for atoms and ions in rows 1-4 of the periodic table have been presented.

Extending conceptual DFT to include external variables: the influence of magnetic fields

An extension of conceptual DFT to include the influence of an external magnetic field is proposed in the context of a program set up to cope with the ever increasing variability of reaction conditions and concomitant reactivity.

Electronic structure behavior of PbO2, IrO2, and SnO2 metal oxide surfaces (110) with dissociatively adsorbed water molecules as a function of the chemical potential

A detailed analysis of the electronic structure of three different electrochemical interfaces as a function of the chemical potential (μ) is performed using the grand canonical density functional theory in the joint density functional theory formulation. Changes in the average number of electrons and the density of states are also described. The evaluation of the global softness, which measures the tendency of the system to gain or lose electrons, is straightforward under this formalism. The observed behavior of these quantities depends on the electronic nature of the electrochemical interfaces.

Extending conceptual DFT to include additional variables: oriented external electric field

The extension of the E = E[N, v] functional for exploring chemical reactivity in a conceptual DFT context to include external electric fields is discussed. Concentrating on the case of a homogeneous field the corresponding response functions are identified and integrated, together with the conventional response functions such as permanent dipole moment and polarizability, in an extended response function tree associated with the E = E[N, v, ε] functional. In a case study on the dihalogens F2, Cl2, Br2, I2 the sensitivity of condensed atomic charges (∂q/∂ε) is linked to the polarizability of the halogen atoms. The non-integrated (∂ρ(r)/∂ε) response function, directly related to the field induced density change, is at the basis of these features. It reveals symmetry breaking for a perpendicular field, not detectable in its atom condensed counterpart, and accounts for the induced dipole moment directly related to the molecular polarizability. The much higher sensitivity of the electronic chemical potential/electronegativity as compared to the chemical hardness is highlighted. The response of the condensed Fukui functions to a parallel electric field increases when going down in the periodic table and is interpreted in terms of the extension of the outer contours in the non-condensed Fukui function. In the case of a perpendicular field the (∂f(r)/∂ε) response function hints at stereoselectivity with a preferential side of attack which is not retrieved in its condensed form. In an application the nucleophilic attack on the carbonyl group in H2CO is discussed. Similar to the dihalogens, stereoselectivity is displayed in the Fukui function for nucleophilic attack (f+) in the case of a perpendicular electric field, and opposite to the one that would arise based on the induced density. Disentangling the expression for the evolution of the Fukui function in the presence of an electric field reveals that this difference can be traced back to local differences in the polarization or induced density between the anionic and the neutral system. This difference may be exploited, e.g. for an appropriately substituted H2CO, to generate enantioselectivity.

Temperature-Dependent Approach to Electronic Charge Transfer

A charge transfer model is developed within the framework of the grand canonical ensemble through the analysis of the behavior of the fractional charge as a function of the chemical potential of the bath when the temperature and the external chemical potential are kept fixed. Departing from the fact that, before the interaction between two species, each one has a zero fractional charge, one can identify two situations after the interaction occurs where the fractional charge of at least one of the species is different from zero, indicating that there has been charge transference. One of them corresponds to the case when one of the species is immersed in a bath conformed by the other one, while the other is related to the case in which both species are present in equal amounts (stoichiometric proportion). Correlations between the fractional charges and average energies, thus obtained with experimental equilibrium constants, kinetic rate constants, hydration constants, and bond enthalpies, indicate that, although at the experimental temperatures, they are very small quantities, they have chemically meaningful information. Additionally, in the stoichiometric case, one also finds a rather good correlation between the equalized chemical potential and the one obtained from experimental information for a test set of diatomic and triatomic molecules.

Establishing the Relationship between Quantum Capacitance and Softness of N-Doped Graphene/Electrolyte Interfaces within the Density Functional Theory Grand Canonical Kohn-Sham Formalism

The joint density functional theory (JDFT) is applied in the context of the grand canonical Kohn-Sham theory to calculate the global and local softness of pristine and N-substituted graphene structures. A comparison is established between the different theoretical approaches to evaluate total capacitance, revealing that the JDFT approach presents the closest result of this property with experimental data. A model of series capacitors is used to determine the quantum and nonquantum contributions of total capacitance, which enables us to determine the limitations of the rigid band approximation for the studied systems. It is found that global chemical softness is proportional to the total capacitance measured in the experiments, when the geometry relaxation is neglected. In this context, it is possible to obtain quantum and total capacitance (and consequently softness) from an average number of electrons vs applied potential plots and the model of series capacitors. Likewise, the relation of capacitance and softness gives rise to a new definition of local capacitance within the JDFT formalism. The evaluation of global and local softness paves the way to analyze electrochemical surface reactivity as a function of applied potential for a solid-electrolyte interface.

An electronic temperature definition for the reactive electronic species: Conciliating practical approaches in conceptual chemical reactivity theory with a rigorous ensemble formulation

By working under the framework of the Helmholtz potential as a functional of the equilibrium density matrix, in this contribution, we provide theoretical evidence about a particular thermodynamic situation, where electronic species display their highest susceptibility to exchange electrons to or from surroundings. This situation is denominated as the electronic temperature condition. Neutral chemical species display their lowest possible hardness value at the electronic temperature condition, and remarkably, under this circumstance, the exchange of any amount of electronic charge will necessarily be translated into a net increase in the corresponding chemical hardness. Chemical response functions defined as partial derivatives of the Helmholtz potential with respect to the (average) number of electrons and evaluated at the electronic temperature condition provide comparable results than those obtained from the coarse quadratic approximation to the exact dependence of the electronic energy vs the number of electrons, including composite quantities as the electrophilicity index. In this context, we show that the exact Helmholtz potential dependence with respect to the number of electrons can accurately be approximated by "temperature dependent" polynomial fits (up to fourth order), evaluated at the electronic temperature condition.

Local electrophilicity

In this work some possibilities for deriving a local electrophilicity are studied. First, we consider the original definition proposed by Chattaraj, Maiti, and Sarkar (J Phys Chem A 107:4973, 2003), in which the local electrophilicity is given by the product of the global electrophilicity, and the Fukui function for charge acceptance is derived by two different approaches, making use of the chain rule for functional derivatives. We also modify the proposals based on the electron density so as to have a definition with the same units of the original definition, which also introduces a dependence in the Fukui function for charge donation. Additionally, we also explore other possibilities using the tools of information theory and the temperature dependent reactivity indices of the density functional theory of chemical reactivity. The poor results obtained from the last two approaches lead us to conjecture that this is due to the fact that the global electrophilicity is not a derivative, like most of the other reactivity indices. The conclusion is that Chattaraj's suggestion seems to be the simplest, but at the same time a very reliable approach to this important property.

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.

A logarithmic informational temperature scale as a criterion to classify molecular systems

In this work, we applied a definition of informational energy and informational temperature using 1865 molecules and showed that such definitions can classify molecules with similar chemical properties such as hardness, softness and chemical potential.

Going beyond the three-state ensemble model: the electronic chemical potential and Fukui function for the general case

The analytical working equations for the chemical potential and the Fukui function for the case of any number of ground and excited states is presented.

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.

Electronegativity and redox reactions

Using the maximum hardness principle, we show that the oxidation potential of a molecule increases as its electronegativity increases and also increases as its electronegativity in its oxidized state increases.

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.