HOMA Index Establishes Similarity to a Reference Molecule

The article shows that the definition of the HOMA index of geometrical aromaticity satisfies the axioms of a similarity function between the examined and benzene ring. Consequently, for purely mathematical reasons, the index works exceptionally well as an index of aromaticity: it expresses a geometric similarity to the archetypal aromatic benzene. Thus, if the molecule is geometrically similar to benzene, then it is also chemically similar, and therefore, it is aromatic. However, the similarity property legitimizes using the HOMA-like indices to express similarity to molecules other than benzene, whether cyclic or linear and existing or hypothetical. The paper demonstrates an example of HOMA-similarity to cyclohexane, which expresses a (relaxed)-saturicity property not accompanied by strong structural strains or steric hindrances. Further, it is also shown that the HOMA index can evaluate the properties of whole molecules, such as 25 unbranched catacondensed isomers of hexacene. The index exhibits a significant quadratic correlation with the total energy differences of planar isomers from which the nonplanar ones deviate. Moreover, the HOMA index of hexacene isomers significantly correlates with the Kekulé count connected to the resonance energy in the Hückel approximation. As a result, the study shows that the HOMA index can be used not only for aromaticity analyses but also as a general chemical descriptor applicable to rings, chains, composed molecular moieties, or even whole molecules.


INTRODUCTION
Aromaticity is not an observable molecular property, although it is essential in modern chemistry.It manifests itself in multiple molecular characteristics and is thus defined in various ways, highlighting its unique features.Different definitions expose the specifics of molecular geometry, reactivity, additional energy stabilization, magnetic currents, and some other features and sometimes disagree.Therefore, it is imperative to understand every aspect of the aromaticity definition.This study is focused on analyzing the definition of the Harmonic Oscillator Model of Aromaticity, HOMA, and geometrical aromaticity index. 1 It is one of the most successful and widely used measures of aromaticity 2−4 despite its foundation being grounded on a simple approximative quantum chemical model.Nevertheless, at the very outset, the authors stated that the index "is quite independent/in principle/of quantum chemical models". 1 Hereafter, we demonstrate why it is the absolutely correct statement.
For hydrocarbon rings, HOMA is defined as follows where R i and R opt (Å) stand for the i-th CC bond length in the analyzed ring and the reference benzene ring (R opt = R B ) for which HOMA = 1, respectively, n is number of CC bonds in the ring, and α = 257.7 Å −2 normalizes the index to be unitless and equal to 0 for a hypothetical perfectly alternating Kekuleć yclohexatriene ring.
−7 = HOMA 1 EN GEO where R av is the averaged bond length in the examined ring, EN is a normalized squared difference of the bond length arithmetic means, and GEO is a normalized bond variance in the examined molecule.
−14 Interestingly, the HOMA index asymptotically converges as the number n of molecular units increases.For example, for alkanes, spiro cyclopropanes, polyynes, and allenes, the index decreases with the increasing n and converges to negative values from −4 to −2. 13 This can suggest an increase in the degree of insulation with an increase in the number of units.For alkenes and polythiophenes, the index increases and converges to values above 0.6, while for polypyrrole, it reaches ca.0.8. 13,15This can suggest the increase of π-electron delocalization with an increase in the number of units. 13,16The HOMA indices were also used to study linear and branched alkanes, and for n-alkanes, a good correlation with the boiling points was found.On the other hand, for n-alkane constitutional isomers, the HOMA index tends to increase with the boiling points from the most branched to the most extended isomer. 17he above data indicate that the HOMA index is more than just a geometrical aromaticity index.However, why it may be so is known. 16First, the HOMA index is connected to a distance in the abstract n-dimensional molecular space and, second, to a statistical value such as a variance. 17Third, the association of HOMA with aromaticity is only through R opt , taken from the aromaticity standard: benzene.The α factor establishes HOMA = 0 for a nonaromatic, hypothetical cyclohexatriene molecule.
This study shows that HOMA expresses a similarity between the studied and reference moiety.If a molecule is geometrically similar to benzene, HOMA-similar, then it also is chemically similar and, thus, aromatic.By replacing benzene with a different reference molecule, the HOMA index with this new reference provides geometrical similarity to the new reference.Consequently, the similarity of the chemical properties of the new reference is also satisfied.Thus, the HOMA index referred to benzene measures aromaticity, but when referred to cyclohexane, it provides similarity to an unstrained, saturated hydrocarbon.We demonstrated that HOMA can also be used to evaluate the properties of entire molecules, such as isomers of hexacene.Ultimately, we argue that the HOMA index can describe a general geometric similarity between chemical moieties and can be shaped according to the unconstrained will of a researcher rather than the form intended by its inventors.

