The unexpected roles of σ and π orbitals in electron donor and acceptor group effects on the 13C NMR chemical shifts in substituted benzenes

The Amsterdam Modeling Suite

In this paper, we present the Amsterdam Modeling Suite (AMS), a comprehensive software platform designed to support advanced molecular and materials simulations across a wide range of chemical and physical systems. AMS integrates cutting-edge quantum chemical methods, including Density Functional Theory (DFT) and time-dependent DFT, with molecular mechanics, fluid thermodynamics, machine learning techniques, and more, to enable multi-scale modeling of complex chemical systems. Its design philosophy allows for seamless coupling between components, facilitating simulations that range from small molecules to complex biomolecular and solid-state systems, making it a versatile tool for tackling interdisciplinary challenges, both in industry and in academia. The suite also emphasizes user accessibility, with an intuitive graphical interface, extensive scripting capabilities, and compatibility with high-performance computing environments.

An investigation of contributors to the spin exchange interactions in organic pentacene-radical dyads using quasi-degenerate perturbation theory

Chromophore-radical (C-R) dyads are a promising class of molecules with potential applications in magnetometry, nuclear magnetic resonance and quantum sensing. Given the vast chemical space that is possible in these systems, computational studies are vital to aid in the rational design of C-R molecules with desired electronic and spin properties. Multireference perturbation theory (MRPT) calculations have been shown to be useful for rationalizing spin correlations in C-R dyads. In this work we apply quasi-degenerate perturbation theory, specifically QD-NEVPT2, for the prediction of vertical transition energies (VTEs) as well as spin-correlation parameters in three-spin-center pentacene-radical dyads containing up to 153 atoms. We find that QD-NEVPT2 performs well in the prediction of JTR, the magnetic coupling parameter between the excited-state triplet and the radical, but underestimates VTEs; this underestimation is attributed to variational averaging over different spin states and active space limitations, and we show that addressing these shortcomings reduces error. The calculated magnitudes and signs of JTR are rationalized through molecular symmetry, coupling distance, and π-structure considerations. The predicted signs of JTR are consistent with and explained via mechanisms of kinetic and potential spin-exchange, allowing for future functional design of magnetic organic molecules. The role of active space choice on VTE accuracy and predicted magnetic coupling is additionally explored.

Global and Local Shielding/Deshielding Variation in C60 and C70 Fullerenes Upon Reduction: Evaluation of Aromaticity, 13C-NMR and Anisotropy Patterns

Fullerenes are statically pleasant species featuring symmetric cages, which can be modified upon reduction. Here, we theoretically account for the variation of 13C-NMR patterns in C60 and C70 upon six-fold reduction and the overall variation of the enabled shielding/deshielding regions induced by π and σ electrons according to different orientations of the external field and the related anisotropy. Our results show a significant modification of the chemical shift given by the main variation of the σ33 (or δ33) shielding component under the principal axis system (PAS) of the chemical shift anisotropy (CSA) at the representative carbon nucleus. For C60 6- a shielding cone property is enabled from any orientation, accounting for a significant spherical aromatic character. In contrast, in C70 6-, a shielding cone is reserved only for an axial-oriented field, with a deshielding cone behavior obtained from the complementary equatorial orientations. The overall anisotropy shows an inner isotropic region for C60 and C60 6-, with a continuous anisotropic outer contour for the latter. In contrast, C70 and C70 6- both show larger anisotropies, given the lesser spherical shape in addition to the modified π-surface. Such information is useful for further rationalizing the implementation of magnetic anisotropic molecular devices into fullerene-based materials.

An Acceptor-Substituted N-Heterocyclic ortho-Quinodimethane: Pushing the Boundaries of Polarization in Donor-Acceptor-Substituted Polyenes

We report the synthesis, isolation, and characterization of a stable donor-acceptor substituted ortho-quinodimethane (oQDM). This system with an imidazolidine scaffold as the donor can also be referred to as acceptor-substituted ortho-N-heterocyclic quinodimethane (oNHQ). We have examined the extent of polarization of the conjugated π-system using single-crystal X-ray diffraction, NMR and UV/vis spectroscopy, cyclic voltammetry, and DFT computations. The bond lengths in the phenyl linker do not exhibit the alternation typical of oQDMs. In addition, the 13C and 15N NMR shifts suggest significant charge separation, an interpretation supported by the diatropic ring current determined by NICSZZ(r) computations, which is characteristic of aromatic compounds. DFT calculations show that polarization is an electronic effect that is amplified by steric influences. More strikingly, the oxidation and reduction potentials of the push-pull substituted oQDM are virtually identical to those of authenticated anionic and cationic derivatives. The results therefore indicate that an aromatic zwitterionic structure represents the electronic structure more accurately than a neutral quinoidal Lewis structure, which indicates that the acceptor-substituted oNHQ is a rare example of an organic zwitterion in which the centers of charge are in conjugation. The ambiphilic reactivity of the acceptor-substituted oNHQ, which is evidenced by the dehydrogenation of ammonia borane and the addition of phenylacetylene via heterolytic C-H bond cleavage, further supports its notation as an organic zwitterion and is reminiscent of frustrated Lewis pairs (FLPs). Thus, the acceptor-substituted oNHQ can be considered to be an intramolecular carbogenic FLP in terms of its reactivity.

