Is Nucleus-Independent Chemical Shift Scan a Reliable Aromaticity Index for Planar Heteroatomic Ring Systems?

How Different are the Diamagnetic and Paramagnetic Contributions to Off-Nucleus Shielding in Aromatic and Antiaromatic Rings?

The spatial variations in the diamagnetic and paramagnetic contributions to the off-nucleus isotropic shielding, σ i s o r = σ i s o d r + σ i s o p r ${\ {{\sigma }_{{\rm i}{\rm s}{\rm o}}\left({\bf r}\right)=\ \sigma }_{{\rm i}{\rm s}{\rm o}}^{{\rm d}}\left({\bf r}\right)+{\sigma }_{{\rm i}{\rm s}{\rm o}}^{{\rm p}}\left({\bf r}\right)}$ , and to the zz component of the off-nucleus shielding tensor, σ z z r = σ z z d r + σ z z p r ${{{\sigma }_{zz}\left({\bf r}\right)=\sigma }_{zz}^{{\rm d}}\left({\bf r}\right)+{\sigma }_{zz}^{{\rm p}}\left({\bf r}\right)}$ , around benzene (C₆ H₆ ) and cyclobutadiene (C₄ H₄ ) are investigated using complete-active-space self-consistent field wavefunctions. Despite the substantial differences between σ i s o r ${{\sigma }_{{\rm i}{\rm s}{\rm o}}\left({\bf r}\right)}$ and σ z z r ${{\sigma }_{zz}\left({\bf r}\right)}$ around the aromatic C₆ H₆ and the antiaromatic C₄ H₄ , the diamagnetic and paramagnetic contributions to these quantities, σ i s o d r ${{\sigma }_{{\rm i}{\rm s}{\rm o}}^{{\rm d}}\left({\bf r}\right)}$ and σ z z d r ${{\sigma }_{zz}^{{\rm d}}\left({\bf r}\right)}$ , and σ i s o p r ${{\sigma }_{{\rm i}{\rm s}{\rm o}}^{{\rm p}}\left({\bf r}\right)}$ and σ z z P r ${{\sigma }_{zz}^{{\rm P}}\left({\bf r}\right)}$ , are found to behave similarly in the two molecules, shielding and deshielding, respectively, each ring and its surroundings. The different signs of the most popular aromaticity criterion, the nucleus-independent chemical shift (NICS), in C₆ H₆ and C₄ H₄ are shown to follow from a change in the balance between the respective diamagnetic and paramagnetic contributions. Thus, the different NICS values for antiaromatic and antiaromatic molecules cannot be attributed to differences in the ease of access to excited states only; differences in the electron density, which determines the overall bonding picture, also play an important role.

New Carbenes and Cyclic Allenes Energetically Comparable to Experimentally Known 1-Azulenylcarbene

1-Azulenylcarbene (18; 0 kJ mol^-1) is experimentally known as the key reactive intermediate for the rearrangement reactions of aryl carbenes in the laboratory. Here, using coupled-cluster methods up to the fc-CCSD(T)/cc-pVTZ//fc-CCSD(T)/cc-pVDZ level, thirteen new carbenes and one new cyclic allene are theoretically identified within the C₁₁H₈ elemental composition that either energetically lie below or very close to 18. While the cyclic allene, bicyclo[5.4.0]undeca-2,3,5,7,9,11-hexene (1; -166 kJ mol^-1), is the experimentally known lowest energy isomer, three other cyclic allenes, bicyclo[5.4.0]undeca-1,2,4,6,8,10-hexene (2; -100 kJ mol^-1), bicyclo[5.4.0]undeca-1,3,4,6,8,10-hexene (3; -97 kJ mol^-1), and bicyclo[6.3.0]undeca-1,2,4,6,8,10-hexene (13; -42 kJ mol^-1), demand new experimental studies. In total, thirty-one isomers are studied in this work (within -166 to +15 kJ mol^-1 from 18) and all are found to be polar (μ ≠ 0). Among these, 1H-benzo[7]annulen-1-ylidene (17; -4 kJ mol^-1; μ = 5.24 D), bicyclo[5.4.0]undeca-2,4,6,8,11-pentaene-10-ylidene (24; 13 kJ mol^-1; μ = 7.59 D), 5-methylene-naphthalen-1-ylidene (26; 15 kJ mol^-1; μ = 5.32 D), 6-methylene-naphthalen-2-ylidene (27; -43 kJ mol^-1; μ = 6.60 D), and 8-methylene-naphthalen-2-ylidene (28; -39 kJ mol^-1; μ = 5.55 D) are competitively polar compared to 18 (μ = 5.39 D). Therefore, these carbene molecules are potential targets for rotational spectroscopists and radioastronomers. Considering the importance of naphthyl and azulenylcarbenes in reactive intermediate chemistry, mechanisms of different rearrangement reactions and plausible formation pathways of some of these new carbenes are studied in this work using density functional theory.