COMPUTATIONS
All structure optimizations were done with the DFT calculations using the B3LYP functional, 18,19 along with the D3 Grimme correction for dispersion forces, 20 the 6-31G** Pople-type basis set, 21 and Gaussian 09 software. 22The 6-31G** basis set was chosen to enable the reproduction of the calculations in any laboratory and because this small basis set yields sufficiently good geometry, energy, and electron density data 23 to build organic physical chemistry indices.It was confirmed that each structure exhibits proper symmetry and is at the energy minimum, owning all positive harmonic frequencies.The XYZ coordinates as well as the original input and output Gaussian 16 files of all molecules calculated are given in the Supporting Information.

EN and GEO Factors.
Let us denote by x 1 ,x 2 ,•••,x n and y 1 ,y 2 ,•••,y n the CC bonds in the analyzed and reference rings x and y, respectively.Then, EN factor (3) can be redrawn as the square of the difference between the CC bond lengths arithmetic means x ̅ = R av and y ̅ = R opt normalized by α (5) The GEO factor ( 4) is then a normalized variance of the CC bond lengths in the examined moiety x So, to calculate the HOMA index, it is enough to know the appropriate arithmetical means and the variation of the CC bond lengths in the single moiety x.Note that only the GEO = α• Var(x) term contains information about the ring size hidden in the summation operation.For the rings other than sixmembered, the summation stops at n ≠ 6, but still, R opt = R B is used.So, it is tacitly assumed that in every reference ring, the CC distance equals that in benzene, even though the matching n-membered hydrocarbon ring does not exist.Hence, HOMA = 1 for every reference ring consisting of k = 3, 4, 5, 7,••• bonds, and it is as ideally aromatic as benzene.However, if arithmetical means of comparing moieties and variation of the ring bond lengths in one of them are sufficient to determine HOMA, then even a single distance can play a role in the reference moiety.It can be present in a selected molecule but can also be equal to the arithmetic mean of the bond lengths in an existing or a hypothetical moiety.

