Polygonal current model: an effective quantifier of aromaticity on the magnetic criterion

On the JAP Method for the Indirect Determination of Delocalized Currents from Experimental Chemical Shifts

The JAP model (after Jirásek, Anderson, and Peeks) to retrieve global current strengths from experimental 1 H chemical shifts has been tested with DFT computations. Both global and local tropicities are correctly predicted in most cases and the quantitative agreement is overall fair. An extension of the model is found to give improvement in an exemplary critical case, where the global delocalized current is negligible and the current density map is dominated by local currents.

Revival of Hückel Aromatic (Poly)benzenoid Subunits in Triplet State Polycyclic Aromatic Hydrocarbons by Silicon Substitution

By employing density functional theory (DFT) calculations we show that mono- and disilicon substitution in polycyclic aromatic hydrocarbons, having two to four benzene units, quenches their triplet state antiaromaticity by creating Hückel aromatic (poly)benzenoid subunit(s) and weakly antiaromatic, or almost nonaromatic silacycle. Therefore, such systems are predicted to be globally aromatic in both the ground state and the first excited triplet state. Putting the silicon atom(s) into various positions of a hydrocarbon provides an opportunity to tune the singlet-triplet energy gaps. They depend on the global aromaticity degree which, in turn, depends on the type of aromatic carbocyclic subunit(s) and the extent of their aromaticity. On the basis of the set of studied compounds, some preliminary rules on how to regulate the extent of global, semiglobal and local aromaticity are proposed. The results of this work extend the importance of Hückel aromaticity concept to excited triplet states which are usually characterized by the Baird type of (anti)aromaticity.

Magnetic Aromaticity of Cycloporphyrin Nanorings

The ascertainment of magnetic aromaticity is not necessarily straightforward, especially for large and bent systems, such as the cycloporphyrin nanorings recently synthesized by the group of Anderson. Six of these cycloporphyrin nanorings were studied here computationally. Indirect methods, based on nuclear shielding and magnetizabilities, and direct methods, based on standard quantum mechanics, were both used effectively to determine their magnetically induced current strength, which mostly confirmed Anderson’s classification. However, in the case of hexanions, and in particular for cyclohexaporphyrin hexacations, a significant cancellation of delocalized diatropic and paratropic flow occurred, showing that the resultant faint aromatic character was a result of competing aromatic and antiaromatic contributions, as also evidenced by the ipsocentric method. A warning is renewed on the use of isotropic shielding to determine the tropicity of the magnetically induced current.

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.

Anisotropy of the vorticity tensor as a magnetic indicator of aromaticity

The vorticity vector of the current density J^B, induced in the electron cloud of a molecule by a magnetic field B, is defined by V^B = ∇ × J^B.

Doubly zwitterionic, di-reduced, highly electron-rich, air-stable naphthalenediimides: redox-switchable islands of aromatic-antiaromatic states

The di-reduced state of the naphthalene moiety and its congeners have long captivated chemists as it is elusive to stabilize these intrinsically reactive electron-rich π-systems and for their emergent multifaceted properties. Herein we report the synthesis and isolation of two-electron (2e-) reduced, highly electron-rich naphthalenediimides (NDIs). A doubly zwitterionic structure is observed for the first time in a naphthalene moiety and validated by single crystal X-ray crystallography and spectroscopic methods. The synthesis avoids hazardous reducing agents and offers an easy, high-yielding route to bench-stable di-reduced NDIs. Notably, we realized high negative first oxidation potentials of up to -0.730 V vs. Fc/Fc+ in these systems, which establish these systems to be one of the strongest ambient stable electron donors. The study also provides the first insights into the NMR spectra of the di-reduced systems revealing a large decrease in diatropicity of the naphthalene ring compared to its 2e- oxidized form. The NICS, NICS-XY global ring current, gauge-including magnetically induced current (GIMIC) and AICD ring current density calculations revealed switching of the antiaromatic and aromatic states at the naphthalene and the imide rings, respectively, in the di-reduced system compared to the 2e- oxidized form. Notably, the substituents at the phosphonium groups significantly tune the antiaromatic-aromatic states and donor ability, and bestow an array of colors to the di-reduced systems by virtue of intramolecular through-space communication with the NDI scaffold. Computational studies showed intramolecular noncovalent interactions to provide additional stability to these unprecedented doubly zwitterionic systems.

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).

Is the pericyclic transition structure of aza-Diels-Alder reaction aromatic?

The notion of aromaticity is an elusive concept, typically delineated as an electron delocalization pattern within a cyclic structure, and is characterized by unusual stability, reactivity, and other properties. Recently, the aromatic concept has been extended to inorganic molecules and saturated systems. In our work, the so-called pericyclic transition state structures (TSs) involved in aza-Diels-Alder reactions are investigated. Interestingly, geometries and energy barriers evidence that these reactions are highly energetic and are extremely asynchronous. Additionally the localized orbital locator, ELF and QTAIM topological analyses have proven that at least one pair of atoms is not bound at the TS configuration. Moreover, HOMA and PDI indexes show that these TSs where aromaticity is not important. Therefore, these results provoke us to rule out the aromatic character of these so-called pericyclic TSs structures.

