Reactivity indicators for degenerate states in the density-functional theoretic chemical reactivity 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.

Connection between nuclear and electronic Fukui functions beyond frontier molecular orbitals

Based on the relationship between average local ionization energy Ī(r) and average local electron affinity Ā(r) with the electronic Fukui functions, i.e., f-(r) and f+(r), respectively, in this paper, we establish a connection between nuclear and electronic Fukui functions beyond frontier molecular orbitals. As a consequence of this connection, we obtain expressions of average nuclear Fukui functions interpreted as a variation of average nucleophilicity or electrophilicity (weighted by the electronic orbital Fukui functions) with respect to nuclear displacements, which goes beyond the highest occupied molecular orbital/or lowest unoccupied molecular orbital consideration. Furthermore, from this connection and considering the frontier molecular orbital approximation, we derive expressions of nuclear Fukui functions in terms of the atom-condensed electronic Fukui functions, which imply a locality in the chemical reactivity and could be used to study the variation of local nucleophilicity or electrophilicity with respect to nuclear displacements. Finally, this new way to interpret the nuclear Fukui function could be useful in the future to study the chemical reactivity related to molecular vibrations, internal rotations, bond dissociation, chemical reaction along the model of reaction coordinate, and so on.

A density fitting scheme for the fast evaluation of molecular electrostatic potential

Molecular electrostatic potential (MEP) is a significant and crucial physical quantity that can be applied to a large number of scenarios, such as the prediction of nucleophilic or electrophilic attacks, fitting atomic charges, σ-hole, and so forth. The computational cost for the MEP has an O(N² ) scaling with the increase of atoms, which is intractable and laborious for macromolecules. Herein, a density fitting molecular electrostatic potential (DF-MEP) is used to reduce the computational costs for the macromolecular MEP. It is found that the accuracy of DF-MEP is almost identical to the conventional molecular electrostatic potential (Conv-MEP), while the computational costs can be reduced to an O(N) scaling, for example, the computational time of 699,200 grids for the Trp-cage molecule (304 atoms) only takes 16.6 s at the B3LYP-D3(BJ)/def2-SVP level of theory with 16 CPU cores compared with 3060.2 s for the Conv-MEP method.

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.

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.

On the Prediction of Lattice Energy with the Fukui Potential: Some Supports on Hardness Maximization in Inorganic Solids

Using perturbation theory within the framework of conceptual density functional theory, we derive a lower bound for the lattice energy of the ionic solids. The main element of the lower bound is the Fukui potential in the nuclei of the molecule corresponding to the unit formula of the solid. Thus, we propose a model to calculate the lattice energy in terms of the Fukui potential. Our method, which is extremely simple, performs well as other methods using the crystal structure information of alkali halide solids. The method proposed here correlates surprisingly well with the experimental data on the lattice energy of a diverse series of solids having even a non-negligible covalent characteristic. Finally, the validity of the maximum hardness principle (MHP) is assessed, showing that in this case, the MHP is limited.

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.

Electronegativity under Confinement

The electronegativity concept was first formulated by Pauling in the first half of the 20th century to explain quantitatively the properties of chemical bonds between different types of atoms. Today, it is widely known that, in high-pressure regimes, the reactivity properties of atoms can change, and, thus, the bond patterns in molecules and solids are affected. In this work, we studied the effects of high pressure modeled by a confining potential on different definitions of electronegativity and, additionally, tested the accuracy of first-order perturbation theory in the context of density functional theory for confined atoms of the second row at the Hartree–Fock level. As expected, the electronegativity of atoms at high confinement is very different than that of their free counterparts since it depends on the electronic configuration of the atom, and, thus, its periodicity is modified at higher pressures.

Proposal of a Fermi-Dirac-Derived Reactivity Descriptor: Beyond the Frontier MO Model

In this paper, we derive a reactivity descriptor stemming from the Fermi-Dirac population scheme, applied to density functional calculations on molecular systems. Assuming that molecular orbitals only marginally change when temperature is slightly increased from 0 K, we study the response of electron density to a change in temperature. Connection with usual conceptual density functional theory descriptors is made, and the T-variation of electron density for some representative examples is given 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.

