Which NICS method is most consistent with ring current analysis? Assessment in simple monocycles

Microhydration of small protonated polyaromatic hydrocarbons: a first principles study

Using first principles methodology, we investigate the microsolvation of protonated benzene (BzH+), protonated coronene (CorH+) and protonated dodecabenzocoronene (DbcH+). Gas phase complexes of these small protonated polyaromatic hydrocarbons (H+PAHs) with mono-, di-, and tri-hydrated water molecules are considered. Their most stable forms are presented, where we discuss their structural, energetic aromaticity and IR and UV spectral features. In particular, we focus on the analysis of the bonding and various non-bonded interactions between these protonated aromatics and water clusters. The strength of non-bonded interactions is quantified and correlated with their electron density profiles. Furthermore, insights into the interfacial interactions and stability of these complexes were obtained through non-covalent index and symmetry-adapted perturbation theory (SAPT0) analyses. We also discuss the effects of the extension of the π aromatic cloud on the water solvation of these protonated aromatics. In particular, we extended our predictions for the S0 → S1 and S0 → T1 wavelength transitions of micro hydrated H+PAHs to deduce those of these species solvated in aqueous solution. The present findings should be useful for understanding, at the microscopic level, the effects of water interacting with H+PAHs, which are relevant for organic chemistry, astrochemistry, atmospheric chemistry, combustion and materials science.

D 4h H©K4H4-: a planar tetracoordinate hydrogen global minimum

Herein we report a square-like D4h H©K4H4- anion with one planar tetracoordinate hydrogen (ptH) center, which is the global minimum (GM) structure and possesses good dynamic stability. The planar structure of the system is preserved by four peripheral K-H-K three-center two-electron (3c-2e) σ bonds together with one 5c-2e σ bond over the HK4 core. The multicenter ionic bonds dominate the stability of ptH, while the contribution of qualitative σ aromaticity is extremely limited.

Aromaticity of six-membered nitro energetic compounds through molecular electrostatic potential, magnetic, electronic delocalization and reactivity-based indices

The electron density depletion associated with π-hole at the ring center typical of energetic compounds was clearly revealed by the molecular electrostatic potential (ESP). In addition, the spatial arrangement of NO2 groups appears to affect the ESP value in the ring center, and therefore sensitivity to detonation. Indeed, for monocyclic nitrobenzene compounds with the same number of NO2 groups, the ESP value in the ring center decreases as the NO2 groups are more closely spaced. As expected, the central rings become less aromatic as NO2 groups are added. The MCI, PDI, PLR, NICSzz(1), FLU indices are all strongly correlated with the ESP values observed in the ring center of the set of nitrobenzenes. Aromaticity indices based on electron delocalization criteria appear to be very sensitive to small variations in aromaticity. Among magnetic-based indices, only NICSzz(1) is capable to predict small changes in aromaticity. The PLR index derived from conceptual DFT is quite relevant for predicting small variations in aromaticity. According to our results, the most suitable aromaticity index is not based on a single criterion, and that selecting it is more subtle. Therefore, it is important to combine information from several criteria to obtain a more complete description of the aromaticity of the studied compounds.

Is azulene's local aromaticity and relative stability driven by the Glidewell-Lloyd rule?

The local aromaticity of azulene has been studied to understand their electronic properties. For this purpose, we have used the magnetic criterion through magnetically induced current density maps, ring current strengths, NICSzz(1), and the bifurcation value of three-dimensional surfaces of NICSzz. On the other hand, the delocalization criterion was used by calculating the MCI and ELFπ. The results show that the five-membered ring (5-MR) is more aromatic than the seven-membered ring (7-MR) and more aromatic than the free C5H5- ring. The opposite case is seen for the seven-membered ring, which is less aromatic than the free C7H7+. The local aromatic rings in azulene are formed due to an intramolecular electron transfer from the 7-MR to the 5-MR. In addition, the proposed resonance structures that allow explaining the properties of azulene, such as the dipole moment or the relative stability (in comparison to other isomers), show a preference for the formation of 5-MRs; for this reason, it is possible to conclude that the aromaticity and relative stability of azulene is driven by the Glidewell-Lloyd rule.

Understanding Aromaticity in [5]Helicene-Bridged Cyclophanes: A Comprehensive Study

This study explores the aromaticity of doubly [5]helicene-bridged (1,4)cyclophane and triply [5]helicene-bridged (1,3,5)cyclophane via calculations of the magnetic response and of electronic aromaticity indices. The primary objective is to assess the π-electron delocalization to determine whether they sustain global ring currents associated with π aromaticity. The molecules show local ring currents in the presence of an external magnetic field. The ring currents flow diatropically in the stacked six-membered rings and in the helicene arms. However, these π currents are not interconnected due to the discontinuity of the π delocalization at the C-C single bonds connecting the central six-membered rings to the helicene arms. Electronic indices suggest that the helicene-arm systems have significantly smaller electron delocalization than benzene. The reduction in the delocalization does not compromise their ability to exhibit ring currents in the presence of an external magnetic field. The analysis provides further evidence that the magnetic criteria yield a different degree of aromaticity for the helicene arms than obtained in the calculation of the electronic aromaticity indices. However, both approaches confirm that the studied molecules are not globally aromatic.

