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 Topo-Geometrical 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 Slater-type-orbital 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 3HC–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 [59].

Several concepts and theories have been proposed in the literature in order to explain the molecular chemical bond [1019]. 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 [2426].

2. Theory and Computational Details

2.1. Similarity Indexes

Reviews about concepts and applications in quantum similarity can be found elsewhere (see, for example, [59, 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 self-similarity 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 self-similarity and calculates only the measures of “shape similarity,” as an other alternative is the of Hodgkim-Richards 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 [3238]. 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 [4046]. This idea has been applied in this study to CH3–F taken as molecule A and CH3–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 self-similarity.

Similar to the contributions of the Euclidean distances Equation (17) has a convergence range [0,), where 0 represents the self-similarity. This means that small distances of CH3–F with respect to CH3–Cl and CH3–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 [4754]. 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, 5557], 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 [5964]:

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., three-dimensional 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 Topo-Geometrical 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 [6770], 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 first-principle electronic structure calculations that make use of (DFT) [50]. All the molecular structures were fully optimized using the gradient generalized approximation (GGA) [7376] local exchange correlations PW91 [7779]. The geometry optimizations were calculated using a standard Slater-type-orbital basis set with triple-quality plus double-polarization 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 14 the (MQS) matrices for the compounds analyzed are shown, taken as an example.

In Tables 1 and 2, is possible to see that the CH3–Cl and CH3–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 sp3 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 CH3–Cl and CH3–Br with respect to CH3–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 CH3–F with the CH3–Cl and 103.18 pm when compared with the CH3–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 (CH3–Cl > CH3–Br).

On the other hand, the predicted values for the Euclidean distance ((b) and (d)) are closest to the reference value (134 pm) of CH3–F with respect to CH3–Cl and CH3–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 [8183].

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 CH3–Cl is the most stable compound with chemical potential of  eV and the most reactive compound is the CH3–F with  eV. While the CH3–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 CH3–F with  eV and  eV to CH3–Br. On the other hand, such trend is also observed with the values of ionization potential of the CH3–Br with a value of eV with the most low energy to remove electrons from the electron cloud and the electron affinity  eV.

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 CH3–F, : 2.10 eV to CH3–Cl, and : 2.02 eV to CH3–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.

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 CH3–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 CH3–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.

Figure 6 shows the nature of the polar interactions in the C–H bonds of CH3–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 CH3–F that only has a single reactive site f+: 0.38.

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 CH3–Br are shown (see Table 8). This value for the Fukui function f+ shows a lower inductive effect with respect to CH3–Cl and to the CH3–F, respectively, according to the electronegativity values.

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 CH3–F have (: 0.02) in the fluoride, to the CH3–Cl have (: 0.88) in the chloride, and to the CH3–Br have (: 0.99) to the Bromide (see Table 5).

These values are quantified by the electrophilic character [8183] 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 CH3–F has the higher inductive effect; therefore, we can see the trend (CH3–F > CH3–Cl > CH3–Br) according to the trends of correlation reported in Figure 3, and these theoretical values are in agreement with experimental results [8587].

From Figures 57, 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 (CH3–F, CH3–Cl, CH3–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 electron-withdrawing (−I). Although this does not imply donation formal of a unit charge, considering that by convection is assumed as electron-withdrawing 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 C3H8I with C3H8F (0.7318), and the lowest values comparing C3H8Br with C3H8Cl (0.9887), forming the overlap quantum dissimilarity trend with respect to C3H8 (C3H8F (0.7427) > C3H8Cl (0.8895) > C3H8Br (0.9025) > C3H8I (0.9252)) in agreement with the trends shown in Figure 3.

In Table 10 we found the higher Euclidean distance of overlap between C3H8I and C3H8F (1.5732) and the lowest is between C3H8Br and C3H8Cl (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 C3H8F and C3H8 (0.9468). However, the lowest is between C3H8Br and C3H8Cl (0.9990). Also the electronic quantum dissimilarity trend with respect to C3H8 C3H8F (0.9468) > C3H8Cl (0.9631) > C3H8Br (0.9643) > C3H8I (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 C3H8F and C3H8 (5.4638), and the lowest is between C3H8Br and C3H8Cl (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 CH3–F, CH3–Cl, and CH3–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 Topo-Geometrical 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 (CH3–F, CH3–Cl, and CH3–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 CH3–F, with the greater hardness (: 5.97 eV) and the most stable the CH3–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 electron-release (+I) and the topological Fukui (f+) with the electron-withdrawing (−I).

Moreover, through this model we can be demonstrated the origin of the polar interactions of the C–H bond (see Figures 57) that can be applied to a much wider range of compounds such as CH3–NH2, CH3–(CH2)n–NH2, CH3–OH, CH3–(CH2)n–OH, CH3–SH, and CH3–(CH2)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 Morales-Bayuelo thanks the Universidad Nacional Andres Bello (Santiago-Chile) for a Ph.D. fellowship (CONICYT (63100003)). Ricardo Vivas-Reyes is indebted to the Universidad de Cartagena (Cartagena, Colombia), for continuous support to the research group.