The Variable S-matrix approach offers a unique way to extract low energy
threshold parameters for a given NN potential. We extract those parameters for
the np system from the NijmII and Reid93 potentials, to all partial waves with
total angular momentum j≤m5, using the generalized effective range
expansion % $$ ({\bf f}^{sj}){l,l’} k^l k^{l’} = -({{\bf a}^{sj}}^{-1}){l,l’}
- \frac12 ({\bf r}^{sj}){l,l’} k^2 + ({\bf v_2}^{sj}){l,l’} k^4 + + ({\bf
v_3}^{sj}){l,l’} k^6 + + ({\bf v_4}^{sj}){l,l’} k^8 + … - {\rm i}
k^{l+l’+1} $ {\bf f}^{sj} = ({\bf S}^{sj} - {\bf 1}) / (2 {\rm i} k)
isthescatteringamplitudeand{\bf S}^{sj}$ is the unitary S-matrix in
coupled channel space with total spin s and total angular momentum j. Our
calculation includes all the relevant contributions of the full amplitude to
order O(k8) in the CM momentum. We also discuss the validity of the
generalized effective range expansion in the region of analyticity $ k \le
m_\pi /2$.