Bond-alternated and bond-equalized hexazine derivatives

Planar cyclic conjugated molecules satisfying the Hückel rule are generally classified as aromatics and π delocalization is the key feature of aromatic compounds. It was shown that the π system of benzene prefers bond-alternation and that the delocalization observed is a consequence of bond-equalization, which is a σ effect. For systems wherein the π bonds are strong, such as those between N atoms, the π-distortivity may outweigh the σ preference for bond-equalization. Thus, one anticipates a bond-alternated structure for N6; however, neither the bond-equalized structure (D6h) nor the bond-alternated one (D3h) is a minimum as evident from the previous literature; calculations have shown that cyclic N6 is non-planar. Herein we show that a Lewis acid coordination strategy can be employed to stabilize the planar structure of N6. We show that the structure can be bond-alternated or bond-equalized depending on the strength of the Lewis acid. Kinetic stability with respect to concerted decomposition to three N-N triple bonded systems was assessed, and a few bond-equalized N6 systems were found to be potential candidates for synthesis.

Unexpected structural preference with metallophilic Ag-Au contacts in silver(I)-N-heterocyclic carbene cluster; experimental and theoretical approach

A novel synthetic donor-atom-selective approach has been adopted for the synthesis of a heterobimetallic cluster of a new NCN-pincer, 1,3-bis-(1-methyl-1H-benzo[d]imidazol-2-yl-methyl)-1H-imidazol-3-ium hexafluorophosphate (1·HPF6). The complex [Ag3(1)3][PF6]3 (2) has been prepared via the Ag2O route; which undergoes transmetallation to yield a cluster that seems to be the first example of the heterobimetallic trinuclear system [Au-Ag2(1)2Cl][PF6]2, 3. Finally, the trinuclear cluster geometries of 2 and 3 were confirmed via SCXRD studies. Interestingly, Au(I) binds preferentially with soft donor Ccarbene, which transmetallated from the cluster of 2. In both the cyclic trinuclear clusters, the M-M interactions were further inspected using gauge independent atomic orbital (GIAO) computations. Both 2 and 3 are luminescent and possess σ-aromaticity; the NICS values indicate that 3 is more aromatic than 2.

On the aromaticity and stability of benzynes in the ground and lowest-lying triplet excited states

In this work, we have revisited the aromaticity of benzyne isomers at the unrestricted density functional theory level (UDFT) using the energetic, magnetic, and delocalization criteria. In addition, this last criterion has also been analyzed employing complete active space (CASSCF) calculations. The results show conservation of aromaticity in these monocycles. Additionally it is observed that this trend is maintained in polycyclic aromatic hydrocarbon derivatives such as biradical didehydrophenanthrenes. Do these results imply a violation of Baird's rule? The answer is No, because this conservation in aromaticity is due to the loss of hydrogen atoms affects only the electronic σ skeleton and exerts a minor influence on the π cloud. Additionally, we have analyzed the relative stability of benzyne isomers and their relationship with experimental ΔES-T values. According to the literature, the stability of the benzynes in the singlet state is due to an effective interaction between the electrons of the biradical centers; however, this effect is completely reversed in the triplet state, which explains why the para isomer has the lowest ΔES-T gap.

Deciphering the origin of the first steps in the degradation of azo dyes: a computational study

Azo dyes find applications across various sectors including food, pharmaceuticals, cosmetics, printing, and textiles. The contaminating effects of dyes on aquatic environments arise from toxic effects caused by their long-term presence in the environment, buildup in sediments, particularly in aquatic species, degradation of pollutants into mutagenic or mutagenic compounds, and low aerobic biodegradability. Therefore, we theoretically propose the first steps of the degradation of azo dyes based on the interaction of hydroperoxyl radical (•OOH) with the dye. This interaction is studied by the OC and ON mechanisms in three azo dyes: azobenzene (AB), disperse orange 3 (DO3), and disperse red 1 (DR1). Rate constants calculated at several temperatures show a preference for the OC mechanism in all the dyes with lower activation energies than the ON mechanism. The optical properties were calculated and because the dye-•OOH systems are open shell, to verify the validity of the results, a study of the spin contamination of the ground [Formula: see text] and excited states [Formula: see text] was previously performed. Most of the excited states calculated are acceptable as doublet states. The absorption spectra of the dye-•OOH systems show a decrease in the intensity of the bands compared to the isolated dyes and the appearance of a new band of the type π → π* at a longer wavelength in the visible region, achieving up to 868 nm. This demonstrates that the reaction with the •OOH radical could be a good alternative for the degradation of the azo dyes.

