Susy Design of Smooth Quantum Rings in Graphene

We develop a suitable technique to design zero-energy graphene models with radial electrostatic potentials capable of achieving electrostatic confinement. Using the Gaussian law for electrostatics, we derive the charge density associated with these potentials that correspond to concentric electrostatic rings. The technique is based on a modified supersymmetric transformation that allows us to design time-reversal invariant interaction terms and to find the corresponding zero-energy bound states in analytical form. Consequently, solutions with the same probability density but different angular momentum are characterized by circular probability currents flowing in opposite directions. The energies of the systems defined in two Dirac valleys (one-valley) have a fourfold (twofold) degeneracy. As an example of the technique, we construct a ring-decorated Coulomb potential that exhibits zero-energy collapse and bound states together.