Elucidating the Electronic Nature of Rh-based Paddlewheel Catalysts from 103 Rh NMR Chemical Shifts: Insights from Quantum Mechanical Calculations

The tremendous importance of dirhodium paddlewheel complexes for asymmetric catalysis is largely the result of an empirical optimization of the chiral ligand sphere about the bimetallic core. It was only recently that a H(C)Rh triple resonance 103 Rh NMR experiment provided the long-awaited opportunity to examine - with previously inconceivable accuracy - how variation of the ligands impacts on the electronic structure of such catalysts. The recorded effects are dramatic: formal replacement of only one out of eight O-atoms surrounding the metal centers in a dirhodium tetracarboxylate by an N-atom results in a shielding of the corresponding Rh-site of no less than 1000 ppm. The current paper provides the theoretical framework that allows this and related experimental observations made with a set of 19 representative rhodium complexes to be interpreted. In line with symmetry considerations, it is shown that the shielding tensor responds only to the donor ability of the equatorial ligands along the perpendicular principal axis. Axial ligands, in contrast, have no direct effect on shielding but may come into play via the electronicc i s ${cis}$ -effect that they exert onto the neighboring equatorial sites. On top of these fundamental interactions, charge redistribution within the core as well as the electronict r a n s ${trans}$ -effect of ligands of different donor strengths is reflected in the recorded 103 Rh NMR shifts.

Synthesis and Characterization of Azido- and Nitratoalkyl Nitropyrazoles as Potential Melt-Cast Explosives

Desirable advancements in the field of explosive materials include the development of novel melt-castable compounds with melting points ranging from 80 to 110 °C. This is particularly important due to the limited performance and high toxicity associated with TNT (trinitrotoluene). In this study, a series of innovative melt-castable explosives featuring nitratoalkyl and azidoalkyl functionalities attached to the 3-nitro-, 4-nitro-, 3,4-dinitropyrazole, or 3-azido-4-nitropyrazole scaffold are introduced. These compounds were synthesized using straightforward methods and thoroughly characterized using various analytical techniques, including single-crystal X-ray diffraction, IR spectroscopy, multinuclear nuclear magnetic resonance (NMR) spectroscopy, mass spectrometry, elemental analysis, and DTA. Furthermore, the energetic properties such as (theoretical) performance data, sensitivities, and compatibilities of the compounds were evaluated and compared among the different structures.

Solid-State 19F NMR Chemical Shift in Square-Planar Nickel-Fluoride Complexes Linked by Halogen Bonds

The halogen bond (XB) is a highly directional class of noncovalent interactions widely explored by experimental and computational studies. However, the NMR signature of the XB has attracted limited attention. The prediction and analysis of the solid-state NMR (SSNMR) chemical shift tensor provide useful strategies to better understand XB interactions. In this work, we employ a computational protocol for modeling and analyzing the 19F SSNMR chemical shifts previously measured in a family of square-planar trans NiII-L2-iodoaryl-fluoride (L = PEt3) complexes capable of forming self-complementary networks held by a NiF···I(C) halogen bond [Thangavadivale, V.; Chem. Sci. 2018, 9, 3767-3781]. To understand how the 19F NMR resonances of the nickel-bonded fluoride are affected by the XB, we investigate the origin of the shielding in trans-[NiF(2,3,5,6-C6F4I)(PEt3)2], trans-[NiF(2,3,4,5-C6F4I)(PEt3)2], and trans-[NiF(C6F5)(PEt3)2] in the solid state, where a XB is present in the two former systems but not in the last. We perform the 19F NMR chemical shift calculations both in periodic and molecular models. The results show that the crystal packing has little influence on the NMR signatures of the XB, and the NMR can be modeled successfully with a pair of molecules interacting via the XB. Thus, the observed difference in chemical shift between solid-state and solution NMR can be essentially attributed to the XB interaction. The very high shielding of the fluoride and its driving contributor, the most shielded component of the chemical shift tensor, are well reproduced at the 2c-ZORA level. Analysis of the factors controlling the shielding shows how the highest occupied Ni/F orbitals shield the fluoride in the directions perpendicular to the Ni-F bond and specifically perpendicular to the coordination plane. This shielding arises from the magnetic coupling of the Ni(3d)/F(2p lone pair) orbitals with the vacant σNi-F* orbital, thereby rationalizing the very highly upfield (shielded) resonance of the component (δ33) along this direction. We show that these features are characteristic of square-planar nickel-fluoride complexes. The deshielding of the fluoride in the halogen-bonded systems is attributed to an increase in the energy gap between the occupied and vacant orbitals that are mostly responsible for the paramagnetic terms, notably along the most shielded direction.