HOMA Establishes a Similarity between Moieties.
A similarity between two elements x and y of a set X is a function × s X X : which for any x,y ∈ X satisfies the following axioms 24 = s x y s x y s y x y X x X s x y s y y Notice that if α = 1, then HOMA becomes a similarity function s HOMA (x,y) The non-negativity is satisfied since for hydrocarbons 0 < (x i − y ̅ ) 2 < 1, so < < x y 0 ( ) 1 ; thus, . The symmetricity is given by the symmetricity of the (x i − y ̅ ) 2 factors.Eventually, if x i = y ̅ for any i, then HOMA = s HOMA (y,y) = 1 is the maximum.Notice that HOMA truncated to a single EN factor, 1 − EN, for α = 1 is the similarity function for analogous reasons.Notably, HOMA is connected also to a distance in abstract molecular space. 17owever, the HOMA non-negativity is not essential in studying the aromaticity of molecules.Notice that α ≈ 250 > 1 provides most figures between 1 and ca.−10, which have been convenient so far.Besides, returning to the strict similarity function is trivial.Thus, for chemical purposes, let us use a similarity-like function satisfying only symmetricity and having maximum axioms (8).Then, the HOMA and HOMA shortened to the single EN factor functions are just such "chemical" similarity functions.
Therefore, for the hexahydroxy and hexafluoro benzene, C 6 (OH) 6 and C 6 F 6 , the rings are geometrically nearly identical to that of benzene, HOMA ≈ 1. 17 The side ring in triphenylene is quite similar to benzene (HOMA = 0.93, Table 1).The offcenter coronene and C 6 Li 6 rings are less similar to benzene.Indeed, the CC bonds in the former vary, and though the latter has D 6h symmetry, the CC bonds are longer than those in benzene.Hence, they have HOMA ≈ 0.8, 17  , equals ca.0.5 (Table 1).Although nothing seems similar between the five-membered pentalene dianion and benzene rings, in terms of similarity established by the HOMA function, they are somehow aromatic, which for C 8 H 6 2 is confirmed by the NICS index. 25Again, according to (7), HOMA-similarity requires similar CC bond arithmetic means in the studied and reference rings and slight variation of the CC bonds in the former.It is difficult to judge whether the central ring in triphenylene (HOMA = 0.28) is still similar to benzene or not.Unsurprisingly, ethene and cyclooctatetraene are HOMA-dissimilar to benzene and are not aromatic: HOMA ≈ 0.1 and −0.1, respectively (Table 1).The CC bond in the former is double and radically different from that in aromatic benzene, whereas the double and single CC bonds in the latter alter, and the molecule has a significant GEO factor.In agreement with the intention of the HOMA index founders, 1 similarity with benzene ends at HOMA = 0, reflecting the nonaromaticity of the perfectly alternating Kekulécyclohexatriene ring.However, according to HOMA, chair cyclohexane and archetypically antiaromatic cyclobutadiene rings have similar HOMA indices around −4.4 (Table 1).Moreover, the HOMA of the ring in antiaromatic pentalene is close to 0. Thus, the HOMA-similarity with benzene does not allow for distinguishing between cycloalkane, nonaromatic, and antiaromatic hydrocarbons.We shall check whether such differentiation could be possible by going beyond the primary definition in Section 3.3.2.
Finally, let us emphasize that the HOMA-similarity criterion is not always intuitive and should be operated deliberately.Equation 7states that a molecule must have similar average bond lengths and as slight bond length variation as possible to be similar to the reference one.An unreflective HOMA analysis of the CC distances in tetracyano, tetrathiocyanato, and 1,1difluoro-2,2-dichloro ethylene, C 2 (CN) 4 , C 2 (NCS) 4 , and F 2 C�CCl 2 , with the CC distances at the B3LYP/D3/6-31G** level equal to 1.3719, 1.3773, and 1.3351 Å, respectively, would lead to the conclusion that the first two molecules are aromatic, while the third is nonaromatic, as HOMA equals 0.85, 0.91, and 0.04, respectively.Therefore, the HOMA analysis must always be placed in a narrowly defined chemical context.1).

Journal of Chemical Information and Modeling
cation (C 7 H 7 + , D 7h ), while the anion has slightly lowered symmetry by the Jahn−Teller effect. 26The persubstituted molecules listed above, like, e.g., hexamethylbenzene, 27 or fragments of a larger symmetric structure, such as the central ring in coronene, also have all identical CC bonds.Nevertheless, one-CC-bond molecules are the limit case eq 7 allows.They would set the HOMA-similarity to the most basic organic molecules such as ethene or acetylene.1. Tricyclic adamantane, acyclic and branched neopentane, and unstrained cycloheptane are exceptionally similar to cyclohexane.Indeed, they are even more HOMA C -similar than the boat cyclohexane.The more strained the cyclic alkane (cyclopentane, cyclobutane, cyclopropane), the less HOMA Csimilar it is.After all, the three-membered ring of bicyclo[1.1.0]butane(HOMA C ≈ 0.36), the central ring in coronene (HOMA C ≈ 0.20), and the cyclobutadiene dianion (HOMA C ≈ 0.10) already seem to be not similar to cyclohexane.
Remark that except for saturated not-strained cyclic adamantane and cycloheptane, HOMA C of neopentane equals 1.00, and HOMA C of isopentane, isohexane, hexane, pentane, and butane exceeds 0.97 (Table 1).Thus, for a moiety to be HOMA C -similar, saturation instead of molecular cyclicity is necessary.Such a similarity allows for a moderate strain as in cyclobutane (HOMA C = 0.96) or a much larger one as in cyclopropane (HOMA C = 0.85) or even in cubane (HOMA C = 0.80).Still, a three-membered ring in bicyclo[1.1.0]butaneexhibits HOMA C = 0.36, and if the six-membered ring is flat as in 12H-coronene, in which each C atom is double-bonded with the external rim, HOMA C = 0.20.The index of very strained, elusive tetrahedrane is equal to −0.81.Thus, the HOMA C index expresses a (relaxed)-saturicity not accompanied by strong structural strains or steric hindrances.
To better see the difference between using the HOMA B and HOMA C indices, let us look at the HOMA C (A) = f(HOMA B (A)) function, where A stands for an arbitrary moiety (Figure 2a).At first glance, the points are spread erratically.Nevertheless, at closer inspection, a regular, nonlinear boundary can be seen (Figure 2a).One boundary branch goes through ethyne, allene, and benzene, and the other through cyclohexane and similar compounds.The former seems composed of primarily unsaturated compounds, while the other is composed of nonstrained structures, irrespectively saturated or not.All interior points seem to correspond to strained and simultaneously unsaturated hydrocarbons.
Notwithstanding partial chaos present in Figure 2a, the HOMA(A) decomposition into the EN(A) and GEO(A) factors reveals where the irregularities come from.The plot of the HOMA(A) function truncated to the 1-EN(A) term forms the border curve where all of the moieties are positioned (black points, Figure 2b).The HOMA(A) function truncated to the 1-GEO(A) term establishes a half straight line (red points, Figure 2b).All of the moieties with GEO = 0 are in the point (1,1) -the beginning of the half straight line.The benzene and cyclohexane molecules and all of the moieties are located at the border curve in Figure 2a.Thus, the border is formed by moieties that have all CC bonds of equal lengths and thus GEO = 0, whereas the interior is formed by the moieties that have (significant) bond length alteration.Seemingly, the curve in Figure 2b is a deformed parabola, but the rotation of this curve by 45°using a pair of parametric equations shows that it is nothing but a parabola (Figure S1 in the Supporting Information).

