“…In light of its potential applicative relevance, in this work, we develop the inverse scattering transform (IST) for the system (2) as a tool to solve the initial-value problem, as well as to obtain explicit soliton solutions. While the IST for "unreduced" matrix NLS systems, and for the "canonical" reductions corresponding to cases 1 and 2 (focusing and defocusing matrix NLS), is well established, both with zero and nonzero boundary conditions (see, for instance, [5,17,[31][32][33][34] and references therein), the IST and the soliton solutions corresponding to the reductions in cases 3 and 4 described above are novel, and present some interesting aspects and additional challenges with respect to the other two cases in that one needs to impose suitable constraints on the norming constants to guarantee that the soliton solutions are smooth for all x, t ∈ R. As a matter of fact, this work also provides several advances as far as the IST for general matrix NLS systems is concerned, including the well-known focusing and defocusing matrix NLS systems corresponding to cases 1 and 2. Specific focus of the work is to: (i) provide a rigorous definition of the norming constants that does not require any unjustified analytic extension of the scattering relations, and clarify the role of the rank of the norming constants in the spectral characterization of the corresponding solutions; (ii) properly account for all the symmetries in the potential matrix, and derive the corresponding symmetries in the scattering data (reflection coefficients and norming constants); (iii) formulate the inverse problem as a Riemann-Hilbert problem (RHP), instead of in terms of Marchenko equations; (iv) obtain novel soliton solutions for the reductions of Eq.…”