Effects of surface acidity on the structure of organometallics supported on oxide surfaces

Well-defined organometallics supported on high surface area oxides are promising heterogeneous catalysts. An important design factor in these materials is how the metal interacts with the functionalities on an oxide support, commonly anionic X-type ligands derived from the reaction of an organometallic M-R with an -OH site on the oxide. The metal can either form a covalent M-O bond or form an electrostatic M⁺⋯^-O ion-pair, which impacts how well-defined organometallics will interact with substrates in catalytic reactions. A less common reaction pathway involves the reaction of a Lewis site on the oxide with the organometallic, resulting in abstraction to form an ion-pair, which is relevant to industrial olefin polymerization catalysts. This Feature Article views the spectrum of reactivity between an organometallic and an oxide through the prism of Brønsted and/or Lewis acidity of surface sites and draws analogies to the molecular frame where Lewis and Brønsted acids are known to form reactive ion-pairs. Applications of the well-defined sites developed in this article are also discussed.

Relativistic Density Functional NMR Tensors Analyzed with Spin-free Localized Molecular Orbitals

The implementation of fast relativistic methods based on density functional theory, in conjunction with localized molecular orbital (LMO) based analysis, allows straightforward interpretations of NMR parameters in terms of contributions from core shells, lone pairs, and bonds, for compounds containing elements from across the periodic table. We present a conceptual review of a frequently used LMO analysis of NMR parameters calculated in the presence of spin-orbit interactions and other relativistic effects. An accompanying example focuses on the ¹⁵ N shielding in a heavy metal complex.

Computational 1 H and 13 C NMR in structural and stereochemical studies

Present review outlines the advances and perspectives of computational ¹ H and ¹³ C NMR applied to the stereochemical studies of inorganic, organic, and bioorganic compounds, involving in particular natural products, carbohydrates, and carbonium ions. The first part of the review briefly outlines theoretical background of the modern computational methods applied to the calculation of chemical shifts and spin-spin coupling constants at the DFT and the non-empirical levels. The second part of the review deals with the achievements of the computational ¹ H and ¹³ C NMR in the stereochemical investigation of a variety of inorganic, organic, and bioorganic compounds, providing in an abridged form the material partly discussed by the author in a series of parent reviews. Major attention is focused herewith on the publications of the recent years, which were not reviewed elsewhere.

What Distinguishes the Strength and the Effect of a Lewis Acid: Analysis of the Gutmann-Beckett Method

IUPAC defines Lewis acidity as the thermodynamic tendency for Lewis pair formation. This strength property was recently specified as global Lewis acidity (gLA), and is gauged for example by the fluoride ion affinity. Experimentally, Lewis acidity is usually evaluated by the effect on a bound molecule, such as the induced 31 P NMR shift of triethylphosphine oxide in the Gutmann-Beckett (GB) method. This type of scaling was called effective Lewis acidity (eLA). Unfortunately, gLA and eLA often correlate poorly, but a reason for this is unknown. Hence, the strength and the effect of a Lewis acid are two distinct properties, but they are often granted interchangeably. The present work analyzes thermodynamic, NMR specific, and London dispersion effects on GB numbers for 130 Lewis acids by theory and experiment. The deformation energy of a Lewis acid is identified as the prime cause for the critical deviation between gLA and eLA but its correction allows a unification for the first time.

Competitive tetrel bond and hydrogen bond in benzaldehyde-CO2: characterization via rotational spectroscopy

