Evaluation of the aromaticity of non-planar and bowl-shaped molecules by NICS criterion

On Aromaticity of the Aromatic α-Amino Acids and Tuning of the NICS Indices to Find the Aromaticity Order

The NICS aromaticity indices of the rings in flexible phenylalanine (Phe), tryptophan (Trp), tyrosine (Tyr), and histidine (His) chiral molecules were analyzed. These molecules have several dozens of conformers, and their rings are slightly non-planar. Therefore, the population-averaged NICS_pav index was defined, and the NICS scans had to be performed with respect to planes found by the least-squares routine. A rule differentiating an obverse and a reverse ring face in aromatic amino acids was formulated. The NICS scan minima corresponding to the obverse and reverse face were unequal, which prompted us to use the term ring face aromaticity/ring face tropicity. It appeared that for Phe, Trp, Tyr, and His, the reverse face has always had higher ring face aromaticity/ring face tropicity than the obverse one. Despite the NICS modifications, uncertainty about the amino acid aromaticity order remained. This motivated us to use the integral INICS index newly proposed by Stanger as well. Then, the following sequence was obtained: Trp(phenyl) > Phe > Trp(pyrrole) > His > Tyr. The juxtaposition of the INICS indices of amino acids with that of some model rings revealed a fair transferability of the values. Finally, analysis of the substituent effect on INICS demonstrated that the aromaticity of Tyr is the lowest due to the strength of the OH group π-electron-donating effect able to perturb enough the ring charge distribution and its magnetic aromaticity. The NICS calculations were executed using the ARONICS program written within the project.

Core-electron contributions to the molecular magnetic response

Orbital contributions to the magnetic response depend on the method used to compute them. Here, we show that dissecting nuclear magnetic shielding tensors using natural localized molecular orbitals (NLMOs) leads to anomalous core contributions. The arbitrariness of the assignment might significantly affect the interpretation of the magnetic response of nonplanar molecules such as C₆₀ or [14]helicene and the assessment of their aromatic character. We solve this problem by computing the core- and σ-components of the induced magnetic field (and NICS) and the magnetically induced current density by removing the valence electrons (RVE). We estimate the core contributions to the magnetic response by performing calculations on the corresponding highly charged molecules, such as C₆H₆³⁰⁺ for benzene, using gauge-including atomic orbitals and canonical molecular orbitals (CMOs). The orbital contributions to nuclear magnetic shielding tensors are usually estimated by employing a natural chemical shielding (NCS) analysis in NLMO or CMO bases. The RVE approach shows that the core contribution to the magnetic response is small and localized at the nuclei, contrary to what NCS calculations suggest. This may lead to a completely incorrect interpretation of the magnetic σ-orbital response of nonplanar structures, which may play a major role in the overall magnetic shielding of the system. The RVE approach is thus a simple and inexpensive way to determine the magnetic response of the core- and σ-electrons.

Theoretical study on the structures, electronic properties, and aromaticity of thia[4]circulenes

The structures, electronic properties, and aromaticity of a series of thia[4]circulenes are predicted using quantum chemical calculations.

Toward a Generalized Hückel Rule: The Electronic Structure of Carbon Nanocones

In this work, we investigate a particular class of carbon nanocones, which we name graphannulenes, and present a generalized Hückel rule (GHR) that predicts the character of their ground state based on simply the three topological indices that uniquely define them. Importantly, this rule applies to both flat and curved systems, encompassing a wide variety of known structures that do not satisfy the "classic" 4n + 2 rule such as coronene, corannulene, and Kekulene. We test this rule at the Hückel level of theory for a large number of systems, including structures that are convex and flat, with a saddle-like geometry, and at the CASSCF level of theory for a selected representative subset. All the performed calculations support the GHR that we propose in this work.

Analysis of Conformational, Structural, Magnetic, and Electronic Properties Related to Antioxidant Activity: Revisiting Flavan, Anthocyanidin, Flavanone, Flavonol, Isoflavone, Flavone, and Flavan-3-ol

Understanding the antioxidant activity of flavonoids is important to investigate their biological activities as well as to design novel molecules with low toxicity and high activity. Aromaticity is a chemical property found in cyclic structures that plays an important role in their stability and reactivity, and its investigation can help us to understand the antioxidant activity of some heterocyclic compounds. In the present study, we applied the density functional theory (DFT) to investigate the properties of seven flavonoid structures with well-reported antioxidant activity: flavan, anthocyanidin, flavanone, flavonol, isoflavone, flavone, and flavan-3-ol. Conformational, structural, magnetic, and electronic analyses were performed using nuclear magnetic resonance, ionization potentials, electron affinity, bond dissociation energy, proton affinity, frontier molecular orbitals (highest occupied molecular orbital (HOMO)/lowest unoccupied molecular orbital (LUMO)), and aromaticity through nucleus-independent chemical shifts to analyze these seven flavonoid structures. We revised the influence of hydroxyl groups on the properties of flavonoids and also investigated the influence of the aromaticity of these seven flavonoids on the antioxidant activity.

Fast prediction of oxygen reduction reaction activity on carbon nanotubes with a localized geometric descriptor

By using the pyramidalization angle as a localized geometric descriptor for oxygen reduction reaction (ORR) activity of carbon nanotubes (CNTs), we show the ORR activity of these systems can be readily predicted with mere structural optimization of CNTs.

