Aromatic Stabilization Energy Research Articles

Comparison of Various Theoretical Measures of Aromaticity within Monosubstituted Benzene.

The effects of monosubstitution on the aromaticity of benzene are assessed using a number of different quantitative schemes. The ability of the mobile π-electrons to respond to an external magnetic field is evaluated using several variants of the NICS scheme which calculate the shielding of points along the axis perpendicular to the molecule. Another class of measures is related to the drive toward the uniformity of C-C bond lengths and strengths. Several energetic quantities are devised to approximate an aromatic stabilization energy and the tendency of the molecule to maintain planarity. There is a lack of consistency in that the various measures of aromaticity lead to differing conclusions as to the effects of substituents on the aromaticity of the ring.

Use of Ferrocenyl Ni(II) and Zn(II) Porphyrins as Active Organic Electrode Materials for Sodium Secondary Batteries

AbstractThis study investigates the potential use of ferrocenyl‐substituted porphyrins as organic electrode materials for sodium secondary batteries. Expressly, 5,15‐di(ferrocenyl)‐10,20‐di(pentafluorophenyl)porphyrins (free base, Fc−H2Por), along with its zinc (Zn[II] complex, Fc−ZnPor) and nickel (Ni[II] complex, Fc−NiPor) derivatives were synthesized and employed as the electrode materials. The charge‐discharge processes disclosed that the ferrocenyl‐substituted porphyrin compound performs electrochemical processes poorly without chelating metals. In contrast, when the ferrocenyl‐substituted porphyrins were metalated with Ni(II) and Zn(II), they exhibited dual ionic charge‐discharge capabilities, indicating their stable electrochemical properties. In particular, long‐term cycle tests emphasized the superior stability and high Coulombic efficiency of Fc−ZnPor electrodes. Computational calculations provided conclusive evidence of variances in aromatic stabilization energy based on the metal center, significantly influencing the chemical stability and electrochemical performance. These findings can deliver valuable insight for further development and optimizations in energy storage applications.

Visible-Light-Promoted Tandem Skeletal Rearrangement/Dearomatization of Heteroaryl Enallenes

Access to complex three-dimensional molecular architectures via dearomatization of ubiquitous aryl rings is a powerful synthetic tool, which faces, however, an inherent challenge to overcome energetic costs due to the loss of aromatic stabilization energy. Photochemical methods that allow one to populate high-energy states can thus be an ideal strategy to accomplish otherwise prohibitive reaction pathways. We present an original dearomative rearrangement of heteroaryl acryloylallenamides that leads to complex fused tricycles. The visible-light-promoted method occurs under mild conditions and tolerates a variety of functional groups. According to DFT modeling used to rationalize the outcome of the cascade, the reaction involves a sequential [2+2] allene–alkene photocycloaddition, which is followed by a selective retro- [2+2] step that paves the way for the dearomatization of the heteroaryl partner. This scenario is original with respect to the reported photochemical reactivity of similar substrates and thus holds promise for ample future developments.

Update for Isomerization Stabilization Energies: The Fulvenization Approach.

An alternative approach for calculating aromatic stabilization energies is proposed based on transforming an (anti)aromatic ring into a fulvene isomer. This fulvenization process gives a value of 34.05 kcal·mol-1 for benzene in the singlet state and a value of -17.85 kcal·mol-1 in the triplet state. Additionally, it is possible to use experimental values (as long as they exist) for the calculation as the gas-phase formation enthalpies of benzene and fulvene, whose difference is 33.72 kcal·mol-1. On the other hand, this same approach has been evaluated on several six-membered rings, including those persubstituted, biradicals, azines, and inorganic analogues, giving results in agreement with those reported in the literature using different criteria. Additionally, it is possible to differentiate the aromaticity of the rings in polycyclic aromatic hydrocarbons according to Clar's rules. Assigning the (anti)aromatic character in various nonbenzenoid rings (neutral and charged), except for five- and seven-membered rings, is also possible. The construction of the fulvene isomers in PAHs is set such that nonaromaticity-related effects are not considered. The results show that the fulvenization approach is an effective and efficient approach that can serve as an alternative or complement to existing tools.

Sufficiently Enriched Dual‐Ion Batteries with Ferrocenyl Substituted Nickel(II) Norcorrole Organic Electrodes

