Guest Editor: Y. Yamakoshi Beilstein J. Org. Chem.2024,20, 3144–3150.https://doi.org/10.3762/bjoc.20.260 Received 15 Aug 2024,
Accepted 19 Nov 2024,
Published 02 Dec 2024
Quantitative assessment of the first acidity constant (pKa) for BFC (27.6 in CH3CN) and FIC (27.8 in CH3CN) shows the methylene protons to be significantly more acidic than those in related cyclopentadiene (32 in CH3CN), indene (34 in CH3CN), or fluorene (37 in CH3CN) and comparable to the methine of 9-perfluorophenylfluorene (28.14 in CH3CN). This work reports quantitative pKa values of BFC and FIC, places those values in a broadened context of CpH-cognate hydrocarbon acidity and presents a Clar–Loschmidt graph perspective to help understand the “surprises”.
A classic textbook tetrad linking hydrocarbon acidity to aromatic stabilization energy comprises cyclopentadiene (CpH), indene (InH), fluorene (FlH), and diphenylmethane (DPMH) [1,2], with pKa values in DMSO equal to 18 [3], 20.1 [3], 22.6 [3], and 32.2 [4], respectively (Scheme 1) [5,6]. The reaction enthalpy for deprotonation of CpH is ca 20 kcal/mol less endothermic than DPMH and ca. 24–27 kcal/mol less endothermic than 1,4-pentadiene (PDH; pKa ≈ 35 in DMSO) [3], values strikingly similar to the resonance stabilization energy estimated for benzene [7,8]. Furthermore, the trend of pKa values for CpH, InH, and FlH correlates with the reduction of the aromatic stabilization energy for the anion across the series [1,2]. At first glance, this model supports the notion that the relative pKa values of cyclopentadienes embedded in polynuclear aromatic hydrocarbons (CpH-PAHs) reflect some general measure of the embedded CpH's aromatic character. In this context, the observation that BFC and FIC (FlH cognates) manifest acidity comparable to CpH, was surprising [9]. This study uses a Clar–Loschmidt graph model to show that the acidity of CpH-PAH is more nuanced [10,11]. Subsequently, literature and computational examples reveal a leitmotif for CpH-PAH-based redox-active carbon-rich materials.
Results and Discussion
From competitive titrations against 9-C6F5-fluorene (pKa = 28.14 in CH3CN) [12,13] and (4-Me-C6F4)(C6H5)CHCN (pKa = 26.98 in CH3CN) [12,13] as upper and lower references, respectively, the acidities of BFC (pKa = 27.6 in CH3CN) and FIC (pKa = 27.8 in CH3CN) could be quantified and ordered relative to a series of organic aromatic acids (Table 1). Rather than falling among the series cyclopentadiene (pKa = 32 in CH3CN), indene (pKa = 34 in CH3CN), or fluorene (pKa = 37 in CH3CN), BFC and FIC rank as stronger acids than CpH in acetonitrile. Extrapolation to DMSO by conversion equations, supports the assertion that BFC and FIC are more acidic than CpH in an “absolute acidity” sense [13].
Table 1:
Relative aromatic hydrocarbon acidities I.
Compound
CH3CN
CH3CN
DMSO
DMSO
Cmptla
Exptl/Estb
Cmptla
Exptl/Estb
fluorene
35.0
37b,c
21.6
22.6d
indene
32.2
34b,c
18.8
20.1d
cyclopentadiene
29.1
32b,c
15.6
18.0d
FIC
27.8
27.8e
14.4
14.4b,e
BFC
28.3
27.6e
14.3
14.3b,e
9-phenylfluorene
30.8
32b,c
17.4
17.9f
fluoradene
22.8
23.9b,c
9.4
10.5f
diphenylmethane
43.4
48b,c
29.9
32.2g
aB97-D/def2-TZVPP(solvent)//B97-D/def2-TZVPP(solvent) referenced to 9-C6F5-fluorene (solvent). bEstimated using experimental data and correlation equations. cRef. [13]. dRef. [3]. eThis work. fRef. [5]. gRef. [4].
How odd is such a pKa value of ca. 14 for a CpH-PAH like BFC or FIC in DMSO? From the initial structural analogy to fluorene there is a substantial 8 pKa unit difference, which is nearly 11 kcal/mol in ambient reaction free energy. Computations (B97-D/def2-TZVPP(THF)//B97-D/def2-TZVPP(THF)) also predict pKa values of 14.4 and 14.3 for FIC and BFC in DMSO, respectively, so to the extent that these values are surprising, the problem must mostly come from a distortion of our perspective, and not from any dispute between theory and experiment. If that is the case, then what might be a better perspective to investigate CpH-PAH structures like BIC and FIC?
Let us consider a perspective based on the teachings of Clar and Loschmidt. Clar held the electron sextet as the paragon of aromaticity [14] and Loschmidt considered benzene a "group element" [15]. Combining these perspectives, one can view PAH networks as graphs comprised of Clar–Loschmidt (CL) "elements" [10,16]. For corannulene, a maximum of two CL elements (sextets) are present in any one resonance form. Therefore, applying the CL perspective to BFC and FIC one sees dibenzofluorene (full name 13H-dibenzo[a,c]fluorene) as a cognate, which highlights these compounds to be like a conformationally planar-locked 2,3-diphenylindene (pKa 17.7 in DMSO [17].
The CL perspective also allows one to create graphs with points representing the aromatic elements, such as, benzene, Cp anion, pyrrole, etc. (Scheme 2, top). Triphenylene is a simple cyclic graph of three benzene elements and their direct connections. A simple modification of the triphenylene graph is the replacement of benzene by Cp anion (Scheme 2, middle). From the CL graph perspective, the anion of cyclopentaphenanthrene is related to that of dibenzofluorene the same way as the Cp anion is related to that of indene. Therefore, cyclopentaphenanthrene should be more acidic than dibenzofluorene, which it is [18]. Monofluorenocorannulene and dibenzofluorene map onto a common CL graph with benzene and indenyl elements (Scheme 2, bottom). These compounds display related electronic structures and are predicted computationally to have similar pKa values (i.e., show similar acidity).
The extreme of this graph example is the five rim benzene elements and a Cp− core, pentabenzocorannulene (Scheme 3, upper) [19,20]. Such a CpH derivative should be a very acidic aromatic hydrocarbon, and computations for this study predict its pKa to be ca 1.1 in DMSO. This high acidity is well seen if one recognizes that, when seen through the CL model, pentabenzocorannulene is more of a well-overlapping pentaphenyl-CpH derivative than a pentabenzocorannulene. Swap out Cp− by the isolelectronic pyrrole and it becomes clear that pentabenzoazacorannulene is also less of an azocorannulene and more a pentaphenylpyrrole, which explains its ease of synthesis and the physical properties [21,22].
The graph with four benzene elements and a CpH core is evident in the commercially available tetrabenzofluorene, TBF (Scheme 3, lower). The parent and the phenyl TBF derivative were also recently prepared by a mechanochemical Scholl reaction [23-25]. Furthermore, the anion of TBF is stable enough to be easily handled and the crystal structures of a number of common ammonium salts of TBF anion, grown from water, have been reported [26].
Neutral highly acidic compounds tend to deprotonate to form highly stable anions. The air stability of TBF conforms to that model. In contrast, although the dianion of BFC can be generated with tert-butoxide under inert atmosphere, in the presence of water it quenches and exposed to oxygen it oxidizes to form the diketone. Access to the dianion of BFC presages an interesting diradical and this was achieved by inclusion of mesityl protecting groups [27]. Extension of the BFC model with thiophene provides further interesting materials [28]. A reasonable corollary to this behavior would assert that derivatives of TBF and pentabenzocorannulene would produce air-stable radical and ionic PAHs, and that coupling such fragments would lead to stable redox-active carbon sheets.
Conclusion
In conclusion, the “surprise” in the surprising pKa for BIC and FIC was in our expectation of the deprotonated forms as poorly delocalized fluorenyl anions. The CL perspective provided us a different way to look at these compounds and interpret them as phenyl-substituted CpHs, thus reconciling their higher acidity. That perspective led us to a CL graph representation that predicted/rationalized additional acidic CpH-PAHs and sharpened our understanding of an “azacorannulene”. Working from graph structures based on chemically stable (“group elementary”) nodes is a useful principle in molecular design and chemical synthesis. Such a perspective is important to understand fundamental physical organic molecular properties as well as to predict desirable unforeseen new material designs (e.g., redox-active electrooptical materials).
Methods
pKa measurements in acetonitrile
The experimental setup and methodology for the pKa determination of BFC and FIC in acetonitrile was essentially the same as described in detail in previous papers [12,13]. A brief description will follow.
The pKa determinations in acetonitrile are based on the determinations of differences of pKa values of two acids. In this case one compound is a reference acid with a previously known pKa value and the other acid is either FIC or BFC. Both compounds, as well as the references are also separately titrated in order to obtain the UV–vis spectra of the acids in neutral as well as in deprotonated forms.
Then, the same titration is done with a mixture of measured acid (FIC or BFC) and a reference acid. Using the spectral data from the titrations of mixtures the dissociation levels α = [A−]/([A−] + [HA]) of both acids in all the mixtures formed during titration are calculated and are then in turn used to calculate the differences of pKa values (∆pKa) of FIC or BFC and the used reference acids according the following equation:
The pKa values of BFC and FIC are estimated as a result of ∆pKa measurements against different reference acids.
A Perkin-Elmer Lambda 40 UV–vis spectrophotometer connected with optical fibre cables to an external cell compartment inside a MBraun Unilab glovebox filled with 99.999% pure argon was used for the spectrophotometric titrations. This setup ensured that during all titrations the moisture and oxygen contents in argon were always under 10 ppm during measurements. Triflic acid (Aldrich, 99+%) and phosphazene base P2-Et (165535-45-5, Aldrich, ≥98%) were used to prepare the acidic and basic titrant solutions, respectively. Acetonitrile (Romil 190 SpS far UV/gradient quality) was used as solvent after drying with molecular sieves (3 Å), which lowered the water content to a range of 2–6 ppm.
Acids with previously published pKa values in acetonitrile were used as reference acids [12,13]. Usually this kind of ∆pKa measurements are done against three reference acids but due to the lack of suitable reference acids in the weakly acidic pKa region in acetonitrile only two reference acids were used. The results are presented in Table 2. We estimate the standard uncertainties of the assigned pKa values of BFC and FIC as 0.15 pKa units.
Table 2:
Results of pKa measurements in acetonitrile.
Acid
Reference acid
pKa (Ref)
ΔpKa
pKa (acid)
assigned pKa
BFC
9-C6F5-fluorene
28.14
0.55
27.59
27.6
(4-Me-C6F4)(C6H5)CHCN
26.98
−0.78
27.76
FIC
9-C6F5-fluorene
28.14
0.38
27.76
27.8
(4-Me-C6F4)(C6H5)CHCN
26.98
−0.94
27.92
pKa predictions using correlation analysis
The pKa values of fluorene, indene, cyclopentadiene, 9-phenylfluorene, and diphenylmethane in acetonitrile were estimated by averaging the values obtained from their experimental pKa values in DMSO [5] and three correlation equations: equations 2.1 and 2.2 in Ref. [13] and correlation equation composed of experimental pKa values of 9-C6F5-fluorene, octafluorofluorene, fluoradene, (4-Me-C6F4)(C6H5)CHCN, (C6F5)(C6H5)CHCN, 9-COOMe-fluorene, and 9-CN-fluorene in acetonitrile [12,13] and DMSO [5]. All these pKa estimations involve significant extrapolation and the agreement between the estimates obtained from different equations is not good. Therefore, the standard uncertainties of these obtained pKa estimates are high: 2 pKa units for fluorene, indene and 9-phenylfluorene, 3 pKa units for cyclopentadiene (as this acid is very different from the acids used for developing the correlation equations), and 3 pKa units for diphenylmethane (as very strong extrapolation is involved with this acid) [12,13].
The pKa values of FIC and BFC in DMSO were estimated from the experimental pKa values in acetonitrile determined in this study using the average values from the same three equations as above, however, using them in reverse. The standard uncertainty estimates of the pKa values of FIC and BFC in DMSO are 0.8 pKa units.
Synthesis
Material used to measure the acidity of compounds BFC and FIC was prepared previously by methods reported in reference [9].
Computational methods
The structural and energetic analyses of the molecular systems for all compounds described in this study were carried out with the B97-D dispersion enabled density functional method [29,30], using an ultrafine grid, together with the def2-TZVPP basis set [31]. Full geometry optimizations were performed and uniquely characterized via second derivatives (Hessian) analysis to establish stationary points and effects of zero point and thermal energy contributions. Effects of solvent employed the COSMO:ab initio continuum method [32,33], using a dielectric as in experiment, and fully optimized using radii of Klamt [34]. Electronic and thermal free energy differences between neutral and anion were compared to a reference pentafluorophenylfluorene, and subsequently converted to pKa values by dividing by 1.36. Visualization and analysis of structural and property results were obtained using Avogadro [35]. Our group develops GAMESS [36] and has also contributed to Gaussian software packages, in this work the G09 ES64L-G09RevE.01 version [37] of the latter was used.
Acknowledgements
High-performance computational work was carried out using the supercomputer at SPST Tianjin University.
Funding
Experimental work at Tartu was carried out using the instrumentation at the Estonian Center of Analytical Chemistry (TT4, https://www.akki.ee) and was supported by the Estonian Research Council grant (PRG690) and by the Estonian Ministry of Education and Research (TK210).
Data Availability Statement
All data that supports the findings of this study is available in the published article and/or the supporting information of this article.
References
Dicks, A. P. J. Chem. Educ.2003,80, 1322–1327. doi:10.1021/ed080p1322
Return to citation in text:
[1]
[2]
Streitwieser, A., Jr. Molecular Orbital Theory for Organic Chemistry; Wiley: New York, NY, USA, 1961.
Return to citation in text:
[1]
[2]
Bordwell, F. G.; Drucker, G. E.; Fried, H. E. J. Org. Chem.1981,46, 632–635. doi:10.1021/jo00316a032
Return to citation in text:
[1]
[2]
[3]
[4]
[5]
Bordwell, F. G.; Matthews, W. S.; Vanier, N. R. J. Am. Chem. Soc.1975,97, 442–443. doi:10.1021/ja00835a048
Return to citation in text:
[1]
[2]
Bordwell, F. G. Acc. Chem. Res.1988,21, 456–463. doi:10.1021/ar00156a004
Return to citation in text:
[1]
[2]
[3]
[4]
Kinoshita, T.; Fujita, M.; Kaneko, H.; Takeuchi, K.; Yoshizawa, K.; Yamabe, T. Bull. Chem. Soc. Jpn.1998,71, 1145–1149. doi:10.1246/bcsj.71.1145
Return to citation in text:
[1]
Baldridge, K. K.; Siegel, J. S. J. Phys. Org. Chem.2004,17, 740–742. doi:10.1002/poc.819
Return to citation in text:
[1]
Liu, S.; Roch, L. M.; Allemann, O.; Xu, J.; Vanthuyne, N.; Baldridge, K. K.; Siegel, J. S. J. Org. Chem.2018,83, 3979–3986. doi:10.1021/acs.joc.8b00320
Return to citation in text:
[1]
[2]
Baldridge, K. K.; Siegel, J. S. Angew. Chem., Int. Ed.2013,52, 5436–5438. doi:10.1002/anie.201300625
Return to citation in text:
[1]
[2]
Sritana-Anant, Y.; Seiders, T. J.; Siegel, J. S. Top. Curr. Chem.1998,196, 1–43. doi:10.1007/3-540-69691-1_1
Return to citation in text:
[1]
Kütt, A.; Leito, I.; Kaljurand, I.; Sooväli, L.; Vlasov, V. M.; Yagupolskii, L. M.; Koppel, I. A. J. Org. Chem.2006,71, 2829–2838. doi:10.1021/jo060031y
Return to citation in text:
[1]
[2]
[3]
[4]
[5]
[6]
Kütt, A.; Tshepelevitsh, S.; Saame, J.; Lõkov, M.; Kaljurand, I.; Selberg, S.; Leito, I. Eur. J. Org. Chem.2021, 1407–1419. doi:10.1002/ejoc.202001649
Return to citation in text:
[1]
[2]
[3]
[4]
[5]
[6]
[7]
[8]
[9]
Clar, E. The Aromatic Sextet; Wiley: New York, NY, USA, 1972.
Return to citation in text:
[1]
Loschmidt, J. Konstitutionsformeln der Organischen Chemie in Graphischer Darstellung; Wilhelm Engelmann: Leipzig, Germany, 1860.
Return to citation in text:
[1]
Rigby, S. S.; Stradiotto, M.; Brydges, S.; Pole, D. L.; Top, S.; Bain, A. D.; McGlinchey, M. J. J. Org. Chem.1998,63, 3735–3740. doi:10.1021/jo9711504
Return to citation in text:
[1]
Bordwell, F. G.; Drucker, G. E. J. Org. Chem.1980,45, 3325–3328. doi:10.1021/jo01304a037
Return to citation in text:
[1]
Vianello, R.; Maksić, Z. B. Eur. J. Org. Chem.2005, 3571–3580. doi:10.1002/ejoc.200500153
Return to citation in text:
[1]
Kanthak, M.; Muth, E.; Dyker, G. Eur. J. Org. Chem.2010, 6084–6091. doi:10.1002/ejoc.201000821
Return to citation in text:
[1]
Yokoi, H.; Hiraoka, Y.; Hiroto, S.; Sakamaki, D.; Seki, S.; Shinokubo, H. Nat. Commun.2015,6, 8215. doi:10.1038/ncomms9215
Return to citation in text:
[1]
Ito, S.; Tokimaru, Y.; Nozaki, K. Angew. Chem., Int. Ed.2015,54, 7256–7260. doi:10.1002/anie.201502599
Return to citation in text:
[1]
Báti, G.; Csókás, D.; Stuparu, M. C. Chem. – Eur. J.2024,30, e202302971. doi:10.1002/chem.202302971
Return to citation in text:
[1]
Hopkinson, A. C.; Lee-Ruff, E.; Maleki, M. Synthesis1986, 366–371. doi:10.1055/s-1986-31637
Return to citation in text:
[1]
Schröder, K.; Haase, D.; Saak, W.; Lützen, A.; Beckhaus, R.; Wichmann, S.; Schellenberg, J. Organometallics2006,25, 3824–3836. doi:10.1021/om0602613
Return to citation in text:
[1]
Frey, O. N.; Stalling, T.; Schlüter, F.; Saak, W.; Schmidtmann, M.; Haase, D.; Beckhaus, R. Dalton Trans.2016,45, 1085–1092. doi:10.1039/c5dt04137f
Return to citation in text:
[1]
Pedersen, V. B. R.; Jasti, R.; Nielsen, M. B. ChemPhotoChem2024,8, e202400122. doi:10.1002/cptc.202400122
Return to citation in text:
[1]
Grimme, S. J. Comput. Chem.2006,27, 1787–1799. doi:10.1002/jcc.20495
Return to citation in text:
[1]
Grimme, S. J. Chem. Phys.2006,124, 034108. doi:10.1063/1.2148954
Return to citation in text:
[1]
Weigend, F.; Ahlrichs, R. Phys. Chem. Chem. Phys.2005,7, 3297–3305. doi:10.1039/b508541a
Return to citation in text:
[1]
Baldridge, K. K.; Biggs, B.; Gilbertson, R. D.; Haley, M. M.; Bläser, D.; Boese, R.; Maulitz, A. H.; Siegel, J. S. Chem. Commun.1998, 1137–1138. doi:10.1039/a801217b
Return to citation in text:
[1]
Baldridge, K.; Klamt, A. J. Chem. Phys.1997,106, 6622–6633. doi:10.1063/1.473662
Return to citation in text:
[1]
Klamt, A.; Jonas, V.; Bürger, T.; Lohrenz, J. C. W. J. Phys. Chem. A1998,102, 5074–5085. doi:10.1021/jp980017s
Return to citation in text:
[1]
Hanwell, M. D.; Curtis, D. E.; Lonie, D. C.; Vandermeersch, T.; Zurek, E.; Hutchison, G. R. J. Cheminf.2012,4, 17. doi:10.1186/1758-2946-4-17
Return to citation in text:
[1]
Schmidt, M. W.; Baldridge, K. K.; Boatz, J. A.; Elbert, S. T.; Gordon, M. S.; Jensen, J. H.; Koseki, S.; Matsunaga, N.; Nguyen, K. A.; Su, S.; Windus, T. L.; Dupuis, M.; Montgomery, J. A., Jr. J. Comput. Chem.1993,14, 1347–1363. doi:10.1002/jcc.540141112
Return to citation in text:
[1]
Gaussian 09, Revision D01; Gaussian, Inc.: Wallingford, CT, 2009.
Return to citation in text:
[1]
Kütt, A.; Leito, I.; Kaljurand, I.; Sooväli, L.; Vlasov, V. M.; Yagupolskii, L. M.; Koppel, I. A. J. Org. Chem.2006,71, 2829–2838. doi:10.1021/jo060031y
13.
Kütt, A.; Tshepelevitsh, S.; Saame, J.; Lõkov, M.; Kaljurand, I.; Selberg, S.; Leito, I. Eur. J. Org. Chem.2021, 1407–1419. doi:10.1002/ejoc.202001649
Liu, S.; Roch, L. M.; Allemann, O.; Xu, J.; Vanthuyne, N.; Baldridge, K. K.; Siegel, J. S. J. Org. Chem.2018,83, 3979–3986. doi:10.1021/acs.joc.8b00320
Kütt, A.; Leito, I.; Kaljurand, I.; Sooväli, L.; Vlasov, V. M.; Yagupolskii, L. M.; Koppel, I. A. J. Org. Chem.2006,71, 2829–2838. doi:10.1021/jo060031y
13.
Kütt, A.; Tshepelevitsh, S.; Saame, J.; Lõkov, M.; Kaljurand, I.; Selberg, S.; Leito, I. Eur. J. Org. Chem.2021, 1407–1419. doi:10.1002/ejoc.202001649
Kütt, A.; Leito, I.; Kaljurand, I.; Sooväli, L.; Vlasov, V. M.; Yagupolskii, L. M.; Koppel, I. A. J. Org. Chem.2006,71, 2829–2838. doi:10.1021/jo060031y
13.
Kütt, A.; Tshepelevitsh, S.; Saame, J.; Lõkov, M.; Kaljurand, I.; Selberg, S.; Leito, I. Eur. J. Org. Chem.2021, 1407–1419. doi:10.1002/ejoc.202001649
Kütt, A.; Leito, I.; Kaljurand, I.; Sooväli, L.; Vlasov, V. M.; Yagupolskii, L. M.; Koppel, I. A. J. Org. Chem.2006,71, 2829–2838. doi:10.1021/jo060031y
13.
Kütt, A.; Tshepelevitsh, S.; Saame, J.; Lõkov, M.; Kaljurand, I.; Selberg, S.; Leito, I. Eur. J. Org. Chem.2021, 1407–1419. doi:10.1002/ejoc.202001649
Kütt, A.; Leito, I.; Kaljurand, I.; Sooväli, L.; Vlasov, V. M.; Yagupolskii, L. M.; Koppel, I. A. J. Org. Chem.2006,71, 2829–2838. doi:10.1021/jo060031y
13.
Kütt, A.; Tshepelevitsh, S.; Saame, J.; Lõkov, M.; Kaljurand, I.; Selberg, S.; Leito, I. Eur. J. Org. Chem.2021, 1407–1419. doi:10.1002/ejoc.202001649
Kütt, A.; Leito, I.; Kaljurand, I.; Sooväli, L.; Vlasov, V. M.; Yagupolskii, L. M.; Koppel, I. A. J. Org. Chem.2006,71, 2829–2838. doi:10.1021/jo060031y
13.
Kütt, A.; Tshepelevitsh, S.; Saame, J.; Lõkov, M.; Kaljurand, I.; Selberg, S.; Leito, I. Eur. J. Org. Chem.2021, 1407–1419. doi:10.1002/ejoc.202001649
Liu, S.; Roch, L. M.; Allemann, O.; Xu, J.; Vanthuyne, N.; Baldridge, K. K.; Siegel, J. S. J. Org. Chem.2018,83, 3979–3986. doi:10.1021/acs.joc.8b00320
Baldridge, K. K.; Biggs, B.; Gilbertson, R. D.; Haley, M. M.; Bläser, D.; Boese, R.; Maulitz, A. H.; Siegel, J. S. Chem. Commun.1998, 1137–1138. doi:10.1039/a801217b
33.
Baldridge, K.; Klamt, A. J. Chem. Phys.1997,106, 6622–6633. doi:10.1063/1.473662
Schmidt, M. W.; Baldridge, K. K.; Boatz, J. A.; Elbert, S. T.; Gordon, M. S.; Jensen, J. H.; Koseki, S.; Matsunaga, N.; Nguyen, K. A.; Su, S.; Windus, T. L.; Dupuis, M.; Montgomery, J. A., Jr. J. Comput. Chem.1993,14, 1347–1363. doi:10.1002/jcc.540141112