  This work analyzes monolayer graphene in external electromagnetic fields,
which is described by the Dirac equation with minimal coupling. Supersymmetric
quantum mechanics allows building new Dirac equations with modified magnetic
fields. Here, we will use complex factorization energies and iterate the method
in order to arrive at Hermitian graphene Hamiltonians. Finally, we compare
these results with the matrix supersymmetric quantum mechanics approach.