AbstractFerrocenyl‐substituted nickel(II) norcorrole and nickel(II) porphyrin (Fc2NiNc and Fc2NiPor, respectively) are prepared as active electrode materials for dual‐ion batteries, which enables the systematically study of aromaticity and antiaromaticity for organic electrode materials. Intrinsically, Fc2NiNc and Fc2NiPor electrodes have excellent electrochemical behaviors with dual ion interaction with Li cations and bis(fluorosulfonyl)amide [FSA]– inorganic anions from the ionic liquid electrolyte, inducing high capacities (259 mAh g−1 for Fc2NiNc and 200 mAh g−1 for Fc2NiPor). The dominance of pseudocapacitive charge storage is exhibited by the Fc2NiNc and Fc2NiPor electrodes giving rise to rapid charge‐discharge processes. Computational calculation conclusively demonstrates the mechanism of multi‐electron redox reaction resulting in high discharge capacities and aromatic stabilization energy toward achieving a stable long‐term cyclability over 1000 cycles at high current density. Furthermore, Li4Ti5O12/IL/Fc2NiNc cells exhibit high power and energy densities (234 Wh kg−1 at 118 W kg−1) and accomplish outstanding cyclability over 10 000 cycles. This work provides a strategy for designing robust organic electrode materials for dual‐ion batteries.

Predicting a New Δ-Proton Sponge-Base of 4,12-Dihydrogen-4,8,12-triazatriangulene through Proton Affinity, Aromatic Stabilization Energy, and Aromatic Magnetism.

Herein, we report designing a new Δ (delta-shaped) proton sponge base of 4,12-dihydrogen-4,8,12-triazatriangulene (compound 1) and calculating its proton affinity (PA), aromatic stabilization, natural bond orbital (NBO), electron density ρ(r), Laplacian of electron density ∇2 ρ(r), (2D-3D) multidimensional off-nucleus magnetic shielding (σzz (r) and σiso (r)), and scanning nucleus-independent chemical shift (NICSzz and NICS). Density functional theory (DFT) at B3LYP/6-311+G(d,p), ωB97XD/6-311+G(d,p), and PW91/def2TZVP were used to compute the magnetic shielding variables. In addition, relevant bases like pyridine, quinoline, and acridine were also studied and compared. The protonation of compound 1 yields a highly symmetric carbocation of three Hückel benzenic rings. Comparing our findings of the studied molecules showed that compound 1 precedes others in PA, aromatic isomerization stabilization energy, and basicity. Therefore, the basicity may be enhanced when a conjugate acid gains higher aromatic features than its unprotonated base. Both multidimensional σzz (r) and σiso (r) off-nucleus magnetic shieldings outperformed electron-based techniques and can visually monitor changes in aromaticity that occur by protonation. The B3LYP/6-311+G(d,p), ωB97XD/6-311+G(d,p), and PW91/def2TZVP levels showed no significant differences in detailing isochemical shielding surfaces.

Analysis of Bonding by Quantum Chemistry─Resolving Delocalization Stabilization in a Mechanistic Basis and New Hückel Model

No general and unique understanding of the mechanism of covalent bonding in physical terms is provided by current computational methods or by a consensus among experts. Bonding is studied by energy decomposition analysis but may also be related to the interatomic motion of valence electrons within the molecule. This dynamical view of the mechanism of bonding is not widely appreciated. The aim here is to make it accessible by translation into a corresponding form of quantum chemical energy analysis. The interatomic electron motion is directly related to the delocalization taking place when atomic basis functions are combined into molecular orbitals. A "tribasis method" is introduced, allowing an atomic basis set to form subsets of (1) strictly localized atomic functions and (2) interatomic bridge functions which allow delocalization. Calculations can then identify ground states without (no bridge functions) and with delocalization. The scheme is based on exact quantum mechanics but demonstrated by a minimal basis treatment of H2+ and H2 in Hartree-Fock and valence bond approximations which show that the bond energy is a sum of repulsive localization and more strongly attractive delocalization energies. The tribasis method is used to reconstruct the Hückel theory of π-electron delocalization in planar hydrocarbon molecules to account for the "overlap problem". In its empirically fitted form, the new theory can accurately resolve both π → π* transition energy and aromatic stabilization energy. The picture of covalent bonding emerging from both hydrogenic and Hückel calculations is that there is a presence of a Pauli repulsion of localization which is overcome by a roughly twice as strong delocalization stabilization to form the bond.