Structures, Energetics, and Spectra of (NH) and (OH) Tautomers of 2-(2-Hydroxyphenyl)-1-azaazulene: A Density Functional Theory/Time-Dependent Density Functional Theory Study

Tautomerization of 2-(2-hydroxyphenyl)-1-azaazulene (2OHPhAZ) in the gas phase and ethanol has been studied using B3LYP, M06-2X, and ωB97XD density functional theory (DFT) with different basis sets. For more accurate data, energies were refined at CCSD(T)/6-311++G(2d,2p) in the gas phase. Nuclear magnetic resonance (NMR), aromaticity, Fukui functions, acidity, and basicity were also calculated and compared with experimental data. Time-dependent density functional theory (TDDFT)-solvation model based on density (TDDFT-SMD) calculations in acetonitrile have been utilized for the simulation of UV-vis electronic spectra. In addition, electronic structures of the investigated system have been discussed. The results reveal that the enol form (2OHPhAZ) is thermodynamically and kinetically stable relative to the keto tautomer (2OPhAZ) and different rotamers (2OHPhAZ-R1:R3) in the gas phase and ethanol. A comparison with the experiment illustrates a good agreement and supports the computational findings.

Six Low-Lying Isomers of C11H8 Are Unidentified in the Laboratory-A Theoretical Study

Isomers of C₁₁H₈ have been theoretically examined using density functional theory and coupled-cluster methods. The current investigation reveals that 2aH-cyclopenta[cd]indene (2), 7-ethynyl-1H-indene (6), 4-ethynyl-1H-indene (7), 6-ethynyl-1H-indene (8), 5-ethynyl-1H-indene (9), and 7bH-cyclopenta[cd]indene (10) remain elusive till date in the laboratory. The puckered low-lying isomer 2 lies at 9 kJ mol^-1 below the experimentally known molecule, cyclobuta[de]naphthalene (3), at the fc-CCSD(T)/cc-pVTZ//fc-CCSD(T)/cc-pVDZ level of theory. 2 lies at 36 kJ mol^-1 above the thermodynamically most stable and experimentally known isomer, 1H-cyclopenta[cd]indene (1), at the same level. It is identified that 1,2-H transfer from 1 yields 2H-cyclopenta[cd]indene (14) and subsequent 1,2-H shift from 14 yields 2. Appropriate transition states have been identified, and intrinsic reaction coordinate calculations have been carried out at the B3LYP/6-311+G(d,p) level of theory. Recently, 1-ethynyl-1H-indene (11) has been detected using synchrotron-based vacuum ultraviolet ionization mass spectrometry. 2-Ethynyl-1H-indene (4) and 3-ethynyl-1H-indene (5) have been synthetically characterized in the past. While the derivatives of 7bH-cyclopenta[cd]indene (10) have been isolated elsewhere, the parent compound remains unidentified till date in the laboratory. Although C₁₁H₈ is a key elemental composition of astronomical interest for the formation of polycyclic aromatic hydrocarbons in the interstellar medium, none of its low-lying isomers have been characterized by rotational spectroscopy though they are having a permanent dipole moment (μ ≠ 0). Therefore, energetic and spectroscopic properties have been computed, and the present investigation necessitates new synthetic studies on C₁₁H₈, in particular 2, 6-10, and also rotational spectroscopic studies on all low-lying isomers.