Similarity to a
Molecule with a Nonzero GEO Factor.Now, the feasibility of discrimination between nonaromatic and antiaromatic hydrocarbons is addressed using the HOMAsimilarity function.The former compounds should exhibit single and double bond alternation in the ring, as in the elusive cyclohexatriene, while the latter should have 4 instead of 6 πelectrons per ring, as Huckel's rule for the aromatic systems predicted.The HOMA B index does not allow for distinguishing these species classes (Section 3.2).However, the HOMA C index enables an additional qualification: the nonaromatic and antiaromatic rings have a negative HOMA C index, and they are dissimilar to cyclohexane (Table 1).Hence, the necessary, though insufficient, condition of being a non-or antiaromatic ring is to be HOMA B and HOMA C dissimilar.Immediately a question appears of whether a sufficient condition cannot be formulated using the HOMA index.The answer would be positive if we could find good references for these hydrocarbons.
The first problem in finding such references is the presence of the GEO factor in the archetypical cyclohexatriene and cyclobutadiene rings.In such a case, the condition of the HOMA-similarity should be supplemented by comparing the GEO factors.Hence, expression 7 has to be rewritten in a new HOMA Δ form as follows where module of variances |•| is necessary to conserve the having a maximum axiom (8) for HOMA Δ .To calculate the HOMA CHT Δ -similarity to cyclohexatriene (CHT), it is necessary first to know the averaged CC distances and their variation in the cyclohexatriene moiety.From Krygowski's original paper addressing the parameters of the HOMA index for cyclohexatriene, 28 the following relationship between the bond distances and the parameter α emerges 17 According to HOMA CHT Δ , the following moieties from Table 1 are the most similar to cyclohexatriene: ring in pentalene, cycloheptatriene anion, cyclooctatetraene, ring in 1,4-dihydropentalene, butadiene, and benzyne: their HOMA CHT Δ is equal to 0.91, 0.81, 0.77, 0.57, 0. 54, and 0.52, respectively.This sequence shows that the HOMA CHT Δ similarity criterion is not discriminative.Indeed, pentalene, 1,4-dihydropentalene, and the cycloheptatriene anion in the singlet state are antiaromatic The calculations were performed at the B3LYP/D3/6-31G** level with the Grimmes D3 empirical dispersion correction.Molecule numbering is as in Figure 3. Moieties 2, 3, 6, 9, and 10 are similar to cyclohexatriene shown in bold.2).