Rupturing aromaticity by periphery overcrowding

The balance between strain relief and aromatic stabilization dictates the form and function of non-planar π-aromatics. Overcrowded systems are known to undergo geometric deformations, but the energetically favourable π-electron delocalization of their aromatic ring(s) is typically preserved. In this study we incremented the strain energy of an aromatic system beyond its aromatic stabilization energy, causing it to rearrange and its aromaticity to be ruptured. We noted that increasing the steric bulk around the periphery of π-extended tropylium rings leads them to deviate from planarity to form contorted conformations in which aromatic stabilization and strain are close in energy. Under increasing strain, the aromatic π-electron delocalization of the system is broken, leading to the formation of a non-aromatic, bicyclic analogue referred to as ‘Dewar tropylium’. The aromatic and non-aromatic isomers have been found to exist in rapid equilibrium with one another. This investigation demarcates the extent of steric deformation tolerated by an aromatic carbocycle and thus provides direct experimental insights into the fundamental nature of aromaticity.

The curious case of the crystalline tri-thorium cluster: cyclic delocalization without aromatic stabilization?

Actinides have been known to form extremely weak homonuclear bonds with their d-type orbitals, and one should therefore expect the superposition of cyclic resonance forms containing such bonds to bring rather marginal aromatic stabilization to the system, if any. It is for this very reason that the discovery of the cyclically delocalized Th3 σ-bonding in the crystalline cluster isolated by Liddle and co-workers has sparked such vigorous discussion on the actual role of molecular aromaticity on the periphery of the periodic table. It has recently been argued that the tri-thorium ring at the heart of the cluster features considerable aromatic stabilization energy comparable to the heterocyclic π-aromatic rings such as thiophene. However, previous investigations involved highly ionized model clusters like Th3Cl64+ or Th310+ in which aromatic stabilization associated with the cyclic delocalization of charge is dramatically exaggerated. In this work we investigate the model tri-thorium clusters at different geometries and ionization states to show that cyclic delocalization of electrons in the isolated crystalline cluster may be associated with rather marginal σ-aromatic stabilization energy thus strongly suggesting its non-aromatic character.

Aromaticity and antiaromaticity in monoheterocyclic three‐membered rings: Application of natural bond orbital theory

AbstractUsing natural bond orbital theory, aromatic stabilization energies (NBO‐ASEs) were calculated for neutral and cationic monoheterocyclic three‐membered rings C2H2X(R)n (X = group 14–17 elements of rows 3–5; R = F, H or SiH3; and n = 0–2). Generally, there was a decrease in the (anti)aromatic nature of the molecule moving down the group as either the stabilizing orbital overlap (aromaticity) or destabilizing orbital interactions (antiaromaticity) were reduced due to larger atoms inducing longer bonds. (Anti)aromatic trends were supported with dissected nucleus independent chemical shift NICS(0)π,zz values (R2 = 0.94) and also by comparison of NBO‐ASEs to literature values for a selection of larger systems including benzene, cyclobutadiene, and the cyclopentadienyl anion. Comparisons are made to other studies that have used a range of approaches to study (anti)aromaticity in some of these molecules, and it was found that non‐dissected NICS values and homodesmotic equations often gave significantly overestimated if not misleading results.

Particle on a Ring Model for Teaching the Origin of the Aromatic Stabilization Energy and the Hückel and Baird Rules.

Simple mathematical models can serve to reveal the essence of experimental phenomena and scientific concepts. The particle in a box (PIB), for example, is widely used in undergraduate programs to teach the quantum mechanical principles behind the UV–vis spectra of conjugated polyenes and polyynes. In this work, the particle on a ring (POR) and the PIB models are used to elucidate the concept of aromaticity in Introductory Chemistry courses. Thus, we explain the origin of the aromatic stabilization energy, Hückel’s rule, and Baird’s rule. Besides applications, the limitations of the POR and PIB models are also discussed.

Investigation of intramolecular hydrogen bonding in naphthoquinone derivatives by quantum chemical calculations

