Research Article  Open Access
Alejandro MoralesBayuelo, Ricardo VivasReyes, "Topological Model on the Inductive Effect in Alkyl Halides Using Local Quantum Similarity and Reactivity Descriptors in the Density Functional Theory", Journal of Quantum Chemistry, vol. 2014, Article ID 850163, 12 pages, 2014. https://doi.org/10.1155/2014/850163
Topological Model on the Inductive Effect in Alkyl Halides Using Local Quantum Similarity and Reactivity Descriptors in the Density Functional Theory
Abstract
We present a topological analysis to the inductive effect through steric and electrostatic scales of quantitative convergence. Using the molecular similarity field based in the local guantum similarity (LQS) with the TopoGeometrical Superposition Algorithm (TGSA) alignment method and the chemical reactivity in the density function theory (DFT) context, all calculations were carried out with Amsterdam Density Functional (ADF) code, using the gradient generalized approximation (GGA) and local exchange correlations PW91, in order to characterize the electronic effect by atomic size in the halogens group using a standard Slatertypeorbital basis set. In addition, in this study we introduced news molecular bonding relationships in the inductive effect and the nature of the polar character in the C–H bond taking into account the global and local reactivity descriptors such as chemical potential, hardness, electrophilicity, and Fukui functions, respectively. These descriptors are used to find new alternative considerations on the inductive effect, unlike to the binding energy and dipole moment performed in the traditional organic chemical.
1. Introduction
The inductive effect is one of the most important electronic effects in chemistry [1]. According to the International Union of Pure and Applied Chemistry (IUPAC) the inductive effect is defined as transmission effect of charge on a chain of atoms via electrostatic induction [2]. The inductive effect is experimentally observable via molecular polarization of polar covalent bonds produced by electrostatic induction, due to differences in electronegativity between the atoms involved [3]. In this contribution, our attention will be focused in the quantification of the molecular polarization from the point of view of the molecular bond in alkyl halides molecules such as _{3}HC–X, (X = F, Cl, Br) taken as example. Using an atomic approach of molecular quantum similarity performed by Cioslowski and Nanayakkara [4] and the Molecular Quantum Similarity Indexes (MQSI) based on the electron density proposed by CarbóDorca and coworkers [5–9].
Several concepts and theories have been proposed in the literature in order to explain the molecular chemical bond [10–19]. Nowadays, the most used theories are the molecular orbital theory [20, 21] and the valence bond method [22, 23]. In this paper was used the concept of chemical bond as minimization of equilibrium interactions determined by the bond distances at the experimental level. We propose to relate these experimental distances with the dis(similarity) electrostatic and steric, using quantum similarity descriptors such as overlap and Coulomb indexes, with their Euclidean distances.
On the other hand, we present in this contribution a theoretical investigation in order to relate the inductive effect in these alkyl halides and analyze its dependence on the bond distance with respect to the variation of electronegativity of substituent atoms (–F, –Cl, –Br). This approach could shed some light onto some molecular aspects that can contribute to the knowledge behind interactions of these molecules.
Another reason for carrying out this study is to relate the dis(similarity) between the inductive effect and the nature of chemical bond with differences of electronegativity between atoms. In addition, this approach suggests a new perspective to the treating of the chemical bond with the inductive effect, taking into account that the search for new interpretations is a field of great interest in organic physical chemistry; in addition to the inductive effect also can be used to determine the stability of molecule with respect to their charge distribution [24–26].
2. Theory and Computational Details
2.1. Similarity Indexes
Reviews about concepts and applications in quantum similarity can be found elsewhere (see, for example, [5–9, 27, 28]); only general definitions will be given at this point. Quantum similarity provides a quantitative measure of the resemblance between two quantum objects A and B.
Carbó and coworkers introduced the similarity indexes [5, 27]; they defined the quantum similarity measure between molecules A and B with electron density and , based on the idea of minimizing the expression for the Euclidean distance: where is the Euclidean distance, involving the overlap integral between the electron density of the molecule A and B and and are the selfsimilarity of molecules A and B [27].
The most common index is the generalized quantum similarity by the cosine, introduced by Carbó et al. [28]; this index can be expressed mathematically as Introducing () elements in the operator () we can write where () is the Carbó index; this index is mathematically defined in the interval (0,1], where 1 is the selfsimilarity and calculates only the measures of “shape similarity,” as an other alternative is the of HodgkimRichards index [29, 30]. This index appears naturally when using the arithmetic mean and can be defined mathematically as And is another alternative of quantification of the quantum similarity.
2.2. Local Similarity Indexes
2.2.1. The Hirshfeld Approach
The Hirshfeld approach [32–38]. This approach is based on partitioning the electron density of a molecule in contributions . In this study, the Hirshfeld approach is used to express the global similarity index (4) to local level.
These contributions are proportional to the weight of the electron density of the isolated molecule in the call promolecular density [39]. The “promolecular density” is obtained as The contribution of atom A in the electron density is The weight () is obtained as where is the electron density of the isolated atom A [40–46]. This idea has been applied in this study to CH_{3}–F taken as molecule A and CH_{3}–X, (X = Cl or Br) designating as B. The atomic contribution of carbon atom (C) in the molecule A is given by with For the atomic contribution of carbon atom (C) in the molecule B is obtained as with so that we can express the numerator in the Carbó index in (4) as where it is possible to get atomic contributions and considerations of similarity.
2.2.2. Description of Molecular Bond Proposed
To study the effects by varying electronegativity by substituents, we propose a particular characterization of the chemical bond which can be considered complementary to the Hirshfeld partition, considering that the inductive effect is dependent on the interatomic bond length [2].
We propose to study the molecular bond based on the physical concepts of bond length, defined as the minimum length (equilibrium) of the interactions according to the chemical bond model proposed by Burrau [24], these equilibrium distances are represented by () and we propose to multiply the equilibrium lengths () with the (MQSI), so that we can express: where is (4).
For the Euclidean distance, we have where is defined in (1).
Taking into account that () is related with the overlap similarity index (OSI)) and with the Coulomb similarity index (CSI)) and can be interpreted as equilibrium distances that minimize the overlap and the coulomb interactions. So we can define the descriptor of chemical bonding (13) as the equilibrium distance that minimizes the structural and electronic indexes; in the sense we can express the total contributions to the molecular bond from the point of view of the quantum similarity indexes as Equation (16) has a range of convergence of (0, 2], where 2 means maximum selfsimilarity.
Similar to the contributions of the Euclidean distances Equation (17) has a convergence range [0,), where 0 represents the selfsimilarity. This means that small distances of CH_{3}–F with respect to CH_{3}–Cl and CH_{3}–Br are associated with strong bonding force and a high inductive character.
3. Global and Local Reactivity Descriptors
The global reactivity descriptors are defined within the (DFT) and were interpreted by Parr and coworkers [47–54]. The chemical potential can be written mathematically in terms of the energy of an electron in the frontier molecular orbitals (higher occupied molecular orbital) HOMO and (lowest unoccupied molecular orbital) LUMO as From (18), one can obtain quantitative expression for the chemical hardness () [48, 55–57], and it is understood as the opposition of the system to distort the electron cloud and mathematically can be written as On the other hand, we have the global electrophilicity () introduced by Parr et al. [58], which is a measure of the stabilization energy of the system when it is saturated by electrons from the external environment and represented mathematically as The Fukui function was one of the descriptors used in this work and defined as as is discontinuous at integer values of . There are three types according to Fukui, f^{+} (22) which contains information on the reactivity local of nucleophilic attack, f^{−} (23) which does the same for an electrophilic attack [59–64]:
4. Molecular Alignment and Computational Details
The principal problem with the calculation of the quantum similarity measures (QSM) is the dependence the relative orientation of the molecules. The problem has been addressed through the use of several methodologies and approaches. Many critera have been used in order to deal with this problem, such as aligning the molecules according to common chemical features (e.g., threedimensional molecular fields such as steric, electrostatic, or hydrophobic fields, among others) [65].
Also these molecules can be aligned based on topological and geometrical features as has been implemented in the TopoGeometrical Superposition Algorithm (TGSA) [66]. In this method the alignment is based on the recognition of a common skeleton, taking only the types of atoms and bond length [67–70], in this paper the recognition of the common carbon atom was used as is shown in Figure 1.
4.1. Density Functional Calculations
In order to characterize the effects by atomic size, all calculations were carried out by using the Amsterdam Density Functional (ADF) program package [71, 72]. ADF is a program for firstprinciple electronic structure calculations that make use of (DFT) [50]. All the molecular structures were fully optimized using the gradient generalized approximation (GGA) [73–76] local exchange correlations PW91 [77–79]. The geometry optimizations were calculated using a standard Slatertypeorbital basis set with triplequality plus doublepolarization functions for all the atoms.
5. Results and Discussion
The results reported in this work are those associated with the molecular quantum similarity indexes, followed by results of the relationship of bond proposed, and finally the results of the reactivity descriptors.
5.1. Molecular Quantum Similarity Analysis
In Tables 1–4 the (MQS) matrices for the compounds analyzed are shown, taken as an example.




