Spatial Contributions to Nuclear Magnetic Shieldings

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.

Magnetic Characterization of the Infinitene Molecule

The origin-independent current density induced by a perpendicular magnetic field in the infinitene molecule has been calculated, confirming the recently presented result by Orozco-Ic et al. (Phys. Chem. Chem. Phys.2022, 24, 6404−640935262148) of two disjointed global current pathways along the edges formed by 24 carbon atoms having the form of the infinity symbol. The current strength has been assessed along the C–C bonds forming the two separate circuits, whose particular shape provides a diamagnetic exaltation which is only 73% of the expected value for this aromatic molecule. Through space currents have been found along the bond paths determined by the electron density gradient, whose strength is 10% that of the aromatic benzene ring current. It is shown that the pair of high-field ¹H NMR experimental signals carry the signature of the two global currents, which are counterrotating inside the fjord regions with respect to the rim of the coronene subunits.

The hidden aromaticity in borazine

The aromaticity of borazine (B₃N₃H₆), also known as the inorganic benzene, is a controversial issue since this compound has several characteristics that could qualify it as an aromatic compound. However, recent studies using magnetic criteria indicate that this compound should be considered as a non-aromatic system. This assignment is mainly due to diatropic currents in the nitrogen atoms without observation of ring currents. The present work shows by means of the magnetic criteria that borazine has a ring current hidden by the local contributions of degenerate orbitals π₁ and π₂. Additionally, the study of borazine's first triplet state antiaromaticity using the magnetic and energetic criteria by means of isomerization stabilization energies (ISEs) together with Baird's and Hückel's rules suggests that borazine is best described as an (weakly) aromatic system.

Dissected magnetically induced current density calculations show that local currents in borazine counteract the ring current and therefore hide its (weak) aromatic character.

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.

Current density and molecular magnetic properties

We give an overview of the molecular response to an external magnetic field perturbing quantum mechanical systems. We present state-of-the-art methods for calculating magnetically-induced current-density susceptibilities. We discuss the essence and properties of current-density susceptibilities and how molecular magnetic properties can be calculated from them. We also review the theory of spin-current densities, how relativity affects current densities and magnetic properties. An overview of the magnetic ring-current criterion for aromaticity is given, which has implications on theoretical and experimental research. The recently reported theory of antiaromaticity and how molecular symmetry affects the magnetic response are discussed and applied to closed-shell paramagnetic molecules. The topology of magnetically induced current densities and its consequences for molecular magnetic properties are also presented with twisted and toroidal molecules as examples.

Magnetic shielding paints an accurate and easy-to-visualize portrait of aromaticity

Chemists are trained to recognize aromaticity semi-intuitively, using pictures of resonance structures and Frost-Musulin diagrams, or simple electron-counting rules such as Hückel's 4n + 2/4n rule. To quantify aromaticity one can use various aromaticity indices, each of which is a number reflecting some experimentally measured or calculated molecular property, or some feature of the molecular wavefunction, which often has no visual interpretation or may not have direct chemical relevance. We show that computed isotropic magnetic shielding isosurfaces and contour plots provide a feature-rich picture of aromaticity and chemical bonding which is both quantitative and easy-to-visualize and interpret. These isosurfaces and contour plots make good chemical sense as at atomic positions they are pinned to the nuclear shieldings which are experimentally measurable through chemical shifts. As examples we discuss the archetypal aromatic and antiaromatic molecules of benzene and square cyclobutadiene, followed by modern visual interpretations of Clar's aromatic sextet theory, the aromaticity of corannulene and heteroaromaticity.

Electronic Currents and Anapolar Response Induced in Molecules by Monochromatic Light

It is shown that the electric dipole- and electric quadrupole–anapole polarizabilities, denoted respectively by fαβ′ and gα,βγ′, and the anapole magnetizability aαβ, are intrinsic properties of the electron cloud of molecules responding to optical fields. fαβ′ is a nonvanishing property of chiral and achiral compounds, whereas aαβ is suitable for enantiomer discrimination of chiral species. They can conveniently be evaluated by numerical integration, employing a formulation complementary to that provided by perturbation theory and relying on the preliminary computation of electronic current density tensors all over the molecular domain. The origin dependence of the dynamic anapolar response is rationalized via related computational techniques employing numerical integration, as well as definitions of molecular property tensors, for example, electric dipole and electric quadrupole polarizabilties and magnetizability. A preliminary application of the theory is reported for the Ra enantiomer of the hydrogen peroxide molecule, evaluating tensor components of electric dipole-anapole polarizability and anapole magnetizability as functions of the dihedral angle ϕ≡∠ H-O-O-H in the range 0∘≤ϕ≤180∘.

Spatial Contributions to 1H NMR Chemical Shifts of Free-Base Porphyrinoids

A recently developed methodology for calculating, analyzing, and visualizing nuclear magnetic shielding densities is used for studying spatial contributions including ring-current contributions to 1H nuclear magnetic resonance (NMR) chemical shifts of aromatic and anti-aromatic free-base porphyrinoids. Our approach allows a visual inspection of the spatial origin of the positive (shielding) and negative (deshielding) contributions to the nuclear magnetic shielding constants. Diatropic and paratropic current-density fluxes yield both shielding and deshielding contributions implying that not merely the tropicity of the current density determines whether the contribution has a shielding or deshielding character. Instead the shielding or deshielding contribution is determined by the direction of the current-density flux with respect to the studied nucleus.

Magnetically induced ring currents in naphthalene-fused heteroporphyrinoids

The magnetically induced current density of an intriguing naphthalene-fused heteroporphyrin has been studied, using the quantum-chemical, gauge-including magnetically induced currents (GIMIC) method. The ring-current strengths and current-density pathways for the heteroporphyrin, its Pd complex, and the analogous quinoline-fused heteroporphyrin provide detailed information about their aromatic properties. The three porphyrinoids have similar current-density pathways and are almost as aromatic as free-base porphyrin. Notably, we show that the global ring current makes a branch at three specific points. Thus, the global current is composed of a total of eight pathways that include 22 π-electrons, with no contributions from 18-electron pathways.

Benchmarking Magnetizabilities with Recent Density Functionals

We have assessed the accuracy of the magnetic properties of a set of 51 density functional approximations, including both recently published and already established functionals. The accuracy assessment considers a series of 27 small molecules and is based on comparing the predicted magnetizabilities to literature reference values calculated using coupled-cluster theory with full singles and doubles and perturbative triples [CCSD(T)] employing large basis sets. The most accurate magnetizabilities, defined as the smallest mean absolute error, are obtained with the BHandHLYP functional. Three of the six studied Berkeley functionals and the three range-separated Florida functionals also yield accurate magnetizabilities. Also, some older functionals like CAM-B3LYP, KT1, BHLYP (BHandH), B3LYP, and PBE0 perform rather well. In contrast, unsatisfactory performance is generally obtained with Minnesota functionals, which are therefore not recommended for calculations of magnetically induced current density susceptibilities and related magnetic properties such as magnetizabilities and nuclear magnetic shieldings. We also demonstrate that magnetizabilities can be calculated by numerical integration of magnetizability density; we have implemented this approach as a new feature in the gauge-including magnetically induced current (GIMIC) method. Magnetizabilities can be calculated from magnetically induced current density susceptibilities within this approach even when analytical approaches for magnetizabilities as the second derivative of the energy have not been implemented. The magnetizability density can also be visualized, providing additional information that is not otherwise easily accessible on the spatial origin of magnetizabilities.