Journal of Chemical Information and Modeling
for their 4n π electrons per ring 29 and the NICS index criterion for the anion. 30Nonplanar cyclooctatetraene with alternating single and double bonds is nonaromatic.Butadiene has bond length alternation but is not cyclic, and we are not inclined to call it nonaromatic in the sense that the cyclohexatriene and cyclooctatetraene rings are.The structural formula of benzyne (cyclohexa-1-yn-3,5-diene) is formally similar to cyclohexatriene's but, according to the NICS criterion, is more aromatic than benzene. 31et us therefore apply the HOMA CHT Δ index to rings very close to cyclohexatriene.A search for rings similar to "Kekuleb enzene" among a series of triply fused hexasubstituted benzenes with the C 3 axis using HOMA B = 0 as a criterion demonstrated that only some molecules fulfilling HOMA B ≈ 0 met another necessary condition such that R(C−C) > R B > R(C�C). 17However, most of the molecules satisfying the two necessary conditions contained the reactive (antiaromatic) pentalene system fused to the central ring.Such compounds are inappropriate as model reference molecules.
Application of the HOMA CHT Δ criterion to the molecules exhibiting HOMA B ≈ 0 considered before 17 shows that HOMA CHT Δ exhibits similarity to cyclohexatriene only for those molecules that satisfied R(C−C) > R B > R(C�C), whereas for the others, HOMA CHT Δ is close to 0 or negative (Table 2).This demonstrates that eq 10 can help search for moieties similar to rings with significant GEO factors and that the search can be more effective than using the classical HOMA index supplemented by additional criteria (Figure 3). 17 moderate success of the HOMA CHT Δ index in examining the similarity to cyclohexatriene prompts the supposition that a good index can also be constructed for antiaromatic compounds.However, the issue of a good reference system is returned.The strain of the ring of the antiaromatic cyclobutadiene determines its HOMA CBD Δ index.Therefore, it is inadequate as a reference for antiaromaticity.On the other hand, the similarity of less strained antiaromatic pentalene to cyclohexatriene makes discrimination between nonaromatic and antiaromatic hydrocarbons using sole HOMA Δ similarity function doubtful.The sole representation of the 4n π-electron requirement in the bond lengths and their variations may not be sufficiently specific.However, further study of this issue goes beyond this project.
At the end of this section, let us mention that the HOMA Δ index can also be constructed based on a reference free of any symmetry constraints.If one assumes the averaged bond length y ̅ and Var(y) in eq 10 are taken from an arbitrary asymmetric reference, then the HOMA Δ value shows the HOMA Δ -similarity to such a symmetryless compound.This can help assess the similarity to asymmetric compounds such as isooctane, on which the octane number is based, but perhaps also to important natural products such as terpenes.Although, for an asymmetric reference, the analysis could become fuzzier, without further detailed research, it is impossible to assess the benefits of such an approach.However, detailed research on this issue again goes beyond this project.
3.3.3.HOMA as an Index of the Entire Molecule.Consider all 25 unbranched catacondensed isomers of hexacene. 32Let us calculate the HOMA B index of the entire structures rather than of the individual rings (Figure 4).This can be done using expressions 1 or 7, in which now n denotes the number of all CC bonds in the structure instead of only in the single ring.The HOMA B values range from 0.76 for benzo[c]picene to 0.60 for hexacene (Figure 4).Notably, the former is the most stable, and the latter is the least stable.Also, the former has the largest, while the latter has the smallest Kekulécount 33 (Figure 4).
The stability of the unbranched catacondensed isomers of hexacene is biased by a steric hindrance, overcrowding, associated with a specific annelation leading to nonplanarity of the structure. 32Their symmetry and planarity are reflected in the PAH 3-digit code. 34The code used here assigns numbers 0, 1, and 2 to the nonterminal rings depending on annelation (number of H atoms in the edge, where the "clockwise" and "counterclockwise" angular annelations are differentiated by numbers 0 and 2, respectively).For instance, hexacene, benzo[c]picene, and hexahelicene are coded by 1111, 0202, and 0000, respectively (Figure 4).Notice that the original Balaban's 3-digit code 34 defines differently: 0 for linear annulations, and 1 or 2 for angular annulation.Still, the two 3digit codes are equivalent.
A large number of nonplanar structures of the C 2 or C 1 point group symmetry (00 or 22 sequence present in the code) cause that correlation between the stabilization energy and molecular descriptors such as the Kekulécount occurs only after the addition of a term associated with the structure deviation from planarity. 32The Kekulécount K is the number of different perfect matchings of the structure with alternating single and double bonds (as in the Kekuléformula of benzene).Within the resonance theory stemming from the Huckel approximation of the molecular orbitals theory, the Kekulécount allows calculating the resonance energy (RE): RE = A•ln(K), where coefficient A = 0.1185 (eV). 33Still, the resonance theory of the conjugated polyhex hydrocarbons refers to planar structures. 35,36Consequently, the correlation of the total energy difference and the logarithm of the Kekulécount occurs only for planar hexacene isomers, while the nonplanar ones deviate substantially (Figure 5a).
However, the HOMA B index calculated for the entire molecules of the unbranched catacondensed isomers of hexacene performs as well as the Kekulécount (Figure 5b).Indeed, there is a significant logarithmic correlation between HOMA B and the total energy differences for planar structures (R 2 = 0.948, Figure 5b).As for the Kekulécount, energies of the nonplanar ones deviate from it.Moreover, the linear correlation between the calculated HOMA B index and the Kekulécount for all 25 structures is satisfactorily significant (R 2 = 0.906, Figure 5c), and if only planar ones are considered, the correlation coefficient increases to R 2 = 0.966 (Figure 5c).Finally, let us suggest that the HOMA B index calculated for the entire molecules of the hexacene isomers can be understood as a similarity to a hypothetical isomer that would consist of six perfectly aromatic benzene rings.Let us also stress that, similarly to the classical HOMA index, the index can be applied to molecules composed of different numbers of rings and fused differently.So, it can be applied to how the topological indices are used in the Chemical Graph Theory.In fact, we have already demonstrated that the HOMA index is a topological index of the Structural Formula version of the Graph Theory. 37

