Vortex patterns for a cubic extreme-type-II superconductor with spherical cavity

The cavity contains a magnetic dipole which is switched off initially. The parameter $\mu$ cotnrols the strength of the dipole.

The left hand sides of the plots show the Gibbs energy of the system as a function of the control parameter. The right hand sides show the complex-valued order parameter $\psi$: $arg(\psi)$ is highlighted at the back sides of the cube, and the isosurface $|\psi|^2 = 0.1$ is shown. As $\mu$ increases, the typical vortex loops are showing.