On the topology of total and diamagnetic induced electronic currents in molecules

An application of the continuous transformation of the origin of the current density (CTOCD) scheme to constrain the diamagnetic induced charge current density (J(d)) to be divergenceless is introduced. This results in a family of J(d) fields perpendicular and proportional to both the gradient of the electron density and the external magnetic field. Since, in the limit of a complete basis set calculation, the paramagnetic component J(p) also becomes divergenceless, we call this scheme CTOCD-DC (CTOCD for Divergenceless Components). CTOCD-DC allows for a topological characterization of both J(d) and J(p) in terms of their stagnation graphs. All stagnation graphs of J(d) from CTOCD-DC contain the zero points of the gradient of the unperturbed electron density (del rho). In this way, an intimate topological relation between rho and the diamagnetic current contribution is revealed. Numerical experiments exemplified by the case of LiNHF in point group symmetry C-1 suggest that the corresponding paramagnetic current contributions J(p) can show tendencies to accumulate pseudo-stagnation lines in proximity of some kind of the zero points of del rho. Common zero points of del rho and the total currents are exactly zero points of the mechanical momentum density.