Unlocking Heteroaromatic Ring Systems through Chalcogen Insertion into Boroles

In this work, we report the reactivity of various annulated borole derivatives toward chalcogen (O, S, and Se) insertion. Among a series of 9-borafluorenes with different boron substituents (Ph, Br, or o-carboranyl) and a mixed thiophene-benzene-fused derivative, only the 9-o-carboranyl-substituted 9-borafluorene yielded the complete set of chalcogen-containing heteroarenes, including the first 1,2-selenaborinine derivative. To evaluate the aromaticity of this heterocyclic analogue of phenanthrene, nucleus-independent chemical shift (NICS) values were computed and compared to those of its lighter group 16 congeners.

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.

Exploration of the potential energy surface in mixed Zintl clusters applying an automatic Johnson polyhedra generator: the case of arachno E6M24- (E = Si, Ge, Sn; M = Sb, Bi)

A new algorithm called Automatic Johnson Cluster Generator (AJCG) is presented, which, as its name indicates, allows the definition of the desired Johnson polyhedron to subsequently carry out all the possible permutations between the atoms that form this polyhedron. This new algorithm allows the exhaustive study of the structures' potential energy surface (PES). In addition, the AJCG algorithm is helpful for the study of three-dimensional compounds such as boranes or Zintl clusters and their structural derivatives with two or more different atoms. The automatic filling of vertices is particularly useful in mixed compounds because of the possibility of taking into account all possible configurations in the structure. As a test system, we investigated the arachno-type E6M24- (E = Si, Ge, Sn; M = Sb, Bi) structure which has eight vertices and complies with Wade-Mingos rules. Initially, we defined a bipyramidal structure (10 vertices), and filled the vertices with the atoms in all possible configurations. Since the selected system has eight atoms, the two remaining vertices were filled with pseudo atoms to complete the structure. After re-optimizing the initial population generated with AJCG, a large number of isomers with energy below 10 kcal mol-1 are identified. These results show that the most stable isomers possess homonuclear M-M bonds, except Sn6Bi24-. Although the overall putative minima differ at the PBE0-D3 and DLPNO-CCSD(T) levels, they are always competitive minima. In addition to using high-precision methodologies to correctly study relative energies, applying solvent effects in highly charged systems becomes mandatory. The aromatic character of these studied systems was demonstrated qualitatively with two- and three-dimensional mapping and quantitatively by calculating the value of the z-component of the induced magnetic field at the cage center, including scalar and spin-orbit correction for relativistic effects. The compounds studied have a high degree of aromaticity, which allows us to establish that despite structural modifications (i.e., from closo to arachno), the aromaticity is preserved.

Tubular Magnetic Shielding Scan (TMSS): A New Technique for Molecular Space Exploration. (i) The Case of Aromaticity of Benzene and [n]Paracyclophanes

Both traditional and novel techniques were employed in this work for magnetic shielding evaluation to shed new light on the magnetic and aromaticity properties of benzene and 12 [n]paracyclophanes with n = 3-14. Density functional theory (DFT) with the B3LYP functional and all-electron Jorge-ATZP and x2c-TZVPPall-s basis sets was utilized for geometry optimization and magnetic shielding calculations, respectively. Additionally, the 6-311+G(d,p) basis set was incorporated for the purpose of comparing the magnetic shielding results. In addition to traditional evaluations such as NICS/NICSzz-Scan, and 2D-3D σiso(r)/σzz(r) maps, two new techniques were implemented: bendable grids (BGs) and cylindrical grids (CGs) of ghost atoms (Bqs). BGs allow for the recording of magnetic shielding from the bent ring levels of [n]pCPs, while CGs provide tubular magnetic shielding scan (TMSS) maps detailing the magnetic shielding from a cylindrical region above and below the ring frame. Our findings suggest that smaller [n]pCPs with n < 6 exhibit deviations in the magnetic shielding above and below the ring, indicating a broken electron delocalization under the ring. In contrast, larger [n]pCPs tend to behave similarly to benzene in terms of magnetic shielding. Moreover, we found that shorter polymethylene chains of [n]pCPs exhibit significantly higher magnetic shielding interactions with the ring. Both of the above techniques offer new and promising tools for characterizing nonplanar aromatic compounds, thereby contributing to a deeper understanding of their magnetic and electronic properties.

Quest for the Most Aromatic Pathway in Charged Expanded Porphyrins