Kick-Fukui: A Fukui Function-Guided Method for Molecular Structure Prediction

Here, we introduce a hybrid method, named Kick-Fukui, to explore the potential energy surface (PES) of clusters and molecules using the Coulombic integral between the Fukui functions in the first screening of the best individuals. In the process, small stable molecules or clusters whose combination has the stoichiometry of the explored species are used as assembly units. First, a small set of candidates has been selected from a large and stochastically generated (Kick) population according to the maximum value of the Coulombic integral between the Fukui functions of both fragments. Subsequently, these few candidates are optimized using a gradient method and density functional theory (DFT) calculations. The performance of the program has been evaluated to explore the PES of various systems, including atomic and molecular clusters. In most cases studied, the global minimum (GM) has been identified with a low computational cost. The strategy does not allow to identify the GM of some silicon clusters; however, it predicts local minima very close in energy to the GM that could be used as the initial population of evolutionary algorithms.

High regioselectivity in the amination reaction of isoquinolinequinone derivatives using conceptual DFT and NCI analysis

DFT-derived reactivity descriptors and Non-Covalent interaction (NCI) analysis were performed to rationalize the regioselectivity in the amination reaction of some isoquinolinequinone derivatives. Statistical analysis was performed to assess robustness of atomic charges to the basis set. Various electronic population schemes including Mulliken population analysis (MPA), electrostatic method (ChelpG), Hirshfeld population analysis (HPA) and Natural population analysis (NPA) have been considered. The results revealed that NPA was the most efficient for this purpose. NCI study using the reduced density gradient (RDG) was performed for revealing weak interactions. Domains defined by isosurfaces of RDG have been integrated to quantitatively study the strength of weak interactions and their stabilities have also been examined. Steric hindrance caused by coordination of ethanol with the neighboring carbonyl prevents the nucleophilic attack on C-6 and therefore leads to preferential C-7 substitution. The quantitative study of NCI clearly demonstrates that the hydrogen bond of carbonyl (2) is more stabilizing. Consequently, the high polarity of hydrogen bond on this carbonyl may explain the high electrophilicity of C-5 compared to C-8. Our work proved that the difference in local reactivity, as well as the steric hindrance are the key elements explaining the high regioselectivity exhibited by the amination reaction.

Unifying Conceptual Density Functional and Valence Bond Theory: The Hardness-Softness Conundrum Associated with Protonation Reactions and Uncovering Complementary Reactivity Modes

In this study, we address the long-standing issue-arising prominently from conceptual density functional theory (CDFT)-of the relative importance of electrostatic, i.e., "hard-hard", versus spin-pairing, i.e., "soft-soft", interactions in determining regiochemical preferences. We do so from a valence bond (VB) perspective and demonstrate that VB theory readily enables a clear-cut resolution of both of these contributions to the bond formation/breaking process. Our calculations indicate that appropriate local reactivity descriptors can be used to gauge the magnitude of both interactions individually, e.g., Fukui functions or HOMO/LUMO orbitals for the spin-pairing/(frontier) orbital interactions and molecular electrostatic potentials (and/or partial charges) for the electrostatic interactions. In contrast to previous reports, we find that protonation reactions cannot generally be classified as either charge- or frontier orbital-controlled; instead, our results indicate that these two bonding contributions generally interplay in more subtle patterns, only giving the impression of a clear-cut dichotomy. Finally, we demonstrate that important covalent, i.e., spin pairing, reactivity modes can be missed when only a single spin-pairing/orbital interaction descriptor is considered. This study constitutes an important step in the unification of CDFT and VB theory.

Predicting Deprotonation Sites Using Alchemical Derivatives

