Energetic and Geometric Characteristics of Substituents, Part 3: The Case of NO 2 and NH 2 Groups in Their Mono-Substituted Derivatives of Six-Membered Heterocycles

: Substituted heterocyclic arenes play important roles in biochemistry, catalysis, and in the design of functional materials. Exemplary six-membered heteroaromatic molecules, that differ from benzene by inclusion of one heteroatom, are pyridine, phosphorine, arsabenzene, and borabenzene. This theoretical study concerns the inﬂuence of the heteroatom present in these molecules on the properties of substituents of two types: electron-donating (ED) NH 2 group and electron-accepting (EA) NO 2 group, attached at the 2-, 3-, or 4-position. The effect is evaluated by the energy of interaction ( E rel ) between the substituent and the substituted system and electronic properties of the substituents described by the charge of the substituent active region (cSAR) index. In addition, several geometric descriptors of the substituent and heteroaromatic ring, as well as changes in the aromaticity, are considered. The latter are assessed using the Electron Density of Delocalized Bonds (EDDBs) property of delocalized π electrons. The obtained results show that the electronegativity (EN) of the heteroatom has a profound effect on the EA/ED properties of the substituents. This effect is also reﬂected in the geometry of studied molecules. The E rel parameter indicates that the relative stability of the molecules is highly related to the electronic interactions between the substituent and the heteroarene. This especially applies to the enhancement or weakening of π -resonance due to the EN of the heteroatom. Additionally, in the 2-heteroarene derivatives, speciﬁc through-space ortho interactions contribute to the heteroatom effects.


Introduction
Substituent effects are among the most important intramolecular interactions in organic chemistry and related fields. Traditionally, for a given series of substituted systems, the relations between the properties of the reaction center and the characteristics of the substituents are used for their analysis. The first quantitative approach to the description of the substituent effect was proposed by Hammett [1,2]. Most often, Hammett substituent constants or alike are used to characterize the substituent effect quantitatively [3]. It is also known that the properties of the substituent itself depend both on its position in the transmitting moiety and on the nature of the reaction site. The first case can be represented by the different substituent constants for the para and meta positions, while the substituent constants σ+ and σ-, characterizing the substituent effects in molecules with positively and negatively charged reaction sites, respectively, are an example of the second case. This is known as the reverse substituent effect [4]. For any substituent in any system, it can be quantified using the charge of the substituent active region (cSAR) concept [5,6].
The aim of this study is to determine to what extent the heteroatom in the sixmembered aromatic ring affects the stability of the substituted system, as well as the electronic and geometric properties of the strongly electron-attracting and electron-donating substituents, the nitro and amino groups, respectively. The following heteroarenes were selected as the studied substituted moieties: pyridine, phosphorine, arsabenzene, and borabenzene, shown in Figure 1. The atomic properties of heteroatoms are summarized in Table 1. It should be emphasized that each of these atoms may interact with the substituent differently, especially when it is in the ortho position. The nitrogen, phosphorus, and arsenic atoms have an electron pair in the plane of the ring, which can cause specific proximity interactions with the substituent. Boron, on the other hand, has an empty orbital in this plane, which is a Lewis acidic site open to interaction with a Lewis base. Additionally, As is the only atom in this group that has a filled 3d subshell.
Symmetry 2022, 14, x FOR PEER REVIEW 2 of 19 and geometric properties of the strongly electron-attracting and electron-donating substituents, the nitro and amino groups, respectively. The following heteroarenes were selected as the studied substituted moieties: pyridine, phosphorine, arsabenzene, and borabenzene, shown in Figure 1. The atomic properties of heteroatoms are summarized in Table 1. It should be emphasized that each of these atoms may interact with the substituent differently, especially when it is in the ortho position. The nitrogen, phosphorus, and arsenic atoms have an electron pair in the plane of the ring, which can cause specific proximity interactions with the substituent. Boron, on the other hand, has an empty orbital in this plane, which is a Lewis acidic site open to interaction with a Lewis base. Additionally, As is the only atom in this group that has a filled 3d subshell.  Nitrogen heteroarenes are extremely important in biochemistry, and they also occur in nucleic acid bases (pyrimidine and imidazole fragments) [12]. Many derivatives of Nheterocycles have found application in medicinal chemistry. For example, various quinoline derivatives are used in the treatment of malaria [13]; Imatinib, which contains pyridine and pyrimidine fragments, is used as a protein-targeted anticancer drug. Another important anticancer drug, Bortezomib, contains a substituted pyrazine ring [14]. Knowing the properties of substituents is an important aspect considered while designing such compounds. Nitro-, amino-, and chloro-derivatives of pyridine, as well as of N-heterocycles with more endocyclic N atoms, i.e., pyrimidine, pyrazine, and triazine, were the subject of Part 2 in this article series [15]. Phosphorine and arsabenzene were synthesized and characterized, for the first time, 50 years ago [16]. Since their discovery, they have been the subject of much experimental and theoretical research [17][18][19][20][21]. These heterocycles (C5H5E, E = element of group 15), like pyridine, can bind to metals via a lone pair at E (η 1 complexes) and/or via a π-electron sextet (η 6 complexes). Of course, the type of complex formed depends on the substituents attached to the ring. However, in the case of unsubstituted phosphabenzene and arsabenzene, the η 6 complexes predominate [22][23][24]. Phosphorine-based complexes and their use in homogeneous catalysis were summarized in a review by Le Floch [25]. Interestingly, it has been suggested that polycyclic aromatic hydrocarbons (PAHs) containing an endocyclic N atom (PANHs) may be responsible for the characteristic 6.22 μm feature in the emission spectra of some galaxies [26,27]. For this reason, astrophysicists predict that about 1 percent of the nitrogen in our Universe may exist in form of PANHs in an interstellar medium of galaxies. Similarly, the existence of P-containing PAHs in the interstellar medium (PAPHs) has been proposed and plausible  Nitrogen heteroarenes are extremely important in biochemistry, and they also occur in nucleic acid bases (pyrimidine and imidazole fragments) [12]. Many derivatives of N-heterocycles have found application in medicinal chemistry. For example, various quinoline derivatives are used in the treatment of malaria [13]; Imatinib, which contains pyridine and pyrimidine fragments, is used as a protein-targeted anticancer drug. Another important anticancer drug, Bortezomib, contains a substituted pyrazine ring [14]. Knowing the properties of substituents is an important aspect considered while designing such compounds. Nitro-, amino-, and chloro-derivatives of pyridine, as well as of N-heterocycles with more endocyclic N atoms, i.e., pyrimidine, pyrazine, and triazine, were the subject of Part 2 in this article series [15]. Phosphorine and arsabenzene were synthesized and characterized, for the first time, 50 years ago [16]. Since their discovery, they have been the subject of much experimental and theoretical research [17][18][19][20][21]. These heterocycles (C5H5E, E = element of group 15), like pyridine, can bind to metals via a lone pair at E (η 1 complexes) and/or via a π-electron sextet (η 6 complexes). Of course, the type of complex formed depends on the substituents attached to the ring. However, in the case of unsubstituted phosphabenzene and arsabenzene, the η 6 complexes predominate [22][23][24]. Phosphorine-based complexes and their use in homogeneous catalysis were summarized in a review by Le Floch [25]. Interestingly, it has been suggested that polycyclic aromatic hydrocarbons (PAHs) containing an endocyclic N atom (PANHs) may be responsible for the characteristic 6.22 µm feature in the emission spectra of some galaxies [26,27]. For this reason, astrophysicists predict that about 1 percent of the nitrogen in our Universe may exist in form of PANHs in an interstellar medium of galaxies. Similarly, the existence of P-containing PAHs in the interstellar medium (PAPHs) has been proposed and plausible mechanisms of their formation have been investigated [28]. However, such PAPHs have not yet been detected [29].
Borabenzene has not been separated in a pure form. However, several of its donoracceptor complexes, in which a Lewis basic site of a ligand donates electrons on the empty orbital of the boron atom, have been synthesized. Already in 1985, the borabenzene-pyridine complex was obtained and characterized using the XRD method [30]. This complex has also been studied by ab initio quantum chemistry methods, with emphasis on the B-N interaction, as well as IR and UV-vis spectra [31]. More recently, complexes of borabenzene with the N-heterocyclic carbenes Ime4 [32] and ImeS [33] were prepared, and their chemical reactivity was investigated [34]. Another example is the 1,3,5-tris (η 6 -borabenzene-η 5cyclopentadienylcobalt) benzene, which was deposited on a Cu(111) surface and studied using STM measurements and DFT calculations, revealing its possible application in molecular spintronics [35]. Additionally, the results of ab initio calculations showed that replacing carbon atoms in PAHs by bonded B and N atoms may modify PAHs' optical properties, leading to interesting nonlinear optical behavior of such compounds [36].
According to the Hückel rule, cyclic compounds are aromatic when they are planar and have 4n + 2 delocalized π electrons (where n is a positive integer). All studied heterocycles (Figure 1), following this rule, have six delocalized π electrons in the ring, just like benzene. However, the presence of a heteroatom in the ring has different consequences [37]. As shown in Table 2, energetic (Aromatic Stabilization Energy, ASE), electronic (electron density in the ring critical point, ρ RCP , and Electron Density of Delocalized Bonds, EDDB), [38,39] magnetic (Nucleus-Independent Chemical Shift, NICS) [40], and geometric (Bird, I6, and Harmonic Oscillator Model of Aromaticity, HOMA) [41,42] aromaticity indices do not speak exactly the one voice. Phosphorus and arsenic analogs of pyridine exhibit lower aromaticity than benzene and pyridine. However, for these two subgroups, there is significant inconsistency between the different aromaticity indices. For example, the ρ RCP , the HOMA, and one of the ASE [43] descriptors suggest higher aromaticity of the pyridine ring than that of benzene. On the other hand, the Bird index indicates that phosphorine is less aromatic than arsabenzene. Therefore, it is also interesting to investigate the effect of the substituent and its position on the aromaticity of the ring.