FINAL REMARKS
4.1.HOMA Index for Molecules Containing Heteroatoms.A significant strength of the HOMA index is its ability to be generalized and parametrized for heteroatoms 3,7−11 as follows where j = CC, CN or NC, CO or ON, ..., NN, etc., and thus k denotes the number of different (unordered) pairs of atoms in the cycle, α j is a normalization factor for the j-th type of bonds, n j is the number of the j-th kind of bonds, R i j is the i-th bond of j-th

Journal of Chemical Information and Modeling
type, and R opt j is the optimal bond for the j-th type.The α j and R opt j parameters are found in the parametrization procedures referred to the experimental technique or computational level at which distances in a given ring are found or calculated.
Expression 12 can be rewritten to a more convenient form as follows where j indexes the kind of bond, k is the number of differentiated bonds, α j is the j-th normalization factor, x ̅ j and y ̅ j are the means of all bonds of the j-th type and the appropriate j-th reference, respectively, and in the variance symbol, the n j number of bonds of j-th-type is hidden.The HOMA index defined in eqs 12 and 13 satisfies symmetricity and has a maximum axioms (8) and thus is a "chemical" similarity function as the primary HOMA index does.
4.2.Use of the HOMA Formula to Express Similarity Other Than Geometrical.The HOMA index was already studied as a function of the electron density properties in bond critical points (BCPs). 12,13In such a case, HOMA expressed similarity between electron density properties (electron density, potential, and kinetic forms of energy) measured in BCPs and was not a geometrical parameter anymore.However, the HOMA expression allows an even further deviation from the original intention.We can assume that x is an atom property such as partial charge, spin, or chemical shift.Then, y would be the appropriate property in the reference atom, and Var(x) would be the variance of this parameter in the examined moiety such as a ring.Since, for years, the 13 C NMR spectra have been used for determining the aromaticity of the compounds, the comparison between the NICS class of indices and the HOMAbased NMR with the benzene C atom chemical shift as a reference would be especially intriguing.Such comparison is even more important because the chemical shift is observable.In contrast, most of the NICS indices are not, and only recently was it proven that the integral NICS index introduced by Stanger, 38 INICS, 39,40 was physically justified through its relation to the (a) Linear correlation between the total energy difference and the natural logarithm of the Kekulécount ln(K) for 25 unbranched catacondensed isomers of hexacene calculated at the B3LYP/D3/6-31G** level.(b) Logarithmic correlation between the total energy difference and the HOMA B index calculated for entire molecules and (c) the linear correlations between the HOMA B index calculated for entire molecules and the Kekulécount for all 25 structures (in blue) and only for the planar ones (in red).Blue empty cycles correspond to the nonplanar molecules (Figure 4).ring current via Ampere−Maxwell's law as demonstrated by Berger et al. 41,42 4.3.Use of the HOMA Formula to Study Acyclic Molecules.The HOMA index revealed both an increase in delocalization in polyene-like structures and an increase in insulation in unsaturated hydrocarbons, 12,13,16 and a correlation with the boiling point in n-alkanes. 17However, juxtaposing the HOMA indices taken against ethane and ethene or acetylene references in analogy to that presented in Section 3.3.1 could uncover interesting new molecular features.