An alchemical transformation is any process, physical or fictitious, that connects two points in the chemical space. A particularly important transformation is the vanishing of a proton, whose energy can be linked to the proton dissociation enthalpy of acids. In this work we assess the reliability of alchemical derivatives in predicting the proton dissociation enthalpy of a diverse series of mono- and polyprotic molecules. Alchemical derivatives perform remarkably well in ranking the proton affinity of all molecules. Additionally, alchemical derivatives could be use also as a predictive tool because their predictions correlate quite well with calculations based on energy differences and experimental values. Although second-order alchemical derivatives underestimate the dissociation enthalpy, the deviation seems to be almost constant. This makes alchemical derivatives extremely accurate to evaluate the difference in proton affinity between two acid sites of polyprotic molecule. Finally, we show that the reason for the underestimation of the dissociation enthalpy is most likely the contribution of higher-order derivatives.

The Role of the Density Response Kernel in the Protonation Process

The interaction between a Lewis base and a proton produces large changes in the electron distribution of the proton-accepting species. The quantum description of the change in the electronic properties can be predicted by perturbative models. Density functional theory-based reactivity indices are useful quantities to predict the changes in the energy and the electron density of electronic systems. However, the perturbation produced by the proton, namely, the electric field of a bare nucleus, usually leads to the use of perturbative terms beyond the first order. In this work, we analyze the effect of the protonation on the electronic structure of different Lewis bases. We also identify the leading term in the perturbative expansion that determines the protonation site and try to relate it with the chemical nature of the Lewis base. Even when the interaction is initially dominated by the electrostatic forces, we found that the electron redistribution effects can play a more relevant role for highly polarizable species.

Orbital-Weighted Dual Descriptor for the Study of Local Reactivity of Systems with (Quasi-) Degenerate States

An alternative response function, based on the dual descriptor in terms of Koopmans' approximation, is hereby proposed for the description of chemical reactivity in systems with (quasi-) degenerate frontier molecular orbitals. This descriptor is constructed from Fukui functions that include contributions from different orbitals, i.e., orbital-weighted Fukui functions. The methodology is applied to three case studies: the first case consists of a series of benchmark organic and inorganic molecules from which the dual descriptor, based only on frontier orbitals, is not appropriate to describe their reactivity. The second case deals with the proper description of chemical reactivity in Diels-Alder reactions between fullerene C₆₀ and cyclopentadiene (CP), revealing the importance of considering secondary orbital interactions for an adequate regioselectivity description. The third, and last case, consists of a series of polycyclic aromatic hydrocarbons (PAHs) possessing molecular orbital degeneracy. By means of analyzing of this descriptor, an alternative approach to the description of aromaticity is proposed. In all cases, the proposed index called "orbital-weighted dual descriptor" has proven to accurately describe the chemical reactivity and aromaticity of the studied systems.

New advances in conceptual-DFT: an alternative way to calculate the Fukui function and dual descriptor

An alternative way of calculating the Fukui function and the partial derivative of second order of the electronic density with respect to the number of electrons N is presented, the new formulas agree with the usual ones but only in cases without degeneracy. The new operative formulas are more general than the previous ones and are the right ones for those problematic cases where one or both of the frontier molecular orbitals are degenerate. Finally, we present a new way of applying the finite difference approximation that leads to more realistic results than the usual formulas. Graphical abstract A new way of calculating the Fukui function is presented that results in a new operative formula of the function. It has also been obtained the partial derivative of second order of the electronic density with respect to the number of electrons N, and it agree with the usual formula of the dual descriptor function but only in cases without degeneration.

How Biochemical Environments Fine-Tune a Redox Process: From Theoretical Models to Practical Applications

