Next Article in Journal
Stabilizing Highly Active Ru Sites by Electron Reservoir in Acidic Oxygen Evolution
Previous Article in Journal
Molecular-Simulation–Inspired Synthesis of [6]-Prismane via Photoisomerisation of Octafluoro[2.2]paracyclophane
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 29
Issue 4
10.3390/molecules29040784
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Substitution Effect of a Single Nitrogen Atom on π-Electronic Systems of Linear Polycyclic Aromatic Hydrocarbons (PAHs): Theoretically Visualized Coexistence of Mono- and Polycyclic π-Electron Delocalization

by
Jong Min Lim
1,†,
Sangdeok Shim
2,†,
Hoa Thi Bui
3,
Jimin Kim
3,
Ho-Joong Kim
4,
Yoon Hwa
5,* and
Sung Cho
3,*
1
Department of Chemistry, Kyungpook National University, Daegu 41566, Republic of Korea
2
Department of Chemistry, Sunchon National University, Suncheon 57922, Republic of Korea
3
Department of Chemistry, Chonnam National University, Gwangju 61186, Republic of Korea
4
Department of Chemistry, Chosun University, Gwangju 61452, Republic of Korea
5
School of Electrical, Computer and Energy Engineering, Arizona State University, Tempe, AZ 85287, USA
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Molecules 2024, 29(4), 784; https://doi.org/10.3390/molecules29040784
Submission received: 12 January 2024 / Revised: 31 January 2024 / Accepted: 5 February 2024 / Published: 8 February 2024
(This article belongs to the Topic Conjugated Polymers: Preparation, Properties and Applications)
Browse Figures
Versions Notes

Abstract

:
We theoretically investigated the nitrogen substitution effect on the molecular structure and π-electron delocalization in linear nitrogen-substituted polycyclic aromatic hydrocarbons (N-PAHs). Based on the optimized molecular structures and magnetic field-induced parameters of fused bi- and tricyclic linear N-PAHs, we found that the local π-electron delocalization of subcycles (e.g., mono- and bicyclic constituent moieties) in linear N-PAHs is preserved, despite deviation from ideal structures of parent monocycles. The introduction of a fused five-membered ring with a pyrrolic N atom (N-5MR) in linear N-PAHs significantly perturbs the π-electronic condition of the neighboring fused six-membered ring (6MR). Monocyclic pyrrole exhibits substantial bond length alternations, strongly influencing the π-electronic systems of both the fused N-5MR and 6MR in linear N-PAHs, depending on the location of shared covalent bonds. A fused six-membered ring with a graphitic N atom in an indolizine moiety cannot generate monocyclic π-electron delocalization but instead contributes to the formation of polycyclic π-electron delocalization. This is evidenced by bifurcated diatropic ring currents induced by an external magnetic field. In conclusion, the satisfaction of Hückel’s 4n + 2 rule for both mono- and polycycles is crucial for understanding the overall π-electron delocalization. It is crucial to consider the unique characteristics of the three types of substituted N atoms and the spatial arrangement of 5MR and 6MR in N-PAHs.

Graphical Abstract

1. Introduction