Materials and Methods
According to previous research [48][49][50], B3LYP functional [51] with the 6-311++G(d,p) [52] basis set was chosen as a computational method. Calculations were performed in the Gaussian 16, Revision C.01 program [53]. For all systems, frequency calculations were performed after geometry optimization. No imaginary frequencies were found. The geometries obtained in the calculations are in good agreement with the experimental data (X-ray diffraction, microwave spectroscopy) for unsubstituted pyridine [54] and arsabenzene [55]. Some differences in C-C bond lengths (~0.015 Å) and C-P-C angles (~1 • ) are observed in phosphorine; however, the uncertainties of experimental parameters are quite large [56]. A comparison of the calculated and experimental geometries is presented in Table S1 (Supplementary Materials).
Substituents were described in terms of the relative energy of the substituent, E rel (X). It is defined as an energetic effect of the homodesmotic reaction shown in Scheme 1.
are quite large [56]. A comparison of the calculated and experimental geometries is presented in Table S1 (Supplementary Materials).
Substituents were described in terms of the relative energy of the substituent, Erel(X). It is defined as an energetic effect of the homodesmotic reaction shown in Scheme 1. Hence, it can be calculated from Equation (1): where E(RX), E(benz), E(benzX), and E(R) are the electronic energies of the substituted R-X system, benzene (the reference system), X-substituted benzene, and unsubstituted R heterocycle, respectively. Hereinafter, Erel(X) is abbreviated as Erel. This parameter evaluates the energetic effect of interactions between the substituent and the substituted system. Negative values correspond to the stabilization of the molecule with respect to the reference system, whereas the positive values indicate destabilizing interactions. The charge of the substituent active region (cSAR) index [5,6] was used to characterize the electron-accepting/donating properties (EA/ED) of the substituents. For the X substituent, it was calculated according to Equation (2): where q(X) is a sum of charges of atoms at substituent X, and q(Cipso) is the atomic charge at the ipso carbon atom. In line with previous studies [4,57,58], atomic charges were evaluated using the Hirshfeld method [59]. Additionally, it should be noted that Hirshfeld charges have been documented to be almost insensitive to the basis set improvement [60]. The interpretation of the cSAR(X) values is as follows: the electron-donating strength of the X group increases with the increase in cSAR(X), whereas the electron-accepting strength decreases. Natural population analysis was performed using NBO 6.0 software (Theoretical Chemistry Institute, University of Wisconsin, Madison) [61]. Quantum theory of atoms in molecules (QTAIM) calculations and analyses were performed using AIMAII software [62].
For analysis of non-covalent interactions, a reduced density gradient (RDG) function was used [63]. This function is a modification of electron density (ρ) and its gradient (∇ρ) (Equation (3)).
Analysis of RDG isosurfaces is frequently used [64,65] to visualize the non-covalent interactions between groups of atoms. Additionally, the character of these interactions can be assessed by looking at the value of sign (λ2)•ρ (where λ2 is the second eigenvalue of the Hessian matrix of ρ) at each point of the isosurface. Negative values indicate strong noncovalent interactions such as hydrogen bonding, values around zero indicate weak ones (van der Waals), and positive values indicate non-bonding ones, e.g., the steric repulsion. Algorithms for the calculation of RDG implemented in Multiwfn software were used [66]. Visualization of isosurfaces was done in the VMD program [67]. Hence, it can be calculated from Equation (1): where E(RX), E(benz), E(benzX), and E(R) are the electronic energies of the substituted R-X system, benzene (the reference system), X-substituted benzene, and unsubstituted R heterocycle, respectively. Hereinafter, E rel (X) is abbreviated as E rel . This parameter evaluates the energetic effect of interactions between the substituent and the substituted system. Negative values correspond to the stabilization of the molecule with respect to the reference system, whereas the positive values indicate destabilizing interactions. The charge of the substituent active region (cSAR) index [5,6] was used to characterize the electron-accepting/donating properties (EA/ED) of the substituents. For the X substituent, it was calculated according to Equation (2): where q(X) is a sum of charges of atoms at substituent X, and q(C ipso ) is the atomic charge at the ipso carbon atom. In line with previous studies [4,57,58], atomic charges were evaluated using the Hirshfeld method [59]. Additionally, it should be noted that Hirshfeld charges have been documented to be almost insensitive to the basis set improvement [60]. The interpretation of the cSAR(X) values is as follows: the electron-donating strength of the X group increases with the increase in cSAR(X), whereas the electron-accepting strength decreases. Natural population analysis was performed using NBO 6.0 software (Theoretical Chemistry Institute, University of Wisconsin, Madison) [61]. Quantum theory of atoms in molecules (QTAIM) calculations and analyses were performed using AIMAII software [62].
For analysis of non-covalent interactions, a reduced density gradient (RDG) function was used [63]. This function is a modification of electron density (ρ) and its gradient (∇ρ) (Equation (3)).
Analysis of RDG isosurfaces is frequently used [64,65] to visualize the non-covalent interactions between groups of atoms. Additionally, the character of these interactions can be assessed by looking at the value of sign (λ 2 )·ρ (where λ 2 is the second eigenvalue of the Hessian matrix of ρ) at each point of the isosurface. Negative values indicate strong noncovalent interactions such as hydrogen bonding, values around zero indicate weak ones (van der Waals), and positive values indicate non-bonding ones, e.g., the steric repulsion. Algorithms for the calculation of RDG implemented in Multiwfn software were used [66]. Visualization of isosurfaces was done in the VMD program [67].
Aromaticity of substituted heteroarenes was evaluated using the Electron Density of Delocalized Bonds (EDDB) [39]. This descriptor derives from partitioning of the oneelectron density matrix in the natural atomic orbital basis into components representing localized bonds (EDLB(r)) and delocalized bonds (EDDB(r)) [38]. EDDB(r) proved to be a good descriptor of local and global aromaticity in several systems, including hydrocarbons [68] and metallacycles [69]. In this study, the component of EDDB(r) coming from π-electron delocalization, EDDB P (π), was used. Additionally, the values obtained for the unsubstituted systems were subtracted from those for the substituted heteroarenes. Such a difference indicates a change in aromaticity due to the substituent effect and is further denoted as ∆EDDB P (π).
The geometry of substituents and heterocyclic rings was evaluated using the parameters presented in Figure 2.
Aromaticity of substituted heteroarenes was evaluated using the Electron Density of Delocalized Bonds (EDDB) [39]. This descriptor derives from partitioning of the one-electron density matrix in the natural atomic orbital basis into components representing localized bonds (EDLB(r)) and delocalized bonds (EDDB(r)) [38]. EDDB(r) proved to be a good descriptor of local and global aromaticity in several systems, including hydrocarbons [68] and metallacycles [69]. In this study, the component of EDDB(r) coming from πelectron delocalization, EDDBP (π), was used. Additionally, the values obtained for the unsubstituted systems were subtracted from those for the substituted heteroarenes. Such a difference indicates a change in aromaticity due to the substituent effect and is further denoted as ΔEDDBP (π).
The geometry of substituents and heterocyclic rings was evaluated using the parameters presented in Figure 2.