CONCLUSIONS
A close inspection of the definition of the HOMA aromaticity index revealed that it has the mathematical property of a similarity function.This property explains why the index, derived based on simple quantum chemical approximation, is such a good measure of geometrical aromaticity.It expresses a similarity to perfectly aromatic benzene.
An expression based on the transformation of the EN and GEO components of HOMA displays that to calculate the index, it is enough to know the arithmetical means of the CC bond lengths in the examined and reference rings and the variation of the CC bond lengths only in the examined moiety.
The similarity property of the HOMA function enables applying the index to evaluate similarity to other molecules, like cyclohexane.The HOMA-similarity to cyclohexane appeared to express a (relaxed)-saturicity not accompanied by strong structural strains or steric hindrances.A slight reformulation of the HOMA definition to allow the reference moiety to have some variation of the bond lengths, denoting a nonzero GEO factor, showed that similarity to the archetypical elusive cyclohexatriene ring can be better discriminated than using the classical HOMA index.
We demonstrated that HOMA can also be used to evaluate the properties of entire molecules, such as isomers of hexacene.The index calculated for all 25 unbranched catacondensed isomers of hexacene shows a significant quadratic correlation with the total energy differences of planar isomers from which the nonplanar ones deviate.Such an index could be interpreted as a similarity to a hypothetical isomer consisting of six perfectly aromatic benzene rings.Moreover, the HOMA index is significantly correlating with the Kekulécount (connected to the resonance energy within the frame of the Huckel approximation of the molecular orbitals theory) for all 25 isomers of hexacene (R 2 > 0.9), but if only planar ones are considered, the correlation is much stronger (R 2 > 0.96).

Data Availability Statement
All quantum chemical DFT calculations were done using commercially available Gaussian 09 software. 22The B3LYP functional, 18,19 along with the D3 Grimme correction for dispersion forces, 20 the 6-31G** Pople-type basis set, 21 are accessible directly in the Gaussian program.Calculations of the HOMA indices and their components were done using commercial Microsoft Excel program, while correlations were performed using commercial SigmaPlot for Windows ver.14. 43

Figure 2 .
Figure 2. Relationships between indices calculated using benzene (B) and chair cyclohexane (C) reference molecules.The HOMA indices (a) and the HOMA indices truncated to single EN (black points) and GEO factors (red points) (b).

Figure 4 .
Figure 4. Structural formulas of all unbranched catacondensed hexacene isomers.Total energy differences vs benzo[c]picene are in blue and parentheses (kcal/mol), HOMA B of the entire structure are in red, and the Kekulécount are in black and brackets.Point group symmetry and the PAH 3-digit code 34 are in the second row in black.

Figure 5 .
Figure 5. (a) Linear correlation between the total energy difference and the natural logarithm of the Kekulécount ln(K) for 25 unbranched catacondensed isomers of hexacene calculated at the B3LYP/D3/6-31G** level.(b) Logarithmic correlation between the total energy difference and the HOMA B index calculated for entire molecules and (c) the linear correlations between the HOMA B index calculated for entire molecules and the Kekulécount for all 25 structures (in blue) and only for the planar ones (in red).Blue empty cycles correspond to the nonplanar molecules (Figure4).

Table 2 .
HOMA B and HOMA CHT Δ Indices Calculated for the Molecules Shown in Figure3, Using, Respectively, Benzene (B) and Cyclohexatriene (CHT) Reference Molecules (See Table1and Eqs 10 and 11, Respectively) a