Kinetic Stability of Si2C5H2 Isomer with a Planar Tetracoordinate Carbon Atom

Dissociation pathways of the global minimum geometry of Si2C5H2 with a planar tetracoordinate carbon (ptC) atom, 2,7-disilatricyclo[4.1.0.01,3]hept-2,4,6-trien-2,7-diyl (1), have been theoretically investigated using density functional theory and coupled-cluster (CC) methods. Dissociation of Si-C bond connected to the ptC atom leads to the formation of 4,7-disilabicyclo[4.1.0]hept-1(6),4(5)-dien-2-yn-7-ylidene (4) through a single transition state. Dissociation of C-C bond connected to the ptC atom leads to an intermediate with two identical transition states and leads back to 1 itself. Simultaneous breaking of both Si-C and C-C bonds leads to an acyclic transition state, which forms an acyclic product, cis-1,7-disilahept-1,2,3,5,6-pentaen-1,7-diylidene (19). Overall, two different products, four transition states, and an intermediate have been identified at the B3LYP/6-311++G(2d,2p) level of theory. Intrinsic reaction coordinate calculations have also been done at the latter level to confirm the isomerization pathways. CC calculations have been done at the CCSD(T)/cc-pVTZ level of theory for all minima. Importantly, all reaction profiles for 1 are found be endothermic in Si2C5H2. These results are in stark contrast compared to the structurally similar and isovalent lowest-energy isomer of C7H2 with a ptC atom as the overall reaction profiles there have been found to be exothermic. The activation energies for Si-C, C-C, and Si-C/C-C breaking are found to be 30.51, 64.05, and 61.85 kcal mol−1, respectively. Thus, it is emphasized here that 1 is a kinetically stable molecule. However, it remains elusive in the laboratory to date. Therefore, energetic and spectroscopic parameters have been documented here, which may be of relevance to molecular spectroscopists in identifying this key anti-van’t-Hoff-Le Bel molecule.

MgC6H2 Isomers: Potential Candidates for Laboratory and Radioastronomical Studies

Eighty three stationary points of MgC₆H₂ isomers spanning from 0 to 215 kcal mol^-1 have been theoretically identified using density functional theory at the B3LYP/6-311++G(2d,2p) level of theory. Among them, four low-lying isomers lying within 23.06 kcal mol^-1 (1 eV) have been further characterized in detail using high-level coupled-cluster (CC) methods. The thermodynamically most stable isomer turns out to be 1-magnesacyclohepta-4-en-2,6-diyne (1). The other three isomers, 3-magnesahepta-1,4,6-triyne (2), 1-magnesacyclohepta-2,3,4-trien-6-yne (3), and 1-magnesahepta-2,4,6-triyne (4) lie 8.24, 19.76, and 21.36 kcal mol^-1, respectively, above 1 at the ae-CCSD(T)/cc-pCVTZ level of theory. All the four isomers are polar with a permanent electric dipole moment (μ ≠ 0). Hence, they are potential candidates for rotational spectroscopic studies. Considering the recent identification of magnesium-bearing hydrocarbons such as, MgC₂H and MgC₄H in IRC+10216, it is believed that the current theoretical data may be of relevance to laboratory molecular spectroscopic and radioastronomical studies on MgC₆H₂ isomers. The energetic and spectroscopic information gathered in this study would aid the detection of low-lying MgC₆H₂ isomers in the laboratory, which are indispensable for radioastronomical studies. It is also noted here that neither the National Institute of Standards and Technology Chemistry WebBook nor the Kinetic Database for Astrochemistry lists any isomer of MgC₆H₂ at the moment. Therefore, these isomers are studied here theoretically for the very first time.

Theoretical Studies of SiC4H2 Isomers Delineate Three Low-Lying Silylidenes Are Missing in the Laboratory