Polycyclic aromatic hydrocarbons (PAHs) have been of great interest in the field of organic electronics for their uses in emitting fluorophores, sensors, transistors, and light-harvesting apparatus [1,2,3,4]. This is due to their unique electronic and optoelectronic properties, which can be effectively tuned by manipulating the number of fused aromatic rings and their geometrical arrangement [5,6,7]. The addition of fused aromatic moieties increases the total number of coupled π-electrons and can modify the π-conjugation pathway and π-electron delocalization, facilitated by the structural homology of the constituent aromatic rings [8]. Generally, with an increasing number of fused aromatic rings, π-electrons become more widely delocalized, causing energetic destabilization and stabilization of the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO), respectively. This leads to a decrease in transition energies of the lowest singlet (S1) and triplet (T1) states. However, it is known that increasing the number of fused aromatic rings has a limit in reducing the HOMO-LUMO gap [9]. In addition, two-dimensional PAHs have a limited number of geometrical configurations due to their inherent nature. Therefore, alternative techniques are needed to control the electronic characteristics of PAHs more effectively, instead of relying on conventional methods such as increasing fused rings or changing geometrical arrangements.
Nowadays, there have been various attempts to replace carbon atoms in PAHs with heteroatoms, such as nitrogen, sulfur, and oxygen atoms. This is performed to precisely adjust their structures and electronic properties similar to the use of dopants in silicon semiconductors [10,11,12,13,14,15]. Among these, nitrogen substitution has gained significant interest, as nitrogen-substituted PAHs (N-PAHs) are prevalent organic compounds, and several effective methods exist to dope nitrogen into graphene and PAHs, such as nitrogen plasma or chemical substitution [11,16,17,18]. There are at least three kinds of nitrogen-substituted six- and five-membered rings used as building blocks for N-PAHs with pyrrolic, pyridinic, and graphitic N atoms, which are widely found in popular organic compounds (Scheme 1) [19,20]. The five valence electrons of the nitrogen atom with the relatively high electron affinity can locally perturb π-electron distribution of the carbon-based PAHs. This raises an important issue of how the heteroatom effect of nitrogen atoms will impact the overall polycyclic (macroscopic) π-electron distribution and delocalization as it relates to the global electronic properties of N-PAHs.
We conducted a theoretical study on the linear N-PAHs that contain a mono-substituted nitrogen atom to investigate the effect of nitrogen atoms in PAHs. We utilized two types of N-substituted aromatic rings, pyridine and pyrrole, which were fused with six-membered aromatic hydrocarbon rings, as shown in Scheme 1. To avoid odd numbers of π-electrons and radicals in N-PAHs, we maintained a fixed number of covalent bonds at individual nitrogen atoms in five- and six-membered rings at three and two covalent bonds, respectively. Additionally, the linear N-PAHs are substituent-free because their π-electronic condition could be largely perturbed by peripheral substituents. As a result, there are three types of indole and carbazole isomers, as well as two types of quinoline and benzoquinoline isomers, respectively. The linear N-PAHs are employed to represent the three different types of local conditions of nitrogen atoms, namely pyridinic, pyrrolic, and graphitic nitrogen atoms. To investigate the π-electronic properties of the linear N-PAHs, we have employed various aromaticity indices based on optimized molecular structures. In indole and quinoline isomers, the π-electrons in the five- and six-membered rings with pyrrolic and pyridinic N atoms simultaneously contribute to both mono- and polycyclic π-electron conjugation, similar to six-membered aromatic hydrocarbons in naphthalene. On the other hand, the six-membered ring with a graphitic N atom in indolizine cannot cause monocyclic π-electron delocalization, and it seems to entirely contribute to the formation of polycyclic π-electron delocalization along the outward conjugation pathway that satisfies Hückel’s 4n + 2 rule. The six-membered rings with a graphitic N atom in indolizine and their expanded carbazole isomers exhibit low NICS(2) values, as well as bifurcated diatropic ring currents in the anisotropy of the induced current density (ACID) maps, indicating the absence of monocyclic π-electron delocalization and aromaticity. Conclusively, the number of π-electrons along the conjugation pathways and the satisfaction of Hückel’s 4n + 2 rule are critical factors in determining both mono- and polycyclic π-electron delocalization and the relative ratio of mono- and polycyclic π-electron distribution in the linear N-PAHs. The six-membered monocycle with a graphitic N atom cannot generate its own effective π-electron delocalization, which is obviously distinguished from the other aromatic polygons with pyrrolic and pyridinic N atoms.

2. Results and Discussion

2.1. Aromaticity Indices of Monocycles

In order to characterize the π-electronic systems of the linear N-PAHs, we have used various aromaticity indices to evaluate energetic, magnetic, and electronic properties, such as resonance energy and magnetic-induced ring current. Since linear N-PAHs consist of individual aromatic rings, we can examine the intricate changes in the π-electron environment by tracking alterations in the aromaticity indices of both individual constituent monocycles and extended polycycles. Comparing the relative aromaticity indices is useful for establishing the order of aromaticity through analog comparisons with similar molecular structures. However, it is worth noting that the aromaticity indices derived from optimized structures and those based on magnetic properties sometimes exhibit discrepancies. When comparing aromaticity indices of molecules with heteroatoms and/or different polygon structures, directly assigning the relative order of aromaticity and π-electron delocalization becomes challenging due to these inconsistencies in aromaticity indices. For instance, the structure-based HOMA and magnetic-based ACID critical isosurface values of monocycles exhibit a remarkably high correlation coefficient of 1.000 (Figure S2 in the Supplementary Materials). Moreover, the values for benzene (0.999, 0.082) and pyridine (1.000, 0.082) surpass those of pyrrole (0.872, 0.066), indicating that six-membered benzene and pyridine are more aromatic than five-membered pyrrole. However, the magnetic-based NICS(0) and NICS(1) values, which are widely employed to evaluate local aromaticity, do not align with the relative order observed in the HOMA and ACID values (Table 1). The highly negative NICS(0) and NICS(1) values are associated with a significant contribution from in-phase components [21], such as σ bonds and atoms around the NICS probe. Occasionally, we solely utilized the out-of-phase component of NICS to compare the strength of the induced diamagnetic ring current. Nevertheless, extracting the out-of-phase component necessitates additional computational effort and still does not align perfectly with the HOMA and ACID values. Despite both ACID and NICS calculations relying on the magnetic property of compounds, it remains challenging to establish the order of aromaticity.
We aimed to find a correlation between NICS and HOMA values based on the distance from the molecular mean plane of monocycles. We observed a significant negative correlation between NICS(2) and NICS(3) values with HOMA and ACID values (Figure 1). This negative correlation implies that the monocycles with higher HOMA values and more negative NICS values are more aromatic. Because the NICS(0) probe on the molecular mean plane in five-membered pyrrole is approximately 15.5% and 21.2% closer to the constituent atoms and covalent bonds, respectively, compared to six-membered benzene and pyridine (Figure S3 in the Supplementary Materials), the magnetic effects of the constituent atoms and covalent bonds in pyrrole are more significant. As Gershoni-Poranne et al. reported that NICSπZZ plots at a height lower than 1.7 Å show σ-effect contaminations [22], the NICS probe in five-membered pyrrole should be positioned further away from the molecular mean plane to make a more accurate comparison of the strength of diatropic ring currents. We conclude that five-membered pyrrole is less aromatic and less effective for π-electron delocalization than six-membered benzene and pyridine. This is likely due to the smaller ring size and reduced symmetry of the five-membered polygon with a pyrrolic N atom, leading to highly concentrated π-electrons and substantially destabilized HOMO and HOMO-1 in comparison to benzene and pyridine. We utilized NICS(2) values to evaluate the strength of diamagnetic-induced ring currents in the linear N-PAHs as a polygon(shape)-independent index of local aromaticity and π-electron delocalization.