Despite the central role of aromaticity in the chemistry of expanded porphyrins, the evaluation of aromaticity remains difficult for these extended macrocycles. The presence of multiple conjugation pathways and different planar and nonplanar π-conjugation topologies makes the quantification of global and local aromaticity even more challenging. In neutral expanded porphyrins, the predominance of the aromatic conjugation pathway passing through the imine-type nitrogens and circumventing the amino NH groups is established. However, for charged macrocycles, the question about the main conjugation circuit remains open. Accordingly, different conjugation pathways in a set of neutral, anionic, and cationic expanded porphyrins were investigated by means of several aromaticity indices rooted in the structural, magnetic, and electronic criteria. Overall, our results reveal the predominance of the conjugation pathway that passes through all nitrogen atoms to describe the aromaticity of deprotonated expanded porphyrins, while the outer pathway through the perimeter carbon atoms becomes the most aromatic in protonated macrocycles. In nonplanar and charged macrocycles, a discrepancy between electronic and magnetic descriptors is observed. Nevertheless, our work demonstrates AV_min remains the best tool to determine the main conjugation pathway of expanded porphyrins.

Electronic Transmutation Concept: Is the Inverse Process Possible? An Evaluation of Main Group Compounds

The electronic transmutation (ET) concept states that when an element with atomic number Z gains an electron, it transmutes into a Z + 1 element, leading to species that possess similar chemical bonding patterns and geometric structures regarding the original (Z + 1) element. In this work, the opposite concept, that is, the inverse ET, is assessed. For this purpose, several main group compounds have been analyzed in terms of the adaptive natural density partitioning. The obtained results suggest that when an atom Z loses an electron, it transmutes into a Z - 1 atom, acquiring its geometrical structure and bonding pattern.

Theoretical Investigation of Geometries and Bonding of Indium Hydrides in the In2Hx and In3Hy (x = 0-4,6; y = 0-5) Series

Boron hydrides have been an object of intensive theoretical and experimental investigation for many decades due to their unusual and somewhat unique bonding patterns. Despite boron being a neighboring element to carbon, boron hydrides almost always form non-classical structures with multi-center bonds. However, we expect indium to form its interesting molecules with non-classical patterns, though such molecules still need to be extensively studied theoretically. In this work, we investigated indium hydrides of In₂H_x (x = 0-4,6) and In₃H_y (y = 0-5) series via DFT and ab initio quantum chemistry methods, performing a global minimum search, chemical bonding analysis, and studies of their thermodynamical stability. We found that the bonding pattern of indium hydrides differs from the classical structures composed of 1c-2e lone pairs and 2c-2e bonds and the bonding pattern of earlier investigated boron hydrides of the B_nH_n+2 series. The studied stoichiometries are characterized by multi-center bonds, aromaticity, and the tendency for indium to preserve the 1c-2e lone pair.

Resonance-assisted intramolecular triel bonds

The possibility that the intramolecular Tr⋯S triel bond is strengthened by resonance is examined by quantum chemical calculations within the planar five-membered ring of TrH₂-CRCR-CRS (Tr = Al, Ga, In; R = NO₂, CH₃). This internal bond is found to be rather short (2.4-2.7 Å) with a large bond energy between 12 and 21 kcal mol^-1. The pattern of bond length alternation and atomic charges within the ring is consistent with resonance involving the conjugated double bonds. This resonance enhances the triel bond strength by some 25%. The electron-withdrawing NO₂ group weakens the bond, but it is strengthened by the electron-donating CH₃ substituent. NICS analysis suggests the presence of a certain degree of aromaticity within the ring.

Aromaticity Concepts Derived from Experiments

Aromaticity, a very important term in organic chemistry, has never been defined unambiguously. Various ways to describe it come from different phenomena that have been experimentally observed. The most important examples related to some theoretical concepts are presented here.

Miscellaneous Passerini Reaction for α-Acyloxy Carboxamide: Synthesis and Process Optimization Study

:

An accelerating effect of “Phase Transfer Catalyst” as additive was exposed for the Passerini three-component reactions and the influence on the reaction rate was studied concerning direct involvement of reactant molecules. The most flexible reaction for the rapid formation of diverse “α-acyloxycarboxamides” using passerini reaction involved multicomponent reactions using miscellaneous 2-(prop-2-ynyloxy)benzaldehyde with various aromatic acid and slightly non-polar fragment i.e. 2-isocyano-2,3,3-trimethylbutane and the representative molecule was characterized with resepct to DEPT135 NMR technique.

The hidden aromaticity in borazine

The aromaticity of borazine (B3N3H6), also known as the inorganic benzene, is a controversial issue since this compound has several characteristics that could qualify it as an aromatic compound. However, recent studies using magnetic criteria indicate that this compound should be considered as a non-aromatic system. This assignment is mainly due to diatropic currents in the nitrogen atoms without observation of ring currents. The present work shows by means of the magnetic criteria that borazine has a ring current hidden by the local contributions of degenerate orbitals π1 and π2. Additionally, the study of borazine's first triplet state antiaromaticity using the magnetic and energetic criteria by means of isomerization stabilization energies (ISEs) together with Baird's and Hückel's rules suggests that borazine is best described as an (weakly) aromatic system.

