Darboux transformations of the one-dimensional stationary Dirac equation

A matricial Darboux operator intertwining two one-dimensional stationary Dirac Hamiltonians is constructed. This operator is such that the potential of the second Dirac Hamiltonian, as well as the corresponding eigenfunctions, are determined through the knowledge of only two eigenfunctions of the first Dirac Hamiltonian. Moreover this operator, together with its adjoint and the two Hamiltonians, generate a quadratic deformation of the superalgebra subtending the usual supersymmetric quantum mechanics. Our developments are illustrated in the free-particle case and the generalized Coulomb interaction. In the latter case, a relativistic counterpart of shape invariance is observed.