2.2. Visualized Diatropic Ring Currents of Bi- and Tricyclic N-PAHs

We investigated the effect of a substituted nitrogen atom in linear PAHs by comparing the aromaticity indices of bicyclic indole and quinoline isomers. To determine the extent of electron delocalization in fused bicyclic PAHs containing a single nitrogen atom, we theoretically visualized magnetic-induced current maps (Figure 2 and Figure S5 in the Supplementary Materials). With the exception of indolizine, the magnetic-induced ring currents of indole and quinoline isomers resembled that of naphthalene, exhibiting outward unidirectional diatropic ring currents (clockwise) characteristic of a conjugated 10-electron system. In contrast, indolizine displayed distinct bifurcated diatropic ring currents, with one flowing around the N-substituted five-membered ring and the other flowing outward (Figure 2). This bifurcated diatropic ring current observed in indolizine indicates the similar magnitudes of the two diatropic ring currents circulating the five-membered monocycle and the nine-membered polycycle, revealing the local- and poly/macrocyclic aromaticity, respectively. Notably, despite the fusion of rings, the N-substituted five-membered ring in indolizine maintained discernible pyrrolic characteristics.
To investigate whether the occurrence of bifurcated diatropic ring currents is common, we expanded the sample series to include linear tricyclic N-PAHs (Scheme 1). By visualizing induced ring current density maps, we confirmed that two carbazole isomers, isocarbazole and pseudocarbazole with an indolizine moiety, exhibit bifurcated diatropic ring currents (Figure 2 and Figure S6 in the Supplementary Materials). This observation suggests that bifurcated diatropic ring currents in linear N-PAHs with an indolizine moiety are not uncommon. We also noted that the presence of a constituent six-membered ring with a graphitic N atom in linear N-PAHs particularly promotes the formation of such bifurcated diatropic ring currents. The role of a graphitic N atom in the π-electron delocalization of PAHs will be discussed in the subsequent section.

2.3. How to Quantify the Bifurcated Diatropic Ring Currents

To evaluate the bifurcated diatropic ring currents, we compared the HOMA and NICS(2) values of indole and carbazole isomers. The HOMA values of the constituent five- and six-membered rings in indole and carbazole isomers were lower than those of their parent monocycles. However, their NICS(2) values generally increased, indicating increased diamagnetic ring currents due to extended π-conjugation and enriched π-electron delocalization caused by the fusion of rings (Table 2 and Table S1 in the Supplementary Materials). The structural deformation of the individual rings in indole and carbazole isomers is directly related to the formation of poly/macrocyclic π-electron conjugation and delocalization, which corresponds well with high HOMA values along outward polycyclic pathways, resembling 10- and 14π-electronic systems.
The π-conjugation pathways of individual constituent rings in the N-PAHs are spatially overlaid with the outward polycyclic conjugation pathways, making it challenging to accurately determine the monocyclic characteristics based on the structure-based HOMA values of the constituent rings alone, even with the selection of monocyclic covalent bonds. However, NICS(2) provides a means to assess the relative diatropic ring currents of the individual constituent rings, providing collective information about both mono- and polycyclic π-electronic systems at the assigned individual rings. To elucidate the distinct nature of the bifurcated diatropic ring currents in indolizine, we refocused on bicyclic indole and quinoline isomers. In bicyclic indole and quinoline isomers, there are only two possible π-conjugation pathways, monocyclic and outward bicyclic. Linear bicyclic N-PAHs form closed loops, resulting in equal inflow and outflow of polycyclic current density in neighboring constituent rings. If the NICS(2) values of individual constituent rings are polygon-independent, we can deduce that the differences of NICS(2) values between neighboring rings in indole isomers are attributable to varying strengths of monocyclic ring currents. This is because the contribution of the polycyclic ring current cancels out when calculating the difference of NICS(2) values. The nearly identical NICS(2) values of constituent rings in quinoline isomers, similar to those of naphthalene, indicate that substituting a carbon atom with a pyridinic nitrogen atom does not disrupt the mono- and polycyclic π-electronic systems. Conversely, the indole isomers exhibit larger differences in NICS(2) values by 6 to 12 times, suggesting heterogeneous π-electronic conditions for the individual constituent rings by the fusion between five- and six-membered rings (Table 2).

2.4. Structural Disparities between the N-5MR and C-6MR

