Current density tensors

Standard Response Calculation of Electric Dipole Polarizability and Specific Optical Rotation Power Densities

A time-dependent method has been developed to solve the standard response equation for the calculation of dynamic molecular property densities, endowed with the characteristic of being origin-invariant, entirely in the atomic orbital basis at both HF and DFT level of theory. The method has been tuned in particular for the calculation of origin-independent electric dipole polarizability density and specific rotation power density. Some demonstrations are given for the hexabenzocoronene molecule and the Tröger's base.

Electronic Current Density Induced by Uniform Magnetic Fields in Clarenes

Some planar and non-planar clarenes have been studied using maps of magnetically induced quantum-mechanical current density and tools from differential topology to assess their magnetic response in connection with recent results by Du and Wang. Bond current strengths have been computed to estimate quantitative measures. Isosurfaces of the divergence of induced Lorentz force density have been shown to provide useful additional criteria, especially in the case of non-planar clarenes.

On the aromaticity of actinide compounds

The chemistry of actinides has flourished since the late 2010s with the synthesis of new actinide complexes and clusters. On the theoretical side, a range of tools is available for the characterization of these heavy element-containing compounds, but discrepancies in the assessment of aromaticity using different tools have led to controversies. In this Perspective, we examine the origin of controversies relating to the aromaticity of metallic compounds, with a focus on actinides. The aromaticity of actinides is important, not because these molecules are numerous or have a special role in catalysis or reactivity, but because this topic pushes theories of aromaticity to their limits. Owing to its reference independence, the magnetic criterion of aromaticity has been the most popular choice for the characterization of the aromaticity of metallic compounds, including actinide compounds. Through examination of several case studies, we show why this criterion might be misleading for metallic species and explain how findings relating to actinide compounds could reshape theories of aromaticity, not just for actinides but perhaps also for well-known hydrocarbons.

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.

Electronic current densities and origin-independent property densities induced by optical fields

The interaction of a molecule with optical fields is customarily interpreted by means of induced time-dependent electric polarizabilities, magnetizabilities and mixed electric-magnetic polarizabilities. In general, these properties can be rationalized by integrals of density functions formulated in terms of induced charge and current densities. In this perspective, we focus on what has been done so far at the theoretical level, and on what can be expected to be unveiled from the topological study of suitable density functions, endowed with the fundamental requirement of origin invariance. Densities characterized by such a property can be integrated all over the configuration space to obtain electric dipole polarizability and optical rotatory power. Corresponding maps visualize domains mainly involved in the molecular response. The diagonal components of origin-independent density tensor functions that, on integration, yield corresponding electric dipole polarizability tensor of benzene, naphthalene, phenanthrene and ovalene, have been computed, confirming the ubiquitous presence of counter-polarization regions in the proximity of the atomic nuclei. They are associated with toroidal electron currents, induced by time derivative of the electric field of impinging radiation. Electron (de)localization in these systems is readily observed and estimated. The optical rotation density of the carbonyl chromophore is studied in detail. Its essential feature is the separation in quadrants of alternating sign of density about the CO bond. The presence of an extrachromophoric perturbation determines asymmetry in the extension of the quadrant distribution, thus causing optical rotation.

Topological data analysis of vortices in the magnetically-induced current density in LiH molecule

A novel strategy for extracting axial (AV) and toroidal (TV) vortices in the magnetically-induced current density (MICD) in molecular systems is introduced, and its pilot application to LiH molecule is demonstrated. It exploits differences in the topologies of AV and TV cores and involves two key steps: selecting a scalar function that can describe vortex cores in MICD and its subsequent topological analysis. The scalar function of choice is Ω^B_α based on the velocity-gradient Ω method known in research on classical flows. The Topological Data Analysis (TDA) is then used to analyze the Ω^B_α scalar field. In particular, TDA robustly assigns distinct topological features of this field to different vortex types in LiH: AV to saddle-maximum separatrices which connect maxima to 2-saddles located on the domain's boundary, TV to a 1-cycle of the super-level sets of the input data. Both are extracted as the most persistent features of the topologically-simplified Ω^B_α scalar field.

Three-Dimensional Fully π-Conjugated Macrocycles: When 3D-Aromatic and When 2D-Aromatic-in-3D?

