Review Article  Open Access
Philicity and Fugality Scales for Organic Reactions
Abstract
Theoretical scales of reactivity and selectivity are important tools to explain and to predict reactivity patterns, including reaction mechanisms. The main achievement of these efforts has been the incorporation of such concepts in advanced texts of organic chemistry. In this way, the modern organic chemistry language has become more quantitative, making the classification of organic reactions an easier task. The reactivity scales are also useful to set up a number of empirical rules that help in rationalizing and in some cases anticipating the possible reaction mechanisms that can be operative in a given organic reaction. In this review, we intend to give a brief but complete account on this matter, introducing the conceptual basis that leads to the definition of reactivity indices amenable to build up quantitative models of reactivity in organic reactions. The emphasis is put on two basic concepts describing electronrich and electrondeficient systems, namely, nucleophile and electrophiles. We then show that the regional nucleophilicity and electrophilicity become the natural descriptors of electrofugality and nucleofugality, respectively. In this way, we obtain a closed body of concepts that suffices to describe electron releasing and electron accepting molecules together with the description of permanent and leaving groups in addition, nucleophilic substitution and elimination reactions.
1. Introduction
The development of reactivity indices to describe organic reactivity has been an active area of research from the dawn of theoretical physical organic chemistry [1–4]. From the earlier semiempirical models proposed in Hückel molecular orbital (HMO) theory, reactivity was described with the aid of static first order reactivity descriptors like atomic charges, free valence index, and bond orders [5, 6]. These reactivity indices were formerly developed around the ground state of reactants. The next generation of reactivity indices began with the elegant theory proposed by Coulson and LonguetHiggins, in a series of papers describing the response functions, including second order quantities like atomic and bond polarizabilities [7–9]. It is worth emphasizing that, at that time, HMO and CoulsonLonguetHiggins theories were conceived not as methods to approximately solve the Schröedinger equation but as models of chemical bond. The third generation of reactivity indices started after the pioneering work of Gilles Klopman, who introduced the concept of charge and frontier controlled reactions, including solvation effects [4]. Nowadays, the treatment of chemical reactivity is mostly performed using the frame of the conceptual density functional theory developed by Parr et al. [10–14], Pearson [15–18] and Pearson and Songstad [19], and other authors [20–22]. This reactivity model converted classical chemical concepts like electronegativity, hardness, and softness into numbers. In this way, atoms, molecules, and charged system were classified into quantitative scales of reactivity. This historical description does not illustrate different and independent approaches to describing chemical reactivity on quantitative basis. As a matter of fact, those response functions defined within the CoulsonLonguetHiggins theory can easily be cast into the form of response function of DFT [23, 24].
This brief review has been organized as follows: in the first section we introduce the basic definitions leading to the quantitative description of electronegativity, hardness, and softness. The key piece to achieve these definitions lies in the concept of electronic chemical potential and its derivatives introduced by Parr and Yang [13, 27] and Pearson and Songstad [19]. The electronic chemical potential is defined therein as the first derivative of the energy with respect to the number of electrons [27]:
In (1), the derivative is taken at constant external potential (r) (i.e., the potential due to the compensating nuclear charges in the system). and are the vertical ionization potential and electron affinity, respectively. The electronic chemical potential as given by (1) becomes the negative of Mulliken electronegativity, and therefore it becomes a natural description of the direction of the electronic flux during a chemical interaction [27]. For this reason, Parr et al. proposed the electronic chemical potential as a quantity measuring the tendency of electrons to escape from the system. This result is relevant for it gives a first appraisal of the global electrondonating/electron accepting pattern of an interacting pair of atoms or molecules. If the chemical potential of a species A, say, is greater than its partner B, then the electronic flux will take place from A to B, thereby suggesting that, during the interaction, A will act as nucleophile and B as electrophile. We will return to this point afterwards to introduce more refined models of electrophilicity and nucleophilicity.
A useful computational definition of the electronic chemical potential may be obtained using Koppmans’s theorem [28] that leaves (1) in terms of the oneelectron energy levels of the frontier molecular orbitals HOMO and LUMO: The electronic chemical potential of stable species is a negative semidefinite quantity. This remark will be of importance later.
Another pertinent quantity is the chemical hardness, defined as the first derivative of the electronic chemical potential with respect to (or the second derivative of the energy with respect to the number of electrons) [27]. Its finite difference version is as follows: Note that chemical hardness is a positive definite quantity. Its associated inverse is the chemical softness .
The functional derivative of the electronic chemical potential with respect to the external potential (r) at constant defines a local response function called the Fukui function (r) [27]: The Fukui function has a more workable operative form obtained after using a Maxwell relationship [27]:
The Fukui function written as in (5) defines it as the change in the electron density at point r in space (r) after the system accepts or releases one electron unit. The Fukui function is a reactivity index itself (a normalized softness) [27]. However, the most relevant role of the Fukui function is to act as a distribution function that may be used to project any global quantity [35]. The most useful form is however its condensed to atom version, , obtained after a regional integration around the atomic center in a molecule [36–38].
The condensed to atom Fukui function may be approached from a threepoint interpolation finite difference formula, or using a singlepoint calculation using a Mullikenlike population analysis. Other approaches have also been proposed [39]. The finite difference approximation leads to the definition of electrophilic, nucleophilic, and radical Fukui functions that will play a key role in what follows [36, 37].
2. Philicity and Fugality Indices
Based on a proposal by Maynard et al. [40], Parr et al. derived global electrophilicity, defined as the stabilization in energy that an electron acceptor atom or molecule undergoes, when it is embedded in an electron bath at constant electronic chemical potential [12]. The global electrophilicity index was given the following working expression:
Equation (6) shows that the best electrophile will be the species displaying a high value of the electronic chemical potential and a low value of chemical hardness or high electronegativity and high softness (or high polarizability). We have implemented Parr’s global electrophilicity index in the form of reactivity scales for a series of classical reactions in organic chemistry [41] that include cycloaddition reactions [42], elimination reactions [43], nucleophilic reactions (both aromatic and aliphatic) [44, 45], addition reactions [46], epoxidation reactions [47], redox, and biochemical processes [46, 48], including the chemistry of carbenes [49].
Nucleophilicity on the other hand cannot be derived within the same model leading to the definition of the electrophilicity index [30, 50, 51]. This problem arises because, for the righthand side of the parabola model used by Parr et al., the electronic chemical potential becomes positive semidefinite. This drawback of the philicity model is important for it is related to an empirical rule stating that big electrophilicity/nucleophilicity differences can be related to stepwise reaction mechanism with a high polar transition state, whereas small differences are related to nonpolar concerted mechanisms [52]. However, it is still possible to propose a nucleophilicity index based on the same energy expansion leading to the electrophilicity index. This index has been derived for the limit case where the charge released by the nucleophile is exactly equal to one electron unit. In this case the nucleophilicity index can be simply represented as the negative of the ionization potential; namely [30], This definition can intuitively be justified for it suggests that the best nucleophile will be the species that destabilizes to a lowest extent in the process of releasing one electron unit of charge. In the following sections, we shall illustrate the reliability and usefulness of the electrophilicity and nucleophilicity indices for a series of classic organic reactions.
The definition of fugality indices requires the introduction of local reactivity indices, where the electrophilic (+) and nucleophilic (−) Fukui functions play a key role [25]. For instance, using (6) together with the additive property of global softness, one of us introduced the concept of semilocal or regional electrophilicity condensed to atom as follows [43, 53]: where is the condensed to atom electrophilic Fukui function. In a similar way, the condensed to atom nucleophilicity index can be expressed as [25] where is the condensed to atom nucleophilic Fukui function.
The relationship with fugality concepts is framed on the regional philicity concepts sketched in Scheme 1 [25].
First of all, fugality quantities are group properties of a molecule [25]. This means that the propensity of a fragment to detach during a heterolytic bond breaking process may be safely described using a regional property of that fragment embedded in the chemical ambient of the remaining moiety of the molecule. It is important to stress this point that fugality quantities are not intrinsic properties of the isolated fragment [25, 43, 53]. With this model in mind, we can readily define nucleofugality, as the group electrophilicity evaluated on the whole molecule, where the highest values of group electrophilicity are expected to be mostly concentrated at the leaving group (LG) moiety. This is a reasonable representation of nucleofugality number, as the nucleofuge is the group that departs from the molecule bearing the bond electron pair during the heterolytic bond cleavage [53]. At the same time, the electrofuge (or permanent group R) is expected to act as an electron releasing fragment, and, therefore, electrofugality may reasonably be described by the group nucleophilicity [25]. The working formulae to quantify nucleofugality and electrofugality numbers are given by respectively.
3. Applications
The phenomenological reactivity theory described in the previous sections may be applied at the ground state of atoms and molecules as well as at the transition state stage of reactions. Note that, in doing so, both global and semilocal (regional) quantities may be used to build up “activation” properties. In what follows we present some of the applications that we and other authors have used to illustrate the reliability and usefulness of this theoretical model of reactivity that embodies reactivity, (regio)selectivity, and site activation [45, 54, 55].
3.1. The Electrophilicity Index
The electrophilicity index developed by Parr et al. (6) has been widely used to explain many organic reactions. The main purpose of these studies was traced to relate them to experimental data that include hydride affinity [26], rate coefficients [44, 56, 57], toxicity indexes [58–60], and many other applications [61–64]. For instance, Campodónico et al. demonstrated that the electrophilicity of quinones may be concisely used to deduce a hydride affinity (HA) scale in the gas phase [26]. The opportunity of having a global electrophilicity hierarchy related to HA is useful because experimental HA data are scarce. Because the measurement of HA is not direct, the electrophilicityHA relationships provide a simple way to establish a sound HA scale. The gas phase hydride affinity may be obtained as the negative enthalpy change for the reaction , where Q is an oxidized molecule and QH^{−} is its reduced form. The hydride affinity has been considered as a descriptor for Lewis’s acidity [65] since it may be a characteristic of the electron accepting capability of the electrophile [29]. Within this model, the ability of quinones to bind an H^{−} ion will be related to their electrophilic response and may be modulated by the presence of a series of substituent that can stabilize an extra negative charge, as illustrated in Scheme 2 [26].
The comparison between the experimentally observed hydride affinity and the electrophilicity index is shown in Figure 1. The regression is the result of the comparison of six experimentally obtained hydride affinities.
The empirical regression equation is and is a suitable way to obtain the hydride affinity of quinones not established up to date. First of all, from the regression equation (12), it is possible to note that the hydride affinities are well correlated with the index. The relationship with the substituent effect is also described by in the sense that electrondonating groups enhance the hydride affinity of quinones, while electroattracting groups diminish the hydride affinity values. When the substituent is the CN or Cl groups (moderate electronwithdrawing groups), results in an electrophilic activation ( = 2.09 eV and = 1.56 eV, resp.) with respect to the reference nonsubstituted quinone ( = 1.24 eV). On the other hand, the presence of the marginal donating methyl group results in an electrophilic deactivation ( = 0.91 eV) [26]. The regression between the electrophilicity index and the predicted hydride affinities [66] is presented in Figure 2.
The comparison in Figure 2 shows two families that cannot be accommodated in a single correlation line. The regression corresponds to the orthoquinone derivatives (dashed line) and the paraquinone derivatives (solid line). These results represent another useful application of the index. The hydride attachment strongly depends on the relative position of the carbonyl groups on the quinine: in the ortholike derivatives, the electrophilicity index well describes the alphalike effect promoted in this system by the presence of an adjacent electronrich atom [67, 68]. These systems are predicted as nucleophiles or marginal electrophiles [69, 70].
3.2. Nucleophilicity
The usefulness of the nucleophilicity index (7) is discussed here for a series of neutral and charged electron donors [30]. The series of neutral nucleophiles relative to water include N_{2}, CO, PH_{3}, H_{2}CO, H_{2}S, furan, (CH_{3})_{3}P, H_{2}O, NH_{3,} and CH_{3}OCH_{3} [71, 72].
This series was experimentally studied by Legon and Millen who proposed a spectroscopic scale of electrophilicity and nucleophilicity based on the intermolecular stretching force constant for the interaction of these neutral nucleophiles [71, 72] towards the hydrogen fluoride HF as probe. In the gas phase, the facile formation of a hydrogen bonded nucleophileHF complex permits the intermolecular stretching of the nucleophile/HF moieties to be evaluated. The nucleophilicity scale is simply given by the ordered hierarchy of the respective force constants. The comparison between the experimental nucleophilicity numbers reported by Legon and Millen with those obtained from our model equation (7) for the whole series of neutral nucleophiles was fair, because ammonia and dimethyl ether strongly deviated from linearity. A deep sight to the available IR data for the complexes H_{3}NHF and (CH_{3})_{2}OHF revealed that the normal mode assignments for these species were not as clean as desired. The calculated intermolecular stretching around the hydrogen bond was hardly contaminated by torsion and bond deformation (wagging) of the sp^{3} groups attached to the heteroatoms N and O. However, when the regression was made with the 7 remaining nucleophiles, a reasonably good correlation was obtained. The comparison is shown in Figure 3.
It is worth mentioning that the model equation (7) also works reasonably well for ionization potential in solution phase [30].
As we will show later, the semilocal or regional nucleophilicity patterns of molecules become of great relevance to describe another property, namely, electrofugality, a group property describing permanent group abilities in heterolytic bond cleavage processes. The model equation describing regional or group nucleophilicity is that quoted as (9). Site reactivity is also useful for predicting regioselectivity, hydrogen bond basicity, and Lewis molecular basicity in polyfunctional species. In order to evaluate the regional nucleophilicity index given by (9), a series of charged and neutral nucleophiles in the gas phase were evaluated. The results are compiled in Table 1.
 