The NICS(2) values of the C-6MR in indole and N-5MR in isoindole are higher than those of the neighboring polygonal N-5MR and C-6MR, respectively. Pino-Rios et al. proposed the Glidewell–Lloyd rule to explain the relatively larger aromaticity indices of 5MR, which prevent the formation of the smallest antiaromatic 4n group [23]. However, when using NICS(2) values to evaluate local aromaticity and π-electron delocalization, the higher NICS(2) value of the C-6MR in indole does not conform to the Glidewell–Lloyd rule. In this context, we need to consider the structural uniqueness of the parent N-5MR (pyrrole) and the location of the shared covalent bond between the N-5MR and 6MR. Although the fused bicyclic N-PAHs are distinct molecules compared to parent monocycles, the constituent 5MR and 6MR in N-PAHs likely retain their monocyclic characteristics, representing the concept of local aromaticity. If we assume that pyrrole is the ideal aromatic structure of N-5MR for the monocyclic π-electron delocalization, then we should compare the structural changes in the N-5MR in PAHs with the structure of pyrrole. Since pyrrole is not a regular pentagon, we can use this comparison to evaluate its monocyclic π-electron delocalization. The different locations of shared bonds in indole isomers result in heterogeneous effects on the combination of mono- and polycyclic π-electron delocalization. By examining the differing bond lengths of individual monocyclic rings in PAHs compared to parent monocycles, it is possible to evaluate the structural effect caused by the fusion of rings. For ease of use, we have chosen to use Pauling bond numbers derived from the lengths of C-C and C-N bonds to measure the structural deviation, as they allow for the comparison of π-electronic contribution for individual covalent bonds [24]. For instance, although the bond length of the shared covalent bond in naphthalene is 2.73% longer than that of benzene, its Pauling bond number is reduced by 15.1%, directly impacting the weakening of monocyclic (local) π-electronic characteristics and the introduction of additional polycyclic π-electron delocalization. Parent pyrrole has three types of covalent bonds with different Pauling bond numbers of 1.740, 1.314, and 1.398, which correspond to the shared covalent bond locations in indole, isoindole, and indolizine, respectively (Figure S7 in the Supplementary Materials). By comparing the Pauling bond numbers of shared covalent bonds in indole isomers with those of parent monocycles, we can evaluate the deviation of Pauling bond numbers individually for shared covalent bonds at the 6MR and 5MR, and the trend of reduced bond numbers aligns well with the trend of NICS(2) values. For instance, the Pauling bond numbers of shared covalent bonds in indole isomers are 1.330, 1.126, and 1.207 (Figure S8 in the Supplementary Materials). The Pauling bond numbers of shared covalent bonds are reduced by 23.6%, 14.3%, and 13.7% for N-5MR, and reduced by 15.2%, 28.2%, and 24.3% for 6MR in indole, isoindole, and indolizine, respectively. The 5MR in indole and the 6MRs in isoindole and indolizine are structurally more deviated from parent monocycles, and this correlates with the relative order of NICS(2) values for the individual constituent rings. The substantial decrease in Pauling bond numbers compared to parent monocycles indicates a discriminatively attenuated monocyclic 6π-electron delocalization and the emergence of newly formed polycyclic 10π-electron delocalization through the fusion of rings. Conclusively, the locations of shared covalent bonds in indole isomers are critical for the selection of relatively weak monocyclic π-electronic characteristics of the constituent N-5MR and 6MR.

2.5. Graphitic N Atom in Linear N-PAHs

The NICS(2) values of the constituent N-5MR and C-6MR in indole and isoindole differ by 0.88 and 0.78 ppm, respectively. In indolizine, the difference is approximately twice as large, 1.59 ppm. It is noteworthy that the NICS(2) value of the N-6MR in indolizine is significantly reduced compared to that in isoindole, despite having a larger HOMA value (Table 2). Interestingly, the N-6MR in indolizine with a graphitic N atom cannot satisfy Hückel’s 4n + 2 rule, which could be confirmed by localized π-electrons in a virtual six-membered monocycle with a graphitic N atom (Figure S7 in the Supplementary Materials). On the other hand, the outward polycyclic conjugation pathway in indolizine satisfies Hückel’s 4n + 2 rule, implying that the N-6MR in indolizine favors the formation of an outward polycyclic 10π-electronic system while minimally contributing to the monocyclic π-electronic system of the N-6MR. Conclusively, the N-5MR in indolizine contributes to both 6π monocyclic delocalization and 10π polycyclic delocalization with N-6MR, which manifests as bifurcated ring currents. The loss of monocyclic π-electronic characteristics of N-6MR in indolizine prevents collision between opposite monocyclic diatropic ring currents at the shared covalent bond. Similarly, the N-6MR with a graphitic N atom in isocarbazole and peudocarbazole cannot generate a monocyclic diatropic ring current, thereby leading to bifurcated diatropic ring currents of the 10- and 14π-electronic systems. The single N atom in an indolizine moiety simultaneously acts as a pyrrolic N atom in the N-5MR and a graphitic N atom in the N-6MR.

2.6. Coexistence of Mono- and Polycyclic π-Electronic Systems in Linear N-PAHs