In Tables 1 and 2, is possible to see that the CH_{3}–Cl and CH_{3}–Br have the higher quantum similarity with an overlap value of 0.98 and a Coulomb value of 0.99. With a Euclidean distance of Overlap of (0.27, see Table 3) and a Euclidean distance of Coulomb of 0.72 (see Table 4). The high values associated with the (MQS) of Coulomb are justified, taking into account that the alkyl halides are taken as examples of hydrocarbon derivatives in which one or more C–H bonds are replaced by C–X bonds (X = F, Cl, Br), and the similarity of C–X bond is the result of dis(similarity) in the carbon sp^{3} hybrid orbital; in this sense, studies based on the molecular orbital theory suggest that the hybrid orbital of the halogen is due to one small proportion p and to the s character [80]; see Figure 2.
To obtain the relations proposed in this paper the bond lengths obtained from reference [31] were used (Table 5).

We can calculate the relations of bonding described in this study, for (MQSI) of the CH_{3}–Cl and CH_{3}–Br with respect to CH_{3}–F, taking into account that the strong differences in electronegativity determine the molecular polarization effects quantified by the MQSI Coulomb (Table 2).
In Table 6 the values associated with the bond relationships proposed in this study are depicted to quantify the inductive effect; with the overlap index, the values associated with the bond relationships are 109.88 pm when comparing the CH_{3}–F with the CH_{3}–Cl and 103.18 pm when compared with the CH_{3}–Br (see Figure 3(a)); these results are consistent with the quantum similarity indexes of overlap (see Table 1). Comparing the values associated with the relationship of bond proposed we can see that the Coulomb descriptor has the same trend (CH_{3}–Cl > CH_{3}–Br).
 