Results
The discussion of the results is presented in four parts. The first is devoted to the energetic and electronic characteristics of substituents. The second part deals with proximity effects. Then, changes in the structural properties of the substituents and the substituted ring are discussed. The last part concerns the influence of substituents on π-electron delocalization. For all studied systems, the obtained values of the considered descriptors of the substituent effect are collected in Table S2 (Supplementary Materials).

Energetic and Electronic Substituent Effects
The energy of a substituent (Erel), Scheme 1, characterizes interactions of a substituent with a given substituted molecule. As mentioned above, its value indicates whether and to what extent these interactions are stabilizing or destabilizing. In order to allow comparisons with well-known systems for substituent effects, its values were calculated relative to benzene as the reference system. Erel can be considered as a global substituent effect descriptor, i.e., an indication of the energetic consequence of all intramolecular interactions induced by substitution. The local descriptor, cSAR, indicates changes induced in the electronic structure of a chosen molecular fragment, e.g., a substituent or a substituted ring.
It should be emphasized that in 2-, 3-, and 4-derivatives, different types of interactions occur ( Figure 3). For the 2-position (ortho interaction), through-bond and throughspace proximity interactions between O (in NO2) or H (in NH2) and the endocyclic heteroatom are expected. In 3-derivatives, the proximity interaction with the heteroatom is replaced by an interaction with the C-H group, thus representing a situation similar to the substituted benzene reference system. Therefore, it can be assumed that in these cases the energetic substituent effect is only a consequence of inductive and resonance interactions with the heteroatom in the meta position. Regarding 4-derivatives, a similar assumption can be made, but in this case, the para interaction occurs.

Results
The discussion of the results is presented in four parts. The first is devoted to the energetic and electronic characteristics of substituents. The second part deals with proximity effects. Then, changes in the structural properties of the substituents and the substituted ring are discussed. The last part concerns the influence of substituents on π-electron delocalization. For all studied systems, the obtained values of the considered descriptors of the substituent effect are collected in Table S2 (Supplementary Materials).