When linear N-PAHs are expanded by fusing two or three constituent rings, a unique type of six-membered monocycle with a graphitic nitrogen atom is produced. This particular feature is present in isomers of indole and carbazole, which highlights local π-electronic characteristics of its neighboring mono- and (sub-)polycycles in linear N-PAHs. This demonstrates the coexistence of mono- and polycyclic electron delocalization in indolizine and the two carbazole isomers, shown as bifurcated diatropic ring currents (Scheme 2). It is important to note that fulfilling Hückel’s 4n + 2 rule is crucial for both constituent rings and fused polycycles with two or three rings to comprehend the overall electronic characteristics. The S1 and T1 states, which are the lowest excited states of linear N-PAHs, are primarily governed by the HOMO-LUMO transition. Therefore, the longest outward conjugation pathway that satisfies Hückel’s 4n + 2 rule is of paramount importance in explaining the representative electronic dynamics in these lowest excited states. Nevertheless, it is still fundamentally significant to understand the local characteristics exhibited by constituent monocycles and sub-polycycles with a smaller number of constituent rings, as overall thermodynamic stability and chemical reactivity are mainly determined by these local characteristics. Moreover, the relative contribution of π-electron delocalization along the longest conjugation path is also determined by these local characteristics. This is because the structural deformation for the longest conjugation path is constrained by the established local energetic benefits and delocalization by constituent mono- and (sub-)polycycles. Therefore, the polarizability and optical properties of the electrons should be considered based on the coexistence of mono- and polycyclic electronic systems in linear N-PAHs.

3. Materials and Methods

The molecular structures of linear PAHs and monocycles were optimized using the density functional theory (DFT) method with Becke’s three-parameter hybrid exchange functionals and the Lee–Yang–Parr correlation functional (B3LYP), employing the 6-311G(2df,p) basis set [25,26,27]. The calculations were performed on a supercomputer (KISTI, Nurion, Daejeon, Korea). Various aromaticity indices were utilized to characterize the π-electrons in linear PAHs. Nuclear-independent chemical shifts (NICS(x)) were calculated to check the isotropic shielding by constituent atoms on the monocyclic π-conjugation pathways at x Å above the molecular planes [21]. We confirmed that the experimental NMR results of indole and indolizine were well matched with the calculated NMR results based on their optimized molecular structures with GIAO-B3LYP/6-311G(2df,p) (Figure S1 in the Supplementary Materials). To plot the anisotropy of the induced current density (ACID) maps, the continuous set of gauge transformations (CSGT) method was employed to calculate the magnetically induced current densities [28,29]. The results were visualized using POVRAY 3.7 for Windows [30]. In order to compare the structure-based aromaticity of the mono- and polycyclic rings in linear PAHs, the harmonic oscillator model of aromaticity (HOMA) values were calculated [31]. To establish a practical nomenclature for monocycles, we employed concise symbols, namely a six-membered ring (6MR), a five-membered ring (5MR), a carbon atom-exclusive cycle (C-), and a mono-nitrogen-atom-substituted cycle (N-). Therefore, benzene, pyridine, and pyrrole are denoted as C-6MR, N-6MR, and N-5MR, respectively.

4. Conclusions

We conducted a theoretical study on the substitution effect of a single nitrogen atom in linear bi- and tricyclic N-PAHs to better understand electron delocalization and bifurcated diatropic ring currents. Using various aromaticity indices for linear N-PAHs, we evaluated the magnetic-induced ring currents in these molecules by using NICS(2) values for individual constituent monocycles. We looked at the substitution effect of three types of N atoms, pyrrolic, pyridinic, and graphitic, taking into consideration both the polygon shape and the number of covalent bonds at the substituted N atom. Despite deviations from the ideal structures of parent monocycles, our analysis showed that the local electron delocalization of subcycles (such as mono- and bicyclic constituent moieties) in linear N-PAHs is preserved, resulting in a less effective electron delocalization within the subcycles, which is related to the emergence of extended global polycyclic π-electron delocalization. In the case of N-6MR with a graphitic N atom, it cannot solely generate monocyclic π-electron delocalization, but contributes to the formation of polycyclic electron delocalization, shown by the bifurcated diatropic ring currents. This work highlights the significance of the satisfaction of Hückel’s 4n + 2 rule for both sub- and macrocycles in comprehending the overall π-electron delocalization, as well as the optimized molecular structures and spatial distribution of π-electrons in linear N-PAHs. To understand the mono- and polycyclic electronic properties, it is crucial to consider the unique characteristics of the three types of substituted N atoms and the spatial arrangement of 5MR and 6MR in N-PAHs. This study of linear bi- and tricyclic N-PAHs will be a milestone for understanding the relationship between the molecular structure of extended N-PAHs and their π-electronic properties.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules29040784/s1, Table S1. Aromaticity indices of the linear tricyclic N-PAHs; Figure S1: H-NMR comparisons between experimental and theoretical (GIAO-B3LYP/6-311G(2df,p)) results of indole and indolizine; Figure S2: Correlation between HOMA and ACID values of monocycles; Figure S3: Averaged distance from atoms and covalent bonds to the center (NICS probe); Figure S4: Magnetic-induced current maps of six- and five-membered monocycles (benzene, pyridine, and pyrrole); Figure S5: Magnetic-induced current maps of fused bicycles with constituent six- and five-membered rings; Figure S6: Magnetic-induced current maps of fused tricycles with constituent six- and five-membered rings; Figure S7: Bond lengths (black), Pauling bond numbers (red), and HOMA parameters of monocycles and N-6MRs with different charge states; Figure S8: Bond lengths (black), Pauling bond numbers (red), and HOMA parameters of bicycles; Figure S9: Bond lengths (black), Pauling bond numbers (red), HOMA parameters of tricycles, and optimized geometry of molecular structures.