is the central atom. 
From the data collected in Table 1, it may be seen that the index consistently distributes the global nucleophilicity values on those atoms that are expected to be more nucleophilic. As expected, the most nucleophilic center is at the heteroatom (N, O, or S) site.
3.3. Nucleofugality
As defined at the beginning of this review, the nucleofugality index as stated in (10) may be used to quantitatively characterize the leaving group abilities of several fragments commonly present as nucleofuges in substitution and elimination reactions in organic chemistry. The basic ideas and concepts are those summarized in Scheme 1. In order to illustrate the reliability and usefulness of the nucleofugality index we have evaluated it for a series of aryl benzoates and discuss the usefulness of the resulting scales to assist in the rationalization of their reaction mechanism in nucleophilic substitution reactions [31]. The acyl group in aryl benzoates may react along stepwise or concerted nucleophilic substitution channels to yield the corresponding amide [32–34]. A general route is sketched in Scheme 3 [31].
In order to assess the effect of the leaving group on the reaction mechanism, we selected a series of reactions where the permanent group is kept fixed. The series is depicted in Table 2 [31].
 
^{a}Data from references [32–34]. 
From the data summarized in Table 2, it is possible to note that the electrophilicity index at the leaving group is mainly driven by the presence of electronwithdrawing groups. The presence of these groups is related to an enhancement of the nucleofugality (leaving group ability) of this group. It is important to emphasize that the trends in nucleofugality coherently compare with the experimental rate coefficients measured by Um and coworkers [32–34]. The addition of electronattracting groups on the leaving group may be responsible for the enhancement of the rate constant, since these groups contribute to the net destabilization of the intermediate from which the leaving group detaches.
3.4. Electrofugality
The electrofugality index given by (11) has been applied to a set of 20 benzhydryl sulfinates experimentally studied by Baidya et al. to establish a quantitative hierarchy of electrofugality [73]. According to Scheme 1, the best electrofuge is the fragment that displays the highest regional nucleophilicity [25]. During the bond cleavage the electrofuge is the fragment that releases the electron density, thereby acting as an electron donor. The results are summarized in Table 3.
 
