“…The separability of Dirac's equation in axially symmetric potentials has been widely studied [38,39], [40, p. 134], even with a nonminimal (or anomalous moment) coupling [41][42][43], driven by applications in, e.g., electron vortex beam experiments [44][45][46][47], Kumakhov (or channeling) radiation studies [48][49][50], the Aharonov-Bohm (-Casher) effect, and neutron confinement [51]. There are only two axisymmetric (electrostatic) potentials for which the Dirac equation is solvable, viz., the constant potential and the −1/ρ potential [52][53][54]. The former opens the possibility of piecewiseconstant potentials capable of modelling a cylindrical WG, e.g., the finite square well [55,Sec.…”