Simple and efficient visualization of aromaticity: bond currents calculated from NICS values

Aromaticity is a fundamental concept in chemistry, underpinning the properties and reactivity of many organic compounds and materials. The ability to easily and accurately discern aromatic behavior is key to leveraging it as a design element, yet most aromaticity metrics struggle to combine accurate quantitative evaluation, intuitive interpretability, and user-friendliness. We introduce a new method, NICS2BC, which uses simple and inexpensive NICS calculations to generate information-rich and easily-interpreted bond-current graphs. We test the quantitative and qualitative characterizations afforded by NICS2BC for a selection of molecules of varying structural and electronic complexity, to demonstrate its accuracy and ease of analysis. Moreover, we show that NICS2BC successfully identifies ring-current patterns in molecules known to be difficult cases to interpret with NICS and enables deeper understanding of local aromaticity trends, demonstrating that our method adds additional insight.

Alkene insertion reactivity of a o-carboranyl-substituted 9-borafluorene

The synthesis of 9-borafluorene with an electron-withdrawing o-carboranyl substituent and its reactions with a series of alkenes are described. The o-carboranyl substituent is bonded via one of the cluster carbon atoms to the boron atom of the 9-borafluorene moiety. In all cases, the reactions afford partly saturated analogs of borepins (i.e. 6,7-dihydroborepins) by unprecedented alkene insertion into the endocyclic B–C bond of the borole ring. Comparative studies with 9-bromo-9-borafluorene illustrate the superior insertion reactivity of the carboranyl-substituted derivative. A suite of experimental and computational techniques disclose the unique properties of the 9-borafluorene and provide insight into how the 9-carboranyl substituent affects its chemical reactivity.

A 9-carboranyl-substituted 9-borafluorene is reported, which is capable of undergoing efficient ring expansion to 6,7-dihydroborepins by a previously unknown alkene insertion.

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.

Analysis of Local and Global Aromaticity in Si3C5 and Si4C8 Clusters. Aromatic Species Containing Planar Tetracoordinate Carbon

The minimum energy structures of the Si3C5 and Si4C8 clusters are planar and contain planar tetracoordinate carbons (ptCs). These species have been classified, qualitatively, as global (π) and local (σ) aromatics according to the adaptive natural density partitioning (AdNDP) method, which is an orbital localization method. This work evaluates these species’ aromaticity, focusing on confirming and quantifying their global and local aromatic character. For this purpose, we use an orbital localization method based on the partitioning of the molecular space according to the topology of the electronic localization function (LOC-ELF). In addition, the magnetically induced current density is analyzed. The LOC-ELF-based analysis coincides with the AdNDP study (double aromaticity, global, and local). Moreover, the current density analysis detects global and local ring currents. The strength of the global and local current circuit is significant, involving 4n + 2 π- and σ-electrons, respectively. The latter implicates the Si-ptC-Si fragment, which would be related to the 3c-2e σ-bond detected by the orbital localization methods in this fragment.

Unraveling the stability of cyclobutadiene complexes using aromaticity markers

Cyclobutadiene (CBD) is the paradigmatic antiaromatic molecule but is known to form highly stable aromatic complexes, e.g. CBD-Fe(CO)3. This intriguing reversal of aromaticity from antiaromatic to aromatic terrain during the complexation process cannot be appropriately handled with single-reference-based theoretical techniques. We explore this aromaticity reversal, for the first time, by a detailed aromaticity analysis using magnetically induced current densities (MICD) and nucleus independent chemical shifts (NICS) using genuine ab initio multi-reference wavefunction-based theory. We trace the dramatic change of aromaticity for a prototypical cyclobutadiene complex, CBD-CH+ (CH+Fe(CO)3), considering a 3D potential energy surface for two independent parameters, namely the approach of CH+ and the automerization cross-section of cyclobutadiene. The 3D potential energy surfaces indicate the presence of a conical intersection/avoided crossing between the ground and the first excited state. The plot of aromaticity indices and the corresponding numerical values show that the change of aromaticity indices is drastic around the conical intersection/avoided crossing and automerization of cyclobutadiene plays a crucial role in the formation of cyclobutadiene complexes. Computations on analogous CBD-Be and CBD-CO systems (Be/COFe(CO)3) emphasize the generality of the conclusions drawn from the CBD-CH+ system.

Acenes and phenacenes in their lowest-lying triplet states. Does kinked remain more stable than straight?

