We choose parameters which reproduce the pseudo-scalar and vector meson masses as done in usual constituent quark models at \(B=0\)  [1]}, [2]}.
Including the first order pertrubative corrections from \(V_{E,M}\) , the meson spectrum takes the form (\(M=m_1+m_2\)  and \(\mu ^{-1} = m_1^{-1} + m_2^{-1}\) )
\(E_{\rm NRQ} (\vec{K} ) = \frac{\, \vec{K} ^2 \,}{\, 2M\,} + E_{\rm in} \,,\) 
