Unraveling the enigma of Craig-type Möbius-aromatic osmium compounds

Nuclear magnetic resonance (NMR) chemical shifts and the magnetically induced current density (MICD) susceptibility of four osmium containing molecules have been calculated at the density functional theory (DFT) level using three relativistic levels of theory. The calculations were performed at the quasi-relativistic level using an effective core potential (ECP) for Os, at the all-electron scalar exact two-component (X2C) relativistic level, and at the relativistic X2C level including spin-orbit coupling (SO-X2C). In earlier studies, the osmapentalene (1) and the osmapentalynes (2 and 3) were considered Craig-type Möbius aromatic and it was suggested that the analogous osmium compound (4) is Craig-type Möbius antiaromatic. Here, the ring-current strengths were obtained with the gauge including magnetically induced currents (GIMIC) method by integrating the MICD susceptibility passing through planes that intersect chemical bonds and by line integration of the induced magnetic field using Ampère-Maxwell's law. The ring-current calculations suggest that 1, 2 and 3 are weakly aromatic and that 4 is nonaromatic. The accuracy of the MICD susceptibility was assessed by comparing calculated NMR chemical shifts to available experimental data. Visualization of the MICD susceptibility shows that the ring current does not pass from one side of the molecular plane to the other, which means that the MICD susceptibility of the studied molecules does not exhibit any Möbius topology as one would expect for Craig-type Möbius aromatic and for Craig-type Möbius antiaromatic molecules. Thus, molecules 1-3 are not Craig-type Möbius aromatic and molecule 4 is not Craig-type Möbius antiaromatic as previously suggested. Calculations of the 1H NMR and 13C NMR chemical shifts of atoms near the Os atom show the importance of including spin-orbit effects. Overall, our study revisits the understanding of the aromaticity of organometallic molecules containing transition metals.

Revisiting Gauge-Independent Kinetic Energy Densities in Meta-GGAs and Local Hybrid Calculations of Magnetizabilities

In a recent study [J. Chem. Theory Comput. 2021, 17, 1457-1468], some of us examined the accuracy of magnetizabilities calculated with density functionals representing the local density approximation (LDA), generalized gradient approximation (GGA), meta-GGA (mGGA), as well as global hybrid (GH) and range-separated (RS) hybrid functionals by assessment against accurate reference values obtained with coupled-cluster theory with singles, doubles, and perturbative triples [CCSD(T)]. Our study was later extended to local hybrid (LH) functionals by Holzer et al. [J. Chem. Theory Comput. 2021, 17, 2928-2947]; in this work, we examine a larger selection of LH functionals, also including range-separated LH (RSLH) functionals and strong-correlation LH (scLH) functionals. Holzer et al. also studied the importance of the physically correct handling of the magnetic gauge dependence of the kinetic energy density (τ) in mGGA calculations by comparing the Maximoff-Scuseria formulation of τ used in our aforementioned study to the more physical current-density extension derived by Dobson. In this work, we also revisit this comparison with a larger selection of mGGA functionals. We find that the newly tested LH, RSLH, and scLH functionals outperform all of the functionals considered in the previous studies. The various LH functionals afford the seven lowest mean absolute errors while also showing remarkably small standard deviations and mean errors. Most strikingly, the best two functionals are scLHs that also perform remarkably well in cases with significant multiconfigurational character, such as the ozone molecule, which is traditionally excluded from statistical error evaluations due to its large errors with common density functionals.

NMR chemical shift of confined 129Xe: coordination number, paramagnetic channels and molecular dynamics in a cryptophane-A biosensor

Advances in hyperpolarisation and indirect detection have enabled the development of xenon nuclear magnetic resonance (NMR) biosensors (XBSs) for molecule-selective sensing in down to picomolar concentration. Cryptophanes (Crs) are popular cages for hosting the Xe "spy". Understanding the microscopic host-guest chemistry has remained a challenge in the XBS field. While early NMR computations of XBSs did not consider the important effects of host dynamics and explicit solvent, here we model the motionally averaged, relativistic NMR chemical shift (CS) of free Xe, Xe in a prototypic CrA cage and Xe in a water-soluble CrA derivative, each in an explicit H2O solvent, over system configurations generated at three different levels of molecular dynamics (MD) simulations. We confirm the "contact-type" character of the Xe CS, arising from the increased availability of paramagnetic channels, magnetic couplings between occupied and virtual orbitals through the short-ranged orbital hyperfine operator, when neighbouring atoms are in contact with Xe. Remarkably, the Xe CS in the present, highly dynamic and conformationally flexible situations is found to depend linearly on the coordination number of the Xe atom. We interpret the high- and low-CS situations in terms of the magnetic absorption spectrum and choose our preference among the used MD methods based on comparison with the experimental CS. We check the role of spin-orbit coupling by comparing with fully relativistic CS calculations. The study outlines the computational workflow required to realistically model the CS of Xe confined in dynamic cavity structures under experimental conditions, and contributes to microscopic understanding of XBSs.

Current-density pathways in figure-eight-shaped octaphyrins

We have calculated the current density induced by an external magnetic field in a set of figure-eight-shaped expanded porphyrinoids. The studied octaphyrins can be divided into three classes (N2, N4, and N6) based on the number of the inner hydrogen atoms of the pyrrole rings. Using the Runge-Kutta method, the current density is split into diatropic and paratropic contributions that are analyzed separately. The calculations show that one common ring current consists of two rather independent pathways. Each of them follows the outer side of the molecular frame of one half of the molecule and passes to the inner side of the frame on the other half. The ring-current pathways are similar to the ones for [12]infinitene. However, the current density of the octaphyrins is more complex having many branching points and pathways. Vertical through-space current-density pathways pass in the middle of the molecules through a plane that is parallel to the figure-eight-shaped view of the molecules when the magnetic field is perpendicular to the plane. The isolectronic N2 and the N4 dication sustain a weak paratropic ring current inside the molecule, which is also observed in the ¹H NMR magnetic shielding constant of the inner hydrogen atoms. The diatropic current-density contribution dominates in the studied molecules. For the N4 and N6 molecules, the global current-density pathways are only diatropic and N6 sustains the strongest global diatropic current-density flux of 13.2 nA T^-1.

Current Density and Spectroscopy—A Themed Issue in Honor of Professor Riccardo Zanasi on the Occasion of His 70th Birthday

It is our great pleasure to introduce the Festschrift of Chemistry to honor Professor Riccardo Zanasi (Figure 1) on the occasion of his 70th birthday and to recognize his important contributions to quantum chemistry, particularly in the field of magnetic response and chiroptical spectroscopies [...]

Integration of global ring currents using the Ampère-Maxwell law

Magnetically induced ring currents are calculated from the magnetic shielding tensor by employing the Ampère-Maxwell law. The feasibility of the method is demonstrated by integrating the zz component of the shielding tensor along the symmetry axis of highly symmetric ring-shaped aromatic, antiaromatic and nonaromatic molecules. The calculated ring-current strengths agree perfectly with the ones obtained by integrating the current-density flux passing through a plane cutting half the molecular ring. The method can be used in combination with all electronic structure codes capable of calculating nuclear magnetic resonance (NMR) shielding tensors in general points in space. We also show that nucleus independent chemical shifts (NICS) along the symmetry axis are related to the spatial derivative of the strength of the global ring-current along the z axis.