A natural scheme for the quantitative analysis of the magnetically induced molecular current density using an...

A new natural scheme is introduced to analyze quantitatively the magnetically induced molecular current density vector field, J. The set of zero points of J, which is called its stagnation graph (SG), has been previously used to study the topological features of the current density of various molecules. Here, the line integrals of the induced magnetic field along edges of the connected subset of the SG are calculated. The edges are oriented such that all weights, i.e., flux values become non-negative, thereby, an oriented flux-weighted (current density) stagnation graph (OFW-SG) is obtained. Since in the exact theoretical limit, J is divergence-free and due to the topological characteristics of such vector fields, the flux of all separate vortices (current density domains) and neighbouring connected vortices can be determined exactly by adding the weights of cyclic subsets of edges (i.e., closed loops) of the OFW-SG. The procedure is exemplified by the minimal example of LiH for a weak homogeneous external magnetic field, B, perpendicular to the chemical bond. The OFW-SG exhibits one closed loop (formally decomposed into two edges), and an open line extending to infinity on both of its ends. The method provides the means of accurately determining the strength of the current density even in molecules with a complicated set of distinct vortices.