Eleven isomers of SiC₄H₂ lying within 50 kcal mol^-1 have been theoretically investigated using density functional theory and high-level coupled-cluster methods. Among them, four isomers, 1-ethynyl-3-silacycloprop-1(2)-en-3-ylidene (1), diethynylsilylidene (2), 1-sila-1,2,3,4-pentatetraenylidene (4), and 1,3-butadiynylsilylidene (5), have already been identified in the laboratory. The current investigation reports three low-lying (<1 eV) silylidenes [2-methylenesilabicyclo[1.1.0]but-1(3)-en-4-ylidene (3), 4-sila-2-methylenebicyclo[1.1.0]but-1(3)-en-4-ylidene (6), and 3-ethynyl-1-silapropadienylidene (7)] and three high-lying (>1 eV) silylidenes [2-sila-(didehydrovinylidene)cyclopropene (8), an isomer with a planar tetracoordinate carbon (ptC) atom (10), and 1-ethynyl-1-silapropadienylidene (11)], which remain elusive in the laboratory to date. Isomer 9 also contains a ptC atom, which turned out to be a transition state at all levels. Though all isomers are polar (μ ≠ 0), rotational spectrum is available only for 4. Using matrix isolation, three isomers (1, 2, and 5) have been trapped in the laboratory at 10 K. Considering the astrochemical relevance of silicon-carbide clusters in the interstellar medium, the current theoretical data demand new molecular spectroscopic studies on SiC₄H₂. Surprisingly, unlike the isovalent C₅H₂ isomers, where the bent carbenes are yet to be identified in the laboratory, the bent silylidenes (2 and 5) have been trapped in the case of SiC₄H₂. In both the cases, molecules with transannular C-C and/or Si-C bonds remain elusive, though they lie in the low-lying region. Using suitable precursors, whether these peculiar geometries (especially 3 and 6) would be identified or not in the laboratory needs to be addressed by molecular spectroscopists. The present investigation documents structural and spectroscopic information of SiC₄H₂ isomers, which may compliment future molecular spectroscopic observations including radioastronomical searches.

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.

Analysis of the Nucleus-Independent Chemical Shifts of [10]Cyclophenacene: Is It an Aromatic or Antiaromatic Molecule?

[10]Cyclophenacene is an important synthetic target that shows a pair of nucleus-independent chemical shift (NICS) values for the center of mass and six-membered rings typical of an aromatic species. This is found in contrast with the global paratropic current density induced by a magnetic field parallel to the main symmetry axis. This apparent contradiction has been analyzed by studying the tensor character of the magnetic response. It turns out that the molecule displays two characters, one paratropic (antiaromatic) and another one diatropic (aromatic), depending on the orientation of the inducing magnetic field. The paratropic response, which cannot be recognized from the NICS values, is associated with a well-defined destabilization of the belt closure, as witnessed by homodesmotic reactions. A scalar measure of magnetic aromaticity, the field-independent current strength, has been introduced, which allows us to reach the conclusion that [10]cyclophenacene is indeed an aromatic molecule, although it is significantly affected by the paratropic response.

CC/B–N substitution in five membered heterocycles. A computational analysis

Novel five-membered azaboroles are aromatic, stable under neutral conditions, isomer stabilization energy is explained using σ-bond and aromatic stabilization energies.

Borazine: spin blocker or not?

Spin-blocker capacity of borazine is investigated formeta-BB,meta-NN andpara-BN structures highlighting the correlation between magnetic properties and aromaticity.

The critical re-evaluation of the aromatic/antiaromatic nature of Ti3(CO)3: a missed opportunity?

The nature of bonding and aromaticity of Ti(3)(CO)(3), a mill-shaped metal-carbonyl complex, is studied carefully. A unique bonding mechanism between metal and carbonyl groups is found in this species. Ti(3)(CO)(3) is an example of a metal-carbonyl complex with prominent metal to carbonyl donation. Moreover, it is proven that not only is Ti(3)(CO)(3) not an antiaromatic complex but also it is the first synthesized example of d-block, σ+π aromatic species. A quick survey among the first row of transition metals in the periodic table shows that other local minima with similar structures and aromaticity are present and Ti(3)(CO)(3) is the first synthesized species of an unknown family.