Abstract
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.
1 Introduction
Nearly 35 years ago, Parr recognized that density-functional theory (DFT) could be used not only as a formal alternative to wavefunction-based quantum chemistry and as a computational tool, but also as an interpretative tool through which chemical reactivity could be elucidated (Parr et al., 1978; Parr and Yang, 1989). The special power of the electron-density perspective arises because the mathematical structure of DFT naturally accommodates fractional numbers of electrons and therefore partial electron transfer (Janak, 1978; Parr and Bartolotti, 1982; Perdew et al., 1982; Yang et al., 2000; Ayers, 2008; Fuentealba and Cardenas, 2013; Miranda-Quintana and Bochicchio, 2014; Miranda-Quintana and Ayers, 2016a; Ayers and Mel, 2018), while in the wavefunction theory, the number of electrons is linked to the dimensionality of the wavefunction, and is inherently an integer. Ironically, the great utility of DFT for conceptual purposes arises from the same feature that is most problematic for computational applications (Merkle et al., 1992; Mori-Sanchez et al., 2006; Ruzsinszky et al., 2007; Cohen et al., 2008a; b;Mori-Sanchez et al., 2008; Yang et al., 2016). However, the difficulty of treating fractional electrons computationally in DFT is not important in the context of this paper, where one can assume that an exact (or otherwise accurate ab initio) density functional is used (Levy, 1979; Lieb, 1983; Bartlett et al., 2005; Ayers, 2006).
The use of density functional theory for chemical reactivity (DFT-CR), often called conceptual DFT or chemical DFT (Parr and Yang, 1989; Chermette, 1999; Geerlings et al., 2003; Ayers et al., 2005; Chattaraj et al., 2006; Gazsquez, 2008; Liu, 2009a; Johnson et al., 2012; De Proft et al., 2014; Fuentealba and Cardenas, 2015; Miranda-Quintana, 2018), is now well established, both mathematically, conceptually, and computationally. The greatest successes of DFT-CR are probably linked to the definition of a quantitative scale for the chemical hardness, (Parr and Pearson, 1983), which has led to increased understanding of the hard/soft acid/base (HSAB) theory (Nalewajski, 1984; Nalewajski et al., 1988; Chattaraj et al., 1991; Gazquez and Mendez, 1994; Mendez and Gazquez, 1994; Chattaraj, 2001; Melin et al., 2004; Ayers, 2005; Ayers et al., 2006; Anderson et al., 2007a; Ayers, 2007; Cardenas and Ayers, 2013; Miranda-Quintana R. A., 2017) and an entirely new principle for chemical reactivity and molecular stability, the maximum hardness principle (MHP) (Pearson, 1987; Zhou and Parr, 1989; Parr and Chattaraj, 1991; Pearson and Palke, 1992; Pearson, 1994; Chattaraj, 1996; Pearson, 1999; Ayers and Parr, 2000; Torrent-Sucarrat et al., 2001, 2002; Miranda-Quintana R. A., 2017). One of the subtlest results in DFT-CR is that the HSAB principle is, in fact, an inexorable consequence of the maximum hardness principle (Chattaraj et al., 1991; Parr and Chattaraj, 1991; Chattaraj and Ayers, 2005; Chattaraj et al., 2007), with a similar result also connecting the HSAB and minimum electrophilicity principles (Chattaraj and Nath, 1994; Chattaraj et al., 2001; Chattaraj, 2007; Noorizadeh, 2007; Liu, 2009b; Morell et al., 2009; Pan et al., 2013; Miranda-Quintana, 2017a; Miranda-Quintana and Ayers, 2018b; 2019). Here it is important to clarify that what is usually understood as a “principle” in chemistry is different from the definition in physics. For example, Heisenberg’s uncertainty principle does not admit violations. The principles of chemical reactivity, on the other hand, may have exceptions. Therefore, instead of speaking of principles, one should speak of “rules of thumb” for chemical reactivity. However, for consistency with the literature, here we will use both terms as synonyms.
Almost all applications and theory in DFT-CR have been based on a one-reagent approach: the response functions of a reactant molecule are computed, and then used to predict its reactivity. Despite the usefulness of this approach, it sometimes fails. That is, sometimes understanding the inherent reactivity of a molecule is insufficient; one must discern how well-matched two reagents are. Early attempts at quantifying well-matched-ness were made by Berkowitz, Geerlings, and then, much later, by one of the present authors (Berkowitz, 1987; Langenaeker et al., 1995; Ayers et al., 2006; Anderson et al., 2007a; Ayers, 2007; Anderson et al., 2007b; Ayers and Cardenas, 2013). In the first paper of this series (Miranda-Quintana et al., 2022b), we derived a general expression for the interaction energy between two reagents using DFT-CR and drew the links to other, more computational, DFT theories like density-functional embedding (Cortona, 1991; Vaidehi et al., 1992; Wesolowski and Warshel, 1993; Govind et al., 1999; Wesolowski, 2004; Wesolowski and Leszczynski, 2006), electronegativity equalization molecular mechanics (Yang and Mortier, 1986; Mortier et al., 1985; Mortier et al., 1986; Yang and Mortier, s1986; Mortier, 1987; Rappe and Goddard, 1991; Bultinck et al., 2002a; Bultinck et al., 2002b; Verstraelen et al., 2013), and density-based energy decomposition analysis (Wu et al., 2009).
While the interaction energy provides the most fundamental perspective on chemical reactivity, sometimes it is simpler to understand chemical reactivity using alternative reactivity rules. For example, the “|dμ| big is good” (DMB) rule of Parr and Yang (Parr, 1994; Miranda-Quintana R. A. and Ayers P. W., 2018; Miranda-Quintana et al., 2018) states that favorable chemical interactions are usually associated with a large change in the chemical potential (Miranda-Quintana et al., 2018; Miranda-Quintana and Ayers, 2019; Miranda-Quintana et al., 2021). Likewise, Sanderson’s electronegativity equalization principle (Mortier et al., 1986) and Pearson’s HSAB (Pearson, 1968a; Pearson, 1968b; Nalewajski, 1984; Nalewajski et al., 1988; Chattaraj et al., 1991; Miranda-Quintana R. A., 2017; Miranda-Quintana R. A. et al., 2017), and the more recent minimum electrophilicity principles (Fuentealba et al., 2000b; Morell et al., 2009; Miranda-Quintana R. A., 2017) put these reactivity descriptors in center stage by telling us how to use them to understand and predict chemical reactivity (Sanderson, 1951; Parr et al., 1978). It is not farfetched to say that the biggest triumph of DFT-CR is not only to provide mathematically precise definitions for the reactivity descriptors, but to also give us a robust framework to derive new ones (Ayers et al., 2018; Geerlings et al., 2020). However, more often than not, these derivations have been solely based on the venerable parabolic model of Parr and Pearson (Parr and Bartolotti, 1982; Parr and Pearson, 1983; Chattaraj et al., 1995; Ayers and Parr, 2008; Alain Miranda-Quintana and Ayers, 2016; Heidar-Zadeh et al., 2016b; Miranda-Quintana and Ayers, 2016b; Cárdenas et al., 2016; Franco-Pérez et al., 2018). While powerful and hugely influential, this model can be seen as the simplest representation of electron transfer during a chemical reaction. This simplicity has often been (rightfully) argued as one of its key advantages, but this also means that elementary proofs of the HSAB (Pearson, 1968a; Chattaraj et al., 1991), minimum electrophilicity, maximum hardness, and DMB principles usually ignore electrostatic and polarization effects, and do not include charge transfer effects beyond second order (Chattaraj et al., 1991; Parr and Chattaraj, 1991; Ayers and Parr, 2000; Ayers, 2005; Chattaraj and Ayers, 2005; Ayers and Cardenas, 2013; Miranda-Quintana et al., 2018; Ayers et al., 2022). Recently, some of these approximations have been relaxed (Alain Miranda-Quintana et al., 2016; Miranda-Quintana and Ayers, 2016a; Miranda-Quintana R. A., 2017; Miranda-Quintana et al., 2021; Miranda-Quintana et al., 2022a), which strengthens the support for these principles, but the more realistic two-reagent picture remains largely unexplored (Ayers et al., 2006; Ayers, 2007; Chattaraj et al., 2007; Miranda-Quintana R. A. et al., 2017).
The chemical potential (Parr et al., 1978),measures the intrinsic Lewis acid/base strength of a molecule and can be considered to be minus one times the electronegativity. The essence of the rule, then, is that favorable molecular interactions between acids and bases “quench” the acidity/basicity of the reagents as much as possible (forming, in the extreme case, nearly inert salts). Alternatively, favorable chemical changes are associated with large changes in molecular electronegativity. Even though the “|dμ| big is good” rule was first formulated more than 30 years ago, its theoretical provenance has only recently begun to be elucidated (Miranda-Quintana et al., 2018; Miranda-Quintana and Ayers, 2019; Miranda-Quintana et al., 2021).
The maximum hardness principle (MHP) is a particularly interesting case, with roots that seem less transparent than many of the aforementioned reactivity rules. (Pearson, 1987; Pearson and Palke, 1992; Pearson, 1993; Chattaraj, 1996; Pearson, 1999). The maximum hardness principle indicates that more stable conformations are associated with a large hardness,
A corollary of this principle is that the harder a molecule is, the more stable it is. The problem is that the mathematical assumptions under which the maximum hardness principle has been proved (fixed molecular geometry and either constant electron number or constant chemical potential) do not match the conditions under which the principle is usually applied, because the MHP is most commonly used to study molecular rearrangements (Torrent-Sucarrat et al., 2001; 2002). For example, there has been substantial recent interest in using the initial hardness response (the change in hardness associated with the initial approach of two reagents) to study pericyclic reactions (De Proft et al., 2006; Ayers et al., 2007; De Proft et al., 2008; Geerlings et al., 2012).
In the remainder of this paper, we differentiate the energy expression in Eq. 33 of the first paper in this series with respect to the number of electrons. This gives the change in chemical potential (related to the first derivative) hardness (second derivative) due to the interactions between two reagents. These expressions are then used to mathematically justify the “|dμ| big is good” and maximum hardness principles.
The key expression [Eq. 33 from (Miranda-Quintana et al., 2022b)] is
In this equation, ΔUAB is the change in the total energy (U = E + Vnn, where Vnn is the nuclear-nuclear repulsion energy) when the reagents A and B come together. ΔN is the change in the number of electrons in the reagent. Superscript “0” indicates that the term is evaluated for the isolated reagent; subscripts index the reagents. As explained in (Miranda-Quintana et al., 2022b), Eq. 3 can be iterated. The only difference between Eq. 3 and the equation Eq. 33 in (Miranda-Quintana et al., 2022b), is that last it was assumed that . That is, every electron that leaves one reagent goes to the other reagent, and the number of electrons in the combined system does not change. In this paper, we will extend this analysis to consider changes in the total number of electrons,
Most of the reactivity indicators that enter into Eq. 3 are well-known in DFT-CR: the Fukui function f(r) (Parr and Yang, 1984; Yang et al., 1984; Ayers P. W. and Levy M., 2000; Heidar-Zadeh et al., 2016a; Fuentealba et al., 2016), the dual descriptor f(2)(r) (Fuentealba and Parr, 1991; Morell et al., 2005, 2006; Ayers et al., 2007; Cardenas et al., 2009b; Geerlings et al., 2012), and the electron density ρ(r). The change in energy and density upon polarization of one reagent by another are defined through,where χ(r,r’) is the linear-response (or polarizability) kernel χ(r,r’). For convenience, Eq. 3 is written in terms of the nuclear charge density instead of the external potential, (Ayers et al., 2009)
The molecular electrostatic potential (Politzer, 1980; Politzer and Truhlar, 1981; Sjoberg and Politzer, 1990; Gadre et al., 1992; Shirsat et al., 1992; Murray et al., 1996; Suresh and Gadre, 1998; Politzer and Murray, 2002),is not traditionally considered a reactivity indicator in DFT-CR, but it can be placed in a DFT context by differentiating the total energy (including Vnn) with respect to the external potential (Ayers and Parr, 2001; Anderson et al., 2007a). The non-additive kinetic and exchange-correlation energies in the first line of Eq. 3 capture electron-pairing and steric effects, (Gordon and Kim, 1972; Wesolowski and Warshel, 1993; Wesolowski and Warshel, 1994; Liu, 2007; Wu et al., 2009)
1.1 The interaction hardness and chemical potential
To elucidate the MHP, we need to compute the change in hardness due to the interactions between A and B. To do this, we define the interaction hardnessand notice that this quantity can be computed by differentiating the interaction energy,
The treatment of the chemical potential is a bit more nuanced, because since it is an intensive property the correct expression is (Miranda-Quintana et al., 2018):with
The total energy can be used in the previous equations because the nuclear-nuclear repulsion energy does not depend on the number of electrons.
As complicated as it is, Eq. 3 already contains assumptions, most notably assumptions about the “effective external potential” that electrons in one reagent feel due to the electrons and nuclei in the second reagent (Ayers and Parr, 2001; Ayers et al., 2005; Cohen and Wasserman, 2007; Cohen et al., 2009; Liu et al., 2009; Elliott et al., 2010; Osorio et al., 2011). We also must assume that the higher-order terms in the Taylor series (which are implicitly neglected or averaged over in a Taylor-series-with-remainder strategy) are negligible. Extension to include higher-order terms can be made, with commensurate increased complexity in Eq. 3 (Senet, 1996; Geerlings and De Proft, 2008; Cardenas et al., 2009a; Heidar-Zadeh et al., 2016b). Finally, we must assume that the derivatives exist, which implicitly requires that the system is not quasi-degenerate for perturbations of the strength relevant for the analysis. Quasi-degeneracy (even exact degeneracy) can be treated, however, if the derivatives are reinterpreted as differentials (Cardenas et al., 2011; Bultinck et al., 2013a; Bultinck et al., 2013b; Pino-Rios et al., 2017; Cerón et al., 2020; Bultinck and Cárdenas, 2022; Cárdenas et al., 2022). Other effects (e.g., temperature-dependence, spin-specificity) can likewise be treated without essential difficulty, merely by an extension of definition and notation (Galvan et al., 1988; Ghanty and Ghosh, 1994; Ayers and Yang, 2006; Garza et al., 2006; Perez et al., 2008; Franco-Perez et al., 2015a; Franco-Perez et al., 2015b; Alain Miranda-Quintana and Ayers, 2016; Miranda-Quintana R. A. and Ayers P. W., 2016; Franco-Perez et al., 2017a; Franco-Perez et al., 2017b; Franco-Pérez et al., 2017; Robles et al., 2018; Gázquez et al., 2019). The following analysis can also be treated at an atom (or functional-group) condensed level: the integrations over space are merely replaced by sums over atom labels (Yang and Mortier, 1986; Fuentealba et al., 2000a; Ayers et al., 2002; Tiznado et al., 2005; Bultinck et al., 2007; Fuentealba et al., 2016; Echegaray et al., 2017). That provides a more computationally practical form for these results and draws the link to electronegativity equalization methods more strongly.
To obtain expressions for µAB and ηAB that are simple enough to be useful, some further assumptions are needed. Before proceeding, let us rewrite Eq. 3 using a shorter, more convenient notation that will help us with the upcoming manipulations:
Here hAB is the Coulomb interaction between the fragments’ Fukui functions, cAB is the Coulomb interaction between the fragments’ dual descriptors, and θA is the Coulomb interaction between the dual descriptor of fragment A and the Fukui function of fragment B. To obtain this expression, we just grouped terms according to the powers of , and neglected the terms that are N-independent (since we are dealing with derivatives with respect to N these terms won’t be relevant). For instance, , denotes the chemical potential of fragment A in the presence of B and includes not only a contribution from the chemical potential of the isolated fragment A, but also contributions from the interaction of A and B (e.g., the term ). From Eq. 14 it is easy to see that the energy will be minimized if the coefficient of the term is as big as possible (since ); while the coefficient of the term is as small as possible (ideally, a negative number). That is:
In order to evaluate Eq. 11 we need to know the interaction energy as a function of the number of electrons. The number of electrons enters the expression for the interaction energy (Eq. 3) through the reactivity indicators and through the extent of electron transfer. We assume that the reactivity indicators in Eq. 3 do not depend on the number of electrons; this is reliable if the expression in Eq. 3 has already been iterated to convergence or, failing that, that the result of the first iteration (where all the reactivity indicators are computed for the isolated reagents) suffices. Typically, we would take , and then solve for the amount of charge transfer that minimizes the interaction energy. However, this leads to expressions that are far too complicated to analyze. Minimizing Eq. 3 requires solving a cubic equation, and even when there is a clear indication of which root should be taken, the resulting expressions are of little help. Hence, we need to introduce a second approximation, assuming that the dependence of ΔNA and ΔNB on the number of electrons is linear,
Here, are just some convenient reference values used as starting points to expand the correct changes in particle numbers in A and B, respectively.
To calculate the chemical potential and hardness of the product, it is convenient to consider an excess charge on the products, namely (Miranda-Quintana R. A., 2017):
So this implies that
Now we can substitute Eq. 19 in Eq. 14, consider an infinitesimal , and truncate at second order:
Notice that we have omitted the terms that do not depend on . The linear and quadratic coefficients are the equations for the chemical potential and (twice the) hardness of the reaction product, respectively. Namely,
These are the fundamental expressions of this manuscript, as they serve as the basis for our forthcoming analyses.
1.2 The maximum hardness principle
A key point in Eq. 17 is how to estimate and . Perhaps the simplest route is to just use the standard parabolic model result, thus:
In the case of the hardness, this leads to a relatively simple expression:where, for the sake of completeness, we have reverted back to the original notation, showing all the contributions to the interaction hardness in terms of both reagents.
Since , and the reaction energy tends to decrease when increases (Eq. 15), it is straightforward to corroborate that the “Fukui function” pairing that minimizes the energy (cf. Eq. 15) also guaranties a maximum hardness value.
Analyzing the term corresponding to the electrostatic interaction of the dual descriptors, namely , is a bit more involved. Since the energy tends to decrease when decreases (Eq. 16), the interaction between the dual descriptors that minimizes the energy will maximize the hardness if
However, this will be true only when:
This might seem like an odd result, since at it (falsely) seems like it introduces an asymmetry between the reactants. However, because , one can rearrange this equation so that the symmetry of the expression with respect to permutation of the reactant labels is clear: . Notice that Eq. 26 means that the MHP will hold when the hardnesses of the reactants are not very different. This implies that there could be cases where minimizing the energy actually implies that the hardness will tend to decrease. Equation 26 is consistent with other results from the literature which indicates that in double-exchange reactions of acids and bases, the HSAB and DMB rules are be fulfilled only when the differences in hardness of the reactants are not too large (Cardenas and Ayers, 2013; Miranda-Quintana et al., 2018). This is not especially concerning as the restriction on the hardness values is rarely implicated. For example, excluding the (very hard) noble gas atoms, Eq. 26 is violated by very few atom pairs within the periodic table (Cárdenas et al., 2016), and the pairs that do violate the constraint (e.g. Cesium and Fluorine) are so extreme that there is little need for additional tools to elucidate their reactivity.
These results provide some support for the MHP, but they rely on the parabolic model (Eq. 23). We can obtain a more realistic estimate of and if we work instead with a simplified version of Eq. 14 where we neglect all cross-terms (i.e., ). Thus, we are still working with a parabolic model, but now the descriptors have some information regarding the perturbation induced by the other reagent. Therefore now we will have
Now the expression for the energy reads:
Hence
In this case we can obtain support for the MHP in the same way as we did before: the coefficient of the term is positive, therefore the “Fukui function pairing” that minimizes the energy also guaranties a maximum hardness value.
As for the dual descriptor interactions, now the MHP will be true if:
As we saw previously, this result also indicates that just maximizing the hardness might not always lead to more favorable interactions between various reagents. Only when the harder reactant is no more than ∼3.7 times higher than the hardness of the softer reagent might the formation of the hardest product will be favored.
1.3 The “|dμ| big is good” principle
Given the generally more complicated nature of the expressions involved in the treatment of the DMB principle, we will only consider the reference and presented in Eq. 27. Substituting these expressions into Eq. 20:
As it was the case in the last part of the discussion on the MHP, we will discard the cross-terms corresponding to the factors, which means that we can write:With this approximation, the energy model reduces to:
Here we have elected to explicitly include the constant (the -independent term):
Equations 32 and 34 can be rewritten as:where the (PP) index indicates that these are the expressions obtained using the Parr-Pearson parabolic model. (Parr and Pearson, 1983)
Without losing any generality, we can assume that A is the acid, namely: . Then, since, the terms in brackets in Eqs 35, 36 are always positive (cf. Eq. 16, proving the DMB is equivalent to showing that:which can be rephrased as:
Closely following the strategy employed in previous approaches to DMB (Miranda-Quintana et al., 2018), we can take (notice that equivalent expressions for the change in reactant B can be obtained by simply exchanging the indices in the following equations)resulting in:where:
It is easy to check that the sign of only depends on the sign of , so we only need to show that . For this we only need to analyze the sign of the expression:which we can rewrite in a form that makes its dependence on the sign of more apparent,
Ignoring, for the moment, the term on the second line, DMB follows if:
This inequality is very likely to hold in most cases (particularly, in the weakly-interacting regime). The last three terms are always positive. The first term is likewise positive if that . This can always be ensured by taking the initial separation of reagents to be sufficiently large (Yañez et al., 2021). For instance, for Fukui functions localized on two atomic sites separated by 5 Å, is less than 3 eV. Hence, it is safe to assume that inequality 46) holds.
It remains to analyze the term , which cannot be factored in terms of .
These results largely point out to the validity of the DMB principle. The only potential incongruity comes in the form of the hard-to-factor terms in the expressions of and . However, these terms appear because we decided to perform a more rigorous mathematical treatment of the foundations of this principle. Should we have chosen to go with more qualitative arguments (as it was the case for the MHP, and some other discussions of reactivity principles), the evidence in favor of DMB would have been even stronger. Just note that, in accordance with Eqs 15, 16, the and conditions that guarantee a minimum interaction energy, are also the ones that, following Eq. 31, will tend to maximize . However, due to the central role of DMB in chemical reactivity, it is illustrative to perform a more detailed analysis of the conditions supporting its validity. Our analysis here indicates that the DMB principle is valid where certain terms become negligible. The confounding terms become negligible when reagents are sufficiently far apart, suggesting that failures of the DMB principle to predict reactivity are most likely to occur in cases where the activated complex in a chemical reaction is tightly bound or does not exist (e.g., barrierless reactions). This is consistent with the (already well-established) reduction of the efficacy of conceptual density functional theory in such cases, due mainly to the importance of higher-order terms in the perturbative expansion.
2 Summary
In this work we have shown that starting from an expression for the interaction energy between reactants deduced in the first part of this series of papers, it is possible to find a theoretical support for the maximum hardness principle. The main difference between this work and other related papers is that here the perturbation of one reactant on another is explicitly accounted for. In summary, the MHP is fulfilled if 1) the electrostatic interaction of the Fukui functions of the reactants is positive, (Berkowitz, 1987; Ayers P. W. and Levy M., 2000; Osorio et al., 2011), 2) the electrostatic interaction of the dual descriptors of the reactants is negative and, (Ayers et al., 2007; Cardenas et al., 2009b), 3) the relation between the hardnesses of the reactants is bounded by the inequality (30).
Similarly, we provided more arguments favoring the DMB principle, which unsurprisingly also seems to hold when we take in to account the full two-reagent picture. Contrary to the MHP case, establishing the validity of the DMB principle requires more caveats, and additional mathematical scrutiny is warranted. However, even simple qualitative discussion of the form of the expression for the change in reagents’ chemical potential support its validity, and also support the favorability of large Coulomb interactions between fragments’ Fukui functions and small Coulomb interactions between fragments’ dual descriptors.
Overall, the two-reagent model discussed in this and the previous contribution provides a more complete picture of chemical reactivity, encompassing several previous approaches, while also strengthening the arguments supporting several reactivity principles. Further applications of this framework are underway and will be presented elsewhere.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Author contributions
PA, SL, CC, and RM contributed to conception and design of the study. PA and RM wrote the first draft of the manuscript. CC, RM, TG wrote sections of the manuscript. All authors contributed to manuscript revision, read, and approved the submitted version.
Acknowledgments
PWA thanks the Canada Research Chairs, NSERC, and Compute Canada for funding. FH-Z was supported by NSERC and Compute Canada. SL was supported as part of the UNC EFRC: Solar Fuels, an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences under Award No. DE-SC0001011. CC acknowledges financial support from Fondecyt through Grant Nos.1220366 and also by Centers of Excellence With Basal-Conicyt Financing, Grant FB0807. RM-Q thanks support from the University of Florida in the form of a start-up grant.
Conflict of interest
The reviewer (WT) declared a past co-authorship with the author (CC) to the handling editor.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
AndersonJ. S. M.MelinJ.AyersP. W. (2007a). Conceptual density-functional theory for general chemical reactions, including those that are neither charge- nor frontier-orbital-controlled. 1. Theory and derivation of a general-purpose reactivity indicator. J. Chem. Theory Comput.3, 358–374. 10.1021/ct600164j
2
AndersonJ. S. M.MelinJ.AyersP. W. (2007b). Conceptual density-functional theory for general chemical reactions, including those that are neither charge- nor frontier-orbital-controlled. 2. Application to molecules where Frontier molecular orbital theory fails. J. Chem. Theory Comput.3, 375–389. 10.1021/ct6001658
3
AyersP. W. (2005). An elementary derivation of the hard/soft-acid/base principle. J. Chem. Phys.122, 141102. 10.1063/1.1897374
4
AyersP. W.AndersonJ. S. M.BartolottiL. J. (2005). Perturbative perspectives on the chemical reaction prediction problem. Int. J. Quantum Chem.101, 520–534. 10.1002/qua.20307
5
AyersP. W. (2006). Axiomatic formulations of the Hohenberg-Kohn functional. Phys. Rev. A73, 012513. 10.1103/physreva.73.012513
6
AyersP. W.CárdenasC. (2013). Communication: A case where the hard/soft acid/base principle holds regardless of acid/base strength. J. Chem. Phys.138, 181106. 10.1063/1.4805083
7
AyersP. W.FiasS.Heidar-ZadehF. (2018). The axiomatic approach to chemical concepts. Comput. Theor. Chem.1142, 83–87. 10.1016/j.comptc.2018.09.006
8
AyersP. W.LevyM. (2000a). Perspective on "Density functional approach to the frontier-electron theory of chemical reactivity" - Parr RG, Yang W (1984. J. Am. Chem. Soc.106, 4049–4050.
9
AyersP. W.LevyM. (2000b). Perspective on "Density functional approach to the frontier-electron theory of chemical reactivity". Theor. Chem. Accounts103, 353–360. 10.1007/978-3-662-10421-7_59
10
AyersP. W.LiuS.LiT. (2009). Chargephilicity and chargephobicity: Two new reactivity indicators for external potential changes from density functional reactivity theory. Chem. Phys. Lett.480, 318–321. 10.1016/j.cplett.2009.08.067
11
AyersP. W.MelL. (2018). Levy constrained search in fock space: An alternative approach to noninteger electron number. 物理化学学报34, 625–630.
12
AyersP. W.MohamedM.Heidar‐ZadehF. (2022). The hard/soft acid/base rule: A perspective from conceptual density‐functional theory. Concept. Density Funct. Theory Towards a New Chem. React. Theory1, 263–279. 10.1002/9783527829941.ch14
13
AyersP. W.MorellC.De ProftF.GeerlingsP. (2007). Understanding the Woodward-Hoffmann rules by using changes in electron density. Chem. Eur. J.13, 8240–8247. 10.1002/chem.200700365
14
AyersP. W.MorrisonR. C.RoyR. K. (2002). Variational principles for describing chemical reactions: Condensed reactivity indices. J. Chem. Phys.116, 8731–8744. 10.1063/1.1467338
15
AyersP. W.ParrR. G. (2008). Local hardness equalization: Exploiting the ambiguity. J. Chem. Phys.128, 184108. 10.1063/1.2918731
16
AyersP. W.ParrR. G.PearsonR. G. (2006). Elucidating the hard/soft acid/base principle: A perspective based on half-reactions. J. Chem. Phys.124, 194107. 10.1063/1.2196882
17
AyersP. W.ParrR. G. (2000). Variational principles for describing chemical reactions: The fukui function and chemical hardness revisited. J. Am. Chem. Soc.122, 2010–2018. 10.1021/ja9924039
18
AyersP. W.ParrR. G. (2001). Variational principles for describing chemical reactions. Reactivity indices based on the external potential. J. Am. Chem. Soc.123, 2007–2017. 10.1021/ja002966g
19
AyersP. W. (2008). The dependence on and continuity of the energy and other molecular properties with respect to the number of electrons. J. Math. Chem.43, 285–303. 10.1007/s10910-006-9195-5
20
AyersP. W. (2007). The physical basis of the hard/soft acid/base principle. Faraday Discuss.135, 161–190. 10.1039/b606877d
21
AyersP. W.YangW. (2006). Legendre-transform functionals for spin-density-functional theory. J. Chem. Phys.124, 224108. 10.1063/1.2200884
22
BartlettR. J.LotrichV. F.SchweigertI. V. (2005). Ab initiodensity functional theory: The best of both worlds?J. Chem. Phys.123, 062205. 10.1063/1.1904585
23
BerkowitzM. (1987). Density functional approach to Frontier controlled reactions. J. Am. Chem. Soc.109, 4823–4825. 10.1021/ja00250a012
24
BultinckP.CárdenasC. (2022). Chemical response functions in (Quasi‐)Degenerate states. Concept. Density Funct. Theory Towards a New Chem. React. Theory1, 93–109. 10.1002/9783527829941.ch6
25
BultinckP.CardenasC.FuentealbaP.JohnsonP. A.AyersP. W. (2013a). Atomic charges and the electrostatic potential are ill-defined in degenerate ground states. J. Chem. Theory Comput.9, 4779–4788. 10.1021/ct4005454
26
BultinckP.CardenasC.FuentealbaP.JohnsonP. A.AyersP. W. (2013b). How to compute the fukui matrix and function for systems with (Quasi-)Degenerate states. J. Chem. Theory Comput.10, 202–210. 10.1021/ct400874d
27
BultinckP.FiasS.Van AlsenoyC.AyersP. W.Carbó-DorcaR. (2007). Critical thoughts on computing atom condensed Fukui functions. J. Chem. Phys.127, 034102. 10.1063/1.2749518
28
BultinckP.LangenaekerW.LahorteP.De ProftF.GeerlingsP.Van AlsenoyC.et al (2002a). The electronegativity equalization method II: Applicability of different atomic charge schemes. J. Phys. Chem. A106, 7895–7901. 10.1021/jp020547v
29
BultinckP.LangenaekerW.LahorteP.De ProftF.GeerlingsP.WaroquierM.et al (2002b). The electronegativity equalization method I: Parametrization and validation for atomic charge calculations. J. Phys. Chem. A106, 7887–7894. 10.1021/jp0205463
30
CárdenasC.AyersP. W.CedilloA. (2011). Reactivity indicators for degenerate states in the density-functional theoretic chemical reactivity theory. J. Chem. Phys.134, 174103–174113. 10.1063/1.3585610
31
CárdenasC.AyersP. W. (2013). How reliable is the hard-soft acid-base principle? An assessment from numerical simulations of electron transfer energies. Phys. Chem. Chem. Phys.15, 13959–13968. 10.1039/c3cp51134k
32
CárdenasC.EchegarayE.ChakrabortyD.AndersonJ. S. M.AyersP. W. (2009a). Relationships between the third-order reactivity indicators in chemical density-functional theory. J. Chem. Phys.130, 244105. 10.1063/1.3151599
33
CárdenasC.EcheverryA.NovoaT.Robles‐NavarroA.GomezT.FuentealbaP. (2022). The fukui function in extended systems: Theory and applications. Concept. Density Funct. Theory Towards a New Chem. React. Theory2, 555–571. 10.1002/9783527829941.ch27
34
CárdenasC.Heidar-ZadehF.AyersP. W. (2016). Benchmark values of chemical potential and chemical hardness for atoms and atomic ions (including unstable anions) from the energies of isoelectronic series. Phys. Chem. Chem. Phys.18, 25721–25734. 10.1039/c6cp04533b
35
CárdenasC.RabiN.AyersP. W.MorellC.JaramilloP.FuentealbaP. (2009b). Chemical reactivity descriptors for ambiphilic reagents: Dual descriptor, local hypersoftness, and electrostatic potential. J. Phys. Chem. A113, 8660–8667. 10.1021/jp902792n
36
CerónM. L.GomezT.CalatayudM.CárdenasC. (2020). Computing the fukui function in solid-state chemistry: Application to alkaline earth oxides bulk and surfaces. J. Phys. Chem. A124, 2826–2833. 10.1021/acs.jpca.0c00950
37
ChattarajP. K. (2007). A minimum electrophilicity perspective of the HSAB principle. Indian J. Phys. Proc. Indian Assoc. Cultiv. Sci.81, 871–879.
38
ChattarajP. K.AyersP. W.MelinJ. (2007). Further links between the maximum hardness principle and the hard/soft acid/base principle: Insights from hard/soft exchange reactions. Phys. Chem. Chem. Phys.9, 3853–3856. 10.1039/b705742c
39
ChattarajP. K.AyersP. W. (2005). The maximum hardness principle implies the hard/soft acid/base rule. J. Chem. Phys.123, 086101. 10.1063/1.2011395
40
ChattarajP. K. (2001). Chemical reactivity and selectivity: Local HSAB principle versus Frontier orbital theory. J. Phys. Chem. A105, 511–513. 10.1021/jp003786w
41
ChattarajP. K.LeeH.ParrR. G. (1991). HSAB principle. J. Am. Chem. Soc.113, 1855–1856. 10.1021/ja00005a073
42
ChattarajP. K.LiuG. H.ParrR. G. (1995). The maximum hardness principle in the Gyftopoulos-Hatsopoulos three-level model for an atomic or molecular species and its positive and negative ions. Chem. Phys. Lett.237, 171–176. 10.1016/0009-2614(95)00280-h
43
ChattarajP. K.NathS. (1994). Maximum hardness and HSAB principles: An ab initio SCF study. Indian J.Chem., Sect.A Inorg. Bio-inorg., Phys. theor.anal.chem.33A, 842–843.
44
ChattarajP. K.PérezP.ZevallosJ.Toro-LabbéA. (2001). Ab initio SCF and DFT studies on solvent effects on intramolecular rearrangement reactions. J. Phys. Chem. A105, 4272–4283. 10.1021/jp0021345
45
ChattarajP. K.SarkarU.RoyD. R. (2006). Electrophilicity index. Chem. Rev.106, 2065–2091. 10.1021/cr040109f
46
ChattarajP. K. (1996). The maximum hardness principle: An overview. Proc. Indian Natl. Sci. Acad. Part A Phys. Sci.62, 513–531.
47
ChermetteH. (1999). Chemical reactivity indexes in density functional theory. J. Comput. Chem.20, 129–154. 10.1002/(sici)1096-987x(19990115)20:1<129::aid-jcc13>3.0.co;2-a
48
CohenA. J.Mori-SánchezP.YangW. (2008a). Fractional charge perspective on the band gap in density-functional theory. Phys. Rev. B77, 115123. 10.1103/physrevb.77.115123
49
CohenA. J.Mori-SánchezP.YangW. (2008b). Insights into current limitations of density functional theory. Science321, 792–794. 10.1126/science.1158722
50
CohenM. H.WassermanA.CarR.BurkeK. (2009). Charge transfer in partition theory. J. Phys. Chem. A113, 2183–2192. 10.1021/jp807967e
51
CohenM. H.WassermanA. (2007). On the foundations of chemical reactivity theory. J. Phys. Chem. A111, 2229–2242. 10.1021/jp066449h
52
CortonaP. (1991). Self-consistently determined properties of solids without band-structure calculations. Phys. Rev. B44, 8454–8458. 10.1103/physrevb.44.8454
53
De ProftF.AyersP. W.FiasS.GeerlingsP. (2006). Woodward-Hoffmann rules in density functional theory: Initial hardness response. J. Chem. Phys.125, 214101. 10.1063/1.2387953
54
De ProftF.ChattarajP. K.AyersP. W.Torrent-SucarratM.ElangoM.SubramanianV.et al (2008). Initial hardness response and hardness profiles in the study of Woodward-Hoffmann rules for electrocyclizations. J. Chem. Theory Comput.4, 595–602. 10.1021/ct700289p
55
EchegarayE.Toro-LabbeA.DikmenliK.Heidar-ZadehF.RabiN.RabiS.et al (2017). “Negative condensed-to-atom fukui functions: A signature of oxidation-induced reduction of functional groups,” in Correlations in condensed matter under extreme conditions (Springer), 269–278. 10.1007/978-3-319-53664-4_19
56
ElliottP.BurkeK.CohenM. H.WassermanA. (2010). Partition density-functional theory. Phys. Rev. A82, 024501. 10.1103/physreva.82.024501
57
Franco-PérezM.AyersP. W.GázquezJ. L.VelaA. (2015a). Local and linear chemical reactivity response functions at finite temperature in density functional theory. J. Chem. Phys.143, 244117. 10.1063/1.4938422
58
Franco-PérezM.AyersP. W.GázquezJ. L.VelaA. (2017a). Thermodynamic responses of electronic systems. J. Chem. Phys.147, 094105. 10.1063/1.4999761
59
Franco-PérezM.GázquezJ. L.AyersP. W.VelaA. (2015b). Revisiting the definition of the electronic chemical potential, chemical hardness, and softness at finite temperatures. J. Chem. Phys.143, 154103. 10.1063/1.4932539
60
Franco-PérezM.GázquezJ. L.AyersP. W.VelaA. (2018). Thermodynamic justification for the parabolic model for reactivity indicators with respect to electron number and a rigorous definition for the electrophilicity: The essential role played by the electronic entropy. J. Chem. Theory Comput.14, 597–606. 10.1021/acs.jctc.7b00940
61
Franco-PérezM.Heidar-ZadehF.AyersP. W.GázquezJ. L.VelaA. (2017b). Going beyond the three-state ensemble model: The electronic chemical potential and fukui function for the general case. Phys. Chem. Chem. Phys.19, 11588–11602. 10.1039/c7cp00224f
62
Franco-PérezM.Polanco-RamírezC.-A.AyersP. W.GázquezJ. L.VelaA. (2017). New Fukui, dual and hyper-dual kernels as bond reactivity descriptors. Phys. Chem. Chem. Phys.19, 16095–16104. 10.1039/c7cp02613g
63
FuentealbaP.CardenasC. (2015). “Density functional theory of chemical reactivity,” in Chemical modelling: Volume 11. Editor Michael SpringborgJ.-O. J. (London: The Royal Society of Chemistry), 151–174.
64
FuentealbaP.CárdenasC. (2013). On the exponential model for energy with respect to number of electrons. J. Mol. Model.19, 2849–2853. 10.1007/s00894-012-1708-5
65
FuentealbaP.CardenasC.Pino-RiosR.TiznadoW. (2016). “Topological analysis of the fukui function,” in Applications of topological methods in molecular chemistry. Editors Remi ChauvinC. L.BernardS.AlikhaniE. (Springer International Publishing), 227–241. 10.1007/978-3-319-29022-5_8
66
FuentealbaP.ParrR. G. (1991). Higher‐order derivatives in density‐functional theory, especially the hardness derivative ∂η/∂N. J. Chem. Phys.94, 5559–5564. 10.1063/1.460491
67
FuentealbaP.PérezP.ContrerasR. (2000a). On the condensed Fukui function. J. Chem. Phys.113, 2544–2551. 10.1063/1.1305879
68
FuentealbaP.Simón-MansoY.ChattarajP. K. (2000b). Molecular electronic excitations and the minimum polarizability principle. J. Phys. Chem. A104, 3185–3187. 10.1021/jp992973v
69
GadreS. R.KulkarniS. A.ShrivastavaI. H. (1992). Molecular electrostatic potentials: A topographical study. J. Chem. Phys.96, 5253–5260. 10.1063/1.462710
70
GalvanM.VelaA.GazquezJ. L. (1988). Chemical reactivity in spin-polarized density functional theory. J. Phys. Chem.92, 6470–6474. 10.1021/j100333a056
71
GarzaJ.VargasR.CedilloA.GalvánM.ChattarajP. K. (2006). Comparison between the frozen core and finite differences approximations for the generalized spin-dependent global and local reactivity descriptors in small molecules. Theor. Chem. Acc.115, 257–265. 10.1007/s00214-005-0002-3
72
GázquezJ. L.Franco‐PérezM.AyersP. W.VelaA. (2019). Temperature‐dependent approach to chemical reactivity concepts in density functional theory. Int. J. Quantum Chem.119, e25797. 10.1002/qua.25797
73
GazquezJ. L.MendezF. (1994). The hard and soft acids and bases principle: An atoms in molecules viewpoint. J. Phys. Chem.98, 4591–4593. 10.1021/j100068a018
74
GazquezJ. L. (2008). Perspectives on the density functional theory of chemical reactivity. J. Mexican Chem. Soc.52, 3–10.
75
GeerlingsP.AyersP. W.Toro-LabbéA.ChattarajP. K.De ProftF. (2012). The woodward-hoffmann rules reinterpreted by conceptual density functional theory. Acc. Chem. Res.45, 683–695. 10.1021/ar200192t
76
GeerlingsP.ChamorroE.ChattarajP. K.De ProftF.GázquezJ. L.LiuS.et al (2020). Conceptual density functional theory: Status, prospects, issues. Theor. Chem. Accounts139, 1–18. 10.1007/s00214-020-2546-7
77
GeerlingsP.De ProftF. (2008). Conceptual DFT: The chemical relevance of higher response functions. Phys. Chem. Chem. Phys.10, 3028–3042. 10.1039/b717671f
78
GeerlingsP.De ProftF.LangenaekerW. (2003). Conceptual density functional theory. Chem. Rev.103, 1793–1874. 10.1021/cr990029p
79
GhantyT. K.GhoshS. K. (1994). Spin-Polarized generalization of the concepts of electronegativity and hardness and the description of chemical binding. J. Am. Chem. Soc.116, 3943–3948. 10.1021/ja00088a033
80
GordonR. G.KimY. S. (1972). Theory for the forces between closed‐shell atoms and molecules. J. Chem. Phys.56, 3122–3133. 10.1063/1.1677649
81
GovindN.WangY. A.CarterE. A. (1999). Electronic-structure calculations by first-principles density-based embedding of explicitly correlated systems. J. Chem. Phys.110, 7677–7688. 10.1063/1.478679
82
Heidar-ZadehF.Miranda-QuintanaR. A.VerstraelenT.BultinckP.AyersP. W. (2016a). When is the fukui function not normalized? The danger of inconsistent energy interpolation models in density functional theory. J. Chem. Theory Comput.12, 5777–5787. 10.1021/acs.jctc.6b00494
83
Heidar-ZadehF.RicherM.FiasS.Miranda-QuintanaR. A.ChanM.Franco-PérezM.et al (2016b). An explicit approach to conceptual density functional theory descriptors of arbitrary order. Chem. Phys. Lett660, 307–312. 10.1016/j.cplett.2016.07.039
84
JanakJ. F. (1978). Proof that∂E∂ni=εin density-functional theory. Phys. Rev. B18, 7165–7168. 10.1103/physrevb.18.7165
85
JohnsonP. A.BartolottiL. J.AyersP. W.FievezT.GeerlingsP. (2012). “Charge density and chemical reactivity: A unified view from conceptual DFT,” in Modern charge density analysis. Editors GattiC.MacchiP. (New York: Springer), 715–764.
86
LangenaekerW.De ProftF.GeerlingsP. (1995). Development of local hardness-related reactivity indices: Their application in a study of the SE at monosubstituted benzenes within the HSAB context. J. Phys. Chem.99, 6424–6431. 10.1021/j100017a022
87
LevyM. (1979). Universal variational functionals of electron densities, first-order density matrices, and natural spin-orbitals and solution of the v -representability problem. Proc. Natl. Acad. Sci. U.S.A.76, 6062–6065. 10.1073/pnas.76.12.6062
88
LiebE. H. (1983). Density functionals for Coulomb systems. Int. J. Quantum Chem.24, 243–277. 10.1002/qua.560240302
89
LiuS. B. (2009a). Conceptual density functional theory and some recent developments. Acta Physico-Chimica Sin.25, 590–600.
90
LiuS. B. (2009b). “Electrophilicity,” in Chemical reactivity theory: A density functional view. Editor ChattarajP. K. (Boca Raton: Taylor & Francis), 179. 10.1201/9781420065442.ch13
91
LiuS.LiT.AyersP. W. (2009). Potentialphilicity and potentialphobicity: Reactivity indicators for external potential changes from density functional reactivity theory. J. Chem. Phys.131, 114106. 10.1063/1.3231687
92
LiuS. (2007). Steric effect: A quantitative description from density functional theory. J. Chem. Phys.126, 244103. 10.1063/1.2747247
93
MelinJ.AparicioF.SubramanianV.GalvánM.ChattarajP. K. (2004). Is the fukui function a right descriptor of Hard−Hard interactions?J. Phys. Chem. A108, 2487–2491. 10.1021/jp037674r
94
MendezF.GazquezJ. L. (1994). Chemical reactivity of enolate ions: The local hard and soft acids and bases principle viewpoint. J. Am. Chem. Soc.116, 9298–9301. 10.1021/ja00099a055
95
MerkleR.SavinA.PreussH. (1992). Singly ionized first‐row dimers and hydrides calculated with the fully‐numerical density‐functional program numol. J. Chem. Phys.97, 9216–9221. 10.1063/1.463297
96
Miranda-QuintanaR. A.ChattarajP. K.AyersP. W. (2017a). Finite temperature grand canonical ensemble study of the minimum electrophilicity principle. J. Chem. Phys.147, 124103. 10.1063/1.4996443
97
Miranda-QuintanaR. A.AyersP. W. (2018a). Dipolar cycloadditions and the “| Δμ| big is good” rule: A computational study. Theor. Chem. Accounts137, 1–7. 10.1007/s00214-018-2391-0
98
Miranda-QuintanaR. A.AyersP. W. (2016a). Fractional electron number, temperature, and perturbations in chemical reactions. Phys. Chem. Chem. Phys.18, 15070–15080. 10.1039/c6cp00939e
99
Miranda-QuintanaR. A.AyersP. W. (2016). Fractional electron number, temperature, and perturbations in chemical reactions. Phys. Chem. Chem. Phys.18, 15070–15080. 10.1039/c6cp00939e
100
Miranda-QuintanaR. A.AyersP. W. (2016b). Interpolation of property-values between electron numbers is inconsistent with ensemble averaging. J. Chem. Phys.144, 244112. 10.1063/1.4953557
101
Miranda-QuintanaR. A.AyersP. W. (2018b). Note: Maximum hardness and minimum electrophilicity principles. J. Chem. Phys.148, 196101. 10.1063/1.5033964
102
Miranda-QuintanaR. A.AyersP. W. (2016c). Systematic treatment of spin-reactivity indicators in conceptual density functional theory. Theor. Chem. Accounts135, 1–18. 10.1007/s00214-016-1995-5
103
Miranda-QuintanaR. A.AyersP. W. (2019). The “|Δμ| big is good” rule, the maximum hardness, and minimum electrophilicity principles. Theor. Chem. Accounts138. 10.1007/s00214-019-2435-0
104
Miranda-QuintanaR. A.BochicchioR. C. (2014). Energy dependence with the number of particles: Density and reduced density matrices functionals. Chem. Phys. Lett.593, 35–39. 10.1016/j.cplett.2013.12.071
105
Miranda-QuintanaR. A. (2018). Density functional theory for chemical reactivity. Toronto: Apple Academic Press.
106
Miranda-QuintanaR. A.DeswalN.RoyR. K. (2022a). Hammett constants from density functional calculations: Charge transfer and perturbations. Theor. Chem. Accounts141, 1–10. 10.1007/s00214-021-02863-5
107
Miranda-QuintanaR. A.Heidar-ZadehF.AyersP. W. (2018). Elementary derivation of the "|Δμ| big is good" rule. J. Phys. Chem. Lett.9, 4344–4348. 10.1021/acs.jpclett.8b01312
108
Miranda-QuintanaR. A.Heidar-ZadehF.FiasS.ChapmanA. E. A.LiuS.MorellC.et al (2022b). Molecular interactions from the density functional theory for chemical reactivity: The interaction energy between two-reagents. Front. Chem.10. 10.3389/fchem.2022.906674
109
Miranda-QuintanaR. A.KimT. D.CárdenasC.AyersP. W. (2017b). The HSAB principle from a finite-temperature grand-canonical perspective. Theor. Chem. Acc.136, 135. 10.1007/s00214-017-2167-y
110
Miranda-QuintanaR. A.Martínez GonzálezM.AyersP. W. (2016). Electronegativity and redox reactions. Phys. Chem. Chem. Phys.18, 22235–22243. 10.1039/c6cp03213c
111
Miranda-QuintanaR. A. (2017a). Note: The minimum electrophilicity and the hard/soft acid/base principles. J. Chem. Phys.146, 046101. 10.1063/1.4974987
112
Miranda-QuintanaR. A. (2017b). Perturbed reactivity descriptors: The chemical hardness. Theor. Chem. Accounts136, 1–8. 10.1007/s00214-017-2109-8
113
Miranda‐QuintanaR. A.AyersP. W.Heidar‐ZadehF. (2021). Reactivity and charge transfer beyond the parabolic model: The “| Δμ| big is good” principle. ChemistrySelect6, 96–100.
114
MorellC.GrandA.Toro-LabbéA. (2005). New dual descriptor for chemical reactivity. J. Phys. Chem. A109, 205–212. 10.1021/jp046577a
115
MorellC.GrandA.Toro-LabbéA. (2006). Theoretical support for using the Δf(r) descriptor. Chem. Phys. Lett.425, 342–346. 10.1016/j.cplett.2006.05.003
116
MorellC.LabetV.GrandA.ChermetteH. (2009). Minimum electrophilicity principle: An analysis based upon the variation of both chemical potential and absolute hardness. Phys. Chem. Chem. Phys.11, 3417–3423. 10.1039/b818534d
117
Mori-SánchezP.CohenA. J.YangW. (2008). Localization and delocalization errors in density functional theory and implications for band-gap prediction. Phys. Rev. Lett.100, 146401. 10.1103/physrevlett.100.146401
118
Mori-SánchezP.CohenA. J.YangW. (2006). Many-electron self-interaction error in approximate density functionals. J. Chem. Phys.125, 201102. 10.1063/1.2403848
119
MortierW. J. (1987). Electronegativity equalization and its applications. Struct. Bond.66, 125–143.
120
MortierW. J.GhoshS. K.ShankarS. (1986). Electronegativity-equalization method for the calculation of atomic charges in molecules. J. Am. Chem. Soc.108, 4315–4320. 10.1021/ja00275a013
121
MortierW. J.Van GenechtenK.GasteigerJ. (1985). Electronegativity equalization: Application and parametrization. J. Am. Chem. Soc.107, 829–835. 10.1021/ja00290a017
122
MurrayJ. S.BrinckT.PolitzerP. (1996). Relationships of molecular surface electrostatic potentials to some macroscopic properties. Chem. Phys.204, 289–299. 10.1016/0301-0104(95)00297-9
123
NalewajskiR. F. (1984). Electrostatic effects in interactions between hard (soft) acids and bases. J. Am. Chem. Soc.106, 944–945. 10.1021/ja00316a020
124
NalewajskiR. F.KorchowiecJ.ZhouZ. (1988). Molecular hardness and softness parameters and their use in chemistry. Int. J. Quantum Chem.34, 349–366. 10.1002/qua.560340840
125
NoorizadehS. (2007). Is there a minimum electrophilicity principle in chemical reactions?Chin. J. Chem.25, 1439–1444. 10.1002/cjoc.200790266
126
OsorioE.FerraroM. B.OñaO. B.CardenasC.FuentealbaP.TiznadoW. (2011). Assembling small silicon clusters using criteria of maximum matching of the Fukui functions. J. Chem. Theory Comput.7, 3995–4001. 10.1021/ct200643z
127
PanS.SolàM.ChattarajP. K. (2013). On the validity of the maximum hardness principle and the minimum electrophilicity principle during chemical reactions. J. Phys. Chem. A117, 1843–1852. 10.1021/jp312750n
128
ParrR. G.BartolottiL. J. (1982). On the geometric mean principle for electronegativity equalization. J. Am. Chem. Soc.104, 3801–3803. 10.1021/ja00378a004
129
ParrR. G.ChattarajP. K. (1991). Principle of maximum hardness. J. Am. Chem. Soc.113, 1854–1855. 10.1021/ja00005a072
130
ParrR. G. (1994). Companions in the search. Int. J. Quantum Chem.49, 739–770. 10.1002/qua.560490515
131
ParrR. G.DonnellyR. A.LevyM.PalkeW. E. (1978). Electronegativity: The density functional viewpoint. J. Chem. Phys.68, 3801–3807. 10.1063/1.436185
132
ParrR. G.PearsonR. G. (1983). Absolute hardness: Companion parameter to absolute electronegativity. J. Am. Chem. Soc.105, 7512–7516. 10.1021/ja00364a005
133
ParrR. G.YangW. (1984). Density functional approach to the frontier-electron theory of chemical reactivity. J. Am. Chem. Soc.106, 4049–4050. 10.1021/ja00326a036
134
ParrR. G.YangW. (1989). Density-functional theory of atoms and molecules. New York: Oxford UP.
135
PearsonR. G. (1968a). Hard and soft acids and bases, HSAB, part 1: Fundamental principles. J. Chem. Educ.45, 581–587. 10.1021/ed045p581
136
PearsonR. G. (1968b). Hard and soft acids and bases, HSAB, part II: Underlying theories. J. Chem. Educ.45, 643–648. 10.1021/ed045p643
137
PearsonR. G. (1999). Maximum chemical and physical hardness. J. Chem. Educ.76, 267–275. 10.1021/ed076p267
138
PearsonR. G.PalkeW. E. (1992). Support for a principle of maximum hardness. J. Phys. Chem.96, 3283–3285. 10.1021/j100187a020
139
PearsonR. G. (1994). Principle of maximum physical hardness. J. Phys. Chem.98, 1989–1992. 10.1021/j100058a044
140
PearsonR. G. (1987). Recent advances in the concept of hard and soft acids and bases. J. Chem. Educ.64, 561–567. 10.1021/ed064p561
141
PearsonR. G. (1993). The principle of maximum hardness. Acc. Chem. Res.26, 250–255. 10.1021/ar00029a004
142
PerdewJ. P.ParrR. G.LevyM.BalduzJ. L.Jr. (1982). Density-functional theory for fractional particle number: Derivative discontinuities of the energy. Phys. Rev. Lett.49, 1691–1694. 10.1103/physrevlett.49.1691
143
PérezP.ChamorroE.AyersP. W. (2008). Universal mathematical identities in density functional theory: Results from three different spin-resolved representations. J. Chem. Phys.128, 204108. 10.1063/1.2916714
144
Pino‐RiosR.YañezO.InostrozaD.RuizL.CardenasC.FuentealbaP.et al (2017). Proposal of a simple and effective local reactivity descriptor through a topological analysis of an orbital‐weighted fukui function. J. Comput. Chem.38, 481–488.
145
PolitzerP. (1980). Electrostatic potential-electronic density relationships in atoms. II. J. Chem. Phys.73, 3264–3267. 10.1063/1.440521
146
PolitzerP.MurrayJ. S. (2002). The fundamental nature and role of the electrostatic potential in atoms and molecules. Theor. Chem. Accounts Theory, Comput. Model. Theor. Chimica Acta)108, 134–142. 10.1007/s00214-002-0363-9
147
PolitzerP.TruhlarD. (1981). Chemical applications of atomic and molecular electrostatic potentials. New York: Plenum.
148
ProftF. D.AyersP. W.GeerlingsP. (2014). “The conceptual density functional theory perspective of bonding,” in The chemical bond: Fundamental aspects of chemical bonding. Editors ShaikS.FrenkingG. (Darmstadt: Wiley), 233–270. 10.1002/9783527664696.ch7
149
RappeA. K.GoddardW. A. (1991). Charge equilibration for molecular dynamics simulations. J. Phys. Chem.95, 3358–3363. 10.1021/j100161a070
150
RoblesA.Franco-PérezM.GázquezJ. L.CárdenasC.FuentealbaP. (2018). Local electrophilicity. J. Mol. Model.24, 245. 10.1007/s00894-018-3785-6
151
RuzsinszkyA.PerdewJ. P.CsonkaG. I.VydrovO. A.ScuseriaG. E. (2007). Density Functionals that are one- and two- are not always many-electron self-interaction-free, as shown for H2+, He2+, LiH+, and Ne2+. J. Chem. Phys.126, 104102. 10.1063/1.2566637
152
SandersonR. T. (1951). An interpretation of bond lengths and a classification of bonds. Science114, 670–672. 10.1126/science.114.2973.670
153
SenetP. (1996). Nonlinear electronic responses, Fukui functions and hardnesses as functionals of the ground‐state electronic density. J. Chem. Phys.105, 6471–6489. 10.1063/1.472498
154
ShirsatR. N.BapatS. V.GadreS. R. (1992). Molecular electrostatics. A comprehensive topographical approach. Chem. Phys. Lett.200, 373–378. 10.1016/0009-2614(92)87006-b
155
SjobergP.PolitzerP. (1990). Use of the electrostatic potential at the molecular surface to interpret and predict nucleophilic processes. J. Phys. Chem.94, 3959–3961. 10.1021/j100373a017
156
SureshC. H.GadreS. R. (1998). A novel electrostatic approach to substituent constants: Doubly substituted benzenes. J. Am. Chem. Soc.120, 7049–7055. 10.1021/ja973105j
157
TiznadoW.ChamorroE.ContrerasR.FuentealbaP. (2005). Comparison among four different ways to condense the fukui function. J. Phys. Chem. A109, 3220–3224. 10.1021/jp0450787
158
Torrent-SucarratM.LuisJ. M.DuranM.SolàM. (2002). Are the maximum hardness and minimum polarizability principles always obeyed in nontotally symmetric vibrations?J. Chem. Phys.117, 10561–10570. 10.1063/1.1517990
159
Torrent-SucarratM.LuisJ. M.DuranM.SolàM. (2001). On the validity of the maximum hardness and minimum polarizability principles for nontotally symmetric vibrations. J. Am. Chem. Soc.123, 7951–7952. 10.1021/ja015737i
160
VaidehiN.WesolowskiT. A.WarshelA. (1992). Quantum‐mechanical calculations of solvation free energies. A combinedabinitiopseudopotential free‐energy perturbation approach. J. Chem. Phys.97, 4264–4271. 10.1063/1.463928
161
VerstraelenT.AyersP. W.Van SpeybroeckV.WaroquierM. (2013). ACKS2: Atom-condensed Kohn-Sham DFT approximated to second order. J. Chem. Phys.138, 074108. 10.1063/1.4791569
162
WesolowskiT. A.LeszczynskiJ. (2006). “One-electron equations for embedded electron density: Challenge for theory and practical payoffs in multi-scale modelling of complex polyatomic molecules,” in Computational chemistry: Reviews of current trends (Singapore: World Scientific), 1–82.
163
WesolowskiT. A. (2004). Quantum chemistry 'without orbitals' - an old idea and recent developments. Chimia58, 311–315. 10.2533/000942904777677885
164
WesolowskiT. A.WarshelA. (1993). Frozen density functional approach for ab initio calculations of solvated molecules. J. Phys. Chem.97, 8050–8053. 10.1021/j100132a040
165
WesolowskiT.WarshelA. (1994). Ab initio free energy perturbation calculations of solvation free energy using the frozen density functional approach. J. Phys. Chem.98, 5183–5187. 10.1021/j100071a003
166
WuQ.AyersP. W.ZhangY. (2009). Density-based energy decomposition analysis for intermolecular interactions with variationally determined intermediate state energies. J. Chem. Phys.131, 164112. 10.1063/1.3253797
167
YañezO.Báez-GrezR.InostrozaD.Pino-RiosR.Rabanal-LeónW. A.Contreras-GarcíaJ.et al (2021). Kick-Fukui: A fukui function-guided method for molecular structure prediction. J. Chem. Inf. Model. 61, 3955–3963. 10.1021/acs.jcim.1c00605
168
YangW.MortierW. J. (1986). The use of global and local molecular parameters for the analysis of the gas-phase basicity of amines. J. Am. Chem. Soc.108, 5708–5711. 10.1021/ja00279a008
169
YangW.ParrR. G.PucciR. (1984). Electron density, Kohn-Sham Frontier orbitals, and Fukui functions. J. Chem. Phys.81, 2862–2863. 10.1063/1.447964
170
YangW.ZhangY.AyersP. W. (2000). Degenerate ground states and a fractional number of electrons in density and reduced density matrix functional theory. Phys. Rev. Lett.84, 5172–5175. 10.1103/physrevlett.84.5172
171
YangX. D.PatelA. H. G.Miranda-QuintanaR. A.Heidar-ZadehF.González-EspinozaC. E.AyersP. W. (2016). Communication: Two types of flat-planes conditions in density functional theory. J. Chem. Phys.145, 031102. 10.1063/1.4958636
172
ZhouZ.ParrR. G. (1989). New measures of aromaticity: Absolute hardness and relative hardness. J. Am. Chem. Soc.111, 7371–7379. 10.1021/ja00201a014
Summary
Keywords
DFT‐density functional theory, chemical reactivity, HSAB (hard-soft-acid-base) concept, chemical potential, variational principle
Citation
Miranda-Quintana RA, Heidar-Zadeh F, Fias S, Chapman AEA, Liu S, Morell C, Gómez T, Cárdenas C and Ayers PW (2022) Molecular interactions from the density functional theory for chemical reactivity: Interaction chemical potential, hardness, and reactivity principles. Front. Chem. 10:929464. doi: 10.3389/fchem.2022.929464
Received
26 April 2022
Accepted
27 June 2022
Published
22 July 2022
Volume
10 - 2022
Edited by
Juan Frau, University of the Balearic Islands, Spain
Reviewed by
William Tiznado, Andres Bello University, Chile
Daniel Glossman-Mitnik, Centro de Investigación de Materiales Avanzados (CIMAV), Mexico
Updates
Copyright
© 2022 Miranda-Quintana, Heidar-Zadeh, Fias, Chapman, Liu, Morell, Gómez, Cárdenas and Ayers.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ramón Alain Miranda-Quintana, quintana@chem.ufl.edu; Tatiana Gómez, tatiana.mgc@gmail.com Carlos Cárdenas, cardena@uchile.cl; Paul W. Ayers, payers@mcmaster.ca
This article was submitted to Theoretical and Computational Chemistry, a section of the journal Frontiers in Chemistry
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.