Inverse Scattering Transform and Solitons for Square Matrix Nonlinear Schrödinger Equations