The larger stability of phenacenes compared to their acene isomers in their ground states is attributed to the larger aromaticity of the former. To our knowledge the relative stability of acenes and phenacenes in their lowest-lying triplet states (T1) has not been discussed yet. Using unrestricted density functional theory calculations, our results show that for the smallest members of the series, acenes in their T1 states are more stable than the corresponding phenacenes. However, when the number of the rings (n) involved increases, the energy difference is reduced and for n > 12, phenacenes become more stable than acenes in their T1 states. To rationalize this trend, we analyze the aromaticity of acenes and phenacenes using a set of aromaticity descriptors. We find that in the T1 states of both acenes and phenacenes, the outer rings form aromatic Clar π-sextets. In acenes, delocalization of spin density in the central rings leads to the preferred formation of the largest antiaromatic diradical. Resonant structures in the form of antiaromatic diradical Baird π-octadectets and π-tetradectets are the major contributors, while the smaller ones, such as π-doublets and π-sextets, contribute the least. In phenacenes, structures with diradical antiaromatic Baird π-sextets in some of the central rings contribute the most. These results are relevant to understand the (anti)aromaticity of larger polycyclic aromatic hydrocarbons in their triplet states.

Assessment of New Expanded Porpholactones as UV/Vis/NIR Chromophores for Dye-Sensitized Solar Cell Applications

Expanded porphyrins arise as an alternative for potential application as chromophores in dye-sensitized solar cells. (DSSCs). The modification of the core of these compounds provides remarkable changes in the photoelectronic behavior. In the present article, the improvement of its properties for a potential application as UV/vis/NIR chromophores in DSSCs has been studied, when an oxazolone moiety has replaced an imine ring in analogy to the porpholactones first synthesized by Crossley et al. ( J. Chem. Soc., Chem. Commun. 1984, 920-922). These expanded porpholactones present a noticeable red shift as well as an increase in the intensity of the Q-bands regarding the parent compounds. The photophysical properties of Sapphyrin have been explored through DFT calculations and vibrationally resolved absorption spectra simulations. Energetic parameters showed favorable electron injection from the chromophore to the TiO₂ semiconductor. In addition, aromaticity was analyzed and rationalized using magnetic and delocalization criteria. Results showed qualitatively similar trends between aromaticity descriptors and Q bands giving a great opportunity to the use this property in the rational design of chromophores. Finally, the nonequilibrium Green's function formalism shows the ability of expanded porpholactones in electron transport.

Relating nucleus independent chemical shifts with integrated current density strengths

Indices based on the nucleus independent chemical shift (NICS) are the most frequently used in analysis of magnetic aromaticity. The magnetically induced current density, on the other hand, is a key concept in defining magnetic aromaticity. The integrated current strength (current strength susceptibility) was found to be a very useful tool in aromaticity studies. There is widely accepted notion that the properly chosen NICS-based index can provide information on the current density strength and direction in a molecule of interest. In this work, a detailed numerical testing of the relationship between the integrated bond current strength and the most employed NICS indices was performed for a set of 43 monocyclic aromatic molecules. Based on the statistical data analysis, the relationship between the bond current strength and its π and σ electron components, on one side, and the isotropic NICS (NICSiso and NICSπ,iso) and zz-component of the NICS tensor (NICSzz and NICSπ,zz), on the other side, was examined. It was found that between the NICSπ,zz(1) and π-electron bond current strenghts there is very good linear correlation. Quite surprisingly, it was revealed that the NICSiso(1) and NICSzz(1) are not correlated with the π electron bond current strengths. On the other hand, a reasonably good linear correlation was found between the NICSzz(1) and total bond current strengths.

On the NICS limitations to predict local and global current pathways in polycyclic systems

Here, we analyze the possibility of predicting local and global current densities in a series of bicyclic hydrocarbons with 4n and 4n + 2 π-electrons from the nucleus-independent chemical shifts (NICS) computations.

Aromatic nature of neutral and dianionic 1,4-diaza-2,3,5,6-tetraborinine derivatives

The aromatically relevant parameters of boron-rich inorganic benzenes-neutral and dianionic 1,4-diaza-2,3,5,6-tetraborinine derivatives (B4N2R6)-have been computationally estimated and evaluated from geometric, electronic, magnetic, and energetic points of view. The majority of the criteria (ASE, NICSzz, ELF, and PDI) indicate that the aromaticity of the neutral B4N2 benzene analogue stabilized by Lewis bases lies in between those of benzene and borazine. On the other hand, the aromaticity of the dianionic B4N2 benzene analogue 4' is controversial. The pronounced aromatic nature of 4' is supported by ELFπ, PDI, and NICSπzz, but ASE, the FiPC-NICS plot, and ACID oppose this. These data confirm that even with the same B4N2-skeletal framework of a 6π-system, the aromatic feature varies depending on the overall charge of the B4N2 systems.

Searching for double σ- and π-aromaticity in borazine derivatives