The rotational spectrum of the 1 : 1 benzaldehyde-CO₂ complex has been investigated using pulsed-jet Fourier transform microwave spectroscopy complemented with quantum chemical calculations. Two isomers, both characterized by one C⋯O tetrel bond (n → π* interaction) and one C-H⋯O hydrogen bond (n → σ* interaction), have been observed in the pulsed jet. Competition between the tetrel bond and the hydrogen bond has been disclosed by natural bond orbital analysis: isomer I is characterized by one dominating OC_CO2⋯O tetrel bond (12.6 kJ mol^-1) and a secondary (C-H)_formyl⋯O hydrogen bond (2.2 kJ mol^-1); by contrast, in isomer II the (C-H)_phenyl⋯O hydrogen bond (7.6 kJ mol^-1) becomes the dominant bond, while the OC_CO2⋯O tetrel bond (5.8 kJ mol^-1) becomes much weaker with respect to that of isomer I. Using intensity measurements the relative population ratio of the two isomers was estimated to be N_I/N_II ≈ 2/1.

Synthesis and electronic structure analysis of the actinide allenylidenes, [{(NR2)3}An(CCCPh2)]- (An = U, Th; R = SiMe3)

The reaction of [AnCl(NR₂)₃] (An = U, Th, R = SiMe₃) with in situ generated lithium-3,3-diphenylcyclopropene results in the formation of [{(NR₂)₃}An(CH[double bond, length as m-dash]C[double bond, length as m-dash]CPh₂)] (An = U, 1; Th, 2) in good yields after work-up. Deprotonation of 1 or 2 with LDA/2.2.2-cryptand results in formation of the anionic allenylidenes, [Li(2.2.2-cryptand)][{(NR₂)₃}An(CCCPh₂)] (An = U, 3; Th, 4). The calculated ¹³C NMR chemical shifts of the C_α, C_β, and C_γ nuclei in 2 and 4 nicely reproduce the experimentally assigned order, and exhibit a characteristic spin-orbit induced downfield shift at C_α due to involvement of the 5f orbitals in the Th-C bonds. Additionally, the bonding analyses for 3 and 4 show a delocalized multi-center character of the ligand π orbitals involving the actinide. While a single-triple-single-bond resonance structure (e.g., An-C[triple bond, length as m-dash]C-CPh₂) predominates, the An[double bond, length as m-dash]C[double bond, length as m-dash]C[double bond, length as m-dash]CPh₂ resonance form contributes, as well, more so for 3 than for 4.

Solid-state 11B NMR studies of coinage metal complexes containing a phosphine substituted diboraanthracene ligand