Several fully π-conjugated macrocycles with puckered or cage-type structures were recently found to exhibit aromatic character according to both experiments and computations. We examine their electronic structures and put them in relation to 3D-aromatic molecules (e.g., closo-boranes) and to 2D-aromatic polycyclic aromatic hydrocarbons. Using qualitative theory combined with quantum chemical calculations, we find that the macrocycles explored hitherto should be described as 2D-aromatic with three-dimensional molecular structures (abbr. 2D-aromatic-in-3D) and not as truly 3D-aromatic. 3D-aromatic molecules have highly symmetric structures (or nearly so), leading to (at least) triply degenerate molecular orbitals, and for tetrahedral or octahedral molecules, an aromatic closed-shell electronic structure with 6n + 2 electrons. Conversely, 2D-aromatic-in-3D structures exhibit aromaticity that results from the fulfillment of Hückel’s 4n + 2 rule for each macrocyclic path, yet their π-electron counts are coincidentally 6n + 2 numbers for macrocycles with three tethers of equal lengths. It is notable that 2D-aromatic-in-3D macrocyclic cages can be aromatic with tethers of different lengths, i.e., with π-electron counts different from 6n + 2, and they are related to naphthalene. Finally, we identify tetrahedral and cubic π-conjugated molecules that fulfill the 6n + 2 rule and exhibit significant electron delocalization. Yet, their properties resemble those of analogous compounds with electron counts that differ from 6n + 2. Thus, despite the fact that these molecules show substantial π-electron delocalization, they cannot be classified as true 3D-aromatics.

Persistent Planar Tetracoordinate Carbon in Global Minima Structures of Silicon-Carbon Clusters

Recently, we reported a series of global minima whose structures consist of carbon rings decorated with heavier group 14 elements. Interestingly, these structures feature planar tetracoordinate carbons (ptCs) and result from the replacement of five or six protons (H+) from the cyclopentadienyl anion (C5H5−) or the pentalene dianion (C8H62−) by three or four E2+ dications (E = Si–Pb), respectively. The silicon derivatives of these series are the Si3C5 and Si4C8 clusters. Here we show that ptC persists in some clusters with an equivalent number of C and Si atoms, i.e., Si5C5, Si8C8, and Si9C9. In all these species, the ptC is embedded in a pentagonal C5 ring and participates in a three-center, two-electron (3c-2e) Si-ptC-Si σ-bond. Furthermore, these clusters are π-aromatic species according to chemical bonding analysis and magnetic criteria.

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.

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.

Cubic aromaticity in ligand-stabilized doped Au superatoms

The magnetic response of valence electrons in doped gold-based M@Au₈L₈ ^q superatoms (M = Pd, Pt, Ag, Au, Cd, Hg, Ir, and Rh; L = PPh₃; and q = 0, +1, +2) is studied by calculating the gauge including magnetically induced currents (GIMIC) in the framework of the auxiliary density functional theory. The studied systems include 24 different combinations of the dopant, total cluster charge, and cluster structure (cubic-like or oblate). The magnetically induced currents (both diatropic and paratropic) are shown to be sensitive to the atomic structure of clusters, the number of superatomic electrons, and the chemical nature of the dopant metal. Among the cubic-like structures, the strongest aromaticity is observed in Pd- and Pt-doped M@Au₈L₈ ⁰ clusters. Interestingly, Pd- and Pt-doping increases the aromaticity as compared to a similar all-gold eight-electron system Au₉L₈ ⁺¹. With the recent implementation of the GIMIC in the deMon2k code, we investigated the aromaticity in the cubic and butterfly-like M@Au₈ core structures, doped with a single M atom from periods 5 and 6 of groups IX-XII. Surprisingly, the doping with Pd and Pt in the cubic structure increases the aromaticity compared to the pure Au case not only near the central atom but encompassing the whole metallic core, following the aromatic trend Pd > Pt > Au. These doped (Pd, Pt)@Au₈ nanoclusters show a closed shell 1S²1P⁶ superatom electronic structure corresponding to the cubic aromaticity rule 6n + 2.

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.

Spatial Contributions to Nuclear Magnetic Shieldings