Inspired by the double-aromatic (σ and π) C6H3 +, C6I6 2+, and C6(SePh)6 2+ ring-shaped compounds, herein we theoretically study their borazine derivative analogues. The systems studied are the cation and dications with formulas B3N3H3 +, B3N3Br6 2+, B3N3I6 2+, B3N3(SeH)6 2+, and B3N3(TeH)6 2+. Our DFT calculations indicate that the ring-shaped planar structures of B3N3H3 +, B3N3I6 2+, and B3N3(TeH)6 2+ are more stable in the singlet state, while those of B3N3Br6 2+ and B3N3(SeH)6 2+ prefer the triplet state. Besides, exploration of the potential energy surface shows that the ring-shaped structure is the putative global minimum only for B3N3I6 2+. According to chemical bonding analysis, B3N3H3 +, B3N3I6 2+, and B3N3(TeH)6 2+ have σ and π delocalized bonds. The number of delocalized σ/π electrons is 2/6 for the first, and 10/6 for the second and third, similar to what their carbon analogs exhibit. Finally, the analysis of the magnetically induced current density allows B3N3H3 +, B3N3I6 2+, and B3N3(TeH)6 2+ to be classified as strongly σ aromatic, and poorly π aromatic compounds.

Aromatic ouroboroi: heterocycles involving a σ-donor-acceptor bond and 4n + 2 π-electrons

The aromaticity and dynamics of a set of recently proposed neutral 5- and 6-membered heterocycles that are closed by dative (donor-acceptor) or multi-center σ bonds, and have resonance forms with a Hückel number of π-electrons, are examined. The donors and acceptors in the rings include N, O, and F, and B, Be, and Mg, respectively. The planar geometry of the rings, coupled with evidence from different measures of aromaticity, namely the NICSzz, and NICSπzz components of the conventional nucleus independent chemical shifts (NICS), and ring current strengths (RCS), indicate non-trivial degrees of aromaticity in certain cases, including the cyclic C3B2OH6 and C3BOH5 isomers, both with three bonds to the O site in the ring. The former is lower in energy by at least 17.6 kcal mol-1 relative to linear alternatives obtained from molecular dynamics simulations in this work. Some of the other systems examined are best described as non-aromatic. Ring opening, closing, and isomerization are observed in molecular dynamics simulations for some of the systems studied. In a few cases, the ring indeed persists.

Borataalkene or boratabenzene? Understanding the aromaticity of 9-borataphenanthrene anions and its central ring

The recently isolated 1 is isoelectronic to 2 and shows reactivity in central ring that resembles both (cyclo)alkene and aromatic ring. Evaluation of local aromaticity reveal that central ring is best described as a non-aromatic following Clar’s rules.

Calculation of magnetic response properties of tetrazines

Magnetic response properties of 1,2,3,5-tetrazine derivatives including the newly synthesized 4,6-diphenyl-1,2,3,5-tetrazine have been studied computationally at the density functional theory (DFT) level. Calculations of magnetically induced current densities and induced magnetic fields show that the unsubstituted 1,2,3,5-tetrazine is almost as aromatic as benzene. Separating the magnetic shielding functions into molecular orbital components provided additional insights into the magnetic response. The aromatic character estimated from magnetically induced current densities and induced magnetic fields shows that NICS^π_zz(0) values and ring-current strengths yield about the same degree of aromaticity, whereas NICS_zz(0) and NICS_zz(1) values are contaminated by σ electron contributions. The studied 1,2,3,5-tetrazine derivatives are less aromatic than the unsubstituted one. Calculations of magnetic response properties of 4,6-diphenyl-1,2,3,5-tetrazine showed that it is the least aromatic among the studied molecules according to the ring-current criterion, while 4,6-[1,2,3,5]-ditetrazinyl-1,2,3,5-tetrazine is as aromatic as 4,6-dimethyl-1,2,3,5-tetrazine and slightly less aromatic than the unsubstituted 1,2,3,5-tetrazine.

Magnetic response properties of 1,2,3,5-tetrazine derivatives including the newly synthesized 4,6-diphenyl-1,2,3,5-tetrazine have been studied computationally at the density functional theory level.

Relativistic effects on the aromaticity of E3M3H3 (E = C–Pb; M = N–Bi) benzene analogues

The relativistic effects on the aromaticity of a set of benzene analogues, E₃M₃H₃ (E = C–Pb; M = N–Bi) heterocycles, using magnetically induced current density (MICD) and the NICS_zz component of the conventional nucleus independent chemical shift (NICS), is hereby examined.

Does induced current density explain the C-H and C-F Perlin effects?

Spin-spin coupling constant (SSCC) data may be useful in providing information on the stereochemistry and intramolecular interactions in molecules. One-bond C-H and C-F SSCCs (1JCH and 1JCF) are amongst the most important NMR parameters used to study the structure of alicyclic six-membered rings, because the Perlin effect, defined as 1JCHeq > 1JCHax, and the fluorine Perlin-like effect, defined as |1JCFeq| > |1JCFax|, are in wide currency to probe the conformations of these compounds. The origin of these effects has been usually attributed either to dipolar interactions or hyperconjugation, while the induced current density (ICD) has been recently correlated to the magnetic shielding tensors in some six-membered ring compounds, and then used to explain the Perlin effect. Accordingly, this work reports an analysis of the ICD as a descriptor of 1JCH and 1JCF in a series of six-membered rings to find out the role of the ICD in the conventional and fluorine Perlin effect. The atoms in molecules (AIM) magnetic responses obtained from density functional theory (DFT) calculations for the studied compounds did not show any relationship of the first-order electronic current density (ΔJ(1)) with the calculated Δ1JCHax-eq and Δ1JCFax-eq values. Consequently, the title effects cannot be precisely explained by ICD. Nevertheless, an interesting relationship between ΔJ(1) and Δdelocalization involving σCH* orbitals is observed.

