Non-Hermitian SUSY hydrogen-like Hamiltonians with real spectra

It is shown that the radial part of the hydrogen Hamiltonian factorizes as the product of two not mutually adjoint first-order differential operators plus a complex constant epsilon. The 1-SUSY approach is used to construct non-Hermitian operators with hydrogen spectra. Other non-Hermitian Hamiltonians are shown to admit an extra 'complex energy' at epsilon. New self-adjoint hydrogen-like Hamiltonians are also derived by using a 2-SUSY transformation with complex conjugate pairs epsilon, bar epsilon.