We look for least energy solutions to the cooperative systems of coupled Schrödinger equations - Δ u i + λ i u i = ∂ i G (u) in R N , N ≥ 3 , u i ∈ H 1 (R N) , ∫ R N | u i | 2 d x ≤ ρ i 2 i ∈ { 1 , ... , K } with G ≥ 0 , where ρ i > 0 is prescribed and (λ i , u i) ∈ R × H 1 (R N) is to be determined, i ∈ { 1 , ⋯ , K } . Our approach is based on the minimization of the energy functional over a linear combination of the Nehari and Pohožaev constraints intersected with the product of the closed balls in L 2 (R N) of radii ρ i , which allows to provide general growth assumptions about G and to know in advance the sign of the corresponding Lagrange multipliers. We assume that G has at least L 2 -critical growth at 0 and admits Sobolev critical growth. The more assumptions we make about G, N, and K, the more can be said about the minimizers of the corresponding energy functional. In particular, if K = 2 , N ∈ { 3 , 4 } , and G satisfies further assumptions, then u = (u 1 , u 2) is normalized, i.e., ∫ R N | u i | 2 d x = ρ i 2 for i ∈ { 1 , 2 } . [ABSTRACT FROM AUTHOR]