Bonding relationship of overlap. ^{ b}Bonding relationship of coulomb. ^{ c}Euclidean distances of overlap. ^{ d}Euclidean distances of coulomb. 
(a)
(b)
(c)
(d)
On the other hand, the predicted values for the Euclidean distance ((b) and (d)) are closest to the reference value (134 pm) of CH_{3}–F with respect to CH_{3}–Cl and CH_{3}–Br, respectively, reflecting the strong correlation between the molecular set from the electrostatic point of view, such that trends were sketched by the proposed scales and ae in agreement with the theoretical and experimental studies about the effect of the halogen substituent [81–83].
5.2. Analysis of Global and Local Reactivity Descriptors
Although the C–H bonds of alkyl halides are covalent and have some polar character, because the halogens are more electronegative than the carbon; in order to characterize from the chemical reactivity point of view, the polar character (C–H) is calculated in the global and local reactivity descriptions based in DFT.
In Table 7 we can see that the CH_{3}–Cl is the most stable compound with chemical potential of eV and the most reactive compound is the CH_{3}–F with eV. While the CH_{3}–Br have a chemical potential intermediate with a value eV, the analysis of the chemical potential leads to some general comments. The difference in chemical potential and electronegativity among the alkyl halides is also shown in the greatest opposition to the electronic flow from the electron cloud, and this conclusion is obtained with the value of chemical hardness to CH_{3}–F with eV and eV to CH_{3}–Br. On the other hand, such trend is also observed with the values of ionization potential of the CH_{3}–Br with a value of eV with the most low energy to remove electrons from the electron cloud and the electron affinity eV.
 
: chemical potential. : chemical hardness. IP^{c}: Ionization potential. EA^{d}: electron affinity. : electrophilicity. 
The electrophilic character is calculated with the largest stabilization energy when it is saturated by electrons from the external environment, with a value of : 0.83 eV to CH_{3}–F, : 2.10 eV to CH_{3}–Cl, and : 2.02 eV to CH_{3}–Br. Shows that these global reactivity indexes are appropriate to study the chemical reactivity of inductive effect from a different way than that has been done in the chemical literature (binding energy). Figure 4 shows the trends of the global reactivity descriptors.
(a)
(b)
(c)
Figure 4 shows the same trend to the hardness () and the ionization potential (IP). According to the trends of overlap and Coulomb similarity indexes (see Figure 3), enabling the relationship of the trends in the quantum similarity descriptors of overlap and Coulomb with the chemical hardness () and the ionization potential (IP).
On the other hand, the trend in the electron affinity (EA) is in agreement with those found for the Euclidean distances of overlap and Coulomb and can be associated with a high electronic similarity determined by the differences in electronegativity in the molecular set and calculated by the Coulomb operator; see Table 2. In order to analyze the type of interactions present in C–X and C–H bonds the condensed Fukui functions were calculated as it shown in Table 8 to develop the topological model presented in this study.

The values of the Fukui function of CH_{3}–F on carbon atom are f^{+}: 0.38 and f^{−}: 0.66 (see Table 8), respectively. Is evident that the molecular polarization due to the delta of electronegativity (Table 5) together with the value of the Fukui function (f^{−}: 0.66) represents the high susceptibility of attacks by electrophilic species. These values found are in agreement with the calculated for the polar character of the hydrogen atoms due to the fluoride in the C–H bond. To observe this type of interactions in Figure 5 the topological Fukui functions f^{+}: 0.38 and f^{−}: 0.66 are shown for the CH_{3}–F, and also we can see the polar character of the C–H bonds and their contributions to the chemical reactivity due to the inductive effect.
(a)
(b)
(c)
Figure 6 shows the nature of the polar interactions in the C–H bonds of CH_{3}–Cl through the contour maps associated with the Fukui functions f^{+}: 0.53 for the carbon atom and f^{+}: 0.42 for the chloride atom, and these two sites are highly reactive to nucleophilic species, unlike to the CH_{3}–F that only has a single reactive site f^{+}: 0.38.
(a)
(b)
(c)
Finally, in Figure 7 the reactive sites f^{+}: 0.45 for carbon and f^{+}: 0.49 and f^{−}: 0.95 for the bromide in the CH_{3}–Br are shown (see Table 8). This value for the Fukui function f^{+} shows a lower inductive effect with respect to CH_{3}–Cl and to the CH_{3}–F, respectively, according to the electronegativity values.
(a)
(b)
(c)
The inductive effect model proposed in this study provides an explanation for the polar character from a different way than that has been done in the chemical literature (through the dipole moment). On the other hand, the () indexes were introduced by Roy et al. [37, 84] and are expressed mathematically as Equation (24) can be understood as a measure of the stabilization energy produced by the flow of electrons from the external environment for an nucleophilic attacks, and we postulate that can be related with the electrophilic character of atom more electronegative, that for the CH_{3}–F have (: 0.02) in the fluoride, to the CH_{3}–Cl have (: 0.88) in the chloride, and to the CH_{3}–Br have (: 0.99) to the Bromide (see Table 5).
These values are quantified by the electrophilic character [81–83] produced by the inductive effect in the compounds studied, taking into account that when decrease the value of increases the stabilization energy of atom most electronegative in the system when it is saturated by electrons from the external medium, and in this since, the CH_{3}–F has the higher inductive effect; therefore, we can see the trend (CH_{3}–F > CH_{3}–Cl > CH_{3}–Br) according to the trends of correlation reported in Figure 3, and these theoretical values are in agreement with experimental results [85–87].
From Figures 5–7, we can identify the polarization of covalent bonds and the prediction of reactivity sites susceptible to attack by electrophilic and nucleophilic species, allowing the analysis of the electronic density distribution in the molecular set (CH_{3}–F, CH_{3}–Cl, CH_{3}–Br). Taking into account that the polarization can induce a dipole in a bond near to the substituent or in one located more away (inductive effect (I)), in this sense the Fukui topological (f^{−}) can interact with the electron release (+I) and the Fukui topological (f^{+}) with the electronwithdrawing (−I). Although this does not imply donation formal of a unit charge, considering that by convection is assumed as electronwithdrawing groups when are most stronger than the hydrogen atom and are designated as negative inductive effects (−I).
Unlike to those groups which have a less capacity of electron withdrawing (high electron release) with respect to hydrogen atom and are designated as positive inductive effect (+I) [88], showing new considerations among the electron flow in the inductive effect and the local chemical reactivity, using the description of molecular bond proposed in this study ((13)–(17)). In this sense we analyze the quantum similarity increasing steric and electronic effects.
5.3. Quantum Dissimilarity Analysis
In order to analyze the quantum dissimilarity increasing the steric and electronic effects for analyze the stability of the relations proposed are analyzed the dissimilarity indexes in a series of isopropane derivatives shown in Figure 8.
With the aim of characterizing the steric and electronic effect in Tables 9 and 10we calculate the overlap quantum dissimilarity.


