Abstract
We calculate for the first time the static properties of the deuteron, within the framework of supersymmetric quantum mechanics, analytically. A new effective potential and its partner are derived from a superpotential so that all parameters are fitted by the experimental data. An analytical expression is obtained for the deuteron wave function and contributions of the orthogonal ¹³S₁ and ¹³D₁ states are determined, explicitly. Compared to one pion exchange, the superpotential produces an electrostatic as well as two pion exchange terms for the potential. The saddle point radius of the potential and the maximum of the wave function are linearly proportional.