Transition metal interactions with Lewis acids (M → Z linkages) are fundamentally interesting and practically important. The most common Z-type ligands contain boron, which contains an NMR active ¹¹B nucleus. We measured solid-state ¹¹B{¹H} NMR spectra of copper, silver, and gold complexes containing a phosphine substituted 9,10-diboraanthracene ligand (B₂P₂) that contain planar boron centers and weak M → BR₃ linkages ([(B₂P₂)M][BAr^F₄] (M = Cu (1), Ag (2), Au (3)) characterized by large quadrupolar coupling (C_Q) values (4.4-4.7 MHz) and large span (Ω) values (93-139 ppm). However, the solid-state ¹¹B{¹H} NMR spectrum of K[Au(B₂P₂)]^- (4), which contains tetrahedral borons, is narrow and characterized by small C_Q and Ω values. DFT analysis of 1-4 shows that C_Q and Ω are expected to be large for planar boron environments and small for tetrahedral boron, and that the presence of a M → BR₃ linkage relates to the reduction in C_Q and ¹¹B NMR shielding properties. Thus solid-state ¹¹B NMR spectroscopy contains valuable information about M → BR₃ linkages in complexes containing the B₂P₂ ligand.

Aromaticity of ortho and meta 8-Cycloparaphenylene and Their Dications: Induced Magnetic Field Analysis with Localized and Delocalized Orbitals in Strained Nanohoops

Dications of cycloparaphenyles ([n]CPPs) are known to exhibit in-plane global aromaticity, contained in a nanobelt structure. Recently synthesized ortho and meta isomers of [n]CPPs break the radial symmetry of π structure incorporating perpendicular oriented π orbitals. Herein we set to explore the aromaticity of neutral and dicationic ortho and meta isomers of [8]CPP by dissecting the induced magnetic field to contributions of the twofold radial/perpendicular π system using delocalized canonical molecular orbitals (CMO), and introducing the natural localized molecular orbitals (NLMO) analysis with DFT methods. The dications sustain a reduced global aromatic character of the radial π system under a perpendicular orientation of the external field which declines from ortho to meta isomer and reinforces local aromaticity of ortho ring while it destroys aromaticity of meta ring. Aromaticity variations are determined by symmetry governed rotational excitations of frontier π orbitals. The parallel orientation reveals a substantial reduction of local aromaticity verified with NICS_π analysis and electron delocalization indices.

Tetryl-Tetrylene Addition to Phenylacetylene

Phenylacetylene adds [Ar*GeH2 -SnAr'], [Ar*GeH2 -PbAr'] and [Ar'SnH2 -PbAr*] at rt in a regioselective and stereoselective reaction. The highest reactivity was found for the stannylene, which reacts immediately upon addition of one equivalent of alkyne. However, the plumbylenes exhibit addition to the alkyne only in reaction with an excess of phenylacetylene. The product of the germylplumbylene addition reacts with a second equivalent of alkyne and the product of a CH-activation, a dimeric lead acetylide, were isolated. In the case of the stannylplumbylene the trans-addition product was characterized as the kinetically controlled product which isomerizes at rt to yield the cis-addition product, which is stabilized by an intramolecular Sn-H-Pb interaction. NMR chemical shifts of the olefins were investigated using two- and four-component relativistic DFT calculations, as spin-orbit effects can be large. Hydride abstraction was carried out by treating [Ar'SnPhC=CHGeH2 Ar*] with the trityl salt [Ph3 C][Al(OC{CF3 })4 ] to yield a four membered ring cation.

Inverse halogen dependence in anion 13C NMR

Halogens cause pronounced and systematic effects on the 13C NMR chemical shift (δ13C) of an adjacent carbon nucleus, usually leading to a decrease in the values across the halogen series. Although this normal halogen dependence (NHD) is known in organic and inorganic compounds containing the carbon atom in its neutral and cationic forms, information about carbanions is scarce. To understand how δ13C changes in molecules with different charges, the shielding mechanisms of CHX3, CX3+, and CX3- (X = Cl, Br, or I) systems are investigated via density functional theory calculations and further analyzed by decomposition into contributions of natural localized molecular orbitals. An inverse halogen dependence (IHD) is determined for the anion series as a result of the negative spin-orbit contribution instead of scalar paramagnetic effects. The presence of a carbon nonbonding orbital in anions allows magnetic couplings that generate a deshielding effect on the nucleus and contradicts the classical association between δ13C and atomic charge.

31P Chemical Shifts in Ru(II) Phosphine Complexes. A Computational Study of the Influence of the Coordination Sphere

The NMR chemical shift has been the most versatile marker of chemical structures, by reflecting global and local electronic structures, and is very sensitive to any change within the chemical species. In this work, Ru(II) complexes with the same five ligands and a variable sixth ligand L (none, H₂O, H₂S, CH₃SH, H₂, N₂, N₂O, NO⁺, C═CHPh, and CO) are studied by using as the NMR reporter the phosphorus P_A of a coordinated bidentate P_A-N ligand (P_A-N = o-diphenylphosphino-N,N'-dimethylaniline). The chemical shift of P_A in RuCl₂(P_A-N)(PR₃)(L) (R = phenyl, p-tolyl, or p-FC₆H₄) was shown to increase as the Ru-P_A bond distance decreases, an observation that was not rationalized. This work, using density functional theory (DFT) calculations, reproduces reasonably well the observed ³¹P chemical shifts for these complexes and the correlation between the shifts and the Ru-P_A bond distance as L varies. An interpretation of this correlation is proposed by using a natural chemical shift (NCS) analysis based on the natural bonding orbital (NBO) method. This analysis of the principal components of the chemical shift tensors shows how the σ-donating properties of L have a particularly high influence on the phosphine chemical shifts.

Electron Attachment Studies with the Potential Radiosensitizer 2-Nitrofuran

Nitrofurans belong to the class of drugs typically used as antibiotics or antimicrobials. The defining structural component is a furan ring with a nitro group attached. In the present investigation, electron attachment to 2-nitrofuran (C₄H₃NO₃), which is considered as a potential radiosensitizer candidate for application in radiotherapy, has been studied in a crossed electron-molecular beams experiment. The present results indicate that low-energy electrons with kinetic energies of about 0-12 eV effectively decompose the molecule. In total, twelve fragment anions were detected within the detection limit of the apparatus, as well as the parent anion of 2-nitrofuran. One major resonance region of ≈0-5 eV is observed in which the most abundant anions NO₂^-, C₄H₃O^-, and C₄H₃NO₃^- are detected. The experimental results are supported by ab initio calculations of electronic states in the resulting anion, thermochemical thresholds, connectivity between electronic states of the anion, and reactivity analysis in the hot ground state.

Reactivity of Substituted Benzenes toward Oxidative Addition Relates to NMR Chemical Shift of the Ipso-Carbon

The oxidative addition of benzene derivatives to Pd⁰ catalysts is a key step in cross-coupling reactions. In this work, we show that the ipso-carbon chemical shift of substituted benzenes, and in particular the δ₂₂ component of the chemical shift tensor, correlates with the free energy barrier for oxidative addition. This correlation is traced back to the electron density in the p_z orbital of the ipso-carbon (perpendicular of the ring-plane), with high electron densities favoring oxidative addition. The correlation between chemical shift and free energy barrier holds true for a variety of substituted benzenes, making chemical shift a useful descriptor for predicting the reactivity of aromatic substrates in oxidative addition.

Counterintuitive deshielding on the 13 C NMR chemical shift for the trifluoromethyl anion

The trifluoromethyl anion (CF₃ ^- ) displays ¹³ C NMR chemical shift (175.0 ppm) surprisingly larger than neutral (CHF₃ , 122.2 ppm) and cation (CF₃ ⁺ , 150.7 ppm) compounds. This unexpected deshielding effect for a carbanion is investigated by density functional theory calculations and decomposition analyses of the ¹³ C shielding tensor into localized molecular orbital contributions. The present work determines the shielding mechanisms involved in the observed behaviour of the fluorinated anion species, shedding light on the experimental NMR data and demystify the classical correlation between electron density and NMR chemical shift. The presence of fluorine atoms induces the carbon lone pair to create a paramagnetic shielding on the carbon nucleus.

Cp2Ti(κ2-tBuNCNtBu): A Complex with an Unusual κ2 Coordination Mode of a Heterocumulene Featuring a Free Carbene

Although there are myriad binding modes of heterocumulenes to metal centers, the monometallic κ2-ECE (E = O, S, NR) coordination mode has not been reported. Herein, the synthesis, isolation, and physical characterization of Cp2Ti(κ2-tBuNCNtBu) (3) (Cp = cyclopentadienyl, tBu = tert-butyl), a strained 4-membered metallacycle bearing a free carbene, is described. Computational (DFT, CASSCF, QT-AIM, ELF) and solid-state CP-MAS 13C NMR spectroscopic analysis indicate that 3 is best described as a free carbene with partial Ti-Cβ bonding that results from Ti-N π-bonding mixing with N-C-N σ-bonding of the bent N-C-N framework. Reactivity studies of 3 corroborate its carbene-like nature: protonation with [LutH]I results in the formation of a Ti-formamidinate (4), while oxidation with S8 yields a Ti-thioureate (5). Additionally, a related bridged dititanamidocarbene, (Cp2Ti)2(μ-η1,η1-CyNCNCy) (10) (Cy = cyclohexyl) is reported. Taken together, this work suggests that the 2-electron reduction of heterocumulene moieties can allow access to unusual free carbene coordination geometries given the proper stabilizing coordination environment from the reducing transition metal fragment.

GIAO versus GIPAW: Comparison of Methods To Calculate 11B NMR Shifts of Icosahedral Closo-Heteroboranes toward Boron-Rich Borides

In this work, we perform first-principle density functional theory calculations with the Perdew-Burke-Ernzerhof (PBE) exchange correlation functional to compare the results of the gauge-including atomic orbital (GIAO) method with the gauge-including projector-augmented wave (GIPAW) approach for isotropic ¹¹B nuclear magnetic resonance shifts. GIPAW had been used successfully for the theoretical calculation of nuclear magnetic parameters of ¹¹B species in strong ionic solid-phase compounds such as borates but had been applied very rarely to structures where boron is mainly involved in complex covalent bonding situations, for example, in icosahedra of boron-rich borides. Thus, we investigate the accuracy of both well-known methods and reliability of the effective treatment of core electrons on a test set containing 16 experimentally known closo-(hetero)dodecaboranes. In general, we find very good agreement between GIAO and GIPAW when compared to experimental observations. However, accidental degeneracies of the shift values are better predicted by GIPAW. The optimized molecular geometries on the PBE level agree well with gaseous electron diffraction data and lead to theoretical isotropic chemical ¹¹B shifts with root-mean-square errors of 2.1 and 1.0 ppm depending on the used model of converting absolute shieldings to chemical shifts. The comparison with results from hybrid functionals (B3LYP, B3LYP-D2, and PBE0) shows a minor improvement in accuracy, which is in agreement with ¹³C shifts of sp³-hybridized species. In order to prove the reliability of the conversion parameters obtained by PBE, we report the calculated ¹¹B shifts of 1,2-, 1,7-, and 1,12-PCB₁₀H₁₁ with GIAO and GIPAW to our knowledge for the first time. Additionally, Bader's analysis is carried out on the converged electron density for all boron species within the molecular test set, yielding no simple direct relation between charge and isotropic shifts.

Understanding 125Te NMR chemical shifts in disymmetric organo-telluride compounds from natural chemical shift analysis

Magnetic coupling of the lone pair: theoretical investigations reveal the origin of ¹²⁵Te chemical shift in disymmetric organotellurides

Origin of the 29Si NMR chemical shift in R3Si–X and relationship to the formation of silylium (R3Si+) ions

The origin in deshielding of ²⁹Si NMR chemical shifts in R₃Si–X, where X = H, OMe, Cl, OTf, [CH₆B₁₁X₆], toluene, and O_X (O_X = surface oxygen), as well as ⁱPr₃Si⁺ and Mes₃Si⁺ were studied using DFT methods.

Carbon-13 NMR Chemical Shift: A Descriptor for Electronic Structure and Reactivity of Organometallic Compounds

Metal-bonded carbon atoms in metal-alkyl, metal-carbene/alkylidene, and metal-carbyne/alkylidyne species often show significantly more deshielded isotropic chemical shifts than their organic counterparts (alkanes, alkenes, and alkynes). While isotropic chemical shift is universally used to characterize a chemical compound in solution, it is an average value of the three principal components of the chemical shift tensor (δ₁₁ > δ₂₂ > δ₃₃). The tensor components, which are accessible by solid-state NMR spectroscopy, can provide detailed information about the electronic structure (frontier molecular orbitals) at the observed nuclei. This information can be accessed in detail by quantum chemical calculations, most notably by an analysis of the paramagnetic contribution to the NMR shielding tensor. The paramagnetic term mainly results from the coupling of occupied and empty molecular orbitals close in energy-the frontier molecular orbitals-under the effect of the external magnetic field (B₀). In organometallic compounds, a large deshielding of the isotropic carbon-13 chemical shift of the metal-bonded carbon atom is commonly related to the coupling between the occupied σ_M-C orbital and low-lying vacant orbitals of π_M═C^* character. The deshielding at the α-carbon hence probes the extent of σ_M-C and π_M═C^* interactions. This molecular orbital view readily explains the strong deshielding and large anisotropy (evidenced by the span Ω = δ₁₁ - δ₃₃) observed in metal alkylidenes and alkylidynes (200 < δ_iso < 400 ppm). Fischer carbenes are generally more deshielded than Schrock or Grubbs alkylidenes due to their low-lying π_M═C^* orbital. Chemical shift hence shows their higher electrophilic character, connecting NMR spectroscopy to reactivity patterns. Similarly, the α-carbon of metal-alkyls display deshielded chemical shifts in specific coordination environments. This deshielding, which is often prominently pronounced for cationic species, indicates the presence of partial π-bond character in the metal-carbon bond, making these bonds topologically equivalent to alkylidene π-bonds. The π-character in metal-alkyl bonds favors (i) α-H abstraction processes in metal bis-alkyl compounds yielding metal alkylidenes, (ii) [2 + 2]-retrocyclization of metallacyclobutanes that participate in olefin metathesis, (iii) olefin insertion in cationic metal alkyls thus explaining polymerization activity trends and the importance of α-H agostic interactions, and (iv) C-H bond activation on metal-alkyls via σ-bond metathesis. The presence of π-character in the metal-carbon bonds involved in these processes rationalizes the parallel reactivity patterns of metal-alkyls toward olefin insertion and σ-bond metathesis and the fact that σ-bond metathesis, olefin insertion, and olefin metathesis are commonly observed with metal atoms in the same ligand field. Because of the similarities in the frontier molecular orbitals involved in these processes, these reactions can be viewed as isolobal. This explains why certain fragments, such as bent metallocenes (d⁰ Cp₂M) or T-shaped L₃M, are ubiquitous in these reactions.

Computational protocols for calculating 13C NMR chemical shifts

The most recent results dealing with the computation of ¹³C NMR chemical shifts in chemistry (small molecules, saturated, unsaturated and aromatic compounds, heterocycles, functional derivatives, coordination complexes, carbocations, and natural products) are reviewed, paying special attention to theoretical background and accuracy, the latter involving solvent effects, vibrational corrections, and relativistic effects.

The 125Te Chemical Shift of Diphenyl Ditelluride: Chasing Conformers over a Flat Energy Surface

The interest in diphenyl ditelluride (Ph₂Te₂) is related to its strict analogy to diphenyl diselenide (Ph₂Se₂), whose capacity to reduce organic peroxides is largely exploited in catalysis and green chemistry. Since the latter is also a promising candidate as an antioxidant drug and mimic of the ubiquitous enzyme glutathione peroxidase (GPx), the use of organotellurides in medicinal chemistry is gaining importance, despite the fact that tellurium has no recognized biological role and its toxicity must be cautiously pondered. Both Ph₂Se₂ and Ph₂Te₂ exhibit significant conformational freedom due to the softness of the inter-chalcogen and carbon–chalcogen bonds, preventing the existence of a unique structure in solution. Therefore, the accurate calculation of the NMR chemical shifts of these flexible molecules is not trivial. In this study, a detailed structural analysis of Ph₂Te₂ is carried out using a computational approach combining classical molecular dynamics and relativistic density functional theory methods. The goal is to establish how structural changes affect the electronic structure of diphenyl ditelluride, particularly the ¹²⁵Te chemical shift.

Substituent effects in the so-called cationπ interaction of benzene and its boron-nitrogen doped analogues: overlooked role of σ-skeleton

Despite massive efforts to pinpoint the substituent effects in the so-called cationπ systems, no consensus has been yet reached on how substituents exercise their effects in the interaction of the aromatic molecule with the metal ion. The π-polarization (the Hunter model) and the direct local effect (the Wheeler-Houk model) are two lines of thought applied to this problem, but the justification of both approaches is based on insufficiently proven assumptions and approximations. In order to shed more light on this issue we propose a new approach which enables us to gauge directly the energetic trends resulting from the interaction of the ring with the cation. In our method we add one more partitioning level to the interacting quantum atoms (IQA) scheme and decompose the IQA interaction energies into contributions resulting from σ and π electron densities of the aromatic ring. The new approach, which is named partitioned-IQA, abbreviated as p-IQA, has been applied to complexes of derivatives of benzene or azaborine interacting with a sodium cation. The p-IQA approach reveals that in these systems both σ and π electronic moieties are polarized. Interestingly, for the majority of cases the σ-polarization outweighs the π one, contrary to the Hunter model. However, the Wheeler-Houk model is not precise, either, since the σ-polarization shows some degree of non-locality. In addition, the substituents are found to have a negligible influence on the ring orbital-overlapping capability, i.e. the covalency. Therefore, the substituent effect in the cationπ interaction is a nonlocal classical effect, indicating that neither Hunter model nor Wheeler-Houk model is able to fully describe all the aspects of the substituent effects. The p-IQA conclusions for the considered systems have been compared with the results from the functional-group SAPT (F-SAPT) method. We believe that the presented partitioning in the IQA framework will provide a deeper insight into the substituent effects in the cationπ interactions, which is beyond the σ-π atomic charge population separation.

"Through-Space" Relativistic Effects on NMR Chemical Shifts of Pyridinium Halide Ionic Liquids

We have investigated, using two-component relativistic density functional theory (DFT) at ZORA-SO-BP86 and ZORA-SO-PBE0 level, the occurrence of relativistic effects on the ¹ H, ¹³ C, and ¹⁵ N NMR chemical shifts of 1-methylpyridinium halides [MP][X] and 1-butyl-3-methylpyridinium trihalides [BMP][X₃ ] ionic liquids (ILs) (X=Cl, Br, I) as a result of a non-covalent interaction with the heavy anions. Our results indicate a sizeable deshielding effect in ion pairs when the anion is I^- and I₃ ^- . A smaller, though nonzero, effect is observed also with bromine while chlorine based anions do not produce an appreciable relativistic shift. The chemical shift of the carbon atoms of the aromatic ring shows an inverse halogen dependence that has been rationalized based on the little C-2s orbital contribution to the σ-type interaction between the cation and anion. This is the first detailed account and systematic theoretical investigation of a relativistic heavy atom effect on the NMR chemical shifts of light atoms in the absence of covalent bonds. Our work paves the way and suggests the direction for an experimental investigation of such elusive signatures of ion pairing in ILs.

Double aromaticity of the B40 fullerene: induced magnetic field analysis of π and σ delocalization in the boron cavernous structure

B₄₀ enables the formation of a strong long range shielding response under different orientations, characterizing the spherical aromatic nature of the cavernous D_2d structure, which was dissected to contributions from π, σ and core electrons.

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.

The halogen effect on the 13C NMR chemical shift in substituted benzenes

X (F, Cl, Br, I) and R (NH₂, NO₂) group effects on ¹³C NMR chemical shifts are explained by π and σ orbitals, respectively.