Can Anti-Aufbau DFT Calculations Estimate Singlet Excited State Aromaticity? Correspondence on “Dibenzoarsepins: Planarization of 8π-Electron System in the Lowest Singlet Excited State”

The simple anti-aufbau DFT approach for estimating singlet excited state aromaticity suggested in a recent Communication published in this journal is shown to produce incorrect results because it targets a linear combination of the singlet and triplet configurations involving the HOMO and LUMO rather than the first singlet excited state. If the S₁ state of a molecule is dominated by the HOMO→LUMO excitation, a comparably simple but theoretically consistent and qualitatively correct approximation to the S₁ wavefunction can be achieved by performing a small “two electrons in two orbitals” CASSCF(2,2) calculation which can be followed by the evaluation of magnetic aromaticity criteria such as NICS.     

Nucleus Independent Chemical Shift (NICS) at Small Distances from the Molecular Plane: The Effect of Electron Density

When close to the molecular plane, the behavior of nucleus independent chemical shift (NICS) as a function of the distance from the molecular plane deviates from its behavior at larger distances. By using a dense grid of NICS-probes (BQs) it is shown that, when close to the molecular plane, maximal (absolute) NICS values are obtained above the atoms. These maxima move towards the center as the grid is elevated until the (absolute) maximum NICS is obtained at the center and stay there when the grid is further elevated. It is shown that this behavior is a result of the current density, which is influenced by the electron density, according to the Biot-Savart law, which, in turn, causes the induced magnetic field measured by the NICS. It is thus concluded that if magnetic aromaticity is studied, the NICS calculations should be carried out at a large enough distance so that only the π-ring current affects the NICS. At distances ≥2 Å, NICS(r)_π,zz=A+B*C^r. Using non-linear correlation for obtaining A, B and C and extrapolate to NICS(1)_π,zz and NICS(1.7)_π,zz is recommended as measures for aromaticity.     

Analysis of Non-innocence of Phosphaquinodimethane Ligands when Charge and Aromaticity Come into Play

Several phosphaquinodimethanes and their M(CO)₅ complexes (M=Cr, Mo, W) and model derivatives have been theoretically investigated regarding the quest of non-innocence. Computed structural and electronic properties of the P-Me/NH₂ substituted phosphaquinodimethanes and tungsten complexes revealed an interesting non-innocent ligand behaviour for the radical anion complexes with distonic ion character and a strong rearomatization of the middle phenyl ring. The latter was further probed taking also geometric aromaticity (HOMA) and quinoid distortion parameters (HOMQc) into account, as well as NICS(1). Furthermore, the effect of the P-substitution was investigated for real (or plausible) complexes and their free ligands focusing on the resulting aromaticity at the middle phenyl ring and vertical one-electron redox processes. The best picture of ligand engagement in redox changes was provided by representing NICS(1) values versus HOMA and the new geometric distortion parameter HOMQc8.     

Bonding,Aromaticity and Isomerization of Furfuraldehyde through Off-Nucleus Isotropic Magnetic Shielding

Off-nucleus isotropic magnetic shielding (σ_iso(r)) and multi-points nucleus independent chemical shift (NICS(0-2 Å)) index were utilized to find the impacts of the isomerization of gas-phase furfuraldehyde (FD) on bonding and aromaticity of FD. Multidimensional (1D to 3D) grids of ghost atoms (bqs) were used as local magnetic probes to evaluate σ_iso(r) through gauge-including atomic orbitals (GIAO) at density functional theory (DFT) and B3LYP functional/6-311+G(d,p) basis set level of theory. 1D σ_iso(r) responses along each bond of FD were examined. Also, a σ_iso(r) 2D-scan was performed to obtain σ_iso(r) behavior at vertical heights of 0–1 Å above the FD plane in its cis, transition state (TS) and trans forms. New techniques for evaluating 2D σ_iso(r) cross-sections are also included. Our findings showed that bonds in cyclic and acyclic parts of FD are dissimilar. Unlike the C−O bond of furanyl, the C=O bond of the formyl group is magnetically different. C−C and C−H bonds in furanyl are found similar to those in aromatic rings. A unique σ_iso(r) trend was observed for the C₂−C₆ bond during FD isomerization. Based on NICS(0-2 Å) values, the degree of aromaticity follows the order of cis FD<trans FD<furan<TS FD.     

The anisotropic effect of functional groups in H NMR spectra is the molecular response property of spatial NICS—the frozen conformational equilibria of 9-arylfluorenes