Local aromaticity mapping in the vicinity of planar and nonplanar molecules

We report on nucleus-independent magnetic shielding (NICS) scans over the centers of six- and five-membered rings in selected metal phthalocyanines (MPc) and fullerene C60 for more accurate characterization of local aromaticity in these compounds. Detailed tests were conducted on model aromatic molecules including benzene, pyrrole, indole, isoindole, and carbazole and subsequently applied to H₂ Pc, ZnPc, Al(OH)Pc, and CuPc. Similar behavior of three selected magnetic probes, Bq, ³ He, and ⁷ Li⁺ , approaching perpendicularly the ring centers, was observed. For better visualization of shielding zone over the centers of aromatic rings, we introduced a simple mathematical procedure: the first and second derivatives of scan curves with respect to magnetic probe position enabled their additional examination. It allowed an easier localization of curve minimum and discrimination between areas in space varying by the magnetic field magnitude and to illustrate local aromaticity of two different kinds of rings in MPc with better resolution. Our results supported earlier reports on very low aromaticity indexes of pyrrole ring incorporated into MPc and significant aromaticity of the central macrocycle. No direct dependence between harmonic oscillator model of aromaticity and NICS was observed. Instead, a correlation between position of scan curve minimum and its magnitude were observed. In addition, the NICS values and ³ He chemical shifts in the middle of neutral C60 and C60^6- anion agreed well with the reported experimental NMR values for He@C60 and He@C60^6- .

Curved π-conjugated corannulene dimer diradicaloids

So far, most reported open-shell singlet diradicaloids are based on planar π-conjugated molecules. Herein, we report the bridged corannulene dimer diradicaloids, Cor-D1 and Cor-D2, both showing a three-dimensional curved π-conjugated structure. Cor-D1 has a small diradical character (y0 = 5.4%) and behaves more like a closed-shell quinoidal compound at room temperature, while Cor-D2 is a typical open-shell diradicaloid with a larger diradical character (y0 = 16.9%). Both compounds exhibited paramagnetic activity at elevated temperatures, with a singlet-triplet energy gap (ΔES-T) of -8.4 and -3.0 kcal mol-1, respectively. X-ray crystallographic analysis revealed that both molecules have a dumbbell-shaped geometry, with the two terminal corannulene bowls bent to opposite directions. The spin is largely delocalized onto the two bowls in Cor-D2 and there are multiple [CH···π] interactions between the neighboring bowls. Chemical oxidation/reduction to their respective dications/dianions results in global aromaticity with [4n + 2] π-electrons delocalized through the periphery of the whole framework.

Comparison of effects of charge delocalization and π-electron delocalization on the stability of monocyclic compounds

Aromaticities and stabilities of ortho, meta, and para isomers of some derivatives of benzene, C₅H₅^- and C₇H₇⁺ were studied and compared basis on the NICS index and their relative energies. For the benzene and C₇H₇⁺ derivatives, the ortho isomers with less stability were more aromatic. This discrepancy was also observed for the molecules with conjugated and non-conjugated π-electrons. However, for the charged conjugated systems, the structures with delocalized charge were more stable. Effect of electron withdrawing (EWGs) and electron donating groups (EDGs) on the electron delocalization and stability of the neutral and charged molecules was investigated. It was observed that the EWDs and EDGs change the stability trend of the neutral systems, while in the case of charged molecules, the isomers with delocalized charge were more stable regardless of the type of substituents. Although both π-electron delocalization and charge delocalization stabilize the aromatic and conjugated systems, effect of charge delocalization on the stability is more than that of π-electron delocalization.

Influence of nitrogen doping in sumanene framework toward hydrogen storage: A computational study

Two conditions are important to obtain appropriate substances for hydrogen storage; high surface area and fitting binding energy (BE). Doping is a key strategy that improves BE. We investigated hydrogen adsorption onto twenty six nitrogen disubstituted isomers of sumanene (C₁₉N₂H₁₂) by MP2/6-311++G(d,p)//B3LYP/6-31+G(d) and M06-2X/6-31+G(d) levels of theory. Effect of nitrogen doping in different positions of sumanene was checked. To obtain better BE, basis set superposition error (BSSE) and zero point energy (ZPE) corrections were used. Anticipating of adsorption sites and extra details about adsorption process was done by molecular electrostatic potential (MEP) surfaces. Various types of density of state (DOS) diagrams such as total DOS (TDOS), projected DOS (PDOS) and overlap population DOS (OPDOS) and natural bond orbital (NBO) analysis were used to find better insight on the adsorption properties. In addition of temperature depending of the BE, HOMO-LUMO gap (HLG), dipole moment, reactivity and stability, bowl depth and natural population analysis (NPA) of the isomers were studied. A physisorption mechanism for adsorption was proposed and a trivial change was seen. Place of nitrogen atoms in sumanene frame causes to binding energy increases or decreases compared with pristine sumanene. The best and the worst isomers and category of isomers were suggested.

On splitting of the NICS(1) magnetic aromaticity index

The NICS(1) magnetic aromaticity index is split into NICS(1) and NICS(−1) indices when the points 1 Å above and below the ring center are inequivalent by symmetry. The two indices characterize the aromaticity of the two ring faces rather than the ring itself.