AbstractIntramolecular hydrogen‐bonded naphthoquinone derivatives are investigated using both DFT (density functional theory) and MP2 (Møller–Plesset perturbation theory) method. Three different sets of naphthoquinone derivatives, each set consisting of syn and anti conformers with respect to the substitution by (–OH/–NH2) at peri and para position of the ring, are considered for investigation. In all cases, the syn conformer is found to be energetically more stable compared with the anti conformer. However, the hydrogen bond strength of the anti conformer is found to be more as obtained from NBO (natural bond orbital) analysis and QTAIM (quantum theory of atoms in molecules) theory. The O···H–O interactions are associated with covalent character while O···H–N interactions are noncovalent in nature. The NCI (noncovalent interaction) isosurface clearly indicates the formation of O···H–O and O···H–N intramolecular hydrogen bond. The energy partition analysis indicates electrostatic interaction as the dominant force in stabilizing these systems. The aromaticity calculation by HOMA (harmonic oscillator model of aromaticity) index and Bird index indicate the syn conformer with higher aromaticity compared to the anti conformer. The overall stability of the syn conformers is thus dominated by the aromatic stabilization energy instead of hydrogen bonding energy. The global reactivity descriptors further support the higher stability of the syn conformers according to maximum hardness principle (MHP) and minimum electrophilicity principle (MEP).

Interplay of Steric Effects and Aromaticity Reversals to Expand the Structural/Electronic Responses of Dihydrophenazines.

To gain insights into the coupling of conformational and electronic variables, we exploited steric hindrance to modulate a polycyclic skeleton with a bent conformation in the S0 state and a twisted conformation in the S1 state under the guidance of photoexcited aromaticity reversals. Polycyclic 5,10-dihydrophenazine (DHP) adopted a bent structure in S0 but involved a bent-to-planar transformation in S1 due to the excited-state aromaticity of the 8π-electron central ring. The N,N'-locations and 1,4,6,9-sites of the DHP skeleton provided a versatile chemical handle for fine-tuning intramolecular steric hindrance. Specifically, N,N'-diphenyl-5,10-dihydrophenazine (DPP-00) and its derivatives DPP-10-DPP-22 were synthesized with different numbers of methyl groups on the 1,4,6,9-sites. X-ray crystal analyses suggested that the DHP skeletons of DPP-00-DPP-22 had more bending configurations along the N···N axis with an increase in the number of methyl groups. Following the bending-promoted interruption of π-conjugation, the absorption spectra of DPP-00-DPP-22 significantly blue-shifted from 416 to 324 nm. By contrast, the emission bands exhibited a reverse shift to longer wavelengths from 459 to 584 nm as the number of methyl substituents increased. Theoretical calculations revealed that introducing methyl groups caused the planar DHP skeleton in S1 to further twist along the N···N axis, resulting in a twisted high-strain conformation. The greater Stokes shift of the more steric-hindered structure can be attributed to the release of larger strain and aromatic stabilization energy. This research highlighted the potential promise associated with the interplay of steric effects and aromaticity reversals in a single fluorophore.

Calculating stabilization energies for aromatic molecules in schools – quantum chemistry, a valuable tool for teaching and education

While aromaticity is one of the most fundamental concepts in chemistry and is taught in many elementary organic chemistry courses, it is hardly possible to activate pupils and students to learn and explore this topic independently. By applying computational chemistry and especially quantum chemistry, we introduce a possibility to activate students, allowing them to test taught concepts and do their own “experiments” cheaply and securely. We have focused in the current report on the structural and energetic aspects of aromaticity, as they are most compatible with other topics of typical chemistry lessons. To evaluate the aromatic stabilization energy, we suggest the simple-to-calculate and easy-to-understand isomerization stabilization energy (ISE) popularized by Paul von Ragué Schleyer and Frank Pühlhofer some 20 years ago. As practical examples, we demonstrate our concept with benzene- and pyridine-based systems and suggest topics for small projects for students and pupils.

Disclosure of Ground-State Zimmerman-Möbius Aromaticity in the Radical Anion of [6]Helicene and Evidence for 4π Periodic Aromatic Ring Currents in a Molecular "Metallic" Möbius Strip.

