Theoretical and computational study of benzenium and toluenium isomers

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 10¹⁴ (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

FIGURE 1. Benzene and benzenium.

FIGURE 2

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

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

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

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

FIGURE 3. SCF-optimized bond lengths of benzene and benzenium (Å).

FIGURE 4

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

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

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

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

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

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

TABLE 8. Overlap populations for benzene and benzenium.

TABLE 9

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 CH₂ 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 d2FsdGVyLmVybWxlckB1dHNhLmVkdQ==.

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

CrossRef Full Text | Google Scholar

Bader, R. F. W. (1985). Atoms in molecules. Acc. Chem. Res. 18, 9–15. doi:10.1021/ar00109a003

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Fox, M. A., and Whitesell, J. K. (2004). Organic chemistry. 3rd ed. Sudbury, MA, Burlington, MA: Janes and Bartlett Publishers.

Google Scholar

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.

Google Scholar

Glendening, E., Clark, L., and Weinhold, F. (2021). Natural bond orbital 7.0 homepage. Available at: https://nbo6.chem.wisc.edu/index.htm.

Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Mulliken, R. S., and Ermler, W. C. (1977). Diatomic molecules: results of ab initio calculations. Amsterdam: Elsevier. (reprinted 2012), Chap. II.

Google Scholar

Mulliken, R. S., and Ermler, W. C. (1981). Polyatomic molecules: results of ab initio calculations. Amsterdam: Elsevier. (reprinted 2012), Chap. XI.

Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Reed, A. B., Weinstock, R. B., and Weinhold, F. (1985). Natural population analysis. J. Chem. Phys. 83, 735–746. doi:10.1063/1.449486

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Rieke, C. A., Brown, W. G., and Mulliken, R. S. (1941). Hyperconjugation. J. Am. Chem. Soc. 63, 41–56. doi:10.1021/ja01846a008

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Zefirov, N. S., and Schechtman, N. M. (1971). Anomeric effect. Usp. Khim 40, 593–624.

CrossRef Full Text | Google Scholar