Author Contributions

Conceptualization, J.M.L., J.K. and S.C.; methodology, S.S. and H.T.B.; software, Y.H. and S.C.; analysis, J.M.L. and H.-J.K.; writing—original draft preparation, J.M.L., J.K. and S.C.; writing—review and editing, H.-J.K., S.S. and Y.H.; visualization, S.C.; project administration, Y.H. and S.C.; funding acquisition, J.M.L. and S.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was financially supported by Chonnam National University (2019-3807) and the Basic Science Research Program through the National Research Foundation of Korea (NRF), funded by the Ministry of Science and ICT (NRF-2022R1F1A1074891 and RS-2023-00279712). The theoretical calculations were carried out with the support of the National Institute of Supercomputing and Network/Korea Institute of Science and Technology Information for supercomputing resources, including technical support (KSC-2020-CRE-0077).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article and Supplementary Materials.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Zhang, L.; Cao, Y.; Colella, N.S.; Liang, Y.; Brédas, J.-L.; Houk, K.N.; Briseno, A.L. Unconventional, Chemically Stable, and Soluble Two-Dimensional Angular Polycyclic Aromatic Hydrocarbons: From Molecular Design to Device Applications. Acc. Chem. Res. 2015, 48, 500–509. [Google Scholar] [CrossRef]
  2. Schulze, M.; Philipp, M.; Waigel, W.; Schmidt, D.; Würthner, F. Library of Azabenz-Annulated Core-Extended Perylene Derivatives with Diverse Substitution Patterns and Tunable Electronic and Optical Properties. J. Org. Chem. 2016, 81, 8394–8405. [Google Scholar] [CrossRef]
  3. Sun, X.; James, T.D.; Anslyn, E.V. Arresting “Loose Bolt” Internal Conversion from −B(OH)2 Groups Is the Mechanism for Emission Turn-on in Ortho-Aminomethylphenylboronic Acid-Based Saccharide Sensors. J. Am. Chem. Soc. 2018, 140, 2348–2354. [Google Scholar] [CrossRef]
  4. Bayn, A.; Feng, X.; Müllen, K.; Haick, H. Field Effect Transistors Based on Polycyclic Aromatic Hydrocarbons for the Detection and Classification of Volatile Organic Compounds. ACS Appl. Mater. Interfaces 2013, 5, 3431–3440. [Google Scholar] [CrossRef]
  5. Dias, J.R. A Periodic Table for Polycyclic Aromatic Hydrocarbons. Acc. Chem. Res. 1985, 18, 241–248. [Google Scholar] [CrossRef]
  6. Nguyen, L.H.; Truong, T.N. Quantitative Structure–Property Relationships for the Electronic Properties of Polycyclic Aromatic Hydrocarbons. ACS Omega 2018, 3, 8913–8922. [Google Scholar] [CrossRef] [PubMed]
  7. Gutzler, R.; Perepichka, D.F. Π-Electron Conjugation in Two Dimensions. J. Am. Chem. Soc. 2013, 135, 16585–16594. [Google Scholar] [CrossRef] [PubMed]
  8. Randić, M. Aromaticity of Polycyclic Conjugated Hydrocarbons. Chem. Rev. 2003, 103, 3449–3606. [Google Scholar] [CrossRef] [PubMed]
  9. Menon, A.; Dreyer, J.A.H.; Martin, J.W.; Akroyd, J.; Robertson, J.; Kraft, M. Optical Band Gap of Cross-Linked, Curved, and Radical Polyaromatic Hydrocarbons. Phys. Chem. Chem. Phys. 2019, 21, 16240–16251. [Google Scholar] [CrossRef] [PubMed]
  10. Borissov, A.; Maurya, Y.K.; Moshniaha, L.; Wong, W.-S.; Żyła-Karwowska, M.; Stępień, M. Recent Advances in Heterocyclic Nanographenes and Other Polycyclic Heteroaromatic Compounds. Chem. Rev. 2022, 122, 565–788. [Google Scholar] [CrossRef]
  11. Alam, K.; Hong, S.W.; Oh, K.H.; Park, J.K. Divergent C–H. Annulation for Multifused N-Heterocycles: Regio- and Stereospecific Cyclizations of N-Alkynylindoles. Angew. Chem. Int. Ed. 2017, 56, 13387–13391. [Google Scholar] [CrossRef]
  12. Stassen, D.; Demitri, N.; Bonifazi, D. Extended O-Doped Polycyclic Aromatic Hydrocarbons. Angew. Chem. Int. Ed. 2016, 55, 5947–5951. [Google Scholar] [CrossRef]
  13. Delaunay, W.; Szűcs, R.; Pascal, S.; Mocanu, A.; Bouit, P.A.; Nyulászi, L.; Hissler, M. Synthesis and Electronic Properties of Polycyclic Aromatic Hydrocarbons Doped with Phosphorus and Sulfur. Dalton Trans. 2016, 45, 1896–1903. [Google Scholar] [CrossRef]
  14. Ashraf, R.S.; Meager, I.; Nikolka, M.; Kirkus, M.; Planells, M.; Schroeder, B.C.; McCulloch, I. Chalcogenophene Comonomer Comparison in Small Band Gap Diketopyrrolopyrrole-Based Conjugated Polymers for High-Performing Field-Effect Transistors and Organic Solar Cells. J. Am. Chem. Soc. 2015, 137, 1314–1321. [Google Scholar] [CrossRef]
  15. Kumar, G.J.; Bogoslavsky, B.; Debnath, S.; Bedi, A. Effect of Chalcogenophenes on Chiroptical Activity of Twisted Tetracenes: Computational Analysis, Synthesis and Crystal Structure Thereof. Molecules 2023, 28, 5074. [Google Scholar] [CrossRef] [PubMed]
  16. Tian, X.; Shoyama, K.; Würthner, F. Nitrogen-Doped Polycyclic Aromatic Hydrocarbons by a One-Pot Suzuki Coupling/Intramolecular Snar Reaction. Chem. Sci. 2023, 14, 284–290. [Google Scholar] [CrossRef] [PubMed]
  17. Boas, C.R.S.V.; Focassio, B.; Marinho, E.; Larrude, D.G.; Salvadori, M.C.; Leão, C.R.; dos Santos, D.J. Characterization of Nitrogen Doped Graphene Bilayers Synthesized by Fast, Low Temperature Microwave Plasma-Enhanced Chemical Vapour Deposition. Sci. Rep. 2019, 9, 13715. [Google Scholar] [CrossRef] [PubMed]
  18. Yen, P.-J.; Ilango, P.R.; Chiang, Y.-C.; Wu, C.-W.; Hsu, Y.-C.; Chueh, Y.-L.; Wei, K.-H. Tunable Nitrogen-Doped Graphene Sheets Produced with in Situ Electrochemical Cathodic Plasma at Room Temperature for Lithium-Ion Batteries. Mater. Today Energy 2019, 12, 336–347. [Google Scholar] [CrossRef]
  19. Sarkar, S.; Sudolská, M.; Dubecký, M.; Reckmeier, C.J.; Rogach, A.L.; Zbořil, R.; Otyepka, M. Graphitic Nitrogen Doping in Carbon Dots Causes Red-Shifted Absorption. J. Phys. Chem. C 2016, 120, 1303–1308. [Google Scholar] [CrossRef]
  20. Zhuang, S.; Lee, E.S.; Lei, L.; Nunna, B.B.; Kuang, L.; Zhang, W. Synthesis of Nitrogen-Doped Graphene Catalyst by High-Energy Wet Ball Milling for Electrochemical Systems. Int. J. Energy Res. 2016, 40, 2136–2149. [Google Scholar] [CrossRef]
  21. Chen, Z.; Wannere, C.S.; Corminboeuf, C.; Puchta, R.; Schleyer, P.v.R. Nucleus-Independent Chemical Shifts (Nics) as an Aromaticity Criterion. Chem. Rev. 2005, 105, 3842–3888. [Google Scholar] [CrossRef] [PubMed]
  22. Gershoni-Poranne, R.; Stanger, A. The Nics-Xy-Scan: Identification of Local and Global Ring Currents in Multi-Ring Systems. Chem. A Eur. J. 2014, 20, 5673–5688. [Google Scholar] [CrossRef] [PubMed]
  23. Pino-Rios, R.; Solà, M. The Relative Stability of Indole Isomers Is a Consequence of the Glidewell-Lloyd Rule. J. Phys. Chem. A 2021, 125, 230–234. [Google Scholar] [CrossRef] [PubMed]
  24. Pauling, L. Atomic Radii and Interatomic Distances in Metals. J. Am. Chem. Soc. 1947, 69, 542–553. [Google Scholar] [CrossRef]
  25. Becke, A.D. A New Mixing of Hartree–Fock and Local Density-Functional Theories. J. Chem. Phys. 1993, 98, 1372–1377. [Google Scholar] [CrossRef]
  26. Lee, C.; Yang, W.; Parr, R.G. Development of the Colle-Salvetti Correlation-Energy Formula into a Functional of the Electron Density. Phys. Rev. B 1988, 37, 785–789. [Google Scholar] [CrossRef]
  27. Cho, S.; Yoon, Z.S.; Kim, K.S.; Yoon, M.-C.; Cho, D.-G.; Sessler, J.L.; Kim, D. Defining Spectroscopic Features of Heteroannulenic Antiaromatic Porphyrinoids. J. Phys. Chem. Lett. 2010, 1, 895–900. [Google Scholar] [CrossRef]
  28. Geuenich, D.; Hess, K.; Köhler, F.; Herges, R. Anisotropy of the Induced Current Density (Acid), a General Method to Quantify and Visualize Electronic Delocalization. Chem. Rev. 2005, 105, 3758–3772. [Google Scholar] [CrossRef]
  29. Stanger, A. Nucleus-Independent Chemical Shifts (Nics):  Distance Dependence and Revised Criteria for Aromaticity and Antiaromaticity. J. Org. Chem. 2006, 71, 883–893. [Google Scholar] [CrossRef]
  30. POVRAY, version 3.7 beta; Persistence of Vision Raytracer Pty: Williamstown, Australia, 2006.
  31. Krygowski, T.M.; Cyrański, M.K. Structural Aspects of Aromaticity. Chem. Rev. 2001, 101, 1385–1420. [Google Scholar] [CrossRef] [PubMed]