We develop a methodology for calculating, analyzing, and visualizing nuclear magnetic shielding densities which are calculated from the current density via the Biot-Savart relation. Atomic contributions to nuclear magnetic shielding constants can be estimated within our framework with a Becke partitioning scheme. The new features have been implemented in the GIMIC program and are applied in this work to the study of the 1H and 13C nuclear magnetic shieldings in benzene (C6H6) and cyclobutadiene (C4H4). The new methodology allows a visual inspection of the spatial origins of the positive (shielding) and negative (deshielding) contributions to the nuclear magnetic shielding constant of a single nucleus, something which has not been hitherto easily accomplished. Analysis of the shielding densities shows that diatropic and paratropic current-density fluxes yield both shielding and deshielding contributions, as the shielding or deshielding is determined by the direction of the current-density flux with respect to the studied nucleus instead of the tropicity. Becke partitioning of the magnetic shieldings shows that the magnetic shielding contributions mainly originate from the studied atom and its nearest neighbors, confirming the localized character of nuclear magnetic shieldings.

Relating nucleus independent chemical shifts with integrated current density strengths

Indices based on the nucleus independent chemical shift (NICS) are the most frequently used in analysis of magnetic aromaticity. The magnetically induced current density, on the other hand, is a key concept in defining magnetic aromaticity. The integrated current strength (current strength susceptibility) was found to be a very useful tool in aromaticity studies. There is widely accepted notion that the properly chosen NICS-based index can provide information on the current density strength and direction in a molecule of interest. In this work, a detailed numerical testing of the relationship between the integrated bond current strength and the most employed NICS indices was performed for a set of 43 monocyclic aromatic molecules. Based on the statistical data analysis, the relationship between the bond current strength and its π and σ electron components, on one side, and the isotropic NICS (NICSiso and NICSπ,iso) and zz-component of the NICS tensor (NICSzz and NICSπ,zz), on the other side, was examined. It was found that between the NICSπ,zz(1) and π-electron bond current strenghts there is very good linear correlation. Quite surprisingly, it was revealed that the NICSiso(1) and NICSzz(1) are not correlated with the π electron bond current strengths. On the other hand, a reasonably good linear correlation was found between the NICSzz(1) and total bond current strengths.

Cubic magnetic response of diamagnetic molecules via third-order electronic current density

An approach to nonlinear magnetic properties of a diamagnetic molecule in the presence of a strong homogeneous magnetic field B is proposed, introducing the notion of third-order electronic current density and associated fourth-rank current density tensors via a perturbation expansion. It is shown that cubic magnetic response properties, e.g., fourth-rank magnetizability and fourth-rank nuclear magnetic shielding, can be expressed by relationships alternative to those arrived through perturbation theory by means of third-order current density. A sum rule, providing a condition in integral form for electron charge conservation to the third order in B, has been obtained.

Program Package for the Calculation of Origin-Independent Electron Current Density and Derived Magnetic Properties in Molecular Systems

We present SYSMOIC, a program package for the calculation of the origin-independent current density induced at first order by an external magnetic field in planar and nonplanar molecular systems. Origin independence is obtained adopting the continuous transformation of the origin of the current density method, implemented at both density functional theory (DFT) and Hartree-Fock (HF) levels. Expansion coefficients for perturbed and unperturbed molecular orbitals, over basis sets containing up to m-type Gaussian functions, can be calculated by the package itself or obtained from a Gaussian calculation. A number of different functionalities presented so far in the literature that are connected to the induced current, such as current density maps for any orientation of the inducing magnetic field, net bond current strengths, stagnation graphs, magnetic shielding densities, vorticities, and anisotropies, are now made available all together in a single multiplatform package installation.

Tests of accuracy for computed magnetic properties via off-diagonal hypervirial relations

Most of the methods presently available to investigate the molecular magnetic response work extremely well for the computation of properties, such as magnetizability and nuclear magnetic shielding, but they provide insufficiently accurate current density maps, in that they do not guarantee exact conservation, leading to unphysical features in maps. The present study starts from the results obtained by Epstein and Sambe and moves forward to generalize them. An off-diagonal hypervirial relationship, connecting the matrix elements of a given differentiable function of position f(r) to its derivatives ∇f(r), via the anticommutator ∇_αf,p^_{α +} with the canonical momentum operator p^, has first been proven. Afterward, this relationship is applied to show that the equations proposed by Sambe to check the quality and conservation of computed electronic current densities can be obtained as particular cases of this general theorem, with a substantial gain in computational efficiency. Connections with previous work by Arrighini, Maestro, and Moccia are outlined, and the implications that hint at future work are discussed.