Energetic and Electronic Substituent Effects
The energy of a substituent (E rel ), Scheme 1, characterizes interactions of a substituent with a given substituted molecule. As mentioned above, its value indicates whether and to what extent these interactions are stabilizing or destabilizing. In order to allow comparisons with well-known systems for substituent effects, its values were calculated relative to benzene as the reference system. E rel can be considered as a global substituent effect descriptor, i.e., an indication of the energetic consequence of all intramolecular interactions induced by substitution. The local descriptor, cSAR, indicates changes induced in the electronic structure of a chosen molecular fragment, e.g., a substituent or a substituted ring.
It should be emphasized that in 2-, 3-, and 4-derivatives, different types of interactions occur ( Figure 3). For the 2-position (ortho interaction), through-bond and through-space proximity interactions between O (in NO 2 ) or H (in NH 2 ) and the endocyclic heteroatom are expected. In 3-derivatives, the proximity interaction with the heteroatom is replaced by an interaction with the C-H group, thus representing a situation similar to the substituted benzene reference system. Therefore, it can be assumed that in these cases the energetic substituent effect is only a consequence of inductive and resonance interactions with the heteroatom in the meta position. Regarding 4-derivatives, a similar assumption can be made, but in this case, the para interaction occurs. As shown in Figure 4, the NO2 and NH2 groups show a significantly different, in most cases opposite, energetic substituent effect. The nitro derivatives in all cases but one have positive Erel values, which indicates lower stability with respect to nitrobenzene.  As shown in Figure 4, the NO 2 and NH 2 groups show a significantly different, in most cases opposite, energetic substituent effect. The nitro derivatives in all cases but one have positive E rel values, which indicates lower stability with respect to nitrobenzene. As shown in Figure 4, the NO2 and NH2 groups show a significantly different, in most cases opposite, energetic substituent effect. The nitro derivatives in all cases but one have positive Erel values, which indicates lower stability with respect to nitrobenzene. In the case of nitro derivatives of pyridine, Erel grows in the sequence 2-> 3-> 4nitropyridine (the abbreviations shown in Figure 4 are further used). This indicates that the resonance-driven para interaction between NO2 and Nendo results in a substantial destabilization. Since the meta resonance effect is negligible [3], the Erel of 2.63 kcal/mol in N-3-NO2 can be associated with a decrease in aromaticity due to substitution, which is discussed in Section 3.4. Although the ortho interaction is associated with the intuitively repulsive NO•••Nendo contact, Erel has the lowest value. Moreover, it has recently been shown by analyzing the RDG function that this interaction is not purely repulsive and has an attractive contribution from van der Waals interaction [15]. Thus, the highest stability of N-2-NO2 can result from the strongest inductive effect due to the close distance between NO2 and Nendo and the proximity interaction. However, it should be remembered that all the Erel values for NO2-substituted N-heterocycles are positive. Thus, exchanging one of the CH groups for N in nitrobenzene leads to destabilization.
Amino derivatives of pyridine are additionally stabilized compared with aniline ( Figure 4). The stability sequence is N-2-NH2 > N-4-NH2 > N-3-NH2. The meta interaction leads to a negligible energetic effect, while the para interaction has roughly the same absolute value as for N-4-NO2. However, with an opposite sign-in this case, stabilization occurs. The ortho interaction in N-2-NH2 leads to the most stabilizing interactions between the substituent and the pyridine ring. A possible cause is the enhancement (with respect In the case of nitro derivatives of pyridine, E rel grows in the sequence 2-> 3-> 4nitropyridine (the abbreviations shown in Figure 4 are further used). This indicates that the resonance-driven para interaction between NO 2 and N endo results in a substantial destabilization. Since the meta resonance effect is negligible [3], the E rel of 2.63 kcal/mol in N-3-NO 2 can be associated with a decrease in aromaticity due to substitution, which is discussed in Section 3.4. Although the ortho interaction is associated with the intuitively repulsive NO···N endo contact, E rel has the lowest value. Moreover, it has recently been shown by analyzing the RDG function that this interaction is not purely repulsive and has an attractive contribution from van der Waals interaction [15]. Thus, the highest stability of N-2-NO 2 can result from the strongest inductive effect due to the close distance between NO 2 and N endo and the proximity interaction. However, it should be remembered that all the E rel values for NO 2 -substituted N-heterocycles are positive. Thus, exchanging one of the CH groups for N in nitrobenzene leads to destabilization.
Amino derivatives of pyridine are additionally stabilized compared with aniline ( Figure 4). The stability sequence is N-2-NH 2 > N-4-NH 2 > N-3-NH 2 . The meta interaction leads to a negligible energetic effect, while the para interaction has roughly the same absolute value as for N-4-NO 2 . However, with an opposite sign-in this case, stabilization occurs. The ortho interaction in N-2-NH 2 leads to the most stabilizing interactions between the substituent and the pyridine ring. A possible cause is the enhancement (with respect to aniline) of inductive and resonance effects by introducing the N atom in the ortho position. Another possible explanation could be the presence of the N···HN proximity interaction. This interaction, however, has been recently proven to be extremely weak and has not been detected by means of RDG and AIM analyses in any six-membered N-heteroarene [15].
Borabenzene derivatives are interesting cases, as the boron atom is the least electronegative in the studied series (Table 1). Additionally, an empty 2p orbital of B in the ring plane allows the formation of a donor-acceptor interaction with electron-rich atoms (Lewis basic sites) of the ortho substituent. Nitro substitution in the 2-position is associated with a small positive energetic effect (1.12 kcal/mol), in the 3-position with a negligible effect (0.13 kcal/mol), and in the 4-position with a moderately negative effect (−2.07 kcal/mol). However, the amino substitution in the 2-and 4-positions leads to a large destabilizing effect (5.00 and 5.38 kcal/mol, respectively). This is due to changes in the electron-donating power of the amino group, shown in Figures 5 and 6, which are discussed below. On the other hand, interactions between the substituent and the substituted system in 3-B-NH 2 are slightly more favorable than in aniline (−0.40 kcal/mol). negative in the studied series (Table 1). Additionally, an empty 2p orbital of B in the ring plane allows the formation of a donor-acceptor interaction with electron-rich atoms (Lewis basic sites) of the ortho substituent. Nitro substitution in the 2-position is associated with a small positive energetic effect (1.12 kcal/mol), in the 3-position with a negligible effect (0.13 kcal/mol), and in the 4-position with a moderately negative effect (−2.07 kcal/mol). However, the amino substitution in the 2-and 4-positions leads to a large destabilizing effect (5.00 and 5.38 kcal/mol, respectively). This is due to changes in the electron-donating power of the amino group, shown in Figures 5 and 6, which are discussed below. On the other hand, interactions between the substituent and the substituted system in 3-B-NH2 are slightly more favorable than in aniline (−0.40 kcal/mol).   Borabenzene derivatives are interesting cases, as the boron atom is the least electronegative in the studied series (Table 1). Additionally, an empty 2p orbital of B in the ring plane allows the formation of a donor-acceptor interaction with electron-rich atoms (Lewis basic sites) of the ortho substituent. Nitro substitution in the 2-position is associated with a small positive energetic effect (1.12 kcal/mol), in the 3-position with a negligible effect (0.13 kcal/mol), and in the 4-position with a moderately negative effect (−2.07 kcal/mol). However, the amino substitution in the 2-and 4-positions leads to a large destabilizing effect (5.00 and 5.38 kcal/mol, respectively). This is due to changes in the electron-donating power of the amino group, shown in Figures 5 and 6, which are discussed below. On the other hand, interactions between the substituent and the substituted system in 3-B-NH2 are slightly more favorable than in aniline (−0.40 kcal/mol).   Phosphorine and arsabenzene contain group 15 heteroatoms of similar electronegativity but differing in other atomic properties (see Table 1). The trends in E rel values are similar for these heterocycles. For both amino and nitro substitution, E rel decreases in the order of 2-, 3-, 4-, meaning that the ortho derivatives are the least stable and the para ones are the most. The difference between NO 2 and NH 2 systems lies in the signs of E rel . In all NO 2 derivatives, the interaction between the substituent and the heteroarene is less, while in all NH 2 derivatives it is more energetically favorable than in their corresponding benzene derivatives.
Electron-accepting and electron-donating properties of the NO 2 and NH 2 groups, respectively, can be assessed using the cSAR(X) approach ( Figure 5). It can be noticed immediately that in meta derivatives the variability in cSAR(X) is small, and the values are close to the values for nitrobenzene and aniline. Although the differences are small, in all cases the values for the NO 2 group in heterocycles are slightly higher than in benzene.
A similar trend occurs in the case of NH 2 derivatives. Thus, the cSAR(X) indicates that in the meta derivatives of all heterocycles, the NH 2 group is slightly more electron-donating than in aniline, and NO 2 is slightly less electron-withdrawing than in nitrobenzene.
In ortho and para systems, where the resonance effect of substituent on heteroatom occurs, the changes in cSAR(X) and E rel in N and B heteroarenes are highly related to the electronegativity of these atoms. This can be illustrated by the resonance structures of the NH 2 and NO 2 -substituted heteroarenes ( Figure 6). The E heteroatoms that are more electronegative than C, i.e., N, will more likely accept electrons donated by the substituent, increasing the negative charge of E ( Figure 6, structures (5)-(8)). Hence, when such a heteroatom is in the ortho or para position, the electron-donating properties of the NH 2 group may be enhanced and the resonance effect strengthened, which, as can be seen from the negative E rel values (Figure 4), also stabilizes the molecule. On the other hand, the electronegative N endo atom is not prone to increase its positive charge ( Figure 6, structures (1)-(4)). Hence, the withdrawal of electrons by resonance in nitro derivatives is weakened, and positive values of E rel are observed ( Figure 4).
Conversely, when the E heteroatom is less electronegative than C, i.e., B, the resonance effect of the electron-accepting group is enhanced while that of the electron-donating group is weakened. This is well illustrated by the energetic and electronic descriptors of the substituent effect in the B-2-NH 2 and B-4-NH 2 systems. The E rel values reveal the strongest destabilizing intramolecular interactions (among the studied systems), which result in the greatest weakening of the electron-donating power of the amino group. The latter is documented by the lowest cSAR(NH 2 ) values ( Figure 5).
In As and P heteroarenes, for X = NO 2 there is a strong enhancement of the EA properties for the 2-substitution. Regarding the 4-derivatives, the NO 2 is more EA than in the 3-derivatives, but only slightly. In the case of X = NH 2 , the amino group in the 2-derivatives is the least electron-donating, and in the 4-derivatives has slightly more ED ability than in the 3-derivatives.
Therefore, the electronegativity of the heteroatom in the arene ring highly influences the interactions within the substituted system. Depending on the type of substituentelectron-withdrawing (NO 2 ) or electron-donating (NH 2 )-interactions can be enhanced or weakened with respect to benzene derivatives. This is reflected in both the energetic (E rel , Figure 4) and electronic (cSAR, Figure 5) properties of the substituents. In general, in the substituted ortho and para nitro and amino series, the order of changes in the energetic and electronic properties of the substituent follows the changes in electronegativity of the heteroatoms (Table 1), as shown in Figure 7. Moreover, cSAR(X) and E rel are well correlated for the series: para NO 2 and NH 2 and ortho NH 2 (R 2 > 0.98). In para systems, the proximity effects are similar to those in benzene derivatives (i.e., with C-H groups), and the resonance between heteroatom and substituent is strong. Consequently, the energetic effect of the intramolecular interactions is almost entirely due to the enhancement or weakening of the resonance effect by the heteroatom, which causes changes in the electronaccepting/donating ability of the substituent. The strongest resonance effects, as expected, occur for B-4-NO 2 and N-4-NH 2 derivatives. The slopes of the linear relation between cSAR(X) and E rel are 0.0191 and −0.0131 for nitro and amino systems, respectively. Their signs are consistent with those of Hammett's substituent constants (σ p = 0.78 and −0.66) [3]. In ortho systems, the resonance effect is also strong, but different proximity effects also occur. The fact that the ortho series of NH 2 is well correlated while that of NO 2 is not can be explained by proximity effects. They are more significant for the bulky NO 2 group than for the small NH 2 group, which is also noticeable from the distance of the benzNH2 and benzNO2 points, in Figure 7, from the NH 2 and NO 2 ortho line series, respectively. Namely, the benzNO2 point is much further from the corresponding series than the benzNH2 point. Moreover, for the NH 2 group, the slopes for the ortho and para series are almost the same. However, in the first case, the ranges of E rel and cSAR(NH 2 ) variability are~1.4 times greater than in the second. This is mainly due to the much greater stabilization of NH 2 interactions with the ring nitrogen atom in the ortho position than in the para position.
Small variability in meta derivatives is consistent with the well-known fact that the meta substituent effects are weaker [3]. by proximity effects. They are more significant for the bulky NO2 group than for the small NH2 group, which is also noticeable from the distance of the benzNH2 and benzNO2 points, in Figure 7, from the NH2 and NO2 ortho line series, respectively. Namely, the benzNO2 point is much further from the corresponding series than the benzNH2 point. Moreover, for the NH2 group, the slopes for the ortho and para series are almost the same. However, in the first case, the ranges of Erel and cSAR(NH2) variability are ~1.4 times greater than in the second. This is mainly due to the much greater stabilization of NH2 interactions with the ring nitrogen atom in the ortho position than in the para position. Small variability in meta derivatives is consistent with the well-known fact that the meta substituent effects are weaker [3]. Figure 7. Relationships between Erel and cSAR(X) for ortho-, meta-, and para-series of amino and nitro derivatives of monosubstituted heteroarenes. Points representing nitro-(benzNO2) and aminobenzene (benzNH2) are shown (including them in the corresponding para series: y = 0.0188x-0.1471, R 2 = 0.962, and y = -0.131x + 0.1041, R 2 = 0.974, respectively).

Proximity Effects
The heteroatom in position 2 provides an opportunity for various through-space interactions. In such a case, the asymmetric proximity of the substituent causes asymmetricity of the two NO (of NO2) or NH (of NH2) bonds. Hence, the difference in length (Figure 8) or the electron density at the bond critical point ( Figure S1) between these two bonds can indicate the strength of the proximity interaction. The bond depicted with (E) is the bond closer to the heteroatom, whereas (H) is the bond further away from the heteroatom.

Proximity Effects
The heteroatom in position 2 provides an opportunity for various through-space interactions. In such a case, the asymmetric proximity of the substituent causes asymmetricity of the two NO (of NO 2 ) or NH (of NH 2 ) bonds. Hence, the difference in length (Figure 8) or the electron density at the bond critical point ( Figure S1) between these two bonds can indicate the strength of the proximity interaction. The bond depicted with (E) is the bond closer to the heteroatom, whereas (H) is the bond further away from the heteroatom. As shown in Figure 8a, in the case of the NO2 substitution, systems deviating from the y = x line are the 2-and 3-NO2 heteroarenes. The systems above this line (marked with green labels) have a longer NO(E) bond than NO(H) bond, which indicates an attraction between the O and E heteroatom. The systems below the line (marked with red labels) represent an opposite case, thus indicating a repulsion between these atoms. The highest As shown in Figure 8a, in the case of the NO 2 substitution, systems deviating from the y = x line are the 2-and 3-NO 2 heteroarenes. The systems above this line (marked with green labels) have a longer NO(E) bond than NO(H) bond, which indicates an attraction between the O and E heteroatom. The systems below the line (marked with red labels) represent an opposite case, thus indicating a repulsion between these atoms. The highest difference occurs in the B-2-NO 2 system. The NO bond pointing towards the B atom is elongated due to a donor-acceptor interaction between the lone pair of O atom and the empty 2p orbital of B. This interaction is detected by the second-order perturbation theory analysis in the NBO theory [70] visible in the results of the NBO calculation as the energy of the intramolecular donor-acceptor interaction O(LP) -> B(LP *), E = 2.76 kcal/mol, and no steric energy. However, the RDG analysis (at the 0.50 isosurface) does not reveal the B···O contact as an interaction (Figure 9). Interestingly, a similar NO bond elongation is observed in As-and P-2-NO 2 systems, but to a much lesser extent, and the RDG analysis shows E···O contacts (E = As or P). Moreover, in these cases, according to the NBO results, no donoracceptor interactions were detected, while the steric energies were 1.61 and 2.55 kcal/mol, respectively. The noticeable difference in the NO(E) bond elongations between these two systems may be related to the difference in the polarizability of As and P atoms (Table 1). In N-2-NO 2 , the bond closer to the N endo atom is much shorter, indicating a steric strain caused by the repulsion between the lone pairs of O and N endo atoms. The estimated (NBO) steric energy is 2.09 kcal/mol, while no donor-acceptor interaction was detected, contrary to the RDG results ( Figure 9). It can therefore be concluded that the results of the RDG and NBO analyses do not speak the same voice. The 3-NO 2 systems lie slightly below the y = x line; therefore, despite the presence of two CH interactions, there is still some asymmetry. It follows that the interaction with the CH group of the ring closer to the As, P, or N heteroatom causes a slight shortening of the NO(E) bond. As regards the NH2-substituted heteroarenes, the differences between the NH bond lengths are an order of magnitude smaller (up to 0.002 Å, full data available in Table S2) compared with the NO of NO2 (up to 0.014 Å). This indicates that the proximity interactions of the NH2 group are relatively weak. The only system in which the NH(E) bond is As regards the NH 2 -substituted heteroarenes, the differences between the NH bond lengths are an order of magnitude smaller (up to 0.002 Å, full data available in Table S2) compared with the NO of NO 2 (up to 0.014 Å). This indicates that the proximity interactions of the NH 2 group are relatively weak. The only system in which the NH(E) bond is elongated is N-2-NH 2 (by 0.0017 Å, Figure 8), suggesting an attractive NH···N interaction. Additionally, according to the literature, it can be assumed that the electronic interactions in the ortho and para positions are alike. Thus, this assumption leads to the conclusion that the energy of NH···N interaction is~3.7 kcal/mol. However, this is not supported by the geometry (angle NH · N equals 71.2 • ) and the results of the RDG analysis (Figure 9). Such a low value of the angle is well below the minimum threshold for the AH···B angle in the H-bond of~110 • suggested by Desiraju [71]. On the other hand, the shortening of NH(E) in relation to NH(H) is observed in B-2-NH 2 and As-2-NH 2 (by 0.0020 and 0.0017 Å, respectively) and negligible in P-2-NH 2 and P-, B-, and As-3-NH 2 (by up to 0.0006 Å).
None of the proximity interactions have a bond critical point. However, the results of the RDG analysis ( Figure 9) indicate that the proximity interactions of the bulky NO 2 group are stronger than in the case of NH 2 . This explains the lack of correlation observed in the ortho NO 2 series in Figure 7. Furthermore, no single proximity interaction of the NH 2 group could be observed in the RDG plots ( Figure 9) and in the NBO analysis results.

Geometric Properties of Substituents and the Substituted Moiety
The changes in structural properties (of both substituents and rings) are shown in Figures 10-12. The values obtained are summarized in Table S2.

Geometric Properties of Substituents and the Substituted Moiety
The changes in structural properties (of both substituents and rings) are shown in Figures 10-12. The values obtained are summarized in Table S2.   The greatest variability in d CN of nitro derivatives is observed in the ortho series (0.073 Å) due to the proximity effects. Smaller variability is observed in the para series (0.02 Å), and the smallest in the meta one (0.012 Å). Moreover, for the N, P, and As series, except for N-3-NO 2 , the CN bonds are longer than in nitrobenzene. However, this bond is significantly shorter in the B-2-and B-4-NO 2 derivatives (by 0.045 and 0.013 Å, respectively). This confirms the strongest resonance effect of the NO 2 group, even resulting in its stabilizing intramolecular interaction in the B-4-NO 2 molecules, i.e., negative E rel values (Figure 4). The CN bond in P-2-NO 2 is slightly longer than in As-2-NO 2 (Figure 10), despite the larger size of the As atom (Table 1). This again indicates that the NO 2 proximity interactions with As may be more favorable than with P.    In the case of the amino derivatives, the greatest variability in d CN occurs in the para series (0.038 Å) and, opposite to the case of the nitro substitution, it is greater than in the ortho series (0.019 Å). This fact suggests that stronger proximity interactions of the NO 2 group than of the NH 2 , as revealed by RDG analysis (Figure 9), result in more significant changes in the geometry of the NO 2 group. The greatest elongation of the CN bond in the B-4-NH 2 system confirms the strongest disturbance of the observed resonance effect. In turn, the shortening of this bond in N-2-and N-4-NH 2 derivatives, in relation to the CN bond in aniline, shows the enhancement of the resonance effect. The meta series has a very small variability of d CN (0.004 Å). In this case, the d CN values are very similar to aniline. This shows that the heteroatom in the meta position scarcely changes the strength of the reverse substituent effect. It can also be noticed in Figure 5, where the electron-donating properties of the amino group show little variability in the meta series.
In the case of ortho-nitro-substituted heteroarenes, changes in the ∠CEC angle ( Figure 11) are much more significant than for the amino substitution. This is another consequence of the strong NO 2 proximity effects. Additionally, ortho-nitro and amino substitution causes opposite changes in this angle. Amino substitution of the studied heterocycles in the para position leads to the highest changes in the ∠CEC. It is also worth noting that they all are negative, meaning that the ∠CEC of all heterocycles decreases upon substitution. Since in the para position no steric effects occur, this can be solely associated with the electronic interactions with the substituted system. Interestingly, the greatest change occurs in B-4-NH 2 , whereas the lowest change occurs in N-4-NH 2 . As mentioned earlier, the B and N heteroatoms have the opposite effect on resonance between the substituent and heterocycle. This can also be noticed by looking at the ∠ONO angle ( Figure 12a). Here, the values are in most cases similar to that of nitrobenzene. However, for B-4-NO 2 , N-2-NO 2 , and N-4-NO 2 they deviate from the others. The smallest ∠ONO occurs in B-4-NO 2 and the highest in N-2 and N-4 derivatives. It follows that the enhancement of π-resonance is associated with a decrease in ∠ONO relative to the nitrobenzene, while an increase indicates a disruption. The greater ∠ONO in N-2 than in N-4 may also result from proximity interactions.
Since the amino group is non-planar, it is worth analyzing its pyramidality instead of the ∠HNH angle. In Figure 12b it can be noticed that in 2-and 4-derivatives of B-heterocycles, the NH 2 group has the highest pyramidality. This can be associated with an increase in the p character of a hybrid orbital at the N atom, again due to the resonance disruption. For the opposite reason, among the para systems, the NH 2 group is the most planar in the pyridine derivative. Moreover, only for this series a very good correlation of pyramidality with E rel and cSAR was found (R 2 > 0.99, Figure S2).
The substituent also affects the angle α in the ring (at the ipso carbon atom, Figure 2). The nitro group increases this angle with respect to the unsubstituted system, while the amino group has the opposite effect. Thus, changes in α (expressed by ∆α) are positive and negative numbers, respectively (see Figure S3). They mainly result from the resonance effect and proximity interactions mentioned above.
Additionally, the correlation coefficients between all pairs of considered parameters for the para series (benzene derivatives were also taken into account) are gathered in Table S3. For example, the pyramidality of the NH 2 group in this series is well correlated with ∆α, ∆∠CEC, d CN , cSAR(X), and E rel .

Substituent Influence on π-Electron Delocalization
For heteroarenes, the aromaticity of a substituted moiety is an important indicator of the changes induced within it upon substitution. In all cases, the substituted heteroarenes are less aromatic than the unsubstituted systems ( Figure S4), as shown by ∆EDDB P (π) in Figure 13. The smallest decreases in aromaticity occur in systems where the interactions between the substituent and the heteroatom are weak (for the reasons discussed earlier; Figure 6). Such systems are all meta and ortho/para, in which the heteroatom hinders the resonance effect of the substituent (N-2-NO 2 , N-4-NO 2 , B-4-NH 2 , B-2-NH 2 ). The strongest decreases in aromaticity occur in the systems where the heteroatom enhances the resonance effect: N-2-NH 2 , N-4-NH 2 , B-2-NO 2 , and B-4-NO 2 . Other highly affected systems are the rings in the 2-and 4-P and As systems, but this is more pronounced for NH 2 substitution than for NO 2 substitution. Moreover, in general, the amino group causes greater aromaticity changes than the nitro one. The only exceptions are the above-mentioned B-2-and B-4 derivatives.
Plots between ∆EDDB P (π) and E rel ( Figure 14) show that for the para-NO 2 series, and for the NH 2 orthoand para-NH 2 series, changes in aromaticity due to substitution correlate with E rel . This shows that the changes in the stability of systems in which the heteroatom and substituent can interact by resonance are, to a high extent, related to changes in the aromaticity of the substituted heterocyclic ring. The poor determination coefficient (R 2 = 0.53) for the ortho-NO 2 series may result from the aforementioned strong proximity interaction of this group, which is another factor affecting the stability of the molecule. For the para series, a higher slope can be observed in the case of NO 2 substitution than in the case of NH 2 substitution (0.088 and 0.051, respectively); similarly for the dependence ∆EDDB P (π) on cSAR(X) (Figure 14c,d). This is consistent with the higher absolute slope value for the NO 2 series of the cSAR(X) vs. E rel relationship (Figure 7). So, we can say that the resonance effect of the nitro group is stronger than that of the amino substituent, just like in aromatic hydrocarbons. Moreover, the dependence of ∆EDDB P (π) on cSAR(X) shows a good correlation between changes in aromaticity and electronic properties of the substituent for the orthoand para-NH 2 series. However, in the case of NO 2 derivatives, in the para series a good correlation is observed (R 2 = 0.979), while a worse correlation is observed for the ortho systems (R 2 = 0.733). This fact can be explained by the proximity effects of the nitro group, as shown above.

Substituent Influence on π-Electron Delocalization
For heteroarenes, the aromaticity of a substituted moiety is an important indicator of the changes induced within it upon substitution. In all cases, the substituted heteroarenes are less aromatic than the unsubstituted systems ( Figure S4), as shown by ΔEDDBP (π) in Figure 13. The smallest decreases in aromaticity occur in systems where the interactions between the substituent and the heteroatom are weak (for the reasons discussed earlier; Figure 6). Such systems are all meta and ortho/para, in which the heteroatom hinders the resonance effect of the substituent (N-2-NO2, N-4-NO2, B-4-NH2, B-2-NH2). The strongest decreases in aromaticity occur in the systems where the heteroatom enhances the resonance effect: N-2-NH2, N-4-NH2, B-2-NO2, and B-4-NO2. Other highly affected systems are the rings in the 2-and 4-P and As systems, but this is more pronounced for NH2 substitution than for NO2 substitution. Moreover, in general, the amino group causes greater aromaticity changes than the nitro one. The only exceptions are the above-mentioned B-2-and B-4 derivatives. Plots between ∆EDDBP (π) and Erel ( Figure 14) show that for the para-NO2 series, and for the NH2 ortho-and para-NH2 series, changes in aromaticity due to substitution correlate with Erel. This shows that the changes in the stability of systems in which the heteroatom and substituent can interact by resonance are, to a high extent, related to changes in the aromaticity of the substituted heterocyclic ring. The poor determination coefficient (R 2 = 0.53) for the ortho-NO2 series may result from the aforementioned strong proximity interaction of this group, which is another factor affecting the stability of the molecule. For the para series, a higher slope can be observed in the case of NO2 substitution than in the case of NH2 substitution (0.088 and 0.051, respectively); similarly for the dependence ∆EDDBP (π) on cSAR(X) (Figure 14c,d). This is consistent with the higher absolute slope value for the NO2 series of the cSAR(X) vs. Erel relationship (Figure 7). So, we can say that the resonance effect of the nitro group is stronger than that of the amino substituent, just like in aromatic hydrocarbons. Moreover, the dependence of ∆EDDBP (π) on cSAR(X) shows a good correlation between changes in aromaticity and electronic properties of the substituent for the ortho-and para-NH2 series. However, in the case of NO2 derivatives, in the para series a good correlation is observed (R 2 = 0.979), while a worse correlation is observed

Conclusions
This study aimed to investigate to what extent the heteroatom in the six-membered aromatic ring affects the energetic and geometric properties of the strongly electron-attracting or electron-donating substituents, i.e., the nitro and amino groups, respectively.

Conclusions
This study aimed to investigate to what extent the heteroatom in the six-membered aromatic ring affects the energetic and geometric properties of the strongly electron-attracting or electron-donating substituents, i.e., the nitro and amino groups, respectively. For this purpose, pyridine, phosphorine, arsabenzene, and borabenzene, substituted at the 2-, 3-, or 4-position, were chosen as the substituted moieties. The energetic effect of interactions between the substituent and the substituted system was estimated using the relative energy parameter E rel . This parameter is based on the energetic effect of an appropriate homodesmotic reaction with benzene as a reference system and enables the comparison of the properties of the substituents. To characterize the electron-accepting/donating properties of substituents, as well as their changes, the charge of the substituent active region (cSAR) model was used. Aromaticity of heteroarenes was evaluated using the Electron Density of Delocalized Bonds (EDDB) descriptor.
The NO 2 and NH 2 groups show a significantly different, in most cases opposite, energetic substituent effect. The nitro derivatives in all cases but one have positive E rel values, which indicate lower stability with respect to nitrobenzene. The opposite is the case with the amine derivatives; negative E rel values show greater stability with respect to aniline. In both cases, the exceptions are borabenzene derivatives, namely B-4-NO 2 , B-2-NH 2 , and B-4-NH 2 . This is due to the strengthening or weakening of the resonance effect by the boron atom (the least electronegative in the studied series), respectively.
In general, in the substituted ortho and para nitro and amine series, the order of changes in the energetic and electronic properties of the substituent is related to the changes in electronegativity of the heteroatoms. These properties of NO 2 and NH 2 in a given heterocycle are almost entirely due to the enhancement or weakening of the resonance effect by the heteroatom. The greatest strengthening occurs for the ortho and para B-NO 2 and N-NH 2 systems, while the greatest weakening occurs for the B-NH 2 and N-NO 2 derivatives.
In ortho systems, in addition to the resonance effect, there are also different proximity effects. They can cause significant changes in the energetic, electronic, and structural properties of the substituent. The easiest way to observe asymmetric proximity interactions is to compare the lengths of both NO (or NH) bonds in the nitro (amino) group. The largest difference occurs in the B-2-NO 2 system. The NO bond towards the B atom is elongated due to a donor-acceptor interaction between the lone O pair and the empty 2p orbital B. In contrast, in N-2-NO 2 , the bond closer to the N endo atom is much shorter, indicating a steric strain caused by the repulsion between the lone pairs of O and N endo atoms. The differences between the NH bond lengths of the NH 2 group are significantly smaller (up to 0.002 Å) compared with the NO bond lengths of the NO 2 group (up to 0.014 Å). Thus, the proximity effects are much greater for the bulky NO 2 group than for the small NH 2 group, which is also confirmed by the RDG analysis results. It is also reflected in the fact that the correlations between cSAR, E rel , and ∆EDDBP (π) for ortho NO 2 systems are worse than for NH 2 derivatives.
The above-mentioned effects are accompanied by changes in the structural properties of both the substituents and the rings. The substitution of the group at different positions of the heteroarene affects the CN bond lengths, the CEC and α (at the ipso carbon) angles in the ring, and the pyramidality of the NH 2 group. The greatest changes in structural parameters in relation to benzene derivatives are generally observed for 2-and/or 4-derivatives.
Substituted heteroarenes are less aromatic than unsubstituted systems, while the latter are less aromatic than benzene. In general, the amino group causes greater aromaticity changes than the nitro one. The only exceptions are the B-2-and B-4 derivatives (B hinders the resonance effect of the NH 2 substituent, hence the smaller changes). The greatest decreases in aromaticity occur in the systems in which the heteroatom enhances the resonance effect: N-2-NH 2 , N-4-NH 2 , B-2-NO 2 , and B-4-NO 2 . Moreover, for para series, higher slopes (or their absolute values) of the linear dependences ∆EDDB P (π) on E rel , ∆EDDB P (π) on cSAR(X), and cSAR(X) on E rel are observed in the case of NO 2 substitution than for the NH 2 one. It can therefore be said that the resonance effect of the nitro group is stronger than that of the amine group, as is the case with aromatic hydrocarbons.
In summary, the ring heteroatom and the substitution position can significantly change the electron-attracting and electron-donating abilities of the substituent. Moreover, the substituent can significantly affect the electronic and structural properties of the ring.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/ 10.3390/sym14010145/s1, Figure S1: Dependence between the electron densities at critical points of two NO (for X = NO 2 ) or NH (for X = NH 2 ) bonds; Figure S2: Dependences of the pyramidality of the amino group on its cSAR and energetic descriptor, E rel , for the studied systems; Figure S3: Obtained ∆α values for heteroarenes substituted in the 2-, 3-, or 4-positions by the nitro or amino groups; Figure S4: Obtained EDDB P (π) values for heteroarenes unsubstituted and substituted in the 2-, 3-, or 4-positions by the nitro or amino groups; Figure S5: Relationships between electron densities at the critical point of CN bonds and their length for the studied nitro and amine derivatives; Table S1: Comparison of experimental and calculated geometries of pyridine, arsabenzene, and phosphorine; Table S2: Values of all calculated parameters for all studied systems; Table S3: Correlation coefficients for all para derivatives, including benzene derivatives (nitrobenzene and aniline, respectively).