Molecules 29 00784 sch001
Scheme 1. Expansion of N-substituted six- and five-membered aromatic rings with CH-hexagonal aromatic rings (the solid and dashed lines denote unidirectional and mixed expansion, respectively, depending on the location of shared bonds).
Scheme 1. Expansion of N-substituted six- and five-membered aromatic rings with CH-hexagonal aromatic rings (the solid and dashed lines denote unidirectional and mixed expansion, respectively, depending on the location of shared bonds).
Molecules 29 00784 sch001
Molecules 29 00784 g001
Figure 1. (a) Correlation between HOMA and NICS(x) values of parent monocycles and (b) correlation coefficients depending on the distance from the molecular mean plane.
Figure 1. (a) Correlation between HOMA and NICS(x) values of parent monocycles and (b) correlation coefficients depending on the distance from the molecular mean plane.
Molecules 29 00784 g001
Molecules 29 00784 g002
Figure 2. Magnetic-induced current maps of indole and carbazole isomers (clockwise vectors: diatropic ring current).
Figure 2. Magnetic-induced current maps of indole and carbazole isomers (clockwise vectors: diatropic ring current).
Molecules 29 00784 g002
Molecules 29 00784 sch002
Scheme 2. Coexistence of mono- and polycyclic π-electronic systems with bifurcated and singular diatropic ring currents.
Scheme 2. Coexistence of mono- and polycyclic π-electronic systems with bifurcated and singular diatropic ring currents.
Molecules 29 00784 sch002
Table 1. Aromaticity indices of aromatic monocycles.
Table 1. Aromaticity indices of aromatic monocycles.
CompoundHOMAACID Critical Isosurface ValueNICS(0)NICS(1)NICS(2)NICS(3)
Benzene
(C-6MR)		0.9990.082−8.44−10.9−5.01−1.96
Pyridine
(N-6MR)		1.0000.082−7.11−10.6−4.93−1.93
Pyrrole
(N-5MR)		0.8720.066−14.0−11.0−4.22−1.51
Table 2. Aromaticity indices of aromatic bicycles.
Table 2. Aromaticity indices of aromatic bicycles.
CompoundLocationHOMAACID Critical Isosurface ValueNICS(2)
NaphthaleneC-6MR0.8120.076−5.50
C-10MR0.8570.085
IndoleN-5MR0.7470.056−4.51
C-6MR0.9370.065−5.39
N-9MR0.8710.069
IsoindoleN-5MR0.7870.078−5.35
C-6MR0.6790.072−4.57
N-9MR0.8650.079
IndolizineN-5MR0.8400.081−5.30
N-6MR0.7370.072−3.71
N-9MR0.8630.072
QuinolineN-6MR0.8450.072−5.43
C-6MR0.8220.070−5.50
N-10MR0.8770.084
IsoquinolineN-6MR0.8500.078−5.42
C-6MR0.8370.078−5.55
N-10MR0.8730.082
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lim, J.M.; Shim, S.; Bui, H.T.; Kim, J.; Kim, H.-J.; Hwa, Y.; Cho, S. Substitution Effect of a Single Nitrogen Atom on π-Electronic Systems of Linear Polycyclic Aromatic Hydrocarbons (PAHs): Theoretically Visualized Coexistence of Mono- and Polycyclic π-Electron Delocalization. Molecules 2024, 29, 784. https://doi.org/10.3390/molecules29040784