In this study, we give a new physical insight into how enzymatic environments influence a redox process. This is particularly important in a biochemical context, in which oxidoreductase enzymes and low-molecular-weight cofactors create a microenvironment, fine-tuning their specific redox potential. We present a new theoretical model, quantitatively backed up by quantum chemically calculated data obtained for key biological sulfur-based model reactions involved in preserving the cellular redox homeostasis during oxidative stress. We show that environmental effects can be quantitatively predicted from the thermodynamic cycle linking ΔΔ G(OX/RED)_ref-ligand values to the differential interaction energy ΔΔ G_int of the reduced and oxidized species with the environment. Our obtained data can be linked to hydrogen-bond patterns found in protein active sites. The thermodynamic model is further understood in the framework of molecular orbital theory. The key insight of this work is that the intrinsic properties of neither a redox couple nor the interacting environment (e.g., ligand) are enough by themselves to uniquely predict reduction potentials. Instead, system-environment interactions need to be considered. This study is of general interest as redox processes are pivotal to empower, protect, or damage organisms. Our presented thermodynamic model allows a pragmatically evaluation on the expected influence of a particular environment on a redox process, necessary to fully understand how redox processes take place in living organisms.

Benchmark values of chemical potential and chemical hardness for atoms and atomic ions (including unstable anions) from the energies of isoelectronic series

We present benchmark values for the electronic chemical potential and chemical hardness from reference data for ionization potentials and electron affinities. In cases where the energies needed to compute these quantities are not available, we estimate the ionization potential of the metastable (di)anions by extrapolation along the isoelectronic series, taking care to ensure that the extrapolated data satisfy reasonable intuitive rules to the maximum possible extent. We also propose suitable values for the chemical potential and chemical hardness of zero-electron species. Because the values we report are faithful to the trends in accurate data on atomic energies, we believe that our proposed values for the chemical potential and chemical hardness are ideally suited to conceptual studies of chemical trends across the periodic table. The critical nuclear charge (Z critical) of the isoelectronic series with 2 < N < 96 has also been reported for the first time.

Analysis of molecular and (di)atomic dual-descriptor functions and matrices

In this work, the dual-descriptor is studied in matrix form [Formula: see text] and both coordinates condensed to atoms, resulting in atomic and diatomic (or where applicable, bond) condensed single values. This double partitioning method of the dual-descriptor matrix is proposed within the Hirshfeld-I atoms-in-molecule framework although it is easily extended to other atoms-in-molecules methods. Diagonalizing the resulting atomic and bond dual-descriptor matrices gives eigenvalues and eigenvectors describing the reactivity of atoms and bonds. The dual-descriptor function is the diagonal element of the underlying matrix. The extra information contained in the atom and bond resolution is highlighted and the effect of choosing either the fragment of molecular response or response of molecular fragment approach is quantified. Graphical Abstract Atom and bond condensed dual descriptor matrices and functions are derived from molecular ones using Hirshfeld-I atoms in molecules weight functions.

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.

Atom and Bond Fukui Functions and Matrices: A Hirshfeld-I Atoms-in-Molecule Approach

The Fukui function is often used in its atom-condensed form by isolating it from the molecular Fukui function using a chosen weight function for the atom in the molecule. Recently, Fukui functions and matrices for both atoms and bonds separately were introduced for semiempirical and ab initio levels of theory using Hückel and Mulliken atoms-in-molecule models. In this work, a double partitioning method of the Fukui matrix is proposed within the Hirshfeld-I atoms-in-molecule framework. Diagonalizing the resulting atomic and bond matrices gives eigenvalues and eigenvectors (Fukui orbitals) describing the reactivity of atoms and bonds. The Fukui function is the diagonal element of the Fukui matrix and may be resolved in atom and bond contributions. The extra information contained in the atom and bond resolution of the Fukui matrices and functions is highlighted. The effect of the choice of weight function arising from the Hirshfeld-I approach to obtain atom- and bond-condensed Fukui functions is studied. A comparison of the results with those generated by using the Mulliken atoms-in-molecule approach shows low correlation between the two partitioning schemes.

Towards the rationalization of catalytic activity values by means of local hyper-softness on the catalytic site: a criticism about the use of net electric charges

