$\phi^+ = (\omega_1 + i \omega_2)/\sqrt{2}$, $\phi^{0}=(\phi_r^0 + i \phi_i^0)/\sqrt{2}$, $\chi^{++}=(\chi_1 +
i \chi_2)/\sqrt{2}$, $\chi^{+}=(\chi_3 + i \chi_4)/\sqrt{2}$, $\chi^{0}=(\chi_r^0 + i \chi_i^0)/\sqrt{2}$ and $\xi^{+}=(\xi_1 + i \xi_2)/\sqrt{2}$. We consider all possible two-body scattering channels ($\Psi_{a} \Psi_{b}
\to \Psi_{c} \Psi_{d}$) not only for the neutral two-body states as initial and final states but also for the singly-, the doubly-, the triply- and the quadruply-charged two-body states. There are totally 91 initial (or final) two-body states, in which 25 are the neutral, 36 are singly-charged, 22 are doubly-charged, 6 are triply-charged, and the last 2 are the quadruply-charged states. We construct the 91 $\times$ 91 transition matrix of high-energy S-wave amplitudes, and then calculate their eigenvalues.

The initial (final) two-body states can be treated separately as $\Psi_D
\Psi_D$, $\Psi_T \Psi_T$ and $\Psi_D \Psi_T$. The high-energy S-wave amplitudes are block-diagonalized by the electric charge and also the discrete $Z_2$ symmetry ($\Phi \to \Phi$ and $\Delta\to -\Delta$). Each submatrix with respect to the $\Psi_D\Psi_D$ or $\Psi_T\Psi_T$ states can also be classified by irreducible decomposition of direct products of the representations for $O(4)$ as &=&(1)\_D\^ (9) (6),\
&=&(1)\_T\^ (44)(36). The only singlet and symmetric representations, (1)\_[a]{}&=& \_[k=1]{} (\^k\_a)\^2,\
(s)\^[ij]{} &=& \^i\_[a]{} \^j\_a - \_[k=1]{} (\^k\_a)\^2, 