In Table 9 we can see the structural effects on the atomic size, the higher overlap quantum dissimilarity values comparing C_{3}H_{8}I with C_{3}H_{8}F (0.7318), and the lowest values comparing C_{3}H_{8}Br with C_{3}H_{8}Cl (0.9887), forming the overlap quantum dissimilarity trend with respect to C_{3}H_{8} (C_{3}H_{8}F (0.7427) > C_{3}H_{8}Cl (0.8895) > C_{3}H_{8}Br (0.9025) > C_{3}H_{8}I (0.9252)) in agreement with the trends shown in Figure 3.
In Table 10 we found the higher Euclidean distance of overlap between C_{3}H_{8}I and C_{3}H_{8}F (1.5732) and the lowest is between C_{3}H_{8}Br and C_{3}H_{8}Cl (0.2608), in agreement with the values of Table 9. With the aim of analyzing the electronic effects on the quantum dissimilarity the MQS of Coulomb are shown in Table 11.

In Table 11 we can see the higher electronic quantum dissimilarity value between C_{3}H_{8}F and C_{3}H_{8} (0.9468). However, the lowest is between C_{3}H_{8}Br and C_{3}H_{8}Cl (0.9990). Also the electronic quantum dissimilarity trend with respect to C_{3}H_{8} C_{3}H_{8}F (0.9468) > C_{3}H_{8}Cl (0.9631) > C_{3}H_{8}Br (0.9643) > C_{3}H_{8}I (0.9675) is in agreement with the outcomes of Table 9 and Figure 3. This result is consistent with the Euclidean distances of Coulomb (Table 12).