Anisotropy of the vorticity tensor as a magnetic indicator of aromaticity

The vorticity vector of the current density J^B, induced in the electron cloud of a molecule by a magnetic field B, is defined by V^B = ∇ × J^B.

Does induced current density explain the C-H and C-F Perlin effects?

Spin-spin coupling constant (SSCC) data may be useful in providing information on the stereochemistry and intramolecular interactions in molecules. One-bond C-H and C-F SSCCs (1JCH and 1JCF) are amongst the most important NMR parameters used to study the structure of alicyclic six-membered rings, because the Perlin effect, defined as 1JCHeq > 1JCHax, and the fluorine Perlin-like effect, defined as |1JCFeq| > |1JCFax|, are in wide currency to probe the conformations of these compounds. The origin of these effects has been usually attributed either to dipolar interactions or hyperconjugation, while the induced current density (ICD) has been recently correlated to the magnetic shielding tensors in some six-membered ring compounds, and then used to explain the Perlin effect. Accordingly, this work reports an analysis of the ICD as a descriptor of 1JCH and 1JCF in a series of six-membered rings to find out the role of the ICD in the conventional and fluorine Perlin effect. The atoms in molecules (AIM) magnetic responses obtained from density functional theory (DFT) calculations for the studied compounds did not show any relationship of the first-order electronic current density (ΔJ(1)) with the calculated Δ1JCHax-eq and Δ1JCFax-eq values. Consequently, the title effects cannot be precisely explained by ICD. Nevertheless, an interesting relationship between ΔJ(1) and Δdelocalization involving σCH* orbitals is observed.

Frequency-dependent current density tensors as density functions of dynamic polarizabilities

Relationships accounting for contributions to the first-order charge density ρ⁽¹⁾ and current density J⁽¹⁾, induced in the electrons of a molecule by a monochromatic plane wave, have been obtained via time-dependent quantum mechanical perturbation theory. Their gauge invariance and invariance in passive translations of the coordinate system have been demonstrated (i) within the (long wavelengths) electric dipole approximation, in which only the electric field contributions to these densities are required, and (ii) within the electric quadrupole approximation, in which nonseparable terms provided by the magnetic field and by the electric field gradient, assumed uniform over the molecular domain, are needed. It is shown that the physical meaning of current density tensors depending on the frequency ω of the monochromatic wave shone on the molecule, and corresponding to derivatives of the current density J⁽¹⁾(r, ω) with respect to components of the perturbing fields, is that of property density. Therefore, frequency-dependent current density tensors can be interpreted as the integrand function in 3d-space integrals defining dynamic molecular response tensors, e.g., electric dipole, mixed electric dipole-magnetic dipole, and electric dipole-quadrupole polarizabilites. Plots of current density tensors are expected to provide important information on the molecular domains giving dominant contributions to these properties.

Delocalization energy retrieved from the current density tensor

The anisotropy of the magnetically induced current density tensor can be computed by its original formulation, known in the literature as ACID, or by its revised definition, AACID, which takes into account the asymmetric nature of the tensor. In polycyclic aromatic hydrocarbons, the excess of the integrated value of AACID over that computed for isolated ethylenes correlate with the Hückel delocalization energy.

Gauge invariance and origin independence of electronic charge density and current density induced by optical fields

Expressions for the first-order polarization charge density ρ ⁽¹⁾ and current density J ⁽¹⁾ induced in a molecule by a monochromatic plane wave, obtained by time-dependent quantum mechanical perturbation theory, have been investigated to assess their gauge invariance and independence of the coordinate system in passive and active translations. The conditions arrived at show that, within the (long wavelengths) dipole approximation, only the electric contributions to these densities are needed to rationalize the phenomenology. To the next higher quadrupole approximation, assuming that the magnetic field and the electric field gradient are uniform over the molecular dimensions, corresponding contributions to ρ ⁽¹⁾ and J ⁽¹⁾ are considered. It has been found that total densities are independent of the origin, whereas the contributions from electric and magnetic fields are not separately invariant. A magnetic contribution to J ⁽¹⁾, which is by itself origin independent, can be defined by means of an approach based on continuous translation of the origin of the coordinate system.