By exploiting the supersymmetry-inspired factorization method through a judicious use of deuteron ground state wave function, higher partial wave nucleon–nucleon potentials, both energy independent and energy dependent, are generated. We adopt the phase function method to deal with the scattering phase shifts and demonstrate the usefulness of our constructed potentials by means of model calculation.

By exploiting supersymmetry-inspired factorization method together with a judiciously chosen deuteron ground-state wave function, approximate higher partial wave nucleon–nucleon potentials are generated. In this context, a minor modification is also introduced to the generated potentials. The n–p scattering phase shifts are computed and analysed via the phase function method.

By exploiting supersymmetry and factorization method, higher partial wave nucleon–nucleon potentials ($\ell = 1,2,3$) for a few selected triplet and singlet states are generated from the ground-state interaction and wave function. The nuclear Hulthen potential and the corresponding wave function with the parameters of Arnold and Mackellar which fit the deuteron binding energy are used as the starting point of our calculation. The scattering phase shifts are computed for the constructed potentials using phase function method to examine the merit of our approach to the problem.