By means of the Spin-Polarized Conceptual Density Functional Theory (SP-CDFT), three 2,6-bis(imino)pyridine catalysts based on iron(II), used for polymerization of ethylene, were studied. The catalysts differed by the substituent group, bearing either -H, -NO2 or -OCH3. To date, catalytic activity, a purely experimental parameter measuring the mass of polyethylene produced per millimole of iron per time and pressure unit at a fixed temperature, has not been explained in terms of local hyper-softness. The latter is a purely theoretical parameter designed for quantifying electronic effects; it is measured using the metal atom responsible for the coordination process with the monomer (ethylene). Because steric effects are not relevant in these kinds of catalysts and only electronic effects drive the catalytic process, an interesting link is found between catalytic activity and the local hyper-softness condensed on the iron atom by means of four functionals (B3LYP, BP86, B97D, and VSXC). This work demonstrates that the use of local hyper-softness, predicted by the SP-CDFT, is a suitable parameter for explaining order relationships among catalytic activity values, thus quantifying the electronic influence of the substituent group inducing this difference; the use of only net electric charges does not lead to clear conclusions. This finding can aid in estimating catalytic activities leading to a more rational design of new catalysts via computational chemistry.

Direct computation of parameters for accurate polarizable force fields

We present an improved electronic linear response model to incorporate polarization and charge-transfer effects in polarizable force fields. This model is a generalization of the Atom-Condensed Kohn-Sham Density Functional Theory (DFT), approximated to second order (ACKS2): it can now be defined with any underlying variational theory (next to KS-DFT) and it can include atomic multipoles and off-center basis functions. Parameters in this model are computed efficiently as expectation values of an electronic wavefunction, obviating the need for their calibration, regularization, and manual tuning. In the limit of a complete density and potential basis set in the ACKS2 model, the linear response properties of the underlying theory for a given molecular geometry are reproduced exactly. A numerical validation with a test set of 110 molecules shows that very accurate models can already be obtained with fluctuating charges and dipoles. These features greatly facilitate the development of polarizable force fields.

An information-theoretic resolution of the ambiguity in the local hardness

A definition of the local hardness, suitable for application in the local hard/soft acid/base principle, is derived by applying information theory.

Computational nanochemistry report on the oxicams--conceptual DFT indices and chemical reactivity

A density functional theory study of eight oxicams was carried out in order to determine their global and local reactivities. These types of reactivities were measured by means of global and local reactivity descriptors coming from the conceptual density functional theory. Net electrophilicity as a global reactivity descriptor and local hypersoftness as a local reactivity descriptor were the used tools to distinguish reactivity and selectivity among these oxicams. Globally, isoxicam presents the highest electron donating capacity; meanwhile, the highest electron accepting capacity is exhibited by droxicam. Locally, two oxicams present neither nucleophilic nor electrophilic relevant reactivity in their peripheral pyridine ring, droxicam and tenoxicam, so that their more reactive zones are found on the respective fused rings. Oxicams have been divided into two subgroups in order to facilitate the local analysis of reactivity. One group is characterized because their most important condensed values for local hypersoftnes are well-separated: 4-meloxicam, lornoxicam, meloxicam, and normeloxicam. Meanwhile, the opposite situation is found in droxicam, isoxicam, piroxicam, and tenoxicam. As a whole, the nucleophilic characteristic noticeably predominates in these eight oxicams instead of an electrophilic behavior, thus meaning a greater tendency to donate electrons rather than withdrawing them; a consequence of this behavior implies a favorable interaction with a hypothetical receptor bearing one or more electron acceptor functional groups rather than electron donor functional groups; this would imply a maximization of this interaction from the covalent point of view.

Computational Nutraceutics: Chemical Reactivity Properties of the Flavonoid Naringin by Means of Conceptual DFT

The M06 family of density functionals has been assessed for the calculation of the molecular structure and properties of the Naringin molecule. The chemical reactivity descriptors have been calculated through Conceptual DFT. The active sites for nucleophilic and electrophilic attacks have been chosen by relating them to the Fukui function indices and the dual descriptorf(2)(r). A comparison between the descriptors calculated through vertical energy values and those arising from the Koopmans' theorem approximation has been performed in order to check for the validity of the last procedure.