|Home | About | Journals | Submit | Contact Us | Français|
A new aromaticity definition is advanced as the compactness formulation through the ratio between atoms-in-molecule and orbital molecular facets of the same chemical reactivity property around the pre- and post-bonding stabilization limit, respectively. Geometrical reactivity index of polarizability was assumed as providing the benchmark aromaticity scale, since due to its observable character; with this occasion new Hydrogenic polarizability quantum formula that recovers the exact value of 4.5 a03 for Hydrogen is provided, where a0 is the Bohr radius; a polarizability based–aromaticity scale enables the introduction of five referential aromatic rules (Aroma 1 to 5 Rules). With the help of these aromatic rules, the aromaticity scales based on energetic reactivity indices of electronegativity and chemical hardness were computed and analyzed within the major semi-empirical and ab initio quantum chemical methods. Results show that chemical hardness based-aromaticity is in better agreement with polarizability based-aromaticity than the electronegativity-based aromaticity scale, while the most favorable computational environment appears to be the quantum semi-empirical for the first and quantum ab initio for the last of them, respectively.
In conceptual chemistry, the aromaticity notion stands as one of the main pillars of structural understanding of molecular stability and reactivity; it spans all relevant periods of modern chemistry, suffering continuous revivals as the development of the physico-mathematical methods allows [1–5].
Historically, the custom definition of aromaticity is that it characterizes the planar molecules with (4n + 2) π-electrons , favoring substitution while resisting addition reactions [7,8], corresponds with high stability respecting the anti-aromatic structure , and ultimately correlates with low diamagnetism or magnetic susceptibility [10,11]. However, through the discovery of the non-planar feature of the most aromatic structure—benzene , along the generalizations of the Hückel rule for conjugated systems leading with the conjugated circuits [13,14], topological conjugated structures [15–17], up to the aromatic zones of molecular fragments [18,19], the aromaticity concept appears as a quite versatile concept that still needs proper quantification.
The main difficulty in capping the observability character of aromaticity may reside in the fact it does not directly relate to the ground state energy of molecules, but rather with their excited and valence orbitals–a condition closely related with the π electrons delocalized in the structure. Such aromaticity paradox of quantifying molecular stability (that usually is conducted by a Variational algorithm towards the ground state) by means of frontier pi electrons, may be considered as the main origin for the most inconsistencies and controversial points relating to the aromaticity concept in general and of its scales’ realization in particular.
Still, advancing aromaticity scales that employ certain physico-chemical molecular properties are useful at least for understanding whether certain molecular properties are inter-related and whether they enter or not the aromaticity sphere of definition. For that, the aromaticity scales should be always judged and validated in a comparative manner through concluding merely qualitatively upon the benchmarking degree of inter-correlation. This way, the aromaticity concept works best within the “transitivity thinking”: if the scale A1 correlates with the scales A2 and A3 with the same degree then the properties underlying the last two scales should be as well correlated and they may be regarded as different faces of the same molecular property. In short, the aromaticity concept and especially through its scales has the role of ordering among the molecular properties in general and of reactivity indices in special.
Returning to the aromaticity definition, it seems that the two main routes for its quantitative evaluation are the energetic and geometric ways; although within a Variational approach they should be closely related, i.e., minimum global energy corresponds with the geometry optimization, since the aromaticity paradox described above the two sides of molecular structures open two different approaches for introducing quantitative indices of aromaticity. For instance, based on geometric (also extended to topological) criteria, the consecrated harmonic oscillator based-molecular aromaticity (HOMA) [20,21] and the recent topological paths and aromatic zones (TOPAZ) , as well as the ultimate topological index of reactivity (TIR)  indices describe in various extent the influence the nuclei motion, the molecular fragment conjugation, or the site with the maximum probability (entropy) in electrophilic substitution have on aromaticity viewed as increasing (for the first two indices) or decreasing (for the last index) as the more delocalized π-electrons in question, respectively. On the other hand, from the energetic perspective, the resonance energy (RE) or its version reported per concerned π-electrons (REPE) [22–24], together with the heat (enthalpy) of formation (as a thermodynamic stabilization criterion of energy) , give another alternative to quantify aromaticity rises with their increase (for the first two indices) and decrease (for the last index), respectively. Moreover, other groups of methods in aromaticity evaluation are developed by employing the molecular magnetic [17,26–30] as well as based on electron delocalization [31–37] properties. It is therefore clear that the aromaticity scale is neither unique in trend nor in quantification and deserves further geometrical-energetic comparative investigation.
In this context, wishing to provide a fresh geometrical vis-à-vis of energetic aromaticity discussion, the present paper introduces the atoms-in-molecules compactness form of aromaticity that is then specialized both at geometrical level though the polarizability information and within the energetic framework though the electronegativity and chemical hardness reactivity indices. Following this, they will be used for ordering ten basic organic compounds, aromatic annulens, amines, hydroxyarenes, and heterocycles with nitrogen, against the corresponding aromaticity within most common semi-empirical and ab initio quantum chemical methods. The present considerations and results aim to further clarify the relationship between the electronegativity and chemical hardness in modeling the molecular stability/reactivity/aromaticity, as well the “computational distance” among their output furnished by various quantum mechanical schemes used in structural chemistry. For all these, the aromaticity is involved both as the motif and the tool having overall the manifestly unifying character among the fundamental concepts and computational schemes of chemistry.
Modeling the chemical bond is certainty key for describing the chemical reactivity and molecular structures’ stability. Yet, since the chemical bond is a dynamic state, for the best assessment of its connection with the stability and reactivity, the pre- and post-bonding stages are naturally considered.
For the pre-bonding stage, the atomic spheres are considered in the atoms-in-molecule (AIM) arrangement, while for the post-bonding stage the molecular orbitals (MOL) of the already formed molecule are employed, see Figure 1; consequently, their ratio would model the compactness degree of a given property of AIM in respect to its counterpart at the MOL level of the chemical bond. Therefore, the actual compactness index of aromaticity and takes the general form
that becomes workable once the property Π is further specified. Note that for the Equation (1) to be properly implemented, the chemical property Π should be equally defined and with the same meaning for the atoms and molecules, for consistency; such that what is compared is the chemical manifestation of the same property of bonding in its pre- or post-stage of formation. In other terms, Equation (1) may be regarded as a kind of “chemical limit” for the chemical bond that may be slightly oriented towards its atom constituents or to its molecular orbitals prescribing therefore the propensity to reactivity or stability, respectively.
It is worth considering also the quantitative difference between AIM and MOL properties of bonding, in which case the result may be regarded as the first kind of absolute aromaticity–for this reason, it is evaluated essentially between pre- and post-bonding stages and not relative to a referential (different) molecule , while the present approach promotes the compactness version of the aromaticity as the measure associated with the molecular stability in analog manner the compactness of rigid spheres in unit cells provides the crystal stability orderings. Yet, the proper scale hierarchy of compactness aromaticity, i.e., the qualitative tendency respecting the quantitative yield of Equation (1), is to be established depending on the implemented chemical property. In what follows, both the geometrically- and energetically-based reactivity indices will be considered, and their associate AIM compactness aromaticity formula and scales formulated. Moreover, once various quantum methods in evaluating the MOL denominator property in Equation (1) are considered, they will become fully quantum.
Since it has been already shown  that the polarizability α of a conducting sphere of radius r is equal to r3, for the atomic dipole systems the induced perturbation on the electronic cloud the actual formula should be corrected as :
with K a dipole related constant that has to be set out.
Historically, while a direct expansion of the Schrödinger differential equation in powers of eigen-energy, firstly done by Waller and Epstein [41,42], gives an analytic value of the polarizability of 4.5 a.u., the same value was found by variational method by Hassé through employing the Hydrogen ground state modified wave-function 
In modern times of quantum mechanics, the exact static dipole polarizabilities for the excited S states of the Hydrogen atom are determined by using the reduced free-particle Green's function method developed by McDowell and Porter  with the general formula for the polarizabilities found to be :
The last two formulas give both the same celebrated Hydrogen static Polarizability of (9 / 2)a30, where a0 is the Bohr radius.
Yet, another Hydrogenic Polarizability formulation can be actually elegantly developed from the first principles of quantum mechanics; it looks like (see Appendix for the complete derivation):
that immediately recovers the Hydrogen exact limit:
Overall, from above discussion we can assume the universal atomic constant K = 4.5 for atoms, while the atoms-in-molecule polarizability may be written as the atoms in molecule superposition providing the contributing atomic radii is known:
The ansatz of summation of the atomic polarizabilities in molecular polarizability relays on the fact they associate with the deformation (or softness) property of electronic frontier distribution that is additive in overlapping phenomena.
On the other hand, for the molecular polarizability in post-bonding stage one may use the volume information:
to advance the molecular (MOL) working expression throughout the normalizing factor involving the number of valence electrons:
that specializes for aromatics to the number of pi-electrons:
With atoms-in-molecule pre-bonding and the post-bonding molecular polarizabilities the related aromaticity geometrically based index may be constructed as their ratio:
The aromaticity scale is set upon the polarizability relation with the deformability power describing the molecular stability; as such, the higher MOL-polarizability respecting the AIM counterpart, the more flexible is the post-bonding molecular system against the external influences; consequently, as APOL decreases molecular stability increases.
In the same line as we proceeded with polarizability, we now set the AIM and MOL versions of energy based electronegativity and chemical hardness reactivity indices.
where the total atoms in molecule nAIM is the sum of the nA atoms of each A-species present in the molecule:
On the other hand, the electronegativity in the post-bonding stage of molecular state may be formulated through employing its relation with the total energy of a given system or state with N (valence) electrons [50–55], followed by the successive transformations:
by means of central difference approximation, the ionization potential and electronic affinity specializations:
while ending with considering the Koopmans’ frozen core approximation 
allowing therefore writing the MOL electronegativity in terms of highest occupied (HOMO) and lowest unoccupied (LUMO) frontier highest occupied and lowest unoccupied molecular orbitals, respectively.
The two forms of electronegativities, given by Equations (13) and (14), are next combined, according with the AIM-MOL aromaticity recipe of Equation (1), to provide the compactness electronegativity-based aromaticity index:
Now, the electronegativity reactivity principle under the form of electronegativity equalization in molecule [57,58] is used for establishing the behavior of the aromaticity scale based on Equation (17). Accordingly, though aromaticity is described as the ratio of AIM to MOL electronegativity, it is clear that the pre-bonding stage of atomic electronegativity equalization (electronic flowing) into the molecular unified orbitals is the dominant phenomena, so that it is expected to prevail. Therefore, the aromaticity index of electronegativity Equation (17) is higher as the molecular stability (formation) is better realized; in short, as AMOL increases, a more aromatic molecular system is assumed.
For chemical hardness description, the AIM pre-bonding formulation happens to have the same analytical form as that found for the AIM electronegativity :
while the MOL version is constructed based on the previous orbital energy prescriptions of Equations (15) and (16) as applied to the second order derivation of the total energy respecting to the total number of electrons in a given state towards the working HOMO-LUMO formulation [60–62]:
It nevertheless resembles the idea that as the frontier gap between the HOMO and the LUMO orbitals is larger as the molecular system is more stable, i.e., less engaged into chemical reactions through its frontier electrons [63–66].
with the scale trend fixed by the maximum hardness principle [59,67], abstracted from above MOL chemical hardness; in other words, the chemical reactivity described by chemical hardness is driven by the post-bonding stage in molecular formation requiring that a molecules is as stable as its HOMO-LUMO gap increases. Consequently, the aromaticity scale based on Equation (20) is arranged from the lowest to highest values that parallels the increasing reactivity and decreasing aromaticity; in short, smaller AHard, bigger aromaticity character for a molecular system.
However, while atomic electronegativity and chemical hardness in evaluating AIM schemes of Equations (13) and (18) may be implemented by appealing various benchmarking scales [53,54,68], the molecular orbital energies in computing MOL counterparts Equations (14) and (19) require dedicated computations for each concerned molecule; as such, for better understanding and interpreting the obtained electronegativity- and chemical hardness-based aromaticity scales it is worth shortly reviewing the main quantum schemes mainly used in computing the (post-bonding) molecular spectra.
Following the Dirac’s quote, once the Schrodinger equation:
was established “The underlying physical laws necessary for the mathematical theory of a large part of physics and the whole of chemistry are thus completely known” .
Unfortunately, the molecular spectra based on the eigen-problem Equation (21) is neither directly nor completely solved without specific atoms-in-molecule and/or symmetry constraints and approximation. As such at the mono-electronic level of approximation the Schrodinger Equation (21) rewrites under the so called independent-electron problem:
with the aid of effective electron Hamiltonian partitioning:
and the correspondent molecular monoelectronic wave-functions (orbitals) fulfilling the conservation rule of probability:
However, when written as a linear combination over the atomic orbitals the resulted MO-LCAO wave-function:
having the diagonal energy-matrix elements as the eigen-solution
to be found in terms of the expansion coefficients matrix (C), the matrix of the Hamiltonian elements:
and the matrix of the (atomic) overlapping integrals:
Yet, the solution of the matrix equation (26) may be unfolded through the Löwdin orthogonalization procedure [70,71], involving the diagonalization of the overlap matrix by means of a given unitary matrix (U), (U)+(U) =(1), by the resumed procedure:
However, the solution given by Equation (27b) is based on the form of effective independent-electron Hamiltonians that can be quite empirically constructed–as in Extended Hückel Theory ; such “arbitrariness” can be nevertheless avoided by the so called self-consistent field (SCF) in which the one-electron effective Hamiltonian is considered such that to depend by the solution of Equation (25) itself, i.e., by the matrix of coefficients (C); the resulted “Hamiltonian” is called Fock operator, while the associated eigen-problem is consecrated as the Hartree-Fock equation:
In matrix representation Equation (31) looks like:
that may be iteratively solved through diagonalization procedure starting from an input (C) matrix or–more physically appealing–from a starting electronic distribution quantified by the density matrix:
with major influence on the Fock matrix elements:
Note that now the one-electron Hamiltonian effective matrix components HμV differ from those of Equation (28) in what they truly represent, this time the kinetic energy plus the interaction of a single electron with the core electrons around all the present nuclei. The other integrals appearing in Equation (34) are generally called the two-electrons-multi-centers integrals and are written as:
From definition (35), there is immediate to recognize the special integral J = (μμ|νν as the Coulomb integral describing repulsion between two electrons with probabilities μ2 and ν2.
Moreover, the Hartree-Fock Equation (32) with implementations given by Equations (33)–(34) are known as Roothaan equations  and constitute the basics for closed-shell (or restricted Hartree-Fock, RHF) molecular orbitals calculations. Their extension to the spin effects provides the equations for the open-shell (or unrestricted Hartree-Fock, UHF) known also as the Pople-Nesbet Unrestricted equations .
The second level of approximation in molecular orbital computations regards the various ways the Fock matrix elements of Equation (34) are considered, namely the approximations of the integrals (35) and of the effective one-electron Hamiltonian matrix elements HμV.
The main route for such an endeavor is undertaken through neglecting at different degrees certain differential overlapping terms (integrals)–as an offset ansatz–although with limited physical justification–while the adjustment with experiment is done (post-factum) by fitting parameters–from where the semi-empirical name of such approximation. Practically, by emphasizing the (nuclear) centers in the electronic overlapping integral (29):
the differential overlap approximation may be considered by two situations.
For both groups of approximations specific methods are outlined below.
The basic NDO approximation was developed by Pople and is known as the Complete Neglect of Differential Overlap CNDO semi-empirical method [75–78]. It employs the approximation (37) such that the molecular rotational invariance is respected through the requirement the integral (37c) depends only on the atoms A or B where the involved orbitals reside–and not by the orbitals themselves. That is the integral γAB in (37b) is seen as the average electrostatic repulsion between an electron in any orbital of A and an electron in any orbital of B:
In these conditions, the working Fock matrix elements of Equation (34) become within RHF scheme:
From Equations (40a&b) follows there appears that the core Hamiltonian has as well the diagonal and off-diagonal components; the diagonal one represents the energy of an electron in an atomic orbital of an atom (say A) written in terms of ionization potential and electron affinity of that atom :
added to the attraction energy respecting the other (B) atoms to produce the one-center-one-electron integrals:
overall expressing the energy an electron in the atomic orbital μ would have if all other valence electrons were removed to infinity. The non-diagonal terms (the resonance integrals) are parameterized in respecting the overlap integral and accounts (through βAB parameter averaged over the atoms involved) on the diatomic bonding involved in overlapping:
The switch to the UHF may be eventually done through implementing the spin equivalence:
although the spin effects are not at all considered since no exchange integral involved. This is in fact the weak point of the CNDO scheme and it is to be slightly improved by the next Semi-empirical methods.
The exchange effect due to the electronic spin accounted within the Intermediate Neglect of Differential Overlap (INDO) method  through considering in Equations (40a) and (41a) the exchange one-center integrals γAA K = (μν|μν) is evaluated as:
in terms of the Slater-Condon parameters G1, F2, … usually used to describe atomic spectra.
The INDO method may be further modified in parameterization of the spin effects as developed by Dewar’s group and led with the Modified Intermediate Neglect of Differential Overlap (MINDO) method [75,81–89] whose basic equations look like:
Apart from specific counting of spin effects, another particularity of MINDO respecting the CNDO/INDO is that all the non-zero two-center Coulomb integrals are set equal and parameterized by the appropriate one center two electrons integrals AA and AB within the Ohno-Klopman expression [90,91]:
The one-center-one-electron integral Hμμ is preserved from the CNDO/INDO scheme of computation, while the resonance integral (41c) is modified as follows:
with the parameter being now dependent on the atoms-in-pair rather than the average of atomic pair involved. As in INDO, the exchange terms, i.e., the one-center-two-electron integrals, are computed employing the atomic spectra and the Gk, Fk, Slater-Condon parameters, see Equation (42) . Finally, it is worth mentioning that the MINDO (also with its MINDO/3 version) improves upon the CNDO and INDO the molecular geometries, heats of formation, being particularly suited for dealing with molecules containing heteroatoms.
This second group of neglecting differential overlaps semi-empirical methods includes along the interaction quantified by the overlap of two orbitals centered on the same atom also the overlap of two orbitals belonging to different atoms. It is manly based on the Modified Neglect of Diatomic Overlap (MNDO) approximation of the Fock matrix, while introducing further types of integrals in the UHF framework [93–99]:
Note that similar expressions can be immediately written within RHF once simply replacing:
in above Fock (46a&b) expressions.
Now, regarding the (Coulombic) two-center-two-electron integrals of type (38c) appearing in Equations (46) there were indentified 22 unique forms for each pair of non-hydrogen atoms, i.e., the rotational invariant 21 integrals (ss|ss), (ss|pσ pσ), (ss|pπ pπ),…,(pσ pσ|pσ pσ), (pπ pπ|pπ pπ),…, (spσ|spσ), (spπ|spπ),…, (pπ pσ|spπ), (pπ pσ|pπ pσ), and the 22nd one that is written as a combination of two of previously ones, namely (pπ pπ′|pπ pπ′) = 0.5[(pπ pπ|pπ pπ) – (pπ pπ|pπ′ pπ′)], with the typical integral approximation relaying on the Equation (44) structure, however slightly modified as:
where additional parameters cA and cB represent the distances of the multipole charges from their respective nuclei. The MNDO one-center one-electron integral has the same form as in NDO methods, i.e., given by Equation (41b) with the average potential of Equation (39) acting on concerned center; still, the resonance integral is modified as:
containing the atomic adjustable parameters and for the orbitals μ and ν of the atoms A and B, respectively. The exchange (one-center-two-electron) integrals are mostly obtained from data on isolated atoms . Basically, MNDO improves MINDO through the additional integrals considered the molecular properties such as the heats of formations, geometries, dipole moments, HOMO and LUMO energies, etc., while problems still remaining with four-member rings (too stable), hypervalent compounds (too unstable) in general, and predicting out-of-plane nitro group in nitrobenzene and too short bond length (~ 0.17 Å) in peroxide–for specific molecules.
while correcting the one-center-two-electron atomic integrals of Equation (44) by the specific (AM) monopole-monopole interaction parameters. In the same line, the nuclei-electronic charges interaction adds an energetic correction within the αAB parameterized form:
The AM1 scheme, while furnishing better results than MNDO for some classes of molecules (e.g., for phosphorous compounds), still provides inaccurate modeling of phosphorous-oxygen bonds, too positive energy of nitro compounds, while the peroxide bond is still too short. In many case the reparameterization of AM1 under the Stewart’s PM3 model [103,104] is helpful since it is based on a wider plethora of experimental data fitting with molecular properties. The best use of PM3 method lays in the organic chemistry applications.
To systematically implement the transition metal orbitals in semi-empirical methods the INDO method is augmented by Zerner’s group either with non-spectroscopic and spectroscopic (i.e., fitting with UV spectra) parameterization [105–107], known as ZINDO/1 and ZINDO/S methods, respectively [108–114]. The working equations are formally the same as those for INDO except for the energy of an atomic electron of Equation (41a) that now uses only the ionization potential instead of electronegativity of the concerned electron. Moreover, for ZINDO/S the core Hamiltonian elements Hμμ is corrected:
by the fr parameterized integrals:
in terms of the one-center-two-electron Coulomb integrals , . Equation (53) conserves nevertheless the molecular rotational invariance through making the difference between the s- and d-Slater orbitals exponents. The same types of integrals correct also the nuclei-electronic interaction energy by quantity:
Since based on fitting with spectroscopic transitions the ZINDO methods are recommended in conjunction with single point calculation and not with geometry optimization, this should be consider by other off-set algorithms.
Beyond either NDO or NDDO methods, the self-consistent computation of molecular orbitals can be made by the so called ab initio approach, directly relaying on the HF equation or on its density functional extension, as will be in next sketched.
The alternative to semi-empirical methods is the full self-consistent calculation or the so called ab initio approach; it is based on computing of all integrals appearing on Equation (34), yet with the atomic Slater type orbitals (STO), exp(−αr), being replaced by the Gaussian type orbitals (GTO) :
in molecular orbitals expansion–a procedure allowing for much simplification in multi-center integrals computation. Nevertheless, at their turn, each GTO may be generalized to a contracted expression constructed upon the primitive expressions of Equation (55a):
where dpμ and αA are called the exponents and the contraction coefficients of the primitives, respectively. Note that the primitive Gaussians involved may be chosen as approximate Slater functions , Hartree-Fock atomic orbitals , or any other set of functions desired so that the computations become faster. In these conditions, a minimal basis set may be constructed with one function for H and He, five functions for Li to Ne, nine functions for Na to Ar, 13 functions for K and Ca, 18 functions for Sc to Kr,..., etc., to describe the core and valence occupancies of atoms [118–120]. Although such basis does not generally provide accurate results (because of its small cardinal), it contains the essential information regarding the chemical bond and may be useful for qualitative studies, as is the present case for aromaticity scales where the comparative trend is studied.
When simple ab initio method is referred it means that the Hartree-Fock equation (31) with full Fock matrix elements [121–124] of Equations (33) and (34) is solved for a Gaussian contracted basis (55). Actually, the method evaluates iteratively the kinetic energy and nuclear-electron attraction energy integrals–for the effective Hamiltonian, along the overlap and electron-electron repulsion energy integrals (for both the Coulomb and exchange terms), respectively written as:
until the consistency in electronic population of Equation (33) between two consecutive steps is achieved.
Note that such calculation assumes the total wave function as a single Slater determinant, while the resultant molecular orbital is described as a linear combination of the atomic orbital basis functions (MO-LCAO). Multiple Slater determinants in MO description projects the configurationally and post-HF methods, and will not be discussed here.
The main weakness of the Hartree-Fock method, namely the lack in correlation energy, is ingeniously restored by the Density Functional method through introducing of the so called effective one-electron exchange-correlation potential, yet with the price of not knowing its analytical form. However, the working equations have the simplicity of the HF ones, while replacing the exchange term in Equation (34) by the exchange-correlation (“XC”) contribution; there resulted the (general) unrestricted matrix form of the Kohn-Sham equations :
in a similar fashion with the Pople-Nesbet equations of Hartree-Fock theory. The restricted (closed-shell) variant is resembled by the density constraint:
in which case the Roothaan analogous equations (for exchange-correlation potential) are obtained.
that can be actually regarded as the solution of the Kohn-Sham equations themselves. The appeared exchange-correlation energy EXC may be at its turn conveniently expressed through the energy density (per unit volume) by the integral formulation:
once the Fock elements of exchange-correlation are recognized to be of density gradient form :
The quest for various approximations for the exchange-correlation energy density f(ρ) had spanned the last decades in quantum chemistry, and was recently reviewed . Here we will thus present the “red line” of its implementation as will be further used for the current aromaticity applications. The benchmark density functional stands the Slater exchange approximation, derived within the so called Xα theory :
with the α taking the values:
Considerably improvement for molecular calculation was given by Becke’s density gradient correction of the local spin density (or Slater exchange) approximation of the exchange energy :
with the parameters a = (3/2)(3/4π)1/3 and b = 0.0042 chosen to fit the experiment. Other exchange functionals were developed along the same line, i.e., having different realization of the gradient function (66), most notable being those of Perdew and collaborators (e.g., Perdew-Wang-91, PW91) .
The correlation contribution was developed on a somewhat different algorithm, namely employing its definition as the difference between the exact and Hartree-Fock (HF) total energy of a polyelectronic system . Without reproducing the results (more detailed are provided in the dedicated review of Ref. ), for the actual purpose we mention only the Lee-Yang-Parr (LYP) correlation functional [131–133] along the Vosko-Wilk-Nusair (VWN) local correlation density functional .
However, the exchange and correlation density functionals combine into the so called hybrid functionals; those used in the present study refer to:
These are the main methods, at both conceptual and computational levels, to be in next used to asses and compare the atoms-in-molecule compactness aromaticity scales for basics organics.
The above reactivity indices-based aromaticity scales are now computed within the presented quantum chemical schemes for a limited yet significant series of benzenoids (see Table 2) containing the “life” atoms of Table 1. The atoms-in-molecule of aromaticity scales of polarizability, electronegativity and chemical hardness of Equations (12), (17), and (20) are directly computed upon the formulas given by Equations (8), (13), and (18), respectively; they are based exclusively on the data of Table 1 with the AIM results listed in Table 2, the 5th, 9th, and 10th columns, respectively.
For the post-bonding evaluations of the same indices, one must note the special case of polarizability that is computed upon the general Equation (11)–thus involving the molecular volume pre-computation. Here it is worth commenting on the fact that one can directly compute the molecular polarizability in various quantum schemes–however, with the deficiency that such procedure does not distinguishes among the stereo-isomers, i.e., molecules VII (1-Naphthol) and VIII (2-Naphthalelon); IX (2-Naphthalenamine) and X (1-Naphthalenamine) in Table 2, since furnishing the same values, respectively; instead the same quantum scheme is able to distinguish between the volumes of two stereo-isomers making the Equation (11) as a more general approach. This way, the molecular volumes are reported in the 6th column of Table 2 as computed within the ab initio–Hartree Fock method of Section 188.8.131.52; note that the HF method was chosen as the reference since it is at the “middle computational distance” between the semi-empirical and density functional methods; it has only the correlation correction missing; however, even the density functional schemes, although encompassing in principle correlation along the exchange–introduces approximations on the last quantum effect. Therefore, the molecular polarizability is computed upon the Equation (11) in the 7th column of Table 2 with the associate polarizability compactness aromaticities displayed in the 8th column of Table 2.
The molecular energetic reactivity indices of electronegativity and chemical hardness are computed upon the Equations (14) and (19) in terms of HOMO and LUMO energies computed within the quantum semi-empirical and ab initio methods presented in Section 2.3; their individual values as well as the resulted quantum compactness aromaticities, when combined with the AIM values of Table 2, in Equations (17) and (20) are systematically communicated in Tables 3 and and4,4, with adequate scaled representations in Figures 2 and and3,3, respectively.
Note that neither the minimal basis set (STO-3G) nor the single point computation frameworks, although both motivated in the present context in which only the bonding and the post-bonding information should be capped in computation, does not affect the foregoing discussion by two main reasons: (i) they have been equally applied for all molecules considered in all quantum methods’ combinations; and (ii) what is envisaged here is the aromaticity trend, i.e., the intra- and inter- scales comparisons rather than the most accurate values since no exact or experimental counterpart available for aromaticity.
Now, because of the observational quantum character of polarizability, one naturally assumes the (geometric) polarizability based- aromaticity scale of Table 2 as that furnishing the actual standard ordering among the considered molecules in accordance with the rule associated with Equation (12); it features the following newly introduced rules along possible generalizations
Under the reserve that these rules and their generalizations should be verified by extra studies upon a larger set of benzenoid aromatics, we will adopt them here in order to analyze their fulfillment with the energetically-based aromaticity scales of electronegativity and chemical hardness, reported in Tables 3 and and44 and drawn in Figures 2 and and3;3; actually, their behavior is analyzed against the aromaticity ordering rules given by Equations (17) and (20), i.e., as being anti-parallel and parallel with the polarizability-based aromaticity trend of Equation (12), with the results systematized in Tables 5 and and6,6, respectively.
From Table 5 there follows that electronegativity based-aromaticity displays the following properties respecting the aromaticity rules derived from polarizability framework:
Overall, it seems electronegativity may be used in modeling compactness of atoms-in-molecules aromaticity–basically without counting on the exchange or correlation effects, or at best within the HF algorithm, while semi-empirical methods seems not well adequate. Yet, for all aromaticity rules formulated, there exists at least one quantum computational environment for which the electronegativity based compactness aromaticity is in agree with each of them.
The situation changes significantly when chemical hardness is considered for compactness aromaticity computation; the specific behavior is abstracted from the analysis of Table 6 and can be summarized as follows:
Overall, when the chemical hardness agrees with one of the above enounced Aroma Rules it does that within more than one computational scheme; however, the best agreement of chemical hardness with polarizability based- aromaticity scales is for the mono- versus bi- (and possible poly-) benzenic rings decreasing of aromaticity orderings, along the manifestly hydroxyl- superior effects in aromaticity than amine- groups substitution within most of the computational quantum schemes, i.e., valid either for semi-empirical and ab initio methods. The stereoisomerism is not covered by chemical hardness modeling aromaticity, and along the electronegativity limited coverage within HF scheme in Table 5, there follows that the energetic reactive indices are not able to prevail over the geometric indices as polarizability or to predict stereoisomerism ordering in aromaticity modeling compactness schemes.
Finally, few words about the output of the various quantum computational schemes respecting the current aromaticity definition given by Equation (1) are worth addressing. As such, one finds that:
The bottom line is rising by the wish to globally combine the ideas of quantum chemical methods used in chemical bonding, reactivity principles, and aromaticity results; upon the above discussions it follows that MINDO3, AM1 (or PM3)–for semi-empirical along Becke hybrid functionals and Hartree-Fock–for ab initio are the suited methods that fulfils most of the reactivity and the present introduced aromaticity bonding rules. However, the best of them overall seems to remain the consecrated HF scheme, since acquiring the highest number of grades summated throughout Tables 5 and and6.6. As such, a new challenge appears since the present results recommend that correlation does not count too much in aromaticity or reactivity modeling. Nevertheless, further studies with larger set of molecules and types of aromatics should be address for testing whether or not the advanced aromaticity (Aroma1–5) rules are preserved or in which degree they may be generalized or modified such that being in accordance with the principles of chemical bonding and reactivity.
Modeling the stability and reactivity of molecules is perhaps the greatest challenge in theoretical and computational chemistry. This is because the main conceptual tools developed as the reactivity indices of electronegativity and chemical hardness along the associate principles are often suspected by the lack of observability character. Therefore, although very useful in formal explanations of chemical bonding and reactivity, it is hard to find their experimental counterpart unless expressed by related measurable quantities as energy, polarizability, refractivity, etc. When the aromaticity concept come into play, it seems no further conceptual clarification is acquired, since no quantum observable or further precise definition can be advanced; in fact, the aromaticity concept associates either with geometrical, energetic, topologic, electronic molecular circuits (currents), or with the less favored entropic site in a molecule, just to name few of its representations.
However, since at the end, the aromaticity appears to describe the stability character of the molecular sample, its connection with a reactivity index seems natural, although systematically ignored so far. In this respect, the present work focuses on how the electronegativity and chemical hardness based-aromaticity scales behave with respect to others constructed on a direct observable quantum quantity–the polarizability in this case. This is because the polarizability quantity is fundamental in quantum mechanics and usually associated with the second order Stark effect that can be computed within the perturbation theory (see Appendix). Then, two ways of seeing a molecular structure were employed in introducing the actual absolute aromaticity definition:
The two molecular perspectives may be associated with the pre- and post-bonding stages of a chemical bond at equilibrium; therefore, the conceptual and computational competition between these two molecular facets should measure the stability or its contrary effect - the reactivity propensity - being therefore the ideal ingredients for an absolute definition of aromaticity. Note that although an AIM-to-MOL difference definition of absolute aromaticity was recently advanced , the actual study of their ratio definition should account for a sort of compactness degree of molecular structure–as described by the specific molecular property used.
In short, for a molecular property to become a candidate for absolute (here with its compactness variant of) aromaticity, it has to fulfill two basic conditions:
In this respect, all the presently considered reactivity indices, i.e., polarizability, electronegativity, and chemical hardness, have equally consecrated quantum definitions as well as atomic and molecular representations [55,139].
At the atomic level, the experimental values based on the ionization potential and electron affinity definitions for electronegativity and chemical hardness were considered, see Equations (14) and (19), respectively, while for the polarizability new Hydrogenic quantum formulation was provided by Equation (6), and in Appendix by Equation (A22), recovering the exact value for the Hydrogen system by Equation (7), in close agreement with other available atomic quantum formulations of Equations (4) and (5). Nevertheless, the AIM level was formed by appropriate averaging of atoms-in-molecule summation for each of the considered reactivity indices, see Equations (8), (13) and (18), and along of their MOL counterparts of Equations (11), (14), and (19) the polarizability-, electronegativity- and chemical hardness- based aromaticity definitions were formulated with the associate qualitative trends established by Equations (12), (17), and (20), respectively. Yet, for MOL level of computations, all major quantum chemical methods for orbital spectra computation were considered, in Section 2.3, and implemented in the current application for some basics aromatics in Section 2. Because of the quantum observable character of polarizability the related aromaticity scale was considered as benchmark for actual study and it offered the possibility of formulating five rules for aromaticity:
These rules are then checked for electronegativity and chemical hardness derived-aromaticity scales with the synopsis of the results in Tables 5 and and6.6. It followed that chemical hardness, although generally in better agreement with these rules for most of the quantum chemical methods considered for its MOL computation, may not be considered infallible against aromaticity, at least for the reason it does not fulfils at all with the Aroms5 rule above. Surprisingly, chemical hardness index is more suited in modeling aromaticity when considered within semi-empirical computational framework, while the electronegativity responds better in conjunction with ab initio methods.
From quantum computational perspective, the consecrated HF method seems to get more marks in fulfillment of above Aroma1-to-5 rules, cumulated for electronegativity and chemical hardness based-aromaticity scales; it leads with the important idea the correlation effects are not determinant in aromaticity phenomenology, an idea confirmed also by the fact the density functional without exchange and correlation produces not-negligible fits with Aroma1, 2, and 4 rules in electronegativity framework.
Overall, few basic ideas in computing aromaticity should be finally emphasized
Future quests should enlarge the basis of the present conclusions by performing comparative aromaticity studies at the level of biomolecules and nanostructures; at the end of the day, the aromaticity concept in general and with its particular specialization should represent just a tool/vehicle in modeling and understanding the chemical bond of atoms in molecules and nanostructures, either in isolated or interacting states.
Useful discussions on historical aspects of aromaticity concept and on the present studied molecules with Dr. Delia Isac, Associate Professor for Didactic Chemistry at West University of Timisoara, now retired, are kindly thanked. As well, the editors and the associate team of International Journal of Molecular Sciences are very much thanked for supporting the open access publishing of highly rated scientific works in general, and of this paper in special. Author is grateful also to Brietta Pike of MDPI for final English editing of the manuscript.
Starting from the consecrated second order perturbation energy :
is specialized for the Stark potential produced by the applied external electric field with the amplitude in the 0x direction:
under the form:
that allows for α– polarizability in (A3) the general hydrogenoid (Z-dependent) formulation
is the reduced squared elementary charge
Now, to evaluate the atomic polarizability in terms of the quantum basic information contained within the atomic quantum numbers (e.g., n, k), one starts recognizing the general operatorial identity over the complete set of quantum (eigen) states:
On the other hand, recalling the basic quantum commutation rule of momentum with space coordinate:
along the companion energy-coordinate commutator:
there can be inferred the quantum relationship:
upon inserting of the above quantum closure relation over the complete set of eigen-states. The first term in the right-hand side of the last expression may be reformulated as:
and along the similar relation that springs out from the second term in (A10) one gets the equation:
Remarkably, the expansion (A13) may be also obtained by requiring that the Kramers-Heisenberg dispersion relation reduce to the Thomas scattering formula at high energies; indeed, through rewriting Equation (A13) in the form:
Now, returning to the evaluation of polarizability given by (A4) one can use the recipe (A13) to facilitate the skipping out of the energy-singularity towards the all-eigen-state summation (A7) with the successive results:
where we recognized the first Bohr radius expression:
Finally, the obtained expression (A15) is unfolded through replacing the coordinate observation with the atomic radius quantum average displacement respecting its instantaneous value:
It allows the immediate formation of the squared coordinate expression:
of which the observed quantum average looks like
The replacement of Equation (A19) in the polarizability (A15) produces its radial averages’ dependency:
Knowing the first and second order quantum averages for the atomic radius of a Hydrogenic system written in terms of the principal and azimuthal quantum numbers n and l, respectively :
the static atomic polarizability (A20) takes the analytical form:
recovering the exact result for the Hydrogen limiting case:
It is worth noting that the present derivation relays on the second order perturbation energy (A1) while the final expression (A22) is assumed to be exact through the Hydrogen checking case (A23), although different by the other reported also as valid formulations, see Equations (4) and (5) in the main text and References [45–47]. Nevertheless, the present atomic polarizability, either under expressions (A15) or (A22) is to be further tested for reliability in modeling of atomic (or ionic) and molecular systems.