Rotation about the single bond adjoining the aryl and fluorene moieties in 9-arylfluorenes can be frozen out on the NMR timescale if methyl groups are located at either one or both of the ortho positions of the aryl substituent. In the ground-state of these rotamers, the planes of the aryl and fluorene moieties are perpendicular to each other and the methyl substituents are consequently positioned either above the fluorene moiety or in-plane with it; thus, the methyl protons are either shielded or deshielded, respectively, due to the ring current effect of the fluorene moiety. This anisotropic effect on the ¹H chemical shifts of the methyl protons has been quantified on the basis of through-space NMR shieldings (TSNMRS) and subsequently Δδ_calcd compared with the experimentally observed chemical shift differences, Δδ_exp. In this context, the experimental anisotropic effects of functional groups in the ¹H NMR have proven to quantitatively be the molecular response property of theoretical spatial nucleus independent chemical shieldings (NICS). Differences between Δδ_calcd and Δδ_exp were, for the first time, also quantified as arising from steric compression.     

Aromaticity of N-heterocyclic carbene and its analogues: Magnetically induced ring current perspective

The N-heterocyclic carbene, imidazole-2-ylidene, and its main group (13-15) analogues contain cyclically conjugated 6π electrons. Experimental ¹H nuclear magnetic resonance (NMR) spectra suggest an increase in aromaticity along a period from left to right. Whereas the order along a group is as follows: period 2 > period 5 > period 4 > period 3 due to change in structure. To understand the order of aromaticity, the magnetically induced ring currents of the molecules are calculated using aromatic ring current shielding, gauge-including magnetically induced currents (GIMIC) method and Stanger's σ-model applying the gauge-including atomic orbitals NMR technique. It is found that GIMIC best describes the order of aromaticity especially along a group where current-profile changes on the bivalent atom down a group due to change in electron density. Moreover, the GIMIC provides the visualization of current by sign modulus and the anisotropy of the induced current density plots.     

Synthesis and Theoretical Investigation of a 1,8‐Bis(bis(diisopropylamino)cyclopropeniminyl)naphthalene Proton Sponge Derivative

We report herein the synthesis and characterization of a new proton sponge derivative, 1,8‐bis(bis(diisopropylamino)cyclopropeniminyl)naphthalene 4 (DACN), as well as its bis‐protonated counterpart 6 . A crystal structure of 6 is presented, along with variable temperature ¹H NMR data on the BF₄^? salt ( 6?BF₄ ). DFT calculations were performed to investigate the structure of the monoprotonated species 7 and to gain insight into the structural and electronic nature of all three species. The proton affinity (PA) of 4 , calculated at the B3LYP/6‐311G++(d,p)//B3LYP/6‐31G(d,p) level, taking into account thermal corrections from the B3LYP/6‐31G(d,p) method, was 282.3 kcal mol^?1, while its pK_a was estimated at 27.0. NICS calculations were performed to examine the changes in aromaticity within these systems upon each successive protonation. Lastly, homodesmotic reaction schemes were used in order to estimate the factors contributing to the strong PA predicted for 4 .     

Substitute Effect on the Structure,Stability of Valence Isomers and Aromaticity of 1‐H‐Boratabenzene

Density functional theory (B3LYP) calculations were performed on the Me and F substituted valence isomeric forms of 1‐H‐boratabenzene. The calculations revealed that the planar benzene analog is the lowest energy isomer. Its aromaticity is analyzed in the light of the nucleus‐independent chemical shift (NICS) and shows that aromaticity increases in F substituted, but decreases in Me substituted. These calculations indicate substitution of BH with BMe and BF doesn't cause significant variation in bond length.     

Structure,stability, and aromaticity of M?SubPc (MB,Al, and Ga): Computational study

A theoretical investigation on the structure, stability, and aromaticity of M‐subphthalocyanine (M? SubPc; M?B, Al, and Ga) was performed at the B3LYP/6‐31+G*//B3LYP/6‐31G* level. The comparison between M? SubPc and the corresponding M? phthalocyanine (M? Pc) was considered. The geometry optimization of the M? SubPc shows that in the Al? SubPc and Ga? SubPc, the steric repulsions among the three azacoupled isoindole moieties increase, as to their macrocycles tend to be far from planarity. The binding energies of Cl? M … aza‐coupled isoindole corrected by the basis set superposition error (BSSE), and the nucleus‐independent chemical shift (NICS) values at the ring center, which are a simple and effective local aromaticity probe, were calculated. The results show that Al? SubPc is less stable than both B? SubPc and Al? Pc for larger steric repulsion, smaller binding energy, and weaker aromaticity. In the same way, Ga? SubPc is less stable than both B? SubPc and Ga? Pc. In addition, the ring expansion reactivity occurring in B? SubPc was confirmed by the global aromaticity mirrored by the electrophilicity index ω values. Therefore, the Al? SubPc and Ga? SubPc remain unknown, while the corresponding compounds Al? Pc and Ga? Pc are known experimentally. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005