^{a}For this compound the algorithm used to evaluate the nucleophilic Fukui function produces negative values. ^{b}Experimental electrofugality from [33]. ^{c}Predicted values using the empirical equation included in Figure 2(b). 
Table 3 reveals that the regional nucleophilicity is mainly concentrated at fragment R which corresponds to the permanent group. It is important to remark that the electrofugality is fairly ordered in terms of inductive effects promoted by the substituents. The usefulness of this model is illustrated in Figure 4, in the perspective of the predictive potential of this tool: with known electrofugality data for a reduced set of molecules at hand, it becomes possible to anticipate the electrofugality number by simply evaluating the group nucleophilicity of the permanent group moiety at the intermediate complexes [25].
3.5. Site Activation: Fragment Reactivity Analysis
In equilibrium thermodynamics, the electronic chemical potential of density functional theory is a global response function, and therefore it is expected to have a uniform distribution within the whole molecular structure. However, intramolecular reactivity put a serious challenge for the application of electronic descriptors of reactivity, and some adaptations are to be introduced in order to account for how two or more fragments, within the same molecule, may interact to produce bond making/breaking processes. We proposed a model framed on nonequilibrium thermodynamics allowing two molecular fragments to be described by different electronic chemical potential, in such a way that they can exchange work, heat, or energy [74, 75]. We have used the intramolecular DielsAlder (IMDA) reaction of Quinone systems to illustrate this model. The benchmark reaction used was the synthesis of the diterpenoid elisabethin A. The reaction is sketched in Scheme 4 [76].
The fragmentation scheme together with the global and regional indices associated with the diene (D) and Dienophile (Dp) fragments is summarized in Table 4.

