Current Plots

Because we are considering a scattering problem, not only the quantum mechanical density but also the quantum mechanical current is of interest. Because the function makeReImAbsContourPlot sets up definitions for the density and gradient, we can define the quantum mechanical current j easily through the following formula.

Here is a 3D plot of the norm of the current density for our example solution. We reuse the function make3DPlot to conveniently adapt to the shape of the region .

We continue with a simple vector field visualization of the current.

We see the characteristic vortices appearing in such structures [7].