In 1966, Zimmerman proposed a type of Möbius aromaticity that involves through‐space electron delocalization; it has since been widely applied to explain reactivity in pericyclic reactions, but is considered to be limited to transition‐state structures. Although the easily accessible hexahelicene radical anion has been known for more than half a century, it was overlooked that it exhibits a ground‐state minimum and robust Zimmerman‐Möbius aromaticity in its central noose‐like opening, becoming, hence, the oldest existing Möbius aromatic system and with smallest Möbius cycle known. Despite its overall aromatic stabilization energy of 13.6 kcal mol−1 (at B3LYP/6‐311+G**), the radical also features a strong, globally induced paramagnetic ring current along its outer edge. Exclusive global paramagnetic currents can also be found in other fully delocalized radical anions of 4N+2 π‐electron aromatic polycyclic benzenoid hydrocarbons (PAH), thus questioning the established magnetic criterion of antiaromaticity. As an example of a PAH with nontrivial topology, we studied a novel Möbius[16]cyclacene that has a non‐orientable surface manifold and a stable closed‐shell singlet ground state at several density functional theory levels. Its metallic monoanion radical (0.0095 eV band gap at HSE06/6‐31G* level) is also wave‐function stable and displays an unusual 4π‐periodic, magnetically induced ring current (reminiscent of the transformation behaviour of spinors under spatial rotation), thus indicating the existence of a new, Hückel‐rule‐evading type of aromaticity.

From Macrocycles to Quantum Rings: Does Aromaticity Have a Size Limit?

ConspectusThe ring currents of aromatic and antiaromatic molecules are remarkable emergent phenomena. A ring current is a quantum-mechanical feature of the whole system, and its existence cannot be inferred from the properties of the individual components of the ring. Hückel's rule states that when an aromatic molecule with a circuit of [4n + 2] π electrons is placed in a magnetic field, the field induces a ring current that creates a magnetic field opposing the external field inside the ring. In contrast, antiaromatic rings with 4n π electrons exhibit ring currents in the opposite direction. This rule bears the name of Erich Hückel, and it grew from his molecular orbital theory, but modern formulations of Hückel's rule incorporate contributions from others, particularly William Doering and Ronald Breslow. It is often assumed that aromaticity is restricted to small molecular rings with up to about 22 π electrons. This Account outlines the discovery of global ring currents in large macrocycles with circuits of up to 162 π electrons. The largest aromatic rings yet investigated are cyclic porphyrin oligomers, which exhibit global ring currents after oxidation, reduction or optical excitation but not in the neutral ground state. The global aromaticity in these porphyrin nanorings leads to experimentally measurable aromatic stabilization energies in addition to magnetic effects that can be studied by NMR spectroscopy. Wheel-like templates can be bound inside these nanorings, providing excellent control over the molecular geometry and allowing the magnetic shielding to be probed inside the nanoring. The ring currents in these systems are well-reproduced by density functional theory (DFT), although the choice of DFT functional often turns out to be critical. Here we review recent contributions to this field and present a simple method for determining the ring current susceptibility (in nA/T) in any aromatic or antiaromatic ring from experimental NMR data by classical Biot-Savart calculations. We use this method to quantify the ring currents in a variety of aromatic rings. This survey confirms that Hückel's rule reliably predicts the direction of the ring current, and it reveals that the ring current susceptibility is surprisingly insensitive to the size of the ring. The investigation of aromaticity in even larger molecular rings is interesting because ring currents are also observed when mesoscopic metal rings are placed in a magnetic field at low temperatures. The striking similarity between the ring currents in molecules and mesoscopic metal rings arises because the effects have a common origin: a field-dependent phase shift in the electronic wave function. The main difference is that the magnetic flux through mesoscopic rings is much greater because of their larger areas, so their persistent currents are nonlinear and oscillatory with the applied field, whereas the flux through aromatic molecules is so small that their response is approximately linear in the applied field. We discuss how nonlinearity is expected to emerge in large molecular nanorings at high magnetic fields. The insights from this work are fundamentally important for understanding aromaticity and for bridging the gap between chemistry and mesoscopic physics, potentially leading to new functions in molecular electronics.

Experimental and Theoretical Evidence for Aromatic Stabilization Energy in Large Macrocycles.

Enhanced thermodynamic stability is a fundamental characteristic of aromatic molecules, yet most previous studies of aromatic stabilization energy (ASE) have been limited to small rings with up to 18 π-electrons. Here we demonstrate that ASE can be detected experimentally in π-conjugated porphyrin nanorings with Hückel circuits of 76-108 π-electrons. This conclusion is supported by analyzing redox potentials to calculate the energy change for isodesmic reactions that convert an aromatic ring to an antiaromatic ring or vice versa. It is also supported by analyzing the energy barriers to conformational equilibria that disrupt aromaticity in the transition state. Both types of experiment indicate that cationic porphyrin nanorings display ASEs of 1-5 kJ mol-1. Density functional theory calculations reproduce the results for both types of experiment and predict ASEs in the range of 1-16 kJ mol-1. The experimental ASEs in porphyrin nanorings are compared with an experimental ASE of [18]annulene of ∼11 kJ mol-1, deduced from analysis of the energy barriers to conformational equilibria in [16], [18], and [20]annulene. Calculated energies of isodesmic reactions give an ASE of ∼37 kJ mol-1 in [18]annulene. This work contributes to a fundamental understanding of aromaticity in large macrocycles.