AMA Style

Lim JM, Shim S, Bui HT, Kim J, Kim H-J, Hwa Y, Cho S. Substitution Effect of a Single Nitrogen Atom on π-Electronic Systems of Linear Polycyclic Aromatic Hydrocarbons (PAHs): Theoretically Visualized Coexistence of Mono- and Polycyclic π-Electron Delocalization. Molecules. 2024; 29(4):784. https://doi.org/10.3390/molecules29040784

Chicago/Turabian Style

Lim, Jong Min, Sangdeok Shim, Hoa Thi Bui, Jimin Kim, Ho-Joong Kim, Yoon Hwa, and Sung Cho. 2024. "Substitution Effect of a Single Nitrogen Atom on π-Electronic Systems of Linear Polycyclic Aromatic Hydrocarbons (PAHs): Theoretically Visualized Coexistence of Mono- and Polycyclic π-Electron Delocalization" Molecules 29, no. 4: 784. https://doi.org/10.3390/molecules29040784

APA Style

Lim, J. M., Shim, S., Bui, H. T., Kim, J., Kim, H. -J., Hwa, Y., & Cho, S. (2024). Substitution Effect of a Single Nitrogen Atom on π-Electronic Systems of Linear Polycyclic Aromatic Hydrocarbons (PAHs): Theoretically Visualized Coexistence of Mono- and Polycyclic π-Electron Delocalization. Molecules, 29(4), 784. https://doi.org/10.3390/molecules29040784

Article Metrics

Back to TopTop