“…The first-order intertwining method presented here generalizes the shape invariant technique that different authors have used previously to describe the behavior of an electron near to the Dirac points in a graphene layer with applied external magnetic fields [23][24][25][26][27][28]. Similarly as Figure 7: Second intertwining for the Morse potential with k = 6α, ν 1 = − 3 2 , ν 2 = − 1 2 , 1 = − 1 2 E − 1 = − 11α 2 2 , 2 = −E − 1 = −11α 2 : (a) generated potential V 2 (x, 2 ) (continuous line) and initial oneṼ 1 (x, 1 ) (dashed line), with energy levels E n (x)| 2 for the ground state (GS, blue) and the excited states n = 1, 2, 3 (red, green, purple); (d) current density for the excited states with the same colors that in (c).…”