Magnetically induced current density in triple-layered beryllium–boron clusters

Magnetically induced current densities reveal the double aromatic character of the examined Be–B clusters.

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.

Calculations of current densities and aromatic pathways in cyclic porphyrin and isoporphyrin arrays

Magnetically induced current density susceptibilities have been studied for a number of cyclic ethyne- and butadiyne-bridged porphyrin and isoporphyrin arrays. The current density susceptibilities have been calculated using the gauge-including magnetically induced current (GIMIC) method, which is interfaced to the TURBOMOLE quantum chemistry code. Aromatic properties and current pathways have been analyzed and discussed by numerical integration of the current density susceptibilities passing selected chemical bonds yielding current strength susceptibilities. Despite the interrupted π-framework, zinc(ii) isoporphyrin sustains a ring current of ca. 10 nA T^-1. Porphyrin and isoporphyrin dimers sustain a significant current strength at the linker, whereas the larger porphyrinoid arrays sustain mainly local ring currents. Isoporphyrin dimers with saturated meso carbons have strong net diatropic ring-current strengths of 20 nA T^-1 fulfilling Hückels aromaticity rule. Porphyrin trimers and tetramers exhibit almost no current strength at the linker. The porphyrin moieties maintain their strong net diatropic ring current.

Can the current density map topology be extracted from the nucleus independent chemical shifts?

Aromatic compounds are characterised by the presence of a ring current when in a magnetic field. As a consequence, current density maps are used to assess (the degree of) aromaticity of a compound. However, often a more discrete set of so-called Nucleus Independent Chemical Shift (NICS) values is used that is derived from the current density. It is shown here that there is no simple one-to-one relationship that allows reconstructing current density maps from only NICS-values. NICS values should therefore not be used as aromaticity indices without analysis of the ab initio computed current density map.

The influence of heteroatoms on the aromatic character and the current pathways of B2N2-dibenzo[a,e]pentalenes

Four polycyclic molecules have been investigated at the DFT level of theory using magnetically induced current densities calculated with the GIMIC program as an indicator of their (anti)aromaticity.

Topological definition of ring currents

A definition of ring currents in a velocity vector field is proposed according to topological criteria: ring currents are axial vortices confined in, or rotating beyond, a separatrix, i.e., the boundary which marks the limits of the vortex.

The making of ring currents

The π-electron diatropic ring current emerging throughout the trimerization of acetylene to benzene.

Analysis of the magnetically induced current density of molecules consisting of annelated aromatic and antiaromatic hydrocarbon rings

The aromatic pathway of molecules with annelated aromatic and antiaromatic rings has been studied by calculating magnetically induced current densities.

New insights into aromatic pathways of carbachlorins and carbaporphyrins based on calculations of magnetically induced current densities

The aromatic pathways of carbaporphyrins and carbachlorins that are based on magnetically induced current density DFT-GIMIC calculations are presented and discussed.

Understanding the molecular switching properties of octaphyrins

Triggering Hückel–Möbius topological and aromaticity switches in octaphyrins by protonation and redox reactions.

“Carbo-aromaticity” and novel carbo-aromatic compounds

Recent advances in experimental and theoretical studies ofcarbo-benzene derivatives, along with the proposition of a generalization of the definition of aromaticity to the two-membered π-rings of triple bonds, suggest relevance for the notion of “carbo-aromaticity”.

The NICS-XY-scan: identification of local and global ring currents in multi-ring systems

Nucleus-independent chemical shift (NICS)-based methods are very popular for the determination of the induced magnetic field under an external magnetic field. These methods are used mostly (but not only) for the determination of the aromaticity and antiaromaticity of molecules and ions, both qualitatively and quantitatively. The ghost atom that serves as the NICS probe senses the induced magnetic field and reports it in the form of an NMR chemical shift. However, the source of the field cannot be determined by NICS. Thus, in a multi-ring system that may contain more than one induced current circuit (and therefore more than one source of the induced magnetic field) the NICS value may represent the sum of many induced magnetic fields. This may lead to wrong assignments of the aromaticity (and antiaromaticity) of the systems under study. In this paper, we present a NICS-based method for the determination of local and global ring currents in conjugated multi-ring systems. The method involves placing the NICS probes along the X axis, and if needed, along the Y axis, at a constant height above the system under study. Following the change in the induced field along these axes allows the identification of global and local induced currents. The best NICS type to use for these scans is NICSπZZ , but it is shown that at a height of 1.7 Å above the molecular plane, NICSZZ provides the same qualitative picture. This method, namely the NICS-XY-scan, gives information equivalent to that obtained through current density analysis methods, and in some cases, provides even more details.