Diels\u2013Alder reactions and electrophilic substitutions with atypical regioselectivity enable functionalization of terminal rings of anthracene

Reversing the regioselectivity of the renowned Diels–Alder reaction by overriding the usual thermodynamic and kinetic governing factors has always been a formidable challenge to synthetic organic chemists. Anthracenes are well-known to undergo [4 + 2]-cycloadditions with dienophiles at their 9,10-positions (central ring) over 1,4-positions (terminal ring) guided by the relative aromatic stabilization energy of the two possible products, and also by harboring the largest orbital coefficients of the highest occupied molecular orbital (HOMO) at the 9,10-positions. We, herein, report a 1,4-selective [4 + 2]-cycloaddition strategy of 9,10-unsubstituted anthracenes by installing electron-donating substituents on the terminal rings which is heretofore unprecedented to the best of our knowledge. The developed synthetic strategy does not require any premeditated engagement of the 9,10-positions either with any sterically bulky or electron-withdrawing substituents and allows delicate calibration of the regioselectivity by modulating the electron-donating strength of the substituents on the terminal rings. Likewise, the regioselective functionalization of the terminal anthracene ring in electrophilic substitution reactions is demonstrated. A mechanistic rationale is offered with the aid of detailed computational studies, and finally, synthetic applications are presented.

Donor-Acceptor Conjugated Macrocycles with Polyradical Character and Global Aromaticity.

SummaryPolyradical character and global aromaticity are fundamental concepts that govern the rational design of cyclic conjugated macromolecules for optoelectronic applications. Here, we report donor-acceptor (D−A) conjugated macromolecules with and without π-spacer derivatives to tune the antiferromagnetic couplings between the unpaired electrons. The macromolecules without π-spacer have a closed-shell electronic configuration and show global nonaromatic character in the singlet and lowest triplet states. However, the derivatives with π-spacer develop a nearly pure open-shell diradical and a very high polyradical character, not reported for D−A type macromolecules. Furthermore, the π-spacer derivatives display global nonaromaticity in the singlet ground state, but global aromaticity in the lowest triplet state, according to Baird's rule. The absorption spectra of the open-shell macromolecules calculated with time-dependent density functional theory indicate intensive light absorption in the near-infrared region and broadening to 2,500 nm, making these materials suitable for numerous optoelectronic applications.

Stability, Aromaticity, and Photophysical Behaviors of Macrocyclic Molecules: A Theoretical Analysis.

The macrocyclic molecules with terthiophene (TTH) isomers unit exhibit intriguing properties in terms of aromaticity, stability, and absorption. In this work, we theoretically designed a series of macrocyclic molecules featured with TTH and dithienothiophene (DTT) π-conjugated building units, which are used to permute pyrrole unit in porphyrin skeleton. Density functional theory and time-dependent DFT methods are used to evaluate the performance of the designed molecules. Our simulations show that molecules 1–3 exhibit excellent optoelectronic performance. Specifically, the molecule with the DTT unit is more stable than the one with TTH unit in terms of aromaticity and aromatic stabilization energy. This is because DTT unit enhances the coplanarity of the molecular material, facilitating electronic communication. Calculation of vertical electronic excitations suggests the absorption feature of these molecules is mainly contributed by the electronic excitations of higher occupied molecular orbital (HOMO) → lowest unoccupied molecular orbital (LUMO)+1 and HOMO-1 → LUMO. Judging from the key parameters determining the overall performance, 3 stands out because of its good planarity, large HOMO–LUMO gap, and strong aromaticity among all molecules. Interestingly, molecule 1 has the current density flow distributes around the outer section of TTH unit; in contrast, molecule 3 with DTT unit has the current density flow located at the inner section of DTT, which is beneficial for stability and aromaticity. Second-order perturbation energies are calculated to rationalize this observation. We expect that these research results can provide valuable insights into the rational design of novel molecular materials for a variety of applications.