The global electrophilicity of quinone 1 is within the range of strong electrophiles [69, 70]. The most electrophilic center of 1 that may react is the C1 carbon. Consequently, the favored IMDA reaction is that involving the diene fragment and the C1C2 double bond. Note the remarkable resolution of the philicity patterns, a job that is performed by the Fukui function (see (8) and (9)). The fragmentation scheme used put over 98% of electrophilicity on the Dp fragment and over 99% of nucleophilicity on the D moiety.
3.6. QuasiStatic Approach: Reactivity Indices along a Reaction Pathway
A final word concerning the reactivity indices in organic reactions is a brief discussion on a nucleophilic aromatic substitution reaction, within a quasistatic scheme, obtained by following the changes in electrophilicity and nucleophilicity along the intrinsic reaction coordinates (IRC). The model reaction is that sketched in Scheme 5 [45, 55].
The model reaction is the involving morpholine and 1X2,4dinitrobenzenes (XDNB, X = F, Cl, Br, and I) [77, 78]. Figures 5 and 6 display the group electrophilicity and nucleophilicity profiles using the arbitrary fragmentation scheme below [25, 45].
(a)
(b)
(c)
(a)
(b)
(c)
As expected, the electrophilicity of the amine moiety (N) is marginal (Figure 5(a)). The electrophilicity condensed at the LG fragment shows a sudden enhancement for the iodine derivative. This result suggests that iodine may detach in the first stage of the reaction depicted in Scheme 5, in agreement with the experimental reports [77, 78]. Figure 5(c) shows the role of the permanent group on the reaction mechanism. In the region of the Meisenheimer complex (MC) all the substrates become electronically saturated. This result may be associated with the end of the charge transfer process, except for iodine which began to detach from the structure as iodide.
Figure 6 shows the nucleophilicity profile within the partitioning Scheme 6. When the index is condensed over the nucleophilic moiety, it seems that it dramatically diminishes downward the MC formation. This result suggests that the charge transfer reaches its minimum after the nucleophilic attack. It is important to note that the property condensed at the LG and PG reaches a maximum value near the MC: the charge is transferred from the nucleophile and accepted for the PG and then redirected to the LG. However, the maximum values for morpholine are different: 7.29, 6.91, 6.56, and 6.11 eV for X = F, Cl, Br, and I, respectively. These values suggest that the substrates interact in different way depending on the LG present. The results reported are in good agreement with the experimental observations, since the leaving group abilities are in the order F > Cl > Br > I when the nucleophile is morpholine [45, 55, 78].
4. Concluding Remarks and Perspectives
In this review, we have shown how conceptual aspects of density functional theory lead to the definition of reactivity indices amenable to build up quantitative models of reactivity in organic reactions. The emphasis has been put on two basic concepts describing electronrich and electrondeficient systems, namely, nucleophile and electrophiles. We then show that the regional nucleophilicity and electrophilicity become the natural descriptors of electrofugality and nucleofugality, respectively. In this way, we have obtained a closed body of concepts that suffices to describe electron releasing and electron accepting molecules together with the ordering of leaving group ability of nucleofuges present in addition, nucleophilic substitution and elimination reactions. A natural perspective of these models is their use along a reaction coordinate following that result in an additional tool to rationalize reaction mechanisms in organic chemistry.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
References
 J. D. Roberts and A. Streitwieser Jr., “Quantum mechanical calculations of orientation in aromatic substitution,” Journal of the American Chemical Society, vol. 74, no. 18, pp. 4723–4725, 1952. View at: Publisher Site  Google Scholar
 A. A. Frost and B. Musulin, “A mnemonic device for molecular orbital energies,” The Journal of Chemical Physics, vol. 21, no. 3, pp. 572–573, 1953. View at: Google Scholar
 L. Salem, “Intermolecular orbital theory of the interaction between conjugated systems. I. General theory,” Journal of the American Chemical Society, vol. 90, no. 3, pp. 543–552, 1968. View at: Publisher Site  Google Scholar
 G. Klopman, “Chemical reactivity and the concept of charge and frontiercontrolled reactions,” Journal of the American Chemical Society, vol. 90, no. 2, pp. 223–234, 1968. View at: Publisher Site  Google Scholar
 E. Z. Hückel, “Quantentheoretische Beiträge zum Benzolproblem,” Zeitschrift für Physik, vol. 70, pp. 204–286, 1931. View at: Publisher Site  Google Scholar
 W. Heitler and F. London, “Wechselwirkung neutraler Atome und homöopolare Bindung nach der Quantenmechanik,” Zeitschrift für Physik, vol. 44, no. 67, pp. 455–472, 1927. View at: Publisher Site  Google Scholar
 C. A. Coulson and H. C. LonguetHiggins, “The electronic structure of conjugated systems. I. General theory,” Proceedings of the Royal Society of London Series A: Mathematical and Physical Sciences, vol. 191, no. 1024, pp. 39–60, 1947. View at: Publisher Site  Google Scholar
 C. A. Coulson and H. C. LonguetHiggins, “The electronic structure of conjugated systems. II. Unsaturated hydrocarbons and their heteroderivatives,” Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences, vol. 192, no. 1028, pp. 16–32, 1947. View at: Publisher Site  Google Scholar
 C. A. Coulson and H. C. LonguetHiggins, “The electronic structure of conjugated systems. III. Bond orders in unsaturated molecules; IV. Force constants and interaction constants in unsaturated hydrocarbons,” Proceedings of the Royal Society of London A, vol. 193, no. 1035, pp. 447–464, 1948. View at: Google Scholar
 R. G. Parr and P. K. Chattaraj, “Principle of maximum hardness,” Journal of the American Chemical Society, vol. 113, no. 5, pp. 1854–1855, 1991. View at: Publisher Site  Google Scholar
 R. G. Parr and R. G. Pearson, “Absolute hardness: companion parameter to absolute electronegativity,” Journal of the American Chemical Society, vol. 105, no. 26, pp. 7512–7516, 1983. View at: Publisher Site  Google Scholar
 R. G. Parr, L. V. Szentpály, and S. Liu, “Electrophilicity index,” Journal of the American Chemical Society, vol. 121, no. 9, pp. 1922–1924, 1999. View at: Publisher Site  Google Scholar
 R. G. Parr and W. Yang, “Density functional approach to the frontierelectron theory of chemical reactivity,” Journal of the American Chemical Society, vol. 106, no. 14, pp. 4049–4050, 1984. View at: Publisher Site  Google Scholar
 J. P. Perdew, R. G. Parr, M. Levy, and J. L. Balduz Jr., “Densityfunctional theory for fractional particle number: derivative discontinuities of the energy,” Physical Review Letters, vol. 49, no. 23, pp. 1691–1694, 1982. View at: Publisher Site  Google Scholar
 R. G. Pearson, “Hard and soft acids and bases,” Journal of the American Chemical Society, vol. 85, no. 22, pp. 3533–3539, 1963. View at: Publisher Site  Google Scholar
 R. G. Pearson, “Absolute electronegativity and hardness: application to inorganic chemistry,” Inorganic Chemistry, vol. 27, no. 4, pp. 734–740, 1988. View at: Publisher Site  Google Scholar
 R. G. Pearson, “The principle of maximum hardness,” Accounts of Chemical Research, vol. 26, no. 5, pp. 250–255, 1993. View at: Publisher Site  Google Scholar
 R. G. Pearson, “Absolute electronegativity and hardness: applications to organic chemistry,” Journal of Organic Chemistry, vol. 54, no. 6, pp. 1423–1430, 1989. View at: Publisher Site  Google Scholar
 R. G. Pearson and J. Songstad, “Application of the principle of hard and soft acids and bases to organic chemistry,” Journal of the American Chemical Society, vol. 89, no. 8, pp. 1827–1836, 1967. View at: Publisher Site  Google Scholar
 P. K. Chattaraj, H. Lee, and R. G. Parr, “HSAB principle,” Journal of the American Chemical Society, vol. 113, no. 5, pp. 1855–1856, 1991. View at: Publisher Site  Google Scholar
 P. W. Ayers and M. Levy, “Perspective on “Density functional approach to the frontierelectron theory of chemical reactivity”,” Theoretical Chemistry Accounts, vol. 103, no. 34, pp. 353–360, 2000. View at: Publisher Site  Google Scholar
 P. W. Ayers and R. G. Parr, “Variational principles for describing chemical reactions: the Fukui function and chemical hardness revisited,” Journal of the American Chemical Society, vol. 122, no. 9, pp. 2010–2018, 2000. View at: Publisher Site  Google Scholar
 M. GonzalezSuarez, A. Aizman, and R. Contreras, “Phenomenological chemical reactivity theory for mobile electrons,” Theoretical Chemistry Accounts, vol. 126, no. 1, pp. 45–54, 2010. View at: Publisher Site  Google Scholar
 M. GonzalezSuarez, A. Aizman, J. SotoDelgado, and R. Contreras, “Bond Fukui functions as descriptor of the electron density reorganization in π conjugated systems,” The Journal of Organic Chemistry, vol. 77, no. 1, pp. 90–95, 2011. View at: Google Scholar
 R. OrmazábalToledo, P. R. Campodónico, and R. Contreras, “Are electrophilicity and electrofugality related concepts? A density functional theory study,” Organic Letters, vol. 13, pp. 822–824, 2011. View at: Publisher Site  Google Scholar
 P. R. Campodónico, A. Aizman, and R. Contreras, “Electrophilicity of quinones and its relationship with hydride affinity,” Chemical Physics Letters, vol. 471, no. 1–3, pp. 168–173, 2009. View at: Publisher Site  Google Scholar
 R. G. Parr and W. Yang, DensityFunctional Theory of Atoms and Molecules, Oxford University Press, New York, NY, USA, 1989.
 T. Koopmans, “Über die Zuordnung von Wellenfunktionen und Eigenwerten zu den Einzelnen Elektronen Eines Atoms,” Physica, vol. 1, no. 1–6, pp. 104–113, 1934. View at: Publisher Site  Google Scholar
 R. Vianello, N. Peran, and Z. B. Maksić, “Hydride affinities of some substituted alkynes: prediction by DFT calculations and rationalization by triadic formula,” Journal of Physical Chemistry A, vol. 110, no. 47, pp. 12870–12881, 2006. View at: Publisher Site  Google Scholar
 R. Contreras, J. Andres, V. S. Safont, P. Campodonico, and J. G. Santos, “A theoretical study on the relationship between nucleophilicity and ionization potentials in solution phase,” Journal of Physical Chemistry A, vol. 107, no. 29, pp. 5588–5593, 2003. View at: Publisher Site  Google Scholar
 P. R. Campodónico, R. OrmazábalToledo, A. Aizman, and R. Contreras, “Permanent group effect on nucleofugality in aryl benzoates,” Chemical Physics Letters, vol. 498, pp. 221–225, 2010. View at: Publisher Site  Google Scholar
 I. Um, H. Han, J. Ahn, S. Kang, and E. Buncel, “Reinterpretation of curved hammett plots in reaction of nucleophiles with aryl benzoates: change in ratedetermining step or mechanism versus groundstate stabilization,” Journal of Organic Chemistry, vol. 67, no. 24, pp. 8475–8480, 2002. View at: Publisher Site  Google Scholar
 I. Um, S. Jeon, and J. Seok, “Aminolysis of 2,4dinitrophenyl Xsubstituted benzoates and Ysubstituted phenyl benzoates in MeCN: Effect of the reaction medium on rate and mechanism,” Chemistry, vol. 12, no. 4, pp. 1237–1243, 2006. View at: Publisher Site  Google Scholar
 I. H. Um, J. Y. Lee, M. Fujio, and Y. Tsuno, “Structurereactivity correlations in nucleophilic substitution reactions of Ysubstituted phenyl Xsubstituted benzoates with anionic and neutral nucleophiles,” Organic & Biomolecular Chemistry, vol. 4, no. 15, pp. 2979–2985, 2006. View at: Publisher Site  Google Scholar
 R. Contreras, J. Andrés, P. Pérez, A. Aizman, and O. Tapia, “Theory of nonlocal (pair site) reactivity from model staticdensity response functions,” Theoretical Chemistry Accounts, vol. 99, pp. 183–191, 1998. View at: Publisher Site  Google Scholar
 W. Yang and W. J. Mortier, “The use of global and local molecular parameters for the analysis of the gasphase basicity of amines,” Journal of the American Chemical Society, vol. 108, no. 19, pp. 5708–5711, 1986. View at: Publisher Site  Google Scholar
 R. R. Contreras, P. Fuentealba, M. Galván, and P. Pérez, “A direct evaluation of regional Fukui functions in molecules,” Chemical Physics Letters, vol. 304, no. 56, pp. 405–413, 1999. View at: Publisher Site  Google Scholar
 P. Fuentealba, P. Pérez, and R. Contreras, “On the condensed Fukui function,” Journal of Chemical Physics, vol. 113, no. 7, pp. 2544–2551, 2000. View at: Publisher Site  Google Scholar
 P. Bultinck, C. Cardenas, P. Fuentealba, P. A. Johnson, and P. W. Ayers, “How to compute the Fukui matrix and function for systems with (Quasi)degenerate states,” Journal of Chemical Theory and Computation, vol. 10, no. 1, pp. 202–210, 2014. View at: Publisher Site  Google Scholar
 A. T. Maynard, M. Huang, W. G. Rice, and D. G. Covell, “Reactivity of the HIV1 nucleocapsid protein p7 zinc finger domains from the perspective of densityfunctional theory,” Proceedings of the National Academy of Sciences of the United States of America, vol. 95, no. 20, pp. 11578–11583, 1998. View at: Publisher Site  Google Scholar
 P. Pérez, A. Aizman, and R. Contreras, “Comparison between experimental and theoretical scales of electrophilicity based on reactivity indexes,” Journal of Physical Chemistry A, vol. 106, no. 15, pp. 3964–3966, 2002. View at: Publisher Site  Google Scholar
 P. Pérez, L. R. Domingo, M. J. Aurell, and R. Contreras, “Quantitative characterization of the global electrophilicity pattern of some reagents involved in 1,3dipolar cycloaddition reactions,” Tetrahedron, vol. 59, no. 17, pp. 3117–3125, 2003. View at: Publisher Site  Google Scholar
 P. R. Campodónico, J. Andrés, A. Aizman, and R. Contreras, “Nucleofugality index in αelimination reactions,” Chemical Physics Letters, vol. 439, pp. 177–182, 2007. View at: Publisher Site  Google Scholar
 P. R. Campodónico, P. Fuentealba, E. A. Castro, J. G. Santos, and R. Contreras, “Relationships between the electrophilicity index and experimental rate coefficients for the aminolysis of thiolcarbonates and dithiocarbonates,” Journal of Organic Chemistry, vol. 70, no. 5, pp. 1754–1760, 2005. View at: Publisher Site  Google Scholar
 R. OrmazábalToledo, R. Contreras, and P. R. Campodónico, “Reactivity indices profile: a companion tool of the potential energy surface for the analysis of reaction mechanisms. Nucleophilic aromatic substitution reactions as test sase,” Journal of Organic Chemistry, vol. 78, pp. 1091–1097, 2013. View at: Google Scholar
 P. Pérez and R. Contreras, “A theoretical analysis of the gasphase protonation of hydroxylamine, methylderivatives and aliphatic amino acids,” Chemical Physics Letters, vol. 293, no. 34, pp. 239–244, 1998. View at: Publisher Site  Google Scholar
 K. Neimann and R. Neumann, “Electrophilic activation of hydrogen peroxide: Selective oxidation reactions in perfluorinated alcohol solvents,” Organic Letters, vol. 2, no. 18, pp. 2861–2863, 2000. View at: Publisher Site  Google Scholar
 R. Das, J. L. Vigneresse, and P. K. Chattaraj, “Redox and Lewis acid–base activities through an electronegativityhardness landscape diagram,” Journal of Molecular Modeling, vol. 19, no. 11, pp. 4857–4864, 2013. View at: Publisher Site  Google Scholar
 W. Zhang, Y. Zhu, D. Wei, Y. Li, and M. Tang, “Theoretical investigations toward the [4 + 2] cycloaddition of ketenes with Nbenzoyldiazenes catalyzed by Nheterocyclic carbenes: mechanism and enantioselectivity,” Journal of Organic Chemistry, vol. 77, no. 23, pp. 10729–10737, 2012. View at: Publisher Site  Google Scholar
 A. Cedillo, R. Contreras, M. Galván, A. Aizman, J. Andrés, and V. S. Safont, “Nucleophilicity index from perturbed electrostatic potentials,” Journal of Physical Chemistry A, vol. 111, no. 12, pp. 2442–2447, 2007. View at: Publisher Site  Google Scholar
 P. Jaramillo, P. Pérez, R. Contreras, W. Tiznado, and P. Fuentealba, “Definition of a nucleophilicity scale,” Journal of Physical Chemistry A, vol. 110, no. 26, pp. 8181–8187, 2006. View at: Publisher Site  Google Scholar
 P. Campodonico, J. G. Santos, J. Andres, and R. Contreras, “Relationship between nucleophilicity/electrophilcity indices and reaction mechanisms for the nucleophilic substitution reactions of carbonyl compounds,” Journal of Physical Organic Chemistry, vol. 17, no. 4, pp. 273–281, 2004. View at: Publisher Site  Google Scholar
 P. R. Campodonico, A. Aizman, and R. Contreras, “Group electrophilicity as a model of nucleofugality in nucleophilic substitution reactions,” Chemical Physics Letters, vol. 422, no. 4–6, pp. 340–344, 2006. View at: Publisher Site  Google Scholar
 S. Jorge, A. Aizman, R. Contreras, and L. R. Domingo, “On the catalytic effect of water in the intramolecular dielsalder reaction of quinone systems: a theoretical study,” Molecules, vol. 17, no. 11, pp. 13687–13703, 2012. View at: Publisher Site  Google Scholar
 R. OrmazábalToledo, R. Contreras, R. A. Tapia, and P. R. Campodónico, “Specific nucleophileelectrophile interactions in nucleophilic aromatic substitutions,” Organic and Biomolecular Chemistry, vol. 11, pp. 2302–2309, 2013. View at: Publisher Site  Google Scholar
 A. Aizman, R. Contreras, and P. Pérez, “Relationship between local electrophilicity and rate coefficients for the hydrolysis of carbenium ions,” Tetrahedron, vol. 61, no. 4, pp. 889–895, 2005. View at: Publisher Site  Google Scholar
 R. Contreras, J. Andrés, L. R. Domingo, R. Castillo, and P. Pérez, “Effect of electronwithdrawing substituents on the electrophilicity of carbonyl carbons,” Tetrahedron, vol. 61, pp. 417–422, 2005. View at: Publisher Site  Google Scholar
 S. Karabunarliev, O. G. Mekenyan, W. Karcher, C. L. Russom, and S. P. Bradbury, “Quantumchemical descriptors for estimating the acute toxicity of electrophiles to the bathed minnow (Pimephales promelas): an analysis based on molecular mechanisms,” Quantitative StructureActivity Relationships, vol. 15, no. 4, pp. 302–310, 1996. View at: Publisher Site  Google Scholar
 P. J. O'Brien, “Molecular mechanisms of quinone cytotoxicity,” ChemicoBiological Interactions, vol. 80, no. 1, pp. 1–41, 1991. View at: Publisher Site  Google Scholar
 D. R. Roy, R. Parthasarathi, B. Maiti, V. Subramanian, and P. K. Chattaraj, “Electrophilicity as a possible descriptor for toxicity prediction,” Bioorganic and Medicinal Chemistry, vol. 13, no. 10, pp. 3405–3412, 2005. View at: Publisher Site  Google Scholar
 P. K. Chattaraj, U. Sarkar, and D. R. Roy, “Electrophilicity index,” Chemical Reviews, vol. 106, no. 6, pp. 2065–2091, 2006. View at: Publisher Site  Google Scholar
 P. Thanikaivelan, V. Subramanian, J. Raghava Rao, and B. Unni Nair, “Application of quantum chemical descriptor in quantitative structure activity and structure property relationship,” Chemical Physics Letters, vol. 323, no. 12, pp. 59–70, 2000. View at: Publisher Site  Google Scholar
 R. Parthasarathi, V. Subramanian, D. R. Roy, and P. K. Chattaraj, “Electrophilicity index as a possible descriptor of biological activity,” Bioorganic & Medicinal Chemistry, vol. 12, no. 21, pp. 5533–5543, 2004. View at: Publisher Site  Google Scholar
 A. CerdaMonje, A. Aizman, R. A. Tapia, C. Chiappe, and R. Contreras, “Solvent effects in ionic liquids: empirical linear energydensity relationships,” Physical Chemistry Chemical Physics, vol. 14, no. 28, pp. 10041–10049, 2012. View at: Publisher Site  Google Scholar
 K. B. Yatsimirskii, “Hydride affinity as a measure of acidity (electrophilicity),” Theoretical and Experimental Chemistry, vol. 17, no. 1, pp. 75–79, 1981. View at: Publisher Site  Google Scholar
 X. Zhu, C. Wang, H. Liang, and J. Cheng, “Theoretical prediction of the hydride affinities of various p and oquinones in DMSO,” Journal of Organic Chemistry, vol. 72, no. 3, pp. 945–956, 2007. View at: Publisher Site  Google Scholar
 E. Buncel and I. Um, “The αeffect and its modulation by solvent,” Tetrahedron, vol. 60, no. 36, pp. 7801–7825, 2004. View at: Publisher Site  Google Scholar
 N. J. Fina and J. O. Edwards, “The alpha effect. A review,” International Journal of Chemical Kinetics, vol. 5, no. 1, pp. 1–26, 1973. View at: Publisher Site  Google Scholar
 L. R. Domingo, M. J. Aurell, P. Pérez, and R. Contreras, “Quantitative characterization of the global electrophilicity power of common diene/dienophile pairs in Diels–Alder reactions,” Tetrahedron, vol. 58, no. 22, pp. 4417–4423, 2002. View at: Publisher Site  Google Scholar
 L. R. Domingo, M. J. Aurell, P. Pérez, and R. Contreras, “Quantitative characterization of the local electrophilicity of organic molecules. Understanding the regioselectivity on DielsAlder reactions,” The Journal of Physical Chemistry A, vol. 106, no. 29, pp. 6871–6875, 2002. View at: Publisher Site  Google Scholar
 A. C. Legon, “Quantitative gasphase electrophilicities of the dihalogen molecules XY = F2, Cl2, Br2, BrCl and ClF,” Chemical Communications, no. 23, pp. 2585–2586, 1998. View at: Google Scholar
 A. C. Legon and D. J. Millen, “Hydrogen bonding as a probe of electron densities: limiting gasphase nucleophilicities and electrophilicities of B and HX,” Journal of the American Chemical Society, vol. 109, no. 2, pp. 356–358, 1987. View at: Publisher Site  Google Scholar
 M. Baidya, S. Kobayashi, and H. Mayr, “Nucleophilicity and nucleofugality of phenylsulfinate (PhSO_{2}^{}): A key to understanding its ambident reactivity,” Journal of the American Chemical Society, vol. 132, no. 13, pp. 4796–4805, 2010. View at: Publisher Site  Google Scholar
 J. SotoDelgado, A. Aizman, L. R. Domingo, and R. Contreras, “Invariance of electrophilicity of independent fragments. Application to intramolecular DielsAlder reactions,” Chemical Physics Letters, vol. 499, no. 4–6, pp. 272–277, 2010. View at: Publisher Site  Google Scholar
 J. SotoDelgado, L. R. Domingo, and R. Contreras, “Quantitative characterization of group electrophilicity and nucleophilicity for intramolecular DielsAlder reactions,” Organic and Biomolecular Chemistry, vol. 8, no. 16, pp. 3678–3683, 2010. View at: Publisher Site  Google Scholar
 T. J. Heckrodt and J. Mulzer, “Total synthesis of elisabethin A: Intramolecular DielsAlder reaction under biomimetic conditions,” Journal of the American Chemical Society, vol. 125, no. 16, pp. 4680–4681, 2003. View at: Publisher Site  Google Scholar
 B. G. Cox and A. J. Parker, “Solvation of ions. XVIII. Proticdipolar aprotic solvent effects on the free energies, enthalpies, and entropies of activation of an SNAr reaction,” Journal of the American Chemical Society, vol. 95, no. 2, pp. 408–410, 1973. View at: Publisher Site  Google Scholar
 I. Um, L. Im, J. Kang, S. S. Bursey, and J. M. Dust, “Mechanistic assessment of S NAr displacement of halides from 1Halo2,4dinitrobenzenes by selected primary and secondary amines: Brønsted and Mayr analyses,” Journal of Organic Chemistry, vol. 77, no. 21, pp. 9738–9746, 2012. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Rodrigo OrmazábalToledo and Renato Contreras. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.