Theoretical and computational study of benzenium and toluenium isomers
Falonne C. Moumbogno Tchodimo
Walter C. Ermler- Department of Chemistry, The University of Texas at San Antonio, San Antonio, TX, United States
Four methods of computational quantum chemistry are used in a study of hyperconjugation in protonated aromatic molecules. Benzene, benzenium, toluene, and four isomeric forms of toluenium are examined using the self-consistent field level of theory followed by configuration interaction and coupled cluster calculations, as well as density functional theory. Results for proton affinities, geometric parameters, atomic populations, dipole moments, and polarizabilities are reported. The calculated results are in good agreement with previous computational studies and with experimental data. The presence of hyperconjugation is evident from the shortened carbon–carbon bond lengths in the aromatic ring and concomitant changes in dipole moments and polarizabilities. The proton affinities of benzene and toluene compare well with experimental values. The examination of all of the toluenium isomers reveals that the position of the methyl group has a minor impact on the strength of hyperconjugation, although the most stable isomer is found to be the para form. Mulliken population analyses indicate that the addition of a proton contributes to aromatic hyperconjugation and increases the strength of π-bonds at the expense of σ-bonds.
1 Introduction
The Baker–Nathan effect was discussed initially in 1935 (Baker and Nathan, 1935). It is rooted in the presence of additional resonance structures in certain molecules and, in general, is described as an electron density delocalization from a σ-bond to a p-orbital or a π-bond. In a series of seminal articles, Mulliken extended and generalized the Baker–Nathan concept and coined the term hyperconjugation (hereafter denoted as HC) (Mulliken, 1939; Rieke et al., 1941). The importance of delocalizing interactions involving σ-bonds has been examined theoretically and computationally (Reed et al., 1988). It has been reported that HC can influence conformational equilibria (Romers et al., 1969; Zefirov and Schechtman, 1971), modify reactivity (Baddeley, 1973; Chang et al., 1987), and determine selectivity (Beckwith and Duggan, 1998). Furthermore, it is has been established that the presence of HC is mainly responsible for the stability of the more highly substituted alkenes (Fox and Whitesell, 2004). The anomeric effect in the reactivity of α-haloglycine esters with various nucleophiles has been attributed to HC, leading to enhanced halogen nucleofugality which facilitates halogen abstraction by hydrogen-bond donor catalysts (Samanta and Stéphane, 2019). For example, the presence of HC has been attributed to improved catalysis in a pyridoxal 5′-phosphate-dependent enzyme (Dajnowicz et al., 2018), as well as being an important factor in stabilizing certain excited, radical, and ionic species (Muller and Mulliken, 1958; Lambert et al., 1987; Raabe et al., 1996). Davies maintained that HC activation in the presence of Sn–C bonds can accelerate a reaction rate by more than a factor of 1014 (Davies, 1999). In organometallic compounds, the hyperconjugative β-proton coupling of an H–C σ-bond in a cyclic radical can be enhanced or decreased according to the symmetry of the π-orbital with which it hyperconjugates according to the so-called Whiffen effect (Davies, 1999). Ermler et al. (1976) carried out accurate ab initio self-consistent field (SCF) calculations to examine the role of HC in the benzenium and p-toluenium cations (Ermler and Mulliken, 1978). A recent study that examines the phenomenon of HC in aromatic systems is consistent with their conclusions (Xiaojuan et al., 2020).
2 Calculations
Benzene and benzenium are shown in Figure 1, and toluene and the four isomeric forms of toulenium are shown in Figure 2. The goal of the present study is to identify and codify the presence of HC in these molecules through a set of high-level calculations of energies, bond lengths, bond angles, electric dipole polarizabilities, and Mulliken populations.
FIGURE 1. Benzene and benzenium.
FIGURE 2. Toluene and meta-, ortho-, para-, and ipso-toluenium.
The calculations reported in this study were accomplished using the Gaussian software suite (Frisch et al., 2016). SCF calculations were carried out followed by single plus double excitation configuration interaction (CISD) and coupled cluster with single and double and perturbative triple excitation (CCSD(T)) calculations. All were conducted using 6-311G** basis sets of contracted Gaussian-type functions (Frisch et al., 2016). In addition, density functional theory (DFT) calculations using B3LYP exchange correlation functionals and 6-311G** basis sets were carried out. At each level of theory, a full geometry optimization was accomplished. The proton affinities of benzene and each toluene isomers were calculated at each level of theory. At each level of theory the proton affinity is defined as the difference between the total energy of the neutral species at its calculated equilibrium geometry less than that of the protonated species at its calculated equilibrium geometry. Finally, dipole moments, defined as expectation values of the dipole moment operator, and dipole polarizabilities, defined as the change in the dipole moment of a molecule when in a static electric field, are reported at all levels of theory (Frisch et al., 2016).
Initially, single-point energy SCF calculations were carried out to determine total energies, molecular orbital (MO) energies, Mulliken populations, and electric dipole polarizabilities. Results for the benzene systems and the toluene systems are compared to those obtained by Ermler et al. (1976) and Ermler and Mulliken (1978), respectively. These preliminary SCF calculations were followed by CISD and CCSD(T) calculations for each of these molecules, and the aforementioned properties were calculated. Additionally, the ortho, meta, and ipso forms of toluenium were also studied at the stated levels of theory. Figure 2 shows the position of the methyl group in the four isomers of protonated toluene. Datasets for the ortho, meta, para, and ipso forms of the toluenium molecules were constructed, and geometry optimizations were completed yielding equilibrium bond lengths and bond angles and the corresponding proton affinities. These structures were then used to calculate electric dipole polarizabilities and populations of each molecule.
3 Results and discussion
Table 1 compares the present results for benzene and benzenium with those of [17]. The results for toluene and toluenium are compared with those of [18] in Table 2.
TABLE 1. Calculated total SCF energy, HOMO energy, and proton affinity for benzene and benzenium compared to the work of Ermler et al. (1976) (a.u.).
TABLE 2. Calculated total SCF energy, HOMO energy, and proton affinity for toluene and p-toluenium compared to the work of Ermler and Mulliken (1978) (a.u.).
The calculated results are in good agreement with those cited in both earlier studies, with differences attributable to the use of the somewhat smaller 6-311G** basis set in the current calculations compared to the σ/π-optimized polarized double zeta basis sets used in the work of Ermler et al. (1976) and Ermler and Mulliken (1978). The main differences in total energies are due to a more rigorous representation of the carbon 1s core electrons by 9 s-type GTOs versus the standard Gaussian 6 s-type GTOs in the 6-31G** basis sets. These core energy differences do not appreciably impact the quality of the MOs that describe the C valence electrons. The calculated highest occupied molecular orbital (HOMO) energy values in Table 1 and Table 2 are in excellent agreement with the work of Ermler et al. (1976) and Ermler and Mulliken (1978). Results for HOMO energies (from natural orbitals) and for proton affinities for the CISD calculations using the same 6-311G** basis set are shown in Table 3. It is noteworthy that the CISD values for the proton affinities for benzenium and toluenium from the present calculations are in close agreement with the work of Ermler et al. (1976) and Ermler and Mulliken (1978) and with the experimentally determined (NIST Chemistry WebBook) values of 178.4 and 187.0 kcal/mol, respectively. The calculated dipole moment of toluene agrees well with the experimentally determined (NIST Chemistry WebBook) value of 0.36 D.
TABLE 3. CISD and CCSD(T) results for benzene and benzenium.
Geometry optimization calculations for each of the molecules were carried out at each of the levels of theory. As expected, in the presence of HC, the bond lengths of certain carbon–carbon bonds in the aromatic ring are expected to be impacted. Figure 3 shows the calculated bond lengths for benzene and benzenium, and Figure 4 shows those for toluene and toluenium.
FIGURE 3. SCF-optimized bond lengths of benzene and benzenium (Å).
FIGURE 4. SCF-optimized bond lengths of toluene and p-toluenium (Å).
The optimum benzene geometric parameters are in close agreement with the experiment. In the presence of HC, the bond lengths of the C1–C6 and C3–C4 bonds are shortened whereas C1–C2, C5–C6, C5–C4, and C3–C2 are lengthened. As seen in Figure 4, which shows the bond lengths in benzenium, all of the carbon–carbon bonds are near 1.39 Å, somewhat shorter than the values of the standard value for an aromatic carbon–carbon bond length (1.47 Å). Additionally, all of the carbon–carbon bond lengths lie near the standard aromatic bond length of 1.39 Å in the toluene molecule as shown in Figure 4, indicative of the presence of HC.
. In addition, the aliphatic carbon–carbon bond length (C7–C1, of 1.511 Å in toluene, shorter that the standard value of 1.54 Å) indicates the presence of HC. Furthermore, the C7–C1 bond length in toluene is calculated as 1.511 Å, but shortens to 1.489 Å in toluenium. Energies of the toluenium isomers and their respective proton affinities are reported in Table 4.
TABLE 4. Total energies, HOMOs, LUMOs, and proton affinities for toluene and the toluenium isomers.
Calculations of electric dipole polarizabilities of benzene, benzenium, toluene, and all toluenium isomers at the optimized geometries were carried out at different levels of theory [SCF, DFT (B3LYP), CISD, and CCSD(T)]. As is shown, the SCF polarizability of benzene is calculated to be 58.45 a.u. compared to that of benzenium, which is smaller by 1.04 a.u. Also shown are the electric dipole polarizabilities calculated at the CISD level of theory, with the benzenium value smaller by 0.8 a.u. The electric dipole polarizability of toluene is 70.36 a.u., compared to that of toluenium, which is smaller by 0.48 a.u. As shown, the calculated values for both benzene and toluenium are all lower than the experimentally determined values (NIST Chemistry WebBook). Furthermore, the differences in polarizabilities between the neutral and protonated molecules indicate a tightening of the electron clouds for the latter, with the effect more pronounced for toluene.
The different structural isomers of toluenium were examined following geometry optimization calculations on each of the isomers. Examination of the methyl group at different positions (ortho, meta, para, or ipso Figure 2) is important because in determining how the position of the methyl group affects HC. The results are presented in Table 5 and Table 6.
TABLE 5. Average bond lengths (Å), polarizabilities (a.u.), and HCH angles of benzene and benzenium at different levels of theory. Experimental values are shown in parentheses.
TABLE 6. Average bond lengths (Å), polarizabilities (a.u.), and HCH angles of toluene and toluenium isomers at different levels of theory. Experimental values are shown in parentheses.
As can be seen in Table 4, the protonation energies are relatively close to one another. That is, there are no major differences between the total energies, HOMOs, and LUMOs of o-, m-, and p-toluenium. This indicates that the methyl group position in the toluenium ion has a small impact on the strength of HC.
Figure 5 shows the bond lengths in the isomers of toluenium. The bond lengths in each molecule are similar to one other with the p-toluenium showing slightly shorter bond lengths, indicating that the most stable isomer is the para form.
FIGURE 5. SCF-optimized bond lengths of o-, m-, and p-toluenium (Å).
The mass spectrometric study of Kuck et al. (1985) on benzenium and toluenium syntheses revealed all of the structures shown in Figure 1 and Figure 2. Their results are consistent with the present calculated ones. Table 7 shows the Mulliken population analyses for benzene, benzenium, toluenium, and p-toluenium. Table 8 reports overlap populations for benzene, and Table 9 shows overlap populations for toluene and the four toluenium isomers.
TABLE 7. Comparison of Mulliken populations analyses as net charges for benzene, benzenium, toluene, and p-toluenium at the SCF and CISD levels of theory.
TABLE 8. Overlap populations for benzene and benzenium.
TABLE 9. Overlap populations for toluene and toluenium isomers.
The data revealed that, as the proton is added to the benzene ring, the overlap population decreases in that bond. For instance, all of the bond lengths in benzene are equal with overlap populations consistently 0.514–0.515. These are in contrast to the values of 0.389–0.391 in benzenium. This indicates that the addition of a proton contributes to the aromatic conjugation. In addition, HC affects the overlap population by increasing the strength of the π-bond and decreasing that of the σ-bond. It is well known that HC is present in benzenium and absent in benzene. With that fact, looking at Table 9, it is seen that the π-bonds in benzene are stronger compared to the σ-bonds in benzenium. That is, the C1–C2 π-bond in benzene is 0.514, in contrast to the σ-bond in benzenium, which is 0.391. While it is expected that σ-bonds should be stronger than π-bonds, due to the presence of HC, this is reversed as seen in benzene and benzenium. The same trend is seen for toluene and toluenium isomers in terms of overlap populations (the results in Table 9 are defined with reference to Figure 2 with the exception that in toluene, C1 is C3 and C4 is C6.). In this instance, it is seen that the presence of HC bonding range is estimated to be 0.370–0.398.
As Mulliken and Ermler noted out in their analyses of other ab intio studies of six-membered ring molecules, all indicate that the structures have the proton attached to one carbon atom to form a CH2 group perpendicular to the plane of the ring (Mulliken and Ermler, 1981). These earlier studies, all at the SCF level of theory, consistently predict proton affinities near 190 kcal/mol, within 10 kcal/mol of experimental values. The present correlated calculations are consistent with these SCF values, in keeping with the fact that all of the neutral and protonated molecules have the same numbers of electrons and the expectation that, in general, electron correlation does strongly impact proton affinities (Chandra and Goursot, 1996).
The Mulliken population analysis procedure was chosen to further analyze all of the SCF, CISD, CCSD(T), and DFT computational results beyond the energetics and other properties. Although it is implemented in Gaussian (Frisch et al., 2016) and nearly all currently available quantum chemistry software suites, Mulliken’s method has well-known shortcomings, especially revealed in his so-called one-half overlap assignment (Mulliken and Ermler, 1977). As was carried out in our earlier calculations on benzenium (Ermler et al., 1976) and p-toluenium (Ermler and Mulliken, 1978), this work compares the results of population analyses described with the same GTO basis set and within the same level of theory, thereby providing a means by which changes in gross and overlap populations from one molecule or isomer to another can be compared in the same context.
The Mulliken population analysis method has been revisited by a host of researchers who have refined it to address its embedded constraints. Noteworthy are the development of natural bond orbitals (NBOs) by (Reed et al. (1985); Reed et al. (1988)), the atoms in molecules (AIMs) approach of Politzer et al. (1970) and Bader (1985), Roby’s method based on electron density operators and projection operators (Roby, 1974), and the one based on electron localization functions (ELFs) (Silvi and Savin, 1994). In terms of the concept of hyperconjugation, the NBO method is especially compelling and can now be implemented in the context of existing software suites (Glendening et al., 2021).
There are also a number of the methods that do not deal with populations, but interpret electronic structure calculations through analyses of the molecular energy and total wave function. Two examples are the energy decomposition analysis (EDA) method (Mo et al., 2000) and a definition of aromaticity indices (Jug et al., 1991). While the analysis of molecular electronic structure calculations reported herein is based on geometric changes, protonation energies, polarizabilities, and Mulliken populations, future work that incorporates one or more of the aforementioned approaches can possibly provide additional insight into the nature of aromatic hyperconjugation.
4 Conclusion
This investigation into the effects of aromatic hyperconjugation in benzene, benzenium, toluene, and its isomers through high-level quantum chemical calculations shows good agreement with the available experimental data. The presence of HC in these molecules is evident from the observed shortening of carbon–carbon bond lengths in the aromatic ring and the diminishing of polarizabilities in the protonated forms. This phenomenon contributes to the stabilization of these molecules and can play a crucial role in determining their structural and electronic properties. The proton affinities of benzene and toluene calculated in this study are in close agreement with experimental values, further supporting the accuracy of the calculations and the role of HC in these systems.
Additionally, the investigation of all of the toluenium isomers reveals that the position of the methyl group has a minor impact on overall strength of the HC. However, the most stable isomer was found to be the para form, consistent with the influence of HC on the molecular stability and structure. The analysis of Mulliken population data further substantiates the presence of HC in the molecules, showing that the addition of a proton contributes to aromatic conjugation and strengthens π-bonds while weakening σ-bonds. This emphasizes the important role of HC in modifying electronic distributions and bonding patterns within these aromatic systems.
Overall, this investigation shed some light on the significance of hyperconjugation in aromatic molecules, highlighting its influence on structural stability, reactivity, and electronic properties. The implications of these results can be extended to various areas, including catalysis, reaction mechanisms, and the stability of aromatic species.
Data availability statement
The datasets presented in this article are not readily available. Requests to access the datasets should be directed to walter.ermler@utsa.edu.
Author contributions
FT and WE: writing – original draft.
Funding
This research was supported in part by the National Science Foundation under grant CHE-1853477. This work received computational resources from the University of Texas at San Antonio’s HPC cluster operated by the Office of Information Technology and the Texas Advanced Computing Center (TACC) at the University of Texas at Austin.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Baddeley, G. (1973). Delocalisation into antibond orbitals. Tetrahedron Lett. 14, 1645–1648. doi:10.1016/s0040-4039(01)96018-6
Baker, J. W., and Nathan, W. S. (1935). The mechanism of aromatic side-chain reactions with special reference to the polar effects of substituents. J. Chem. Soc., 1844–1847. doi:10.1039/jr9350001844
Beckwith, A. L. J., and Duggan, P. J. (1998). The quasi-HOMO-anomeric interaction in substituted tetrahydropyranyl radicals: diastereoselectivity. Tetrahedron 54, 6919–6928. doi:10.1016/s0040-4020(98)00373-1
Chandra, A. K., and Goursot, A. (1996). Calculation of proton affinities using density functional procedures: a critical study. J. Phys. Chem. 100, 11596–11599. doi:10.1021/jp9603750
Chang, J. W. A., Taira, K., Urano, S., and Sorenstein, D. G. (1987). Stereoelectronic effects on the nucleophilic addition of phosphite to the carbonyl double bond. ab initio molecular orbital calculations on reaction surfaces and the σ-effect. Tetrahedron 43, 3863–3874. doi:10.1016/s0040-4020(01)81668-9
Dajnowicz, S., Parks, J. M. P., Hu, X., Johnston, R. C., Kovalesky, A. Y., and Mueser, T. C. (2018). Hyperconjugation promotes catalysis in a pyridoxal 5′-phosphate-dependent enzyme. ACS Catal. 8 (7), 6733–6737. doi:10.1021/acscatal.8b01911
Davies, A. G. (1999). The Whiffen effect of symmetry-enhanced and symmetry-forbidden hyperconjugation in spin-paired molecules. J. Chem. Soc. Perkin Trans., 2461–2467. doi:10.1039/a905970i
Ermler, W. C., and Mulliken, R. S. (1978). Ab initio SCF computations on toluene and the toluenium ion using a large contracted Gaussian basis set. J. Am. Chem. Soc. 100, 1647–1653. doi:10.1021/ja00474a001
Ermler, W. C., Mulliken, R. S., and Clementi, E. (1976). Ab initio SCF computations on benzene and the benzenium ion using a large contracted Gaussian basis set. J. Am. Chem. Soc. 98, 388–394. doi:10.1021/ja00418a011
Fox, M. A., and Whitesell, J. K. (2004). Organic chemistry. 3rd ed. Sudbury, MA, Burlington, MA: Janes and Bartlett Publishers.
Frisch, M. J., Trucks, G. W., Schlegel, H. B., Scuseria, G. E., Robb, M. A., Cheeseman, J. R., et al. (2016). Gaussian 16, revision C.01. Wallingford, CT: Gaussian, Inc.
Glendening, E., Clark, L., and Weinhold, F. (2021). Natural bond orbital 7.0 homepage. Available at: https://nbo6.chem.wisc.edu/index.htm.
Jug, K., and Köster, A. M. (1991). Aromaticity as a multi-dimensional phenomenon. J. Phys. Org. Chem. 4, 163–169. doi:10.1002/poc.610040307
Kuck, D., Schneider, J., and Grutzmacher, H.-F. (1985). A study of gaseous benzenium and toluenium ions generated from 1,4- dihydro- and 1-methyl-1,4-dihydro-benzoic acids. J. Chem. Soc. Perkin Trans. II, 689–696. doi:10.1039/p29850000689
Lambert, J. B., Wang, G.-T., Finzel, R. B., and Teramura, D. H. (1987). Stabilization of positive charge by.beta.-silicon. J. Am. Chem. Soc. 109, 7838–7845. doi:10.1021/ja00259a036
Mo, Y., Gao, J., and Peyerimhoff, S. D. (2000). Energy decomposition analysis of intermolecular interactions using a block-localized wave function approach. J. Chem. Phys. 112, 5530–5538. doi:10.1063/1.481185
Muller, N., and Mulliken, R. S. (1958). Strong or isovalent hyperconjugation in some alkyl radicals and their positive ions. J. Am. Chem. Soc. 80, 3489–3497. doi:10.1021/ja01547a001
Mulliken, R. S., and Ermler, W. C. (1977). Diatomic molecules: results of ab initio calculations. Amsterdam: Elsevier. (reprinted 2012), Chap. II.
Mulliken, R. S., and Ermler, W. C. (1981). Polyatomic molecules: results of ab initio calculations. Amsterdam: Elsevier. (reprinted 2012), Chap. XI.
Mulliken, R. S. (1939). Intensities of electronic transitions in molecular spectra. IV. Cyclic dienes and hyperconjugation. J. Chem. Phys. 7, 339–352. doi:10.1063/1.1750446
Politzer, P., and Harris, R. R. (1970). Properties of atoms in molecules. I. Proposed definition of the charge on an atom in a molecule. J. Am. Chem.Soc. 92, 6451–6454. doi:10.1021/ja00725a011
Raabe, G., Gais, H. J., and Fleischhauer, J. (1996). Ab initio study of the effect of fluorination upon the structure and configurational stability of σ-sulfonyl carbanions: the role of negative hyperconjugation. J. Am. Chem. Soc. 11, 4622–4630. doi:10.1021/ja953034t
Reed, A. B., Weinstock, R. B., and Weinhold, F. (1985). Natural population analysis. J. Chem. Phys. 83, 735–746. doi:10.1063/1.449486
Reed, A. E., Curtiss, L. A., and Weinhold, F. (1988). Intermolecular interactions from a natural bond orbital, donor–acceptor viewpoint. Chem. Rev. 88, 899–926. doi:10.1021/cr00088a005
Rieke, C. A., Brown, W. G., and Mulliken, R. S. (1941). Hyperconjugation. J. Am. Chem. Soc. 63, 41–56. doi:10.1021/ja01846a008
Roby, K. R. (1974). Quantum theory of chemical valence concepts I. Definition of the charge on an atom in a molecule and of occupation numbers for electron density shared between atoms. Mol. Phys. 27, 81–104. doi:10.1080/00268977400100071
Romers, C., Altona, C., Buys, H. R., and Havinga, E. (1969). Geometry and conformational properties of some five- and six membered heterocyclic compounds containing oxygen or sulfur. Top. Stereochem 4, 39–97. doi:10.1002/9780470147139.ch2
Samanta, S. S., and Stéphane, P. R. (2019). Synthesis and reactivity of α-haloglycine esters: hyperconjugation in action. Eur. J. Org. Chem. 2019 (39), 6597–6605. doi:10.1002/ejoc.201901033
Silvi, B., and Savin, A. (1994). Classification of chemical bonds based on topological analysis of electron localization functions. Nature 371, 683–686. doi:10.1038/371683a0
Xiaojuan, N., Yinong, Z., Gurjyot, S., and Liu, F. (2020). π-Orbital Yin–Yang Kagome bands in anilato-based metal–organic frameworks. Phys. Chem. Chem. Phys. 44, 25827–25832. doi:10.1039/D0CP03941A
Keywords: benzenium, toluenium, proton affinity, hyperconjugation, DFT
Citation: Tchodimo FCM and Ermler WC (2023) Theoretical and computational study of benzenium and toluenium isomers. Front. Chem. 11:1253599. doi: 10.3389/fchem.2023.1253599
Received: 05 July 2023; Accepted: 19 September 2023;
Published: 01 November 2023.
Edited by:
Sam P. De Visser, The University of Manchester, United KingdomReviewed by:
Slavko Radenković, University of Kragujevac, SerbiaGrant Hill, The University of Sheffield, United Kingdom
Copyright © 2023 Tchodimo and Ermler. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Walter C. Ermler, Walter.Ermler@utsa.edu