The Oriented and Flux-Weighted Current Density Stagnation Graph of LiH

A scheme is introduced to quantitatively analyze the magnetically induced molecular current density vector field $\mathbf{J}$. After determining the set of zero points of $\mathbf{J}$, which is called its {\em stagnation graph} (SG), the line integrals $\Phi_{\ell_i}=-\frac{1}{\mu_0} \int_{\ell_i} \mathbf{B}_\mathrm{ind}\cdot\mathrm{d}\mathbf{l}$ along all edges $\ell_i$ of the connected subset of the SG are determined. The edges $\ell_i$ are oriented such that all $\Phi_{\ell_i}$ are non-negative and they are weighted with $\Phi_{\ell_i}$. An oriented flux-weighted (current density) stagnation graph (OFW-SG) is obtained. Since $\mathbf{J}$ is in the exact theoretical limit divergence free and due to the topological characteristics of such vector fields the flux of all separate vortices and neighbouring vortex combinations can be determined by adding the weights of cyclic subsets of edges of the OFW-SG. The procedure is exemplified by the case of LiH for a perpendicular and weak homogeneous external magnetic field $\mathbf{B}$}