In Table 12 the higher Euclidean distance of Coulomb is between C_{3}H_{8}F and C_{3}H_{8} (5.4638), and the lowest is between C_{3}H_{8}Br and C_{3}H_{8}Cl (0.6427), in agreement with the results of Table 11. This result shows the electronegativity character and the atomic size on the steric and electronic effects related to the inductive effects and extrapolating the result of Figure 3.
6. Conclusions
We present a topological analysis on the inductive effect in CH_{3}–F, CH_{3}–Cl, and CH_{3}–Br alkyl halides taken as an example; using local molecular quantum similarity and chemical reactivity in DFT context, this quantification was performed through scales of quantitative convergence using descriptors based in MQSI of overlap and Coulomb (see Figure 3), presenting a good alternative relationship and the chemical bonds, and it also explains the polar character associated with the C–H bonds from the point of view of the global and local chemical reactivity, in a different way than that has been done in the chemical literature (dipole moments and energy binding); in addition these results are consistent with the effects produced by halogen substituents according to experimental results.
The alignment method used is the TopoGeometrical Superposition Algorithm (TGSA), allowing calculating appropriate measures of the asymmetry in the central carbon atom of the alkyl halides. This TGSA method is a straightforward procedure to cope with the problem of relative orientation of the molecules to evaluate MQS; also this work shows the Hirschfield partitioning of the electronic density in the proposed model of covalent bond (see (13)–(17)) can describe and relate the inductive effect in the molecular set (CH_{3}–F, CH_{3}–Cl, and CH_{3}–Br).
The reactivity global indexes can characterize the reactivity of the compounds tested, the higher value of the chemical potential (: −3.14 eV) to the CH_{3}–F, with the greater hardness (: 5.97 eV) and the most stable the CH_{3}–Cl with a chemical potential (: −4.09 eV) and (: 3.97 eV) (see Table 7). The characterization of the reactivity sites calculated by the condensed Fukui functions (see Table 8) helped identify areas susceptible to attack by electrophilic and nucleophilic species, associating topological Fukui functions (f^{−}) with the electronrelease (+I) and the topological Fukui (f^{+}) with the electronwithdrawing (−I).
Moreover, through this model we can be demonstrated the origin of the polar interactions of the C–H bond (see Figures 5–7) that can be applied to a much wider range of compounds such as CH_{3}–NH_{2}, CH_{3}–(CH_{2})_{n}–NH_{2}, CH_{3}–OH, CH_{3}–(CH_{2})_{n}–OH, CH_{3}–SH, and CH_{3}–(CH_{2})_{n}–SH, with (n) large, due to the nature of the TGSA alignment method and in general systems where the inductive effect is marked by the high delta in electronegativity and polarizability, considering that the characterization of inductive effect is an area of growing interest in organic physical chemistry.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
Alejandro MoralesBayuelo thanks the Universidad Nacional Andres Bello (SantiagoChile) for a Ph.D. fellowship (CONICYT (63100003)). Ricardo VivasReyes is indebted to the Universidad de Cartagena (Cartagena, Colombia), for continuous support to the research group.
References
 L. G. Wade Jr., Organic Chemistry, Pearson Prentice Hall, New Jersey, NJ, USA, 6th edition, 2006.
 P. Muller, “Glossary of terms used in physical organic chemistry,” Pure and Applied Chemistry, vol. 66, no. 5, pp. 1077–1148, 1994. View at: Publisher Site  Google Scholar
 UPAC, Compendium of Chemical Terminology, Gold Book, 2nd edition, 1997, Online corrected version (2006).
 J. Cioslowski and A. Nanayakkara, “Similarity of atoms in molecules,” Journal of the American Chemical Society, vol. 115, no. 24, pp. 11213–11215, 1993. View at: Google Scholar
 R. Carbó, M. Arnau, and L. J. Leyda, “How similar is a molecule to another? An electron density measure of similarity between two molecular structures,” International Journal of Quantum Chemistry, vol. 17, pp. 1185–1189, 1980. View at: Publisher Site  Google Scholar
 R. CarbóDorca and L. D. J. Mercado, “Density gradient quantum similarity,” Journal of Computational Chemistry, vol. 31, pp. 2195–2212, 2010. View at: Google Scholar
 X. Gironés and R. CarbóDorca, “Modelling toxicity using molecular quantum similarity measures,” QSAR & Combinatorial Science, vol. 25, no. 7, pp. 579–589, 2006. View at: Google Scholar
 R. CarbóDorca, E. Besalú, and L. D. Mercado, “Communications on quantum similarity. Part 3: a geometricquantum similarity molecular superposition algorithm,” Journal of Computational Chemistry, vol. 32, no. 4, pp. 582–599, 2011. View at: Publisher Site  Google Scholar
 R. CarbóDorca and X. Gironés, “Foundation of quantum similarity measures and their relationship to QSPR: density function structure, approximations, and application examples,” International Journal of Quantum Chemistry, vol. 101, no. 1, pp. 8–20, 2005. View at: Publisher Site  Google Scholar
 J. D. Watson and F. Crick, “Molecular structure of nucleic acids: a structure for deoxyribose nucleic acid,” Nature, vol. 171, pp. 737–738, 1953. View at: Publisher Site  Google Scholar
 P. Su and H. J. Li, “Energy decomposition analysis of covalent bonds and intermolecular interactions,” The Journal of Chemical Physics, vol. 131, no. 1, Article ID 014102, 2009. View at: Publisher Site  Google Scholar
 F. M. Bickelhaupt and E. J. Baerends, “Reviews in computational chemistry,” in WileyVCH, K. B. Lipkowitz and D. B. Boyd, Eds., vol. 15, pp. 1–86, John Wiley and Sons, New York, NY, USA, 2000. View at: Google Scholar
 P. Bao and Z. Yu, “Theoretical studies on the role of πelectron delocalization in determining the conformation of Nbenzylideneaniline with three types of LMO basis sets,” Journal of Computational Chemistry, vol. 27, no. 7, pp. 809–824, 2006. View at: Publisher Site  Google Scholar
 R. S. Berry, S. A. Rice, and J. Ross, Physical Chemistry, Oxford University Press, Oxford, UK, 1st edition, 1980.
 R. S. Berry, S. A. Rice, and J. Ross, Physical Chemistry, Oxford University Press, Oxford, UK, 2nd edition, 2000.
 F. Albert Cotton and G. Wilkinson, Advanced Inorganic Chemistry, John Wiley & Sons, New York, NY, USA, 1st edition, 1962.
 F. Albert Cotton and G. Wilkinson, Advanced Inorganic Chemistry, John Wiley & Sons, New York, NY, USA, 6th edition, 1999.
 W. Stumm and J. J. Morgan, Aquatic Chemistry: Chemical Equilibria and Rates in Natural Waters, John Wiley & Sons, New York, NY, USA, 1st edition, 1970.
 W. Stumm and J. J. Morgan, Aquatic Chemistry: Chemical Equilibria and Rates in Natural Waters, John Wiley & Sons, New York, NY, USA, 3rd edition, 1996.
 I. Fleming, Frontier Orbitals and Organic Chemical Reactions, John Wiley & Sons, New York, NY, USA, 1977.
 K. N. Houk, in Pericyclic Reactions, A. P. Marchand and R. E. Lehr, Eds., vol. 2, pp. 181–271, Academic Press, New York, NY, USA, 1977.
 L. Pauling, The Nature of the Chemical Bond, Cornell University Press, New York, NY, USA, 1960.
 F. Weinhold and C. Landis, Valency and Bonding: A Natural Bond Orbital DonorAcceptor Perspective, Cambridge University Press, Cambridge, UK, 2005.
 Ø. Burrau, “Berechnung des energiewertes des wasserstoffmolekelions (H^{+}_{2}) im normalzustand,” Naturwissenschaften, vol. 15, no. 1, pp. 16–17, 1927. View at: Publisher Site  Google Scholar
 E. Dumont and P. Chaquin, “Inductive effects on proton affinity of benzene derivatives: analysis using fictitious hydrogen atoms,” Journal of Physical Chemistry A, vol. 113, no. 12, pp. 2990–2994, 2009. View at: Publisher Site  Google Scholar
 A. P. Smith, A. E. McKercher, and R. C. Mawhinney, “Inductive effect: a quantum theory of atoms in molecules perspective,” Journal of Physical Chemistry A, vol. 115, no. 45, pp. 12544–12554, 2011. View at: Publisher Site  Google Scholar
 P. Bultinck, X. Gironés, and R. CarbóDorcaz, “Molecular quantum similarity: theory and applications,” Reviews in Computational Chemistry, vol. 21, p. 127, 2005. View at: Publisher Site  Google Scholar
 R. Carbó, J. Arnau, and L. Leyda, “How similar is a molecule to another? An electron density measure of similarity between two molecular structures,” International Journal of Quantum Chemistry, vol. 17, no. 6, pp. 1185–1189, 1980. View at: Publisher Site  Google Scholar
 E. E. Hodgkin and W. G. Richards, “Molecular similarity based on electrostatic potential and electric field,” International Journal of Quantum Chemistry, vol. 14, pp. 105–110, 1987. View at: Publisher Site  Google Scholar
 E. E. Hodgkin and W. G. Richards, “Molecular similarity,” Chemische Berichte, vol. 24, p. 1141, 1988. View at: Google Scholar
 R. C. Weast, Handbook of Chemistry & Physics, CRC Press, New York, NY, USA, 65th edition.
 F. L. Hirshfeld, “Bondedatom fragments for describing molecular charge densities,” Theoretica Chimica Acta, vol. 44, no. 2, pp. 129–138, 1977. View at: Publisher Site  Google Scholar
 R. F. Nalewajski and R. G. Parr, “Information theory, atoms in molecules, and molecular similarity,” Proceedings of the National Academy of Sciences, vol. 97, no. 12, pp. 8879–8882, 2000. View at: Publisher Site  Google Scholar
 R. G. Parr, P. W. Ayers, and R. F. Nalewajski, The Journal of Physical Chemistry A, vol. 109, no. 17, pp. 3957–3959, 2005. View at: Publisher Site
 F. de Proft, C. van Alsenoy, A. Peeters, W. Langenaeker, and P. Geerlings, “Atomic charges, dipole moments, and Fukui functions using the Hirshfeld partitioning of the electron density,” Journal of Computational Chemistry, vol. 23, no. 12, pp. 1198–1209, 2002. View at: Publisher Site  Google Scholar
 P. W. Ayers, “Information theory, the shape function, and the hirshfeld Atom,” Theoretical Chemistry Accounts, vol. 115, no. 5, pp. 370–378, 2006. View at: Publisher Site  Google Scholar
 R. K. Roy, S. Pal, and K. Hirao, “On nonnegativity of Fukui function indices,” The Journal of Chemical Physics, vol. 110, no. 17, article 8236, 1999. View at: Publisher Site  Google Scholar
 F. de Proft, R. VivasReyes, A. Peeters, C. van Alsenoy, and P. Geerlings, “Hirshfeld partitioning of the electron density: atomic dipoles and their relation with functional group properties,” Journal of Computational Chemistry, vol. 24, no. 4, pp. 463–470, 2003. View at: Publisher Site  Google Scholar
 M. Randić, “Design of molecules with desired properties: molecular similarity approach to property optimization,” in Concept and Applications of Molecular Similarity, M. A. Johnson and G. M. Maggiora, Eds., pp. 77–145, John Wiley & Sons, New York, NY, USA, 1990. View at: Google Scholar
 R. Carbó, L. Leyda, and M. Arnau, “How similar is a molecule to another? An electron density measure of similarity between two molecular structures,” International Journal of Quantum Chemistry, vol. 17, no. 6, pp. 1185–1189, 1980. View at: Publisher Site  Google Scholar
 R. Carbó and B. Calabuig, “Molsimil88: molecular similarity calculations using a CNDOlike approximation,” Computer Physics Communications, vol. 55, no. 1, pp. 117–126, 1989. View at: Publisher Site  Google Scholar
 R. Carbó and B. Calabuig, “Molecular quantum similarity measures and Ndimensional representation of quantum objects. I. Theoretical foundations,” International Journal of Quantum Chemistry, vol. 42, no. 6, pp. 1681–1693, 1992. View at: Publisher Site  Google Scholar
 R. Carbó, B. Calabuig, L. Vera, and E. Besalú, “Molecular quantum similarity: theoretical framework, ordering principles, and visualization techniques,” Advances in Quantum Chemistry, vol. 25, pp. 253–313, 1994. View at: Publisher Site  Google Scholar
 M. Sola, J. Mestres, R. Carbó, and M. Duran, “Use of ab initio quantum molecular similarities as an interpretative tool for the study of chemical reactions,” Journal of the American Chemical Society, vol. 116, no. 13, pp. 5909–5915, 1994. View at: Publisher Site  Google Scholar
 G. Boon, C. van Alsenoy, F. de Proft, P. Bultinck, and P. Geerlings, “Molecular quantum similarity of enantiomers of amino acids: a case study,” Journal of Molecular Structure, vol. 727, no. 1–3, pp. 49–56, 2005. View at: Publisher Site  Google Scholar
 G. Boon, C. van Alsenoy, F. de Proft, P. Bultinck, and P. Geerlings, “Similarity and chirality: quantum chemical study of dissimilarity of enantiomers,” Journal of Physical Chemistry A, vol. 107, no. 50, pp. 11120–11127, 2003. 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 and W. Yang, Density Functional Theory of Atoms and Molecules, Oxford University Press, New York, NY, USA, 1989.
 H. Chermette, “Chemical reactivity indexes in density functional theory,” Journal of Computational Chemistry, vol. 20, no. 1, pp. 129–154, 1999. View at: Google Scholar
 P. Geerlings, F. de Proft, and W. Langenaeker, “Conceptual density functional theory,” Chemical Reviews, vol. 103, no. 5, pp. 1793–1874, 2003. View at: Publisher Site  Google Scholar
 S. B. Liu, “Conceptual density functional theory and some recent developments,” Acta PhysicoChimica Sinica, vol. 25, no. 3, pp. 590–600, 2009. View at: Google Scholar
 P. W. Ayers, J. S. M. Anderson, and L. J. Bartolotti, “Perturbative perspectives on the chemical reaction prediction problem,” International Journal of Quantum Chemistry, vol. 101, no. 5, pp. 520–534, 2005. View at: Publisher Site  Google Scholar
 R. G. Pearson, Chemical Hardness, Applications from Molecules to Solids, John Wiley and Sons, Weinheim, Germany, 1997.
 P. W. Ayers, “The physical basis of the hard/soft acid/base principle,” Faraday Discuss, vol. 135, pp. 161–190, 2007. View at: Publisher Site  Google Scholar
 T. Mineva, E. Sicilia, and N. Russo, “Densityfunctional approach to hardness evaluation and its use in the study of the maximum hardness principle,” Journal of the American Chemical Society, vol. 120, no. 35, pp. 9053–9058, 1998. View at: Publisher Site  Google Scholar
 G. de Luca, E. Sicilia, N. Russo, and T. Mineva, “On the hardness evaluation in solvent for neutral and charged systems,” Journal of the American Chemical Society, vol. 124, no. 7, pp. 1494–1499, 2002. View at: Publisher Site  Google Scholar
 T. Mineva, V. Parvanov, I. Petro, N. Neshev, and Russo, “Fukui indices from perturbed Kohn−Sham orbitals and regional softness from mayer atomic valences,” The Journal of Physical Chemistry A, vol. 105, no. 10, pp. 1959–1967, 2001. View at: Publisher Site  Google Scholar
 R. G. Parr, L. von Szentpaly, 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, pp. 1691–1694, 1982. View at: Publisher Site  Google Scholar
 W. T. Yang, Y. K. Zhang, and P. W. Ayers, “Degenerate ground states and a fractional number of electrons in density and reduced density matrix functional theory,” Physical Review Letters, vol. 84, pp. 5172–5175, 2000. View at: Publisher Site  Google Scholar
 W. T. Yang, R. G. Parr, and R. Pucci, The Journal of Chemical Physics, vol. 81, no. 6, article 2862, 1984. View at: Publisher Site
 P. W. Ayers, W. T. Yang, and L. J. Bartolotti, “Fukui function,” in Chemical Reactivity Theory: A Density Functional View, P. K. Chattaraj, Ed., pp. 255–267, CRC Press, Boca Raton, Fla, USA, 2009. View at: 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
 A. MoralesBayuelo, H. Ayazo, and R. VivasReyes, “Threedimensional quantitative structureactivity relationship CoMSIA/CoMFA and LeapFrog studies on novel series of bicyclo [4.1.0] heptanes derivatives as melaninconcentrating hormone receptor R1 antagonists,” European Journal of Medicinal Chemistry, vol. 45, no. 10, pp. 4509–4522, 2010. View at: Publisher Site  Google Scholar
 X. Gironés, D. Robert, and R. CarbóDorca, “TGSA: a molecular superposition program based on topogeometrical considerations,” Journal of Computational Chemistry, vol. 22, no. 2, pp. 255–263, 2001. View at: Google Scholar
 A. MoralesBayuelo, J. Torres, and R. VivasReyes, “Hückel treatment of pyrrole and pentalene as a function of cyclopentadienyl using local quantum similarity index (LQSI) and the topogeometrical superposition approach (TGSA),” Journal of Theoretical and Computational Chemistry, vol. 11, no. 1, pp. 223–239, 2012. View at: Publisher Site  Google Scholar
 A. MoralesBayuelo, R. Baldiris, J. Torres, J. E. Torres, and R. VivasReyes, “Theoretical study of the chemical reactivity and molecular quantum similarity in a series of derivatives of 2adamantylthiazolidine4one using density functional theory and the topogeometrical superposition approach,” International Journal of Quantum Chemistry, vol. 112, no. 14, pp. 2681–2687, 2012. View at: Publisher Site  Google Scholar
 A. MoralesBayuelo, J. Torres, and R. VivasReyes, “Quantum molecular similarity analysis and quantitative definition of catecholamines with respect to biogenic monoamines associated: scale alpha and beta of quantitative convergence,” International Journal of Quantum Chemistry, vol. 112, no. 14, pp. 2637–2642, 2012. View at: Publisher Site  Google Scholar
 A. MoralesBayuelo and R. VivasReyes, “Topological model to quantify the global reactivity indexes as local in Diels—Alder reactions, using density function theory (DFT) and local quantum similarity (LQS),” Journal of Mathematical Chemistry, vol. 51, no. 1, pp. 125–143, 2013. View at: Publisher Site  Google Scholar
 Amsterdam Density Functional (ADF) Code, Vrije Universiteit: Amsterdam, The Netherlands, 2007.
 D. Young, Computational Chemistry, John Wiley & Sons, New York, NY, USA, 2001.
 D. C. Langreth and J. P. Perdew, “Theory of nonuniform electronic systems. I. Analysis of the gradient approximation and a generalization that works,” Physical Review B, vol. 21, no. 12, pp. 5469–5493, 1980. View at: Publisher Site  Google Scholar
 D. C. Langreth and J. P. Perdew, “Theory of nonuniform electronic systems. I. Analysis of the gradient approximation and a generalization that works,” Physical Review B, vol. 21, no. 12, pp. 5469–5493, 1980. View at: Publisher Site  Google Scholar
 J. P. Perdew and Y. Wang, “Accurate and simple density functional for the electronic exchange energy: generalized gradient approximation,” Physical Review B, vol. 33, no. 12, pp. 8800–8802, 1986. View at: Publisher Site  Google Scholar
 J. P. Perdew, “Densityfunctional approximation for the correlation energy of the inhomogeneous electron gas,” Physical Review B, vol. 33, no. 12, pp. 8822–8824, 1986. View at: Publisher Site  Google Scholar
 J. P. Perdew and Y. Wang, “Accurate and simple analytic representation of the electrongas correlation energy,” Physical Review B, vol. 45, no. 23, pp. 13244–13249, 1992. View at: Publisher Site  Google Scholar
 J. P. Perdew, K. Burke, and M. Ernzerhof, “Generalized gradient approximation made simple,” Physical Review Letters, vol. 77, no. 18, pp. 3865–3868, 1996. View at: Publisher Site  Google Scholar
 J. P. Perdew, J. A. Chevary, S. H. Vosko et al., “Atoms, molecules, solids, and surfaces: applications of the generalized gradient approximation for exchange and correlation,” Physical Review B, vol. 46, no. 11, pp. 6671–6687, 1992. View at: Publisher Site  Google Scholar
 A. R. Katritzky and R. D. Topsom, “The σ and πinductive effects,” Journal of Chemical Education, vol. 48, no. 7, pp. 427–431, 1971. View at: Google Scholar
 S. B. Liu, Chemical Reactivity Theory: A Density Functional View, Taylor and Francis, Boca Raton, Fla, USA, 2009.
 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. K. Chattaraj, B. Maiti, and U. Sarkar, “Philicity: a unified treatment of chemical reactivity and selectivity,” Journal of Physical Chemistry A, vol. 107, no. 25, pp. 4973–4975, 2003. View at: Publisher Site  Google Scholar
 R. K. Roy, S. Krishnamurti, P. Geerlings, and S. Pal, “Local softness and hardness based reactivity descriptors for predicting intra and intermolecular reactivity sequences: carbonyl compounds,” The Journal of Physical Chemistry A, vol. 102, no. 21, pp. 3746–3755, 1998. View at: Publisher Site  Google Scholar
 J. L. Moncada and G. S. Morán, “Caracterización de la reactividad intrínseca de los halobencenos en el modelo conceptual de la teoría de funcionales de la densidad (TFD),” Anales de Química, vol. 103, no. 3, pp. 23–27, 2007. View at: Google Scholar
 J. McMurry, Organic Chemistry, Brooks/Cole Publishing, Albany, NY, USA, 4th edition, 1996.
 S. I. Neils, Physical Organic Chemistry, Longman Scientific and Technical Singapore Publishers, Singapore, 2nd edition, 1995.
 B. Tchoubar, Mecanismos de Reacción en Química Orgánica, Limusa, México City, Mexico, 1965.
Copyright
Copyright © 2014 Alejandro MoralesBayuelo and Ricardo VivasReyes. 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.