Magnetically-induced current density investigation in carbohelicenes and azahelicenes

A computational study of the magnetically-induced current (MIC) density has been carried out for a variety of ortho fused polycyclic aromatic molecules at the density functional theory level with the gauge including magnetically induced current (GIMIC) method. With this method, the aromatic character of each ring in a homologous series of carbohelicenes with an increasing number of fused benzene rings is assessed and compared with other aromaticity criteria such as the Nucleus Independent Chemical Shift [NICS(0), NICSzz(0)] and Bond Length Alternation (BLA) parameters. All criteria indicate that the two outer rings are the most aromatic ones [i.e. higher induced current, more negative NICS(0) and NICSzz(0) values, and smaller BLA values]. For the large helicenes (n > 10), the current drops along the following four rings and then rises again. Additionally, we have proven that this behavior is not due to a difference of the local magnetic field coming from a difference of orientation of the ring with respect to the external magnetic field (oriented along the helical axis). Upon fusing additional benzene rings to form hexa-peri-hexabenzo[7]helicene, some rings (B, D, and F) are a lot less aromatic (even non-aromatic) than the others. The NICS(0) and NICSzz(0) values exaggerate this behavior because they are all positive values, which is a signature of antiaromaticity. Then, when substituting one, three, or four benzene rings with pyrrole ones to form mono-aza-[7]-helicene, tri-aza-[7]-helicene, and tetra-aza-[7]-helicene, remarkable changes in the electronic structures of the helicenes are observed. Indeed, the induced currents are always smaller in the pyrrole rings than in the benzene ones. This has been further investigated using the streamline and the color map representations, which indicate that the diatropic current density passing through the plane cutting the C-N bonds in the pyrrole rings is stronger but is more localized than the current density passing through the plane cutting the C-C bonds in the benzene rings. This gives a positive but smaller total induced current for the pyrrole rings than for the benzene ones. For these systems, the NICS(0) and even more the NICSzz(0) values are not fully reliable to probe the local aromaticity, contrary to the induced current. Indeed, the NICSzz(0) values for the tri-aza-[7]-helicene molecule range from -14.23 to 1.14 ppm, which cannot lead to the same conclusion as the induced current values (10.03 to 12.87 nA T-1).

Interpreting Aromaticity and Antiaromaticity through Bifurcation Analysis of the Induced Magnetic Field

In all molecules, a current density is induced when the molecule is subjected to an external magnetic field. In turn, this current density creates a particular magnetic field. In this work, the bifurcation value of the induced magnetic field is analyzed in a representative set of aromatic, non-aromatic and antiaromatic monocycles, as well as a set of polycyclic hydrocarbons. The results show that the bifurcation value of the ring-shaped domain adequately classifies the studied molecules according to their aromatic character. For aromatic and nonaromatic molecules, it is possible to analyze two ring-shaped domains, one diatropic (inside the molecular ring) and one paratropic (outside the molecular ring). Meanwhile, for antiaromatic rings, only a diatropic ring-shaped domain (outside the molecular ring) is possible to analyze, since the paratropic domain (inside the molecular ring) is irreducible with the maximum value (attractor) at the center of the molecular ring. In some of the studied cases, i. e., in heteroatomic species, bifurcation values do not follow aromaticity trends and present some inconsistencies in comparison to ring currents strengths, showing that this approximation provides only a qualitative estimation about (anti)aromaticity.

Carbon rings decorated with group 14 elements: new aromatic clusters containing planar tetracoordinate carbon

A simple and chemically intuitive approach is used to design ptC-containing E–C clusters (E = Si–Pb).

Metal–ligand ring aromaticity in a 2D coordination polymer used as a photosensitive electronic device

14-Membered metal–ligand ring present in a photosensitive 2D coordination polymer, [Zn₂(fum)₂(4-phpy)₄(H₂O)₂] (H₂fum = fumaric acid and 4-phpy = 4-phenyl pyridine) shows aromatic character as evident by the nucleus-independent chemical shifts (NICS) calculation.

The pseudo-π model of the induced magnetic field: fast and accurate visualization of shielding and deshielding cones in planar conjugated hydrocarbons and spherical fullerenes

Hydrogen skeletal models accurately reproduce the π-induced magnetic field of planar PAHs and spherical fullerenes.