## Factorization of the R-matrix. II.

### November 27, 2006

Department of Mathematics, St Petersburg Technology Institute St.Petersburg, Russia. E-mail: S.Derkachov@pobox.spbu.ru Abstract. We study the general rational solution of the Yang-Baxter equation with the supersymmetry algebra $s\ell \left(2|1\right)$  . The R-operator acting in the tensor product of two arbitrary representations of the supersymmetry algebra can be represented as the product of the simpler ”building blocks” – $\mathcal{ℛ}$  -operators.
Contents

1 Introduction

In the previous paper [1we have shown that the general R-matrix can be represented as the product of the much more simpler $\mathcal{ℛ}$  -operators. In this paper we shall consider the general rational solution of the Yang-Baxter equation with the supersymmetry algebra $s\ell \left(2|1\right)$  and show that there exists the same factorization. In fact all calculations are very similar to the $s\ell \left(3\right)$  -example and modifications due to supersymmetry are simple. The generalization of the previous results [1to the algebra of supersymmetry is mainly motivated by the possible applications to the super Yang-Mills theory [2, 3, 12, 5.
The presentation is organized as follows. In Section 2 we collect the standard facts about the algebra $s\ell \left(2|1\right)$  and its representations. We represent the lowest weight modules by polynomials in one even variable ( $z$  ) and two odd variables ( $\theta ,\overline{\theta }$  ) and the $s\ell \left(2|1\right)$  -generators as first order differential operators. We use the notation and formulae from the paper [8. In Section 3 we derive the defining relation for the general R-matrix, i.e. the solution of the Yang-Baxter equation acting on tensor products of two arbitrary representations, the elements of which are polynomials in variables ${z}_{1},{\theta }_{1},{\overline{\theta }}_{1}$  and ${z}_{2},{\theta }_{2},{\overline{\theta }}_{2}$  . In Section 4 we introduce the natural defining equations for the $\mathcal{ℛ}$  -operators and show that the general R-matrix can be represented as the product of such much more simple operators.
Finally, in Section 4 we summarize. In Appendix we calculate the matrix elements of the $\mathcal{ℛ}$  -operators and as consequence obtain the matrix elements of R-matrix in full agreement with the results of the paper [8.

2 $sl\left(2|1\right)$  lowest weight modules

The superalgebra $s\ell \left(2|1\right)$  has eight generators: four odd ${\mathbf{V}}^{±},{\mathbf{W}}^{±}$  and four even $\mathbf{S},{\mathbf{S}}^{±}$  and $\mathbf{B}$  .
Using the natural notations ${\mathbf{E}}_{31}={\mathbf{S}}^{-};{\mathbf{E}}_{21}=-{\mathbf{W}}^{-};{\mathbf{E}}_{32}={\mathbf{V}}^{-};{\mathbf{E}}_{13}={\mathbf{S}}^{+};{\mathbf{E}}_{23}={\mathbf{W}}^{+};{\mathbf{E}}_{12}={\mathbf{V}}^{+}$  ${\mathbf{E}}_{11}=\mathbf{B}-\mathbf{S};{\mathbf{E}}_{22}=-2\mathbf{B};{\mathbf{E}}_{33}=\mathbf{B}+\mathbf{S}.$  the commutation relations for the generators of $s\ell \left(2|1\right)$  can be written compactly in the form [6, 7$\left[{\mathbf{E}}_{AB},{\mathbf{E}}_{CD}\right]={\delta }_{CB}{\mathbf{E}}_{AD}-\left(-{\right)}^{\left(\overline{A}+\overline{B}\right)\left(\overline{C}+\overline{D}\right)}{\delta }_{AD}{\mathbf{E}}_{CB};A,B,C,D=1,2,3$  where the graded commutator is defined as (we choose the grading $\overline{1}=\overline{3}=0,\overline{2}=1$  ) $\left[{\mathbf{E}}_{AB},{\mathbf{E}}_{CD}\right]\equiv {\mathbf{E}}_{AB}\cdot {\mathbf{E}}_{CD}-\left(-{\right)}^{\left(\overline{A}+\overline{B}\right)\left(\overline{C}+\overline{D}\right)}{\mathbf{E}}_{CD}\cdot {\mathbf{E}}_{AB}.$  There are two central elements [6, 7, 10 ${\mathbf{C}}_{2}=\frac{1}{2}{\sum }_{AB}\left(-{\right)}^{\overline{B}}{\mathbf{E}}_{AB}{\mathbf{E}}_{BA}={\mathbf{S}}^{2}-{\mathbf{B}}^{2}+{\mathbf{S}}^{+}{\mathbf{S}}^{-}+{\mathbf{V}}^{+}{\mathbf{W}}^{-}+{\mathbf{W}}^{+}{\mathbf{V}}^{-};{\mathbf{C}}_{3}=\frac{1}{6}{\sum }_{ABC}\left(-{\right)}^{\overline{B}+\overline{C}}{\mathbf{E}}_{AB}{\mathbf{E}}_{BC}{\mathbf{E}}_{CA}$  The Verma module is the generic lowest weight $s\ell \left(2|1\right)$  -module ${\mathbf{V}}_{\Lambda };\Lambda =\left(\ell ,b\right)$  . As a linear space ${\mathbf{V}}_{\Lambda }$  is spanned by the basis with even elements ${\mathbf{a}}_{k},{\mathbf{b}}_{k}$  ${\mathbf{a}}_{k}={\mathbf{S}}_{+}^{k}{\mathbf{a}}_{0},k=0,1,2\cdots ;{\mathbf{b}}_{k}={\mathbf{S}}_{+}^{k-1}{\mathbf{W}}_{+}{\mathbf{V}}_{+}{\mathbf{a}}_{0},k=1,2\cdots$  and odd elements ${\mathbf{v}}_{k},{\mathbf{w}}_{k}$  ${\mathbf{v}}_{k}={\mathbf{S}}_{+}^{k}{\mathbf{V}}_{+}{\mathbf{a}}_{0};{\mathbf{w}}_{k}={\mathbf{S}}_{+}^{k}{\mathbf{W}}_{+}{\mathbf{a}}_{0},k=0,1,2\cdots$  The vector ${\mathbf{a}}_{0}$  is the lowest weight vector:
${\mathbf{S}}_{-}{\mathbf{a}}_{0}={\mathbf{V}}_{-}{\mathbf{a}}_{0}={\mathbf{W}}_{-}{\mathbf{a}}_{0}=0;\mathbf{S}{\mathbf{a}}_{0}=\ell \cdot {\mathbf{a}}_{0};\mathbf{B}{\mathbf{a}}_{0}=b\cdot {\mathbf{a}}_{0}$  We shall use the representation ${V}_{\Lambda }$  of $s\ell \left(2|1\right)$  in the infinite-dimensional space $\mathbb{C}\left[Z\right]$  where $Z=\left(z,\theta ,\overline{\theta }\right)$  of polynomials in even variable $z$  and odd variables $\theta ,\overline{\theta }$  with the monomial basis ${\left\{{z}^{k},\theta \overline{\theta }{z}^{k};\theta {z}^{k},\overline{\theta }{z}^{k}\right\}}_{k=0}^{\infty }$  and lowest weight vector ${a}_{0}=1$  [8. The action of $s\ell \left(2|1\right)$  in ${V}_{\Lambda }$  is given by the first-order differential operators
 $\begin{array}{c}{S}^{-}=-\partial ;{V}^{-}={\partial }_{\theta }+\frac{1}{2}\overline{\theta }\partial ;{W}^{-}={\partial }_{\overline{\theta }}+\frac{1}{2}\theta \partial \end{array}$ (2.1)
${V}^{+}=-\left[z{\partial }_{\theta }+\frac{1}{2}\overline{\theta }z\partial +\frac{1}{2}\overline{\theta }\theta {\partial }_{\theta }\right]-\left(\ell -b\right)\overline{\theta };{W}^{+}=-\left[z{\partial }_{\overline{\theta }}+\frac{1}{2}\theta z\partial +\frac{1}{2}\theta \overline{\theta }{\partial }_{\overline{\theta }}\right]-\left(\ell +b\right)\theta$  ${S}^{+}={z}^{2}\partial +z\theta {\partial }_{\theta }+z\overline{\theta }{\partial }_{\overline{\theta }}+2\ell z-b\theta \overline{\theta };S=z\partial +\frac{1}{2}\theta {\partial }_{\theta }+\frac{1}{2}\overline{\theta }{\partial }_{\overline{\theta }}+\ell ;B=\frac{1}{2}\overline{\theta }{\partial }_{\overline{\theta }}-\frac{1}{2}\theta {\partial }_{\theta }+b$  It is possible to derive the closed expressions for the elements of the basis ${a}_{k}={S}_{+}^{k}\cdot 1=\left(2\ell {\right)}_{k}\left[z-\frac{kb}{2\ell }\cdot \theta \overline{\theta }\right]\cdot {z}^{k-1};{b}_{k}={S}_{+}^{k-1}{W}_{+}{V}_{+}\cdot 1=\frac{\ell -b}{2\ell }\left(2\ell {\right)}_{k}\left[z+\left(b+\ell +\frac{k}{2}\right)\cdot \theta \overline{\theta }\right]\cdot {z}^{k-1}$  ${v}_{k}={S}_{+}^{k}{V}_{+}\cdot 1=-\left(\ell -b\right)\left(2\ell +1{\right)}_{k}{z}^{k}\overline{\theta };{w}_{k}={S}_{+}^{k}{W}_{+}\cdot 1=-\left(\ell +b\right)\left(2\ell +1{\right)}_{k}{z}^{k}\theta ;\left(2\ell {\right)}_{k}\equiv \frac{\Gamma \left(2\ell +k\right)}{\Gamma \left(2\ell \right)}$  It is evident that for the generic $\ell \ne -\frac{n}{2}$  the module ${V}_{\Lambda }$  is an irreducible lowest weight $s\ell \left(2|1\right)$  -module isomorphic to ${\mathbf{V}}_{\Lambda }$  but for the special values of the spin $\ell =-\frac{n}{2}$  there exists the finite dimensional invariant subspace. There are three cases depending on the relation between $b$  and $n$  [6, 9, 10. The first case is for generic $b\ne ±\frac{n}{2}$  (typical representations) and there exists the $4n$  -dimensional invariant subspace. The second and third cases appear for $b=±\frac{n}{2}$  (atypical representations). For the chiral representation $b=-\frac{n}{2}$  the $\left(2n+1\right)$  -dimensional invariant subspace is spanned on the vectors ${\Phi }_{k}^{+}={\left(z-\frac{\theta \overline{\theta }}{2}\right)}^{k},k=0...n;{W}_{k}=\theta {z}^{k},k=0...n-1$  and for the antichiral representation $b=\frac{n}{2}$  the $\left(2n+1\right)$  -dimensional invariant subspace is spanned on the vectors ${\Phi }_{k}^{-}={\left(z+\frac{\theta \overline{\theta }}{2}\right)}^{k},k=0...n;{V}_{k}=\overline{\theta }{z}^{k},k=0...n-1$  We shall use the three-dimensional chiral representation $V$  .In the basis ${\mathbf{e}}_{1}={S}_{+}\cdot 1=-z+\frac{\theta \overline{\theta }}{2},{\mathbf{e}}_{2}={W}_{+}\cdot 1=\theta ,{\mathbf{e}}_{3}=1$  the $s\ell \left(2|1\right)$  -generators take the form
 $\begin{array}{c}{\mathbf{s}}_{-}=\left(\begin{array}{ccc}0& 0& 0\\ 0& 0& 0\\ 1& 0& 0\end{array}\right);{\mathbf{w}}_{-}=\left(\begin{array}{ccc}0& 0& 0\\ -1& 0& 0\\ 0& 0& 0\end{array}\right);{\mathbf{v}}_{-}=\left(\begin{array}{ccc}0& 0& 0\\ 0& 0& 0\\ 0& 1& 0\end{array}\right);\mathbf{s}=\left(\begin{array}{ccc}\frac{1}{2}& 0& 0\\ 0& 0& 0\\ 0& 0& -\frac{1}{2}\end{array}\right)\end{array}$ (2.2)
${\mathbf{s}}_{+}=\left(\begin{array}{ccc}0& 0& 1\\ 0& 0& 0\\ 0& 0& 0\end{array}\right);{\mathbf{w}}_{+}=\left(\begin{array}{ccc}0& 0& 0\\ 0& 0& 1\\ 0& 0& 0\end{array}\right);{\mathbf{v}}_{+}=\left(\begin{array}{ccc}0& 1& 0\\ 0& 0& 0\\ 0& 0& 0\end{array}\right);\mathbf{b}=\left(\begin{array}{ccc}-\frac{1}{2}& 0& 0\\ 0& -1& 0\\ 0& 0& -\frac{1}{2}\end{array}\right)$  There exists the second three-dimensional representation antichiral representation $\overline{V}$  . In the basis ${\mathbf{e}}_{1}={S}_{+}\cdot 1=-z-\frac{\theta \overline{\theta }}{2},{\mathbf{e}}_{2}={V}_{+}\cdot 1=\overline{\theta },{\mathbf{e}}_{3}=1$  the $s\ell \left(2|1\right)$  -generators take the similar form as in chiral representation but ${\mathbf{v}}^{±}↔{\mathbf{w}}^{±};\mathbf{b}\to -\mathbf{b}$  .
We use the standard definition for the matrix ${A}_{ik}$  of the linear operator $\mathbf{A}$  in the basis $\left\{{\mathbf{e}}_{k}\right\}$  $\mathbf{A}{\mathbf{e}}_{k}={\sum }_{i}{\mathbf{e}}_{i}{A}_{ik}$

3 Yang-Baxter equation and Lax operator

The Yang-Baxter equation is the following three term relation [11, 12 ${\mathbb{R}}_{{\Lambda }_{1}{\Lambda }_{2}}\left(u-v\right){\mathbb{R}}_{{\Lambda }_{1}{\Lambda }_{3}}\left(u\right){\mathbb{R}}_{{\Lambda }_{2}{\Lambda }_{3}}\left(v\right)={\mathbb{R}}_{{\Lambda }_{2}{\Lambda }_{3}}\left(v\right){\mathbb{R}}_{{\Lambda }_{1}{\Lambda }_{3}}\left(u\right){\mathbb{R}}_{{\Lambda }_{1}{\Lambda }_{2}}\left(u-v\right)$  for the operators ${\mathbb{R}}_{{\Lambda }_{i}{\Lambda }_{j}}\left(u\right):{V}_{{\Lambda }_{i}}\otimes {V}_{{\Lambda }_{j}}\to {V}_{{\Lambda }_{i}}\otimes {V}_{{\Lambda }_{j}}$  . We start from the simplest solutions of Yang-Baxter equation and derive the defining equation for the general $\mathbb{R}$  -operator [8, 15. First we put ${\Lambda }_{1}={\Lambda }_{2}={\Lambda }_{3}=\left(-\frac{1}{2},-\frac{1}{2}\right)$  in Yang-Baxter equation and consider the restriction on the invariant subspace $V\otimes V\otimes V$  . We obtain the equation ${\mathbb{R}}_{12}\left(u-v\right){\mathbb{R}}_{13}\left(u\right){\mathbb{R}}_{23}\left(v\right)={\mathbb{R}}_{23}\left(v\right){\mathbb{R}}_{13}\left(u\right){\mathbb{R}}_{12}\left(u-v\right)$  where the operator ${\mathbb{R}}_{12}\left(u\right)$  acts on the first and second copy of $V$  in the tensor product $V\otimes V\otimes V$  and similarly for the other $\mathbb{R}$  -operators. The solution is well known [11, 12 ${\mathbb{R}}_{12}\left(u\right)=u+{P}_{12}$  where ${P}_{12}$  is the (graded)permutation operator in $V\otimes V$  . We choose the basis $\left\{{\mathbf{e}}_{1},{\mathbf{e}}_{2},{\mathbf{e}}_{3}\right\}$  in $V$  so that ${\mathbf{e}}_{1},{\mathbf{e}}_{3}$  are even elements and ${\mathbf{e}}_{2}$  is odd element and our grading is: $\overline{1}=\overline{3}=0,\overline{2}=1$  .
The permutation operator acts on the basis as follows ${P}_{12}{\mathbf{e}}_{i}\otimes {\mathbf{e}}_{k}=\left(-1{\right)}^{\overline{i}\overline{k}}{\mathbf{e}}_{k}\otimes {\mathbf{e}}_{i}$  and additional sign arises for ${\mathbf{e}}_{2}\otimes {\mathbf{e}}_{2}$  only. Secondly we choose ${\Lambda }_{1}={\Lambda }_{2}=\left(-\frac{1}{2},-\frac{1}{2}\right);{\Lambda }_{3}=\Lambda =\left(\ell ,b\right)$  and consider the restriction on the invariant subspace $V\otimes V\otimes {V}_{\Lambda }$  . The restriction of the operator ${\mathbb{R}}_{{\Lambda }_{1}\Lambda }\left(u\right)$  to the space $V\otimes {V}_{\Lambda }$  coincides up to normalization and shift of spectral parameter with the Lax-operator [11, 12, 13$L\left(u\right):V\otimes {V}_{\Lambda }\to V\otimes {V}_{\Lambda }$  and the Yang-Baxter equation coincides with the defining equation for the Lax-operator ${\mathbb{R}}_{12}\left(u-v\right){L}^{\left(1\right)}\left(u\right){L}^{\left(2\right)}\left(v\right)={L}^{\left(2\right)}\left(v\right){L}^{\left(1\right)}\left(u\right){\mathbb{R}}_{12}\left(u-v\right)$  where ${L}^{\left(1\right)}\left(u\right)$  is the operator which acts nontrivially on the first copy of $V$  and ${V}_{\Lambda }$  in the tensor product $V\otimes V\otimes {V}_{\Lambda }$  and ${L}^{\left(2\right)}\left(u\right)$  is the operator which acts nontrivially on the second copy of $V$  and ${V}_{\Lambda }$  . The solution coincides up to additive constant with the Casimir operator ${\mathbf{C}}_{2}$  for the representation $V\otimes {V}_{\Lambda }$  [12, 14$L\left(u\right)\equiv u+2\mathbf{s}\otimes S-2\mathbf{b}\otimes B+{\mathbf{v}}_{+}\otimes {W}_{-}+{\mathbf{s}}_{+}\otimes {S}_{-}-{\mathbf{w}}_{-}\otimes {V}_{+}+{\mathbf{w}}_{+}\otimes {V}_{-}+{\mathbf{s}}_{-}\otimes {S}_{+}-{\mathbf{v}}_{-}\otimes {W}_{+}$  where $\mathbf{s},\mathbf{b},{\mathbf{s}}_{±},{\mathbf{v}}_{±},{\mathbf{w}}_{±}$  are $s\ell \left(2|1\right)$  -generators in the chiral representation ( 2.2 ) and $S,B,{S}_{±},{V}_{±},{W}_{±}$  are generators in the generic representation. The algebra $s\ell \left(2|1\right)$  has two three-dimensional representations – chiral $V$  and antichiral $\overline{V}$  so that there exists the second Lax-operator $\overline{L}\left(u\right):\overline{V}\otimes {V}_{\Lambda }\to \overline{V}\otimes {V}_{\Lambda };{\mathbb{R}}_{12}\left(u-v\right){\overline{L}}^{\left(1\right)}\left(u\right){\overline{L}}^{\left(2\right)}\left(v\right)={\overline{L}}^{\left(2\right)}\left(v\right){\overline{L}}^{\left(1\right)}\left(u\right){\mathbb{R}}_{12}\left(u-v\right)$  The explicit expression for the second Lax-operator is the same but now $\mathbf{s},\mathbf{b},{\mathbf{s}}_{±},{\mathbf{v}}_{±},{\mathbf{w}}_{±}$  are $s\ell \left(2|1\right)$  -generators in the antichiral representation $\overline{V}$  . In matrix form we obtain [14, 8 $L\left(u\right)=\left(\begin{array}{ccc}S+B+u& -{W}_{-}& {S}_{-}\\ {V}_{+}& 2B+u& {V}_{-}\\ {S}_{+}& {W}_{+}& B-S+u\end{array}\right);\overline{L}\left(u\right)=\left(\begin{array}{ccc}S-B+u& -{V}_{-}& {S}_{-}\\ {W}_{+}& -2B+u& {W}_{-}\\ {S}_{+}& {V}_{+}& -B-S+u\end{array}\right)$  Note the arising some additional signs due to grading. For example, we have ${\mathbf{v}}_{+}{\mathbf{e}}_{2}={\mathbf{e}}_{1}$  but ${\mathbf{v}}_{+}\otimes {W}_{-}{\mathbf{e}}_{2}=-{\mathbf{v}}_{+}{\mathbf{e}}_{2}{W}_{-}=-{\mathbf{e}}_{1}{W}_{-}$  so that one obtains ${\mathbf{v}}_{+}{\mathbf{e}}_{2}={\mathbf{e}}_{1}⇒{\mathbf{v}}_{+}=\left(\begin{array}{ccc}0& 1& 0\\ 0& 0& 0\\ 0& 0& 0\end{array}\right);{\mathbf{v}}_{+}\otimes {W}_{-}{\mathbf{e}}_{2}=-{\mathbf{e}}_{1}{W}_{-}⇒{\mathbf{v}}_{+}\otimes {W}_{-}=\left(\begin{array}{ccc}0& -{W}_{-}& 0\\ 0& 0& 0\\ 0& 0& 0\end{array}\right)$  We shall use the chiral Lax-operator $L\left(u\right)$  in defining equation for the general $\mathbb{R}$  -operator. The Lax-operator depends on three parameters $u,\ell ,b$  and we shall use the parametrization ${u}_{1}=u+b+\ell ,{u}_{2}=u+2b,{u}_{3}=u+b-\ell ;\ell =\frac{{u}_{1}-{u}_{3}}{2},b={u}_{2}-\frac{{u}_{1}+{u}_{3}}{2}$  The explicit form of the Lax operator $L\left({u}_{1},{u}_{2},{u}_{3}\right)$  in the functional representation ${V}_{\Lambda }$  is
 $\begin{array}{c}L\left({u}_{1},{u}_{2},{u}_{3}\right)=\left(\begin{array}{ccc}z\partial +\overline{\theta }{\partial }_{\overline{\theta }}+{u}_{1}& -\left({\partial }_{\overline{\theta }}+\frac{1}{2}\theta \partial \right)& -\partial \\ {L}_{21}& \overline{\theta }{\partial }_{\overline{\theta }}-\theta {\partial }_{\theta }+{u}_{2}& {\partial }_{\theta }+\frac{1}{2}\overline{\theta }\partial \\ {L}_{31}& {L}_{32}& -z\partial -\theta {\partial }_{\theta }+{u}_{3}\end{array}\right)\end{array}$ (3.1)
${L}_{21}=-\left(z-\frac{\theta \overline{\theta }}{2}\right){\partial }_{\theta }-\frac{1}{2}\overline{\theta }z\partial +\left({u}_{2}-{u}_{1}\right)\overline{\theta };{L}_{32}=-\left(z+\frac{\theta \overline{\theta }}{2}\right){\partial }_{\overline{\theta }}-\frac{1}{2}\theta z\partial +\left({u}_{3}-{u}_{2}\right)\theta$  ${L}_{31}={z}^{2}\partial +z\left(\theta {\partial }_{\theta }+\overline{\theta }{\partial }_{\overline{\theta }}\right)+\left({u}_{1}-{u}_{3}\right)z+\left({u}_{1}+{u}_{3}-2{u}_{2}\right)\frac{\theta \overline{\theta }}{2}$  There exists the useful factorized representation for the Lax-operator
 $\begin{array}{c}L\left({u}_{1},{u}_{2},{u}_{3}\right)\equiv \left(\begin{array}{ccc}1& 0& 0\\ -\overline{\theta }& 1& 0\\ z+\frac{\theta \overline{\theta }}{2}& -\theta & 1\end{array}\right)\left(\begin{array}{ccc}{u}_{1}& {D}^{-}& -\partial \\ 0& {u}_{2}-1& -{D}^{+}\\ 0& 0& {u}_{3}\end{array}\right)\left(\begin{array}{ccc}1& 0& 0\\ \overline{\theta }& 1& 0\\ -z+\frac{\theta \overline{\theta }}{2}& \theta & 1\end{array}\right)\end{array}$ (3.2)
where ${D}^{±}$  are covariant derivatives ${D}^{-}=-{\partial }_{\overline{\theta }}+\frac{1}{2}\theta \partial ,{D}^{+}=-{\partial }_{\theta }+\frac{1}{2}\overline{\theta }\partial$  The $L$  -operator is $s\ell \left(2|1\right)$  -invariant by construction and as consequence one obtains the equality
 $\begin{array}{c}{\mathbf{M}}^{-1}\cdot L\left(u\right)\cdot \mathbf{M}={\mathbb{S}}^{-1}\cdot L\left(u\right)\cdot \mathbb{S};\mathbb{S}={e}^{\alpha {V}_{-}}\cdot {e}^{\overline{\alpha }{W}_{-}}\cdot {e}^{\left(\lambda +\frac{\alpha \overline{\alpha }}{2}\right){S}_{-}};\mathbf{M}=\left(\begin{array}{ccc}1& 0& 0\\ -\overline{\alpha }& 1& 0\\ \lambda +\frac{\alpha \overline{\alpha }}{2}& -\alpha & 1\end{array}\right)\end{array}$ (3.3)
Finally we put ${\Lambda }_{1}=\left(-\frac{1}{2},-\frac{1}{2}\right)$  in Yang-Baxter equation ${\mathbb{R}}_{{\Lambda }_{1}{\Lambda }_{2}}\left(u-v\right){\mathbb{R}}_{{\Lambda }_{1}{\Lambda }_{3}}\left(u\right){\mathbb{R}}_{{\Lambda }_{2}{\Lambda }_{3}}\left(v\right)={\mathbb{R}}_{{\Lambda }_{2}{\Lambda }_{3}}\left(v\right){\mathbb{R}}_{{\Lambda }_{1}{\Lambda }_{3}}\left(u\right){\mathbb{R}}_{{\Lambda }_{1}{\Lambda }_{2}}\left(u-v\right)$  change the numeration of the representation spaces ${\Lambda }_{2}\to {\Lambda }_{1}=\left({\ell }_{1},{b}_{1}\right);{\Lambda }_{3}\to {\Lambda }_{2}=\left({\ell }_{2},{b}_{2}\right)$  and consider the restriction on the invariant subspace $V\otimes {V}_{{\Lambda }_{1}}\otimes {V}_{{\Lambda }_{2}}$  . In this way one obtains the defining equation for the $\mathbb{R}$  -operator ${L}_{1}\left(u-v\right){L}_{2}\left(u\right){\mathbb{R}}_{{\Lambda }_{1}{\Lambda }_{2}}\left(v\right)={\mathbb{R}}_{{\Lambda }_{1}{\Lambda }_{2}}\left(v\right){L}_{2}\left(u\right){L}_{1}\left(u-v\right)$  The operator ${L}_{k}$  acts nontrivially on the tensor product $V\otimes {V}_{{\Lambda }_{k}}$  which is isomorphic to $V\otimes \mathbb{C}\left[{Z}_{k}\right]$  where ${Z}_{k}=\left({z}_{k},{\theta }_{k},{\overline{\theta }}_{k}\right)$  and the operator ${\mathbb{R}}_{{\Lambda }_{1}{\Lambda }_{2}}\left(u\right)$  acts nontrivially on the tensor product ${V}_{{\Lambda }_{1}}\otimes {V}_{{\Lambda }_{2}}$  which is isomorphic to $\mathbb{C}\left[{Z}_{1},{Z}_{2}\right]=\mathbb{C}\left[{Z}_{1}\right]\otimes \mathbb{C}\left[{Z}_{2}\right]$  . Note that obtained defining equation is slightly different from the ones which was used in [1and [8. The defining equation which is similar to [1, 8is
 $\begin{array}{c}{\mathbb{R}}_{{\Lambda }_{1}{\Lambda }_{2}}^{-1}\left(v-u\right){L}_{1}\left(u\right){L}_{2}\left(v\right)={L}_{2}\left(v\right){L}_{1}\left(u\right){\mathbb{R}}_{{\Lambda }_{1}{\Lambda }_{2}}^{-1}\left(v-u\right)\end{array}$ (3.4)
There exists the well known automorphism of the Yang-Baxter equation ${\mathbb{R}}_{{\Lambda }_{1}{\Lambda }_{2}}\left(u\right)\to {\mathbb{R}}_{{\Lambda }_{1}{\Lambda }_{2}}^{-1}\left(-u\right)$  .
In the simplest $s\ell \left(2\right)$  case we have ${\mathbb{R}}_{{\ell }_{1}{\ell }_{2}}\left(u\right)\sim {\mathbb{R}}_{{\ell }_{1}{\ell }_{2}}^{-1}\left(-u\right)$  but for the more complicated algebras the action of this automorphism is nontrivial. To proceed in close analogy with [1, 8we shall use the defining equation ( 3.4 ) so that we derive the expression for the operator ${\mathbb{R}}_{{\Lambda }_{1}{\Lambda }_{2}}^{-1}\left(v-u\right)$  .

4 The general R-matrix

It is useful to extract the operator of permutation ${\mathbb{P}}_{12}:\mathbb{C}\left[{Z}_{1}\right]\otimes \mathbb{C}\left[{Z}_{2}\right]\to \mathbb{C}\left[{Z}_{2}\right]\otimes \mathbb{C}\left[{Z}_{1}\right];{\mathbb{P}}_{12}\Psi \left({Z}_{1},{Z}_{2}\right)=\Psi \left({Z}_{2},{Z}_{1}\right)$  from the $\mathbb{R}$  -operator ${\mathbb{R}}_{{\Lambda }_{1}{\Lambda }_{2}}^{-1}\left(v-u\right)={\mathbb{P}}_{12}{\stackrel{ˇ}{\mathbb{R}}}_{{\Lambda }_{1}{\Lambda }_{2}}\left(u;v\right)$  and solve the defining equation for the $\stackrel{ˇ}{\mathbb{R}}$  -operator. The main defining equation for the $\stackrel{ˇ}{\mathbb{R}}$  -operator is $\stackrel{ˇ}{\mathbb{R}}\left(u;v\right){L}_{1}\left({u}_{1},{u}_{2},{u}_{3}\right){L}_{2}\left({v}_{1},{v}_{2},{v}_{3}\right)={L}_{1}\left({v}_{1},{v}_{2},{v}_{3}\right){L}_{2}\left({u}_{1},{u}_{2},{u}_{3}\right)\stackrel{ˇ}{\mathbb{R}}\left(u;v\right)$  ${u}_{1}=u+{b}_{1}+{\ell }_{1},{u}_{2}=u+2{b}_{1},{u}_{3}=u+{b}_{1}-{\ell }_{1};{v}_{1}=v+{b}_{2}+{\ell }_{2},{v}_{2}=v+2{b}_{2},{v}_{3}=v+{b}_{2}-{\ell }_{2}$  The operator $\stackrel{ˇ}{\mathbb{R}}$  interchanges all parameters in the product of two $L$  -operators and similar to the $s\ell \left(3\right)$  -case $\stackrel{ˇ}{\mathbb{R}}$  -operator can be represented as the product of the simpler ”elementary building blocks” - $\mathcal{ℛ}$  -operators.
Proposition 1 There exists operator ${\mathcal{ℛ}}_{1}$  which is the solution of the defining equations
 $\begin{array}{c}{\mathcal{ℛ}}_{1}{L}_{1}\left({u}_{1},{u}_{2},{u}_{3}\right){L}_{2}\left({v}_{1},{v}_{2},{v}_{3}\right)={L}_{1}\left({v}_{1},{u}_{2},{u}_{3}\right){L}_{2}\left({u}_{1},{v}_{2},{v}_{3}\right){\mathcal{ℛ}}_{1}\end{array}$ (4.1)
${\mathcal{ℛ}}_{1}={\mathcal{ℛ}}_{1}\left({u}_{1}|{v}_{1},{v}_{2},{v}_{3}\right);{\mathcal{ℛ}}_{1}\left({u}_{1}|{v}_{1},{v}_{2},{v}_{3}\right)={\mathcal{ℛ}}_{1}\left({u}_{1}+\lambda |{v}_{1}+\lambda ,{v}_{2}+\lambda ,{v}_{3}+\lambda \right)$  and these requirements fix the operator ${\mathcal{ℛ}}_{1}$  up to overall normalization constant ${\mathcal{ℛ}}_{1}\sim {\mathbb{S}}_{1}^{-1}\cdot \left[\frac{\Gamma \left({z}_{2}{\partial }_{2}+{u}_{1}-{v}_{3}+1\right)}{\Gamma \left({z}_{2}{\partial }_{2}+{v}_{1}-{v}_{3}+1\right)}\left({f}_{1}+{\overline{\theta }}_{2}{\partial }_{{\overline{\theta }}_{2}}\right)-\frac{\Gamma \left({z}_{2}{\partial }_{2}+{u}_{1}-{v}_{3}\right)}{\Gamma \left({z}_{2}{\partial }_{2}+{v}_{1}-{v}_{3}+1\right)}{z}_{2}{\partial }_{{\theta }_{2}}{\partial }_{{\overline{\theta }}_{2}}\right]\cdot {\mathbb{S}}_{1}$  ${\mathbb{S}}_{1}={e}^{\frac{{\theta }_{2}{\overline{\theta }}_{2}}{2}{\partial }_{2}}\cdot {e}^{{\theta }_{1}{V}_{2}^{-}}\cdot {e}^{{\overline{\theta }}_{1}{W}_{2}^{-}}\cdot {e}^{-\left({z}_{1}+\frac{{\theta }_{1}{\overline{\theta }}_{1}}{2}\right){S}_{2}^{-}};{f}_{1}=\frac{{v}_{1}-{v}_{2}}{{u}_{1}-{v}_{1}}$
Proposition 2 There exists operator ${\mathcal{ℛ}}_{2}$  which is the solution of the defining equations
 $\begin{array}{c}{\mathcal{ℛ}}_{2}{L}_{1}\left({u}_{1},{u}_{2},{u}_{3}\right){L}_{2}\left({v}_{1},{v}_{2},{v}_{3}\right)={L}_{1}\left({u}_{1},{v}_{2},{u}_{3}\right){L}_{2}\left({v}_{1},{u}_{2},{v}_{3}\right){\mathcal{ℛ}}_{2}\end{array}$ (4.2)
${\mathcal{ℛ}}_{2}={\mathcal{ℛ}}_{2}\left({u}_{1},{u}_{2}|{v}_{2},{v}_{3}\right);{\mathcal{ℛ}}_{2}\left({u}_{1},{u}_{2}|{v}_{2},{v}_{3}\right)={\mathcal{ℛ}}_{1}\left({u}_{1}+\lambda ,{u}_{2}+\lambda |{v}_{2}+\lambda ,{v}_{3}+\lambda \right)$  and these requirements fix the operator ${\mathcal{ℛ}}_{2}$  up to overall normalization constant ${\mathcal{ℛ}}_{2}\sim {\mathbb{S}}_{2}^{-1}\cdot \left[{f}_{2}+{u}_{12}\cdot {\theta }_{2}{\partial }_{{\theta }_{2}}+{v}_{23}\cdot {\overline{\theta }}_{1}{\partial }_{{\overline{\theta }}_{1}}+\left({z}_{12}+{\theta }_{1}{\overline{\theta }}_{2}\right){\partial }_{{\overline{\theta }}_{1}}{\partial }_{{\theta }_{2}}+\left({u}_{2}-{v}_{2}\right){\theta }_{2}{\overline{\theta }}_{1}{\partial }_{{\overline{\theta }}_{1}}{\partial }_{{\theta }_{2}}\right]\cdot {\mathbb{S}}_{2}$  ${\mathbb{S}}_{2}={e}^{{\theta }_{1}{\partial }_{{\theta }_{2}}}\cdot {e}^{{\overline{\theta }}_{2}{\partial }_{{\overline{\theta }}_{1}}}\cdot {e}^{\frac{{\theta }_{1}{\overline{\theta }}_{1}}{2}{\partial }_{1}}\cdot {e}^{-\frac{{\theta }_{2}{\overline{\theta }}_{2}}{2}{\partial }_{2}};{f}_{2}=\frac{{u}_{21}{v}_{23}}{{v}_{2}-{u}_{2}},{u}_{12}={u}_{1}-{u}_{2},{v}_{23}={v}_{2}-{v}_{3}$
Proposition 3 There exists operator ${\mathcal{ℛ}}_{3}$  which is the solution of the defining equations
 $\begin{array}{c}{\mathcal{ℛ}}_{3}{L}_{1}\left({u}_{1},{u}_{2},{u}_{3}\right){L}_{2}\left({v}_{1},{v}_{2},{v}_{3}\right)={L}_{1}\left({u}_{1},{u}_{2},{v}_{3}\right){L}_{2}\left({v}_{1},{v}_{2},{u}_{3}\right){\mathcal{ℛ}}_{3}\end{array}$ (4.3)
${\mathcal{ℛ}}_{3}={\mathcal{ℛ}}_{3}\left({u}_{1},{u}_{2},{u}_{3}|{v}_{3}\right);{\mathcal{ℛ}}_{3}\left({u}_{1},{u}_{2},{u}_{3}|{v}_{3}\right)={\mathcal{ℛ}}_{3}\left({u}_{1}+\lambda ,{u}_{2}+\lambda ,{u}_{3}+\lambda |{v}_{3}+\lambda \right)$  and these requirements fix the operator ${\mathcal{ℛ}}_{3}$  up to overall normalization constant ${\mathcal{ℛ}}_{3}\sim {\mathbb{S}}_{3}^{-1}\cdot \left[\frac{\Gamma \left({z}_{1}{\partial }_{1}+{u}_{1}-{v}_{3}+1\right)}{\Gamma \left({z}_{1}{\partial }_{1}+{u}_{1}-{u}_{3}+1\right)}\left({f}_{3}+{\theta }_{1}{\partial }_{{\theta }_{1}}\right)+\frac{\Gamma \left({z}_{1}{\partial }_{1}+{u}_{1}-{v}_{3}\right)}{\Gamma \left({z}_{1}{\partial }_{1}+{u}_{1}-{u}_{3}+1\right)}{z}_{1}{\partial }_{{\theta }_{1}}{\partial }_{{\overline{\theta }}_{1}}\right]\cdot {\mathbb{S}}_{3}$  ${\mathbb{S}}_{3}={e}^{-\frac{{\theta }_{1}{\overline{\theta }}_{1}}{2}{\partial }_{1}}\cdot {e}^{{\theta }_{2}{V}_{1}^{-}}\cdot {e}^{{\overline{\theta }}_{2}{W}_{1}^{-}}\cdot {e}^{-\left({z}_{2}+\frac{{\theta }_{2}{\overline{\theta }}_{2}}{2}\right){S}_{1}^{-}};{f}_{3}=\frac{{u}_{2}-{u}_{3}}{{u}_{3}-{v}_{3}}$
Proposition 4 The $\stackrel{ˇ}{\mathbb{R}}$  -operator can be factorized as follows $\stackrel{ˇ}{\mathbb{R}}\left(u;v\right)={\mathcal{ℛ}}_{1}\left({u}_{1};{v}_{1},{u}_{2},{u}_{3}\right){\mathcal{ℛ}}_{2}\left({u}_{1},{u}_{2};{v}_{2},{u}_{3}\right){\mathcal{ℛ}}_{3}\left({u}_{1},{u}_{2},{u}_{3};{v}_{3}\right)$
There exist six equivalent ways to represent $\stackrel{ˇ}{\mathbb{R}}$  in an factorized form which differ by the order of $\mathcal{ℛ}$  -operators and their parameters. All these expressions and the proof of the factorization of the $\stackrel{ˇ}{\mathbb{R}}$  -operator can be obtained using the pictures similar to [1.
The defining system of equations for the $\mathcal{ℛ}$  -operator can be reduced to the simpler system which clearly shows the property of $s\ell \left(2|1\right)$  -covariance of the $\mathcal{ℛ}$  -operator.
Lemma 1 The defining equation ( 4.1 ) for the operator ${\mathcal{ℛ}}_{1}$  is equivalent to the system of equations
 $\begin{array}{c}{\mathcal{ℛ}}_{1}\left[{L}_{1}\left({u}_{1},{u}_{2},{u}_{3}\right)+{L}_{2}\left({v}_{1},{v}_{2},{v}_{3}\right)\right]=\left[{L}_{1}\left({v}_{1},{u}_{2},{u}_{3}\right)+{L}_{2}\left({u}_{1},{v}_{2},{v}_{3}\right)\right]{\mathcal{ℛ}}_{1}\end{array}$ (4.4)
${\mathcal{ℛ}}_{1}{z}_{1}={z}_{1}{\mathcal{ℛ}}_{1},{\mathcal{ℛ}}_{1}{\theta }_{1}={\theta }_{1}{\mathcal{ℛ}}_{1},{\mathcal{ℛ}}_{1}{\overline{\theta }}_{1}={\overline{\theta }}_{1}{\mathcal{ℛ}}_{1}$
 $\begin{array}{c}{\mathcal{ℛ}}_{1}\cdot \left({V}_{2}^{-}+{\overline{\theta }}_{1}{S}_{2}^{-}\right)=\left({V}_{2}^{-}+{\overline{\theta }}_{1}{S}_{2}^{-}\right)\cdot {\mathcal{ℛ}}_{1}\end{array}$ (4.5)
Lemma 2 The defining equation ( 4.2 ) for the operator ${\mathcal{ℛ}}_{2}$  is equivalent to the system of equations
 $\begin{array}{c}{\mathcal{ℛ}}_{2}\left[{L}_{1}\left({u}_{1},{u}_{2},{u}_{3}\right)+{L}_{2}\left({v}_{1},{v}_{2},{v}_{3}\right)\right]=\left[{L}_{1}\left({u}_{1},{v}_{2},{u}_{3}\right)+{L}_{2}\left({v}_{1},{u}_{2},{v}_{3}\right)\right]{\mathcal{ℛ}}_{2}\end{array}$ (4.6)
$\left[{\mathcal{ℛ}}_{2},{z}_{1}-\frac{{\theta }_{1}{\overline{\theta }}_{1}}{2}\right]=0,{\mathcal{ℛ}}_{2}{\theta }_{1}={\theta }_{1}{\mathcal{ℛ}}_{2};\left[{\mathcal{ℛ}}_{2},{z}_{2}+\frac{{\theta }_{2}{\overline{\theta }}_{2}}{2}\right]=0,{\mathcal{ℛ}}_{2}{\overline{\theta }}_{2}={\overline{\theta }}_{2}{\mathcal{ℛ}}_{2}$
Lemma 3 The defining equation ( 4.3 ) for the operator ${\mathcal{ℛ}}_{3}$  is equivalent to the system of equations
 $\begin{array}{c}{\mathcal{ℛ}}_{3}\left[{L}_{1}\left({u}_{1},{u}_{2},{u}_{3}\right)+{L}_{2}\left({v}_{1},{v}_{2},{v}_{3}\right)\right]=\left[{L}_{1}\left({u}_{1},{u}_{2},{v}_{3}\right)+{L}_{2}\left({v}_{1},{v}_{2},{u}_{3}\right)\right]{\mathcal{ℛ}}_{3}\end{array}$ (4.7)
${\mathcal{ℛ}}_{3}{z}_{2}={z}_{2}{\mathcal{ℛ}}_{3},{\mathcal{ℛ}}_{3}{\theta }_{2}={\theta }_{2}{\mathcal{ℛ}}_{3},{\mathcal{ℛ}}_{3}{\overline{\theta }}_{2}={\overline{\theta }}_{2}{\mathcal{ℛ}}_{3}$
 $\begin{array}{c}{\mathcal{ℛ}}_{3}\cdot \left({W}_{1}^{-}+{\theta }_{2}{S}_{1}^{-}\right)=\left({W}_{1}^{-}+{\theta }_{2}{S}_{1}^{-}\right)\cdot {\mathcal{ℛ}}_{3}\end{array}$ (4.8)
The relations in the first line are simply the rules of commutation of $\mathcal{ℛ}$  -operators with $s\ell \left(2|1\right)$  -generators written in a compact form. In explicit notations we have for ${\Lambda }_{1}=\left({\ell }_{1},{b}_{1}\right)$  and ${\Lambda }_{2}=\left({\ell }_{2},{b}_{2}\right)$  $\mathcal{ℛ}:{V}_{{\Lambda }_{1}}\otimes {V}_{{\Lambda }_{2}}\to {V}_{{\Lambda }_{1}^{\prime }}\otimes {V}_{{\Lambda }_{2}^{\prime }}$  ${\mathcal{ℛ}}_{1}:{\Lambda }_{1}^{\prime }=\left({\ell }_{1}-{\xi }_{1},{b}_{1}+{\xi }_{1}\right);{\Lambda }_{2}^{\prime }=\left({\ell }_{2}+{\xi }_{1},{b}_{2}-{\xi }_{1}\right);{\xi }_{1}=\frac{{u}_{1}-{v}_{1}}{2}$  ${\mathcal{ℛ}}_{2}:{\Lambda }_{1}^{\prime }=\left({\ell }_{1},{b}_{1}-{\xi }_{2}\right);{\Lambda }_{2}^{\prime }=\left({\ell }_{2},{b}_{2}+{\xi }_{2}\right);{\xi }_{2}={u}_{2}-{v}_{2}$  ${\mathcal{ℛ}}_{3}:{\Lambda }_{1}^{\prime }=\left({\ell }_{1}+{\xi }_{3},{b}_{1}+{\xi }_{3}\right);{\Lambda }_{2}^{\prime }=\left({\ell }_{2}-{\xi }_{3},{b}_{2}-{\xi }_{3}\right);{\xi }_{3}=\frac{{u}_{3}-{v}_{3}}{2}$  The $s\ell \left(2|1\right)$  -invariance of $\mathbb{R}$  -matrix follows directly from the properties of $\mathcal{ℛ}$  -operators so that the general R-matrix ${\mathbb{R}}_{{\Lambda }_{1}{\Lambda }_{2}}^{-1}\left(v-u\right)={\mathbb{P}}_{12}{\stackrel{ˇ}{\mathbb{R}}}_{{\Lambda }_{1}{\Lambda }_{2}}\left(u;v\right)$  is automatically $s\ell \left(2|1\right)$  -invariant.
Proof Now we are going to the proof of equivalence of defining equation ( 4.3 ) to the system ( 4.7 ) and derivation of explicit formula for the operator ${\mathcal{ℛ}}_{3}$  . First we show that the system ( 4.7 ) is the direct consequence of the eq. ( 4.3 ). Let us make the shift ${u}_{k}\to {u}_{k}+\lambda ,{v}_{1}\to {v}_{1}+\mu ,{v}_{2}\to {v}_{2}+\nu ,{v}_{3}\to {v}_{3}+\lambda$  in the defining equation ( 4.3 ).The $\mathcal{ℛ}$  -operator is invariant under this shift and $L$  -operators transform as follows ${L}_{1}\to {L}_{1}+\lambda \cdot 1l;{L}_{2}\to {L}_{2}+\lambda \cdot 1l+\left(\mu -\lambda \right)\left(\begin{array}{ccc}1& 0& 0\\ -{\overline{\theta }}_{2}& 0& 0\\ {z}_{2}+\frac{{\theta }_{2}{\overline{\theta }}_{2}}{2}& 0& 0\end{array}\right)+\left(\nu -\lambda \right)\left(\begin{array}{ccc}0& 0& 0\\ {\overline{\theta }}_{2}& 1& 0\\ -{\theta }_{2}{\overline{\theta }}_{2}& -{\theta }_{2}& 0\end{array}\right)$  After all one obtains the equation which contains the arbitrary parameters $\lambda$  , $\mu$  and $\nu$  and as consequence we derive the system ( 4.7 ) and equation ( 4.8 ). Next we show that from the systems of equations ( 4.7 ), ( 4.8 ) follows eq. ( 4.3 ). This will be almost evident if we rewrite these equations in equivalent form using the $s\ell \left(2|1\right)$  -invariance of the $L$  -operator and the commutativity of ${\mathcal{ℛ}}_{3}$  and ${z}_{2},{\theta }_{2},{\overline{\theta }}_{2}$  . We substitute the factorized representation ( 3.2 ) for the operator ${L}_{2}$  in the defining equation for the operator ${\mathcal{ℛ}}_{3}$  ${\mathcal{ℛ}}_{3}{L}_{1}\left({u}_{1},{u}_{2},{u}_{3}\right)\mathbf{M}\left(\begin{array}{ccc}{v}_{1}& {D}_{2}^{-}& -{\partial }_{2}\\ 0& {v}_{2}-1& -{D}_{2}^{+}\\ 0& 0& {v}_{3}\end{array}\right){\mathbf{M}}^{-1}={L}_{1}\left({u}_{1},{u}_{2},{v}_{3}\right)\mathbf{M}\left(\begin{array}{ccc}{v}_{1}& {D}_{2}^{-}& -{\partial }_{2}\\ 0& {v}_{2}-1& -{D}_{2}^{+}\\ 0& 0& {u}_{3}\end{array}\right){\mathbf{M}}^{-1}{\mathcal{ℛ}}_{3}$  and perform the similarity transformation ${\mathbf{M}}^{-1}\cdots \mathbf{M}$  of this matrix equation using the commutativity ${\mathcal{ℛ}}_{3}$  and ${z}_{2},{\theta }_{2},{\overline{\theta }}_{2}$  . Then using the $s\ell \left(2|1\right)$  -invariance of $L$  -operator ( 3.3 ) ${\mathbf{M}}^{-1}\cdot {L}_{1}\cdot \mathbf{M}={\mathbb{S}}^{-1}\cdot {L}_{1}\cdot \mathbb{S};\mathbb{S}={e}^{{\theta }_{2}{V}_{1}^{-}}\cdot {e}^{{\overline{\theta }}_{2}{W}_{1}^{-}}\cdot {e}^{-\left({z}_{2}+\frac{{\theta }_{2}{\overline{\theta }}_{2}}{2}\right){S}_{1}^{-}};\mathbf{M}=\left(\begin{array}{ccc}1& 0& 0\\ -{\overline{\theta }}_{2}& 1& 0\\ {z}_{2}+\frac{{\theta }_{2}{\overline{\theta }}_{2}}{2}& -{\theta }_{2}& 1\end{array}\right)$  we derive the equation for the transformed operator $\mathbf{R}=\mathbb{S}\cdot {\mathcal{ℛ}}_{3}\cdot {\mathbb{S}}^{-1}$
 $\begin{array}{c}\mathbf{R}\cdot {L}_{1}\left({u}_{1},{u}_{2},{u}_{3}\right)\mathbf{L}\left({v}_{1},{v}_{2},{v}_{3}\right)={L}_{1}\left({u}_{1},{u}_{2},{v}_{3}\right)\mathbf{L}\left({v}_{1},{v}_{2},{u}_{3}\right)\cdot \mathbf{R}\end{array}$ (4.9)
where $\mathbf{L}\left({v}_{1},{v}_{2},{v}_{3}\right)\equiv \mathbb{S}\cdot \left(\begin{array}{ccc}{v}_{1}& {D}_{2}^{-}& -{\partial }_{2}\\ 0& {v}_{2}-1& -{D}_{2}^{+}\\ 0& 0& {v}_{3}\end{array}\right)\cdot {\mathbb{S}}^{-1}=\left(\begin{array}{ccc}{v}_{1}& {D}_{2}^{-}+{W}_{1}^{-}& -{\partial }_{2}+{\partial }_{1}\\ 0& {v}_{2}-1& -{D}_{2}^{+}-{V}_{1}^{-}\\ 0& 0& {v}_{3}\end{array}\right)$  To derive the system of equations which is equivalent to the system ( 4.7 ), ( 4.8 ) written in terms of $\mathbf{R}$  we repeat the same trick with the shift of parameters and obtain the system of equations
 $\begin{array}{c}\mathbf{R}\cdot \left[{L}_{1}\left({u}_{1},{u}_{2},{u}_{3}\right)+\mathbf{L}\left({v}_{1},{v}_{2},{v}_{3}\right)\right]=\left[{L}_{1}\left({u}_{1},{u}_{2},{v}_{3}\right)+\mathbf{L}\left({v}_{1},{v}_{2},{u}_{3}\right)\right]\cdot \mathbf{R}\end{array}$ (4.10)
 $\begin{array}{c}\mathbf{R}\cdot {L}_{1}\left({u}_{1},{u}_{2},{u}_{3}\right)\left(\begin{array}{cc}1& 0\\ 0& 1\\ 0& 0\end{array}\right)={L}_{1}\left({u}_{1},{u}_{2},{v}_{3}\right)\left(\begin{array}{cc}1& 0\\ 0& 1\\ 0& 0\end{array}\right)\mathbf{R}\end{array}$ (4.11)
It is evident that all equations of the system ( 4.11 ) contained in the equation ( 4.10 ) except only one $\left(12\right)$  -equation $\mathbf{R}{W}_{1}^{-}={W}_{1}^{-}\mathbf{R}$  . We use the system of equation
 $\begin{array}{c}\mathbf{R}\cdot \left[{L}_{1}\left({u}_{1},{u}_{2},{u}_{3}\right)+\mathbf{L}\left({v}_{1},{v}_{2},{v}_{3}\right)\right]=\left[{L}_{1}\left({u}_{1},{u}_{2},{v}_{3}\right)+\mathbf{L}\left({v}_{1},{v}_{2},{u}_{3}\right)\right]\cdot \mathbf{R};\mathbf{R}{W}_{1}^{-}={W}_{1}^{-}\mathbf{R}\end{array}$ (4.12)
as defining system for operator $\mathbf{R}$  . This system is equivalent to the system ( 4.7 ), ( 4.8 ). Returning to the system ( 4.9 ) (it is the system ( 4.3 ) written in terms of $\mathbf{R}$  ) we note that it is possible to factorize the matrix $diag\left({v}_{1};{v}_{2}-1;1\right)$  from the right $\mathbf{R}\cdot {L}_{1}\left({u}_{1},{u}_{2},{u}_{3}\right)\left(\begin{array}{ccc}1& \frac{{D}_{2}^{-}+{W}_{1}^{-}}{{v}_{2}-1}& -{\partial }_{2}+{\partial }_{1}\\ 0& 1& -{D}_{2}^{+}-{V}_{1}^{-}\\ 0& 0& {v}_{3}\end{array}\right)={L}_{1}\left({u}_{1},{u}_{2},{v}_{3}\right)\left(\begin{array}{ccc}1& \frac{{D}_{2}^{-}+{W}_{1}^{-}}{{v}_{2}-1}& -{\partial }_{2}+{\partial }_{1}\\ 0& 1& -{D}_{2}^{+}-{V}_{1}^{-}\\ 0& 0& {u}_{3}\end{array}\right)\cdot \mathbf{R}$  In comparison with ( 4.12 ) there are three new equations only
 $\begin{array}{c}\mathbf{R}\cdot {L}_{1}\left({u}_{1},{u}_{2},{u}_{3}\right)\left(\begin{array}{c}{\partial }_{1}\\ -{V}_{1}^{-}\\ {v}_{3}\end{array}\right)={L}_{1}\left({u}_{1},{u}_{2},{v}_{3}\right)\left(\begin{array}{c}{\partial }_{1}\\ -{V}_{1}^{-}\\ {u}_{3}\end{array}\right)\mathbf{R}\end{array}$ (4.13)
Indeed the system ( 4.12 ) contains the equations $\left[\mathbf{R},{D}_{2}^{±}\right]=\left[\mathbf{R},{\partial }_{2}\right]=\left[\mathbf{R},{W}_{1}^{-}\right]=0$  and by conditions ( 4.11 ) we obtain the three new equations. It is easy to check that these equations follow from the system ( 4.12 ). Finally the systems of equations ( 4.12 ) is defining and it remains to find the solution. First of all $\left[\mathbf{R},{z}_{2}\right]=\left[\mathbf{R},{\theta }_{2}\right]=\left[\mathbf{R},{\overline{\theta }}_{2}\right]=\left[\mathbf{R},{D}_{2}^{±}\right]=\left[\mathbf{R},{\partial }_{2}\right]=0$  and therefore the operator $\mathbf{R}$  depends on the variables ${z}_{1},{\theta }_{1},{\overline{\theta }}_{1}$  only. For simplicity we use the natural transformation $\mathbf{R}={e}^{\frac{1}{2}{\theta }_{1}{\overline{\theta }}_{1}{\partial }_{1}}\mathbf{r}{e}^{-\frac{1}{2}{\theta }_{1}{\overline{\theta }}_{1}{\partial }_{1}}$  change ${z}_{1},{\theta }_{1},{\overline{\theta }}_{1}\to z,\theta ,\overline{\theta }$  and obtain the system of equations
 $\begin{array}{c}\mathbf{r}{\partial }_{\overline{\theta }}={\partial }_{\overline{\theta }}\mathbf{r};\mathbf{r}\left(z\partial +\overline{\theta }{\partial }_{\overline{\theta }}\right)=\left(z\partial +\overline{\theta }{\partial }_{\overline{\theta }}\right)\mathbf{r};\mathbf{r}\left(\theta {\partial }_{\theta }-\overline{\theta }{\partial }_{\overline{\theta }}\right)=\left(\theta {\partial }_{\theta }-\overline{\theta }{\partial }_{\overline{\theta }}\right)\mathbf{r}\end{array}$ (4.14)
$\mathbf{r}\left({z}^{2}\partial +z\left(\theta {\partial }_{\theta }+\overline{\theta }{\partial }_{\overline{\theta }}\right)+z\left({u}_{1}-{u}_{3}\right)+\theta \overline{\theta }\left({u}_{3}-{u}_{2}\right)\right)=$
 $\begin{array}{c}=\left({z}^{2}\partial +z\left(\theta {\partial }_{\theta }+\overline{\theta }{\partial }_{\overline{\theta }}\right)+z\left({u}_{1}-{v}_{3}\right)+\theta \overline{\theta }\left({v}_{3}-{u}_{2}\right)\right)\mathbf{r}\end{array}$ (4.15)
 $\begin{array}{c}\mathbf{r}\left(z{\partial }_{\overline{\theta }}+\left({u}_{2}-{u}_{3}\right)\theta \right)=\left(z{\partial }_{\overline{\theta }}+\left({u}_{2}-{v}_{3}\right)\theta \right)\mathbf{r}\end{array}$ (4.16)
 $\begin{array}{c}\mathbf{r}\left(-z\left({\partial }_{\theta }+\overline{\theta }\partial \right)+\theta \overline{\theta }{\partial }_{\theta }+\left({u}_{2}-{u}_{1}\right)\overline{\theta }\right)=\left(-z\left({\partial }_{\theta }+\overline{\theta }\partial \right)+\theta \overline{\theta }{\partial }_{\theta }+\left({u}_{2}-{u}_{1}\right)\overline{\theta }\right)\mathbf{r}\end{array}$ (4.17)
First of all the equation ( 4.17 ) is not independent. It is the consequence of equations ( 4.15 ) and $\mathbf{r}{\partial }_{\overline{\theta }}={\partial }_{\overline{\theta }}\mathbf{r}$  due to commutation relation $\left[{\mathbf{S}}^{+},{\mathbf{W}}^{-}\right]={\mathbf{W}}^{+}$  . The general solution of the equations ( 4.14 ) is $\mathbf{r}=\mathbf{a}\left[z\partial \right]+\mathbf{b}\left[z\partial \right]\cdot \theta {\partial }_{\theta }+\mathbf{c}\left[z\partial \right]\cdot z{\partial }_{\theta }{\partial }_{\overline{\theta }}$  The equations ( 4.15 ) and ( 4.16 ) results in recurrence relations $\mathbf{a}\left[z\partial \right]-\mathbf{a}\left[z\partial -1\right]=\left({u}_{2}-{u}_{3}\right)\cdot \mathbf{c}\left[z\partial \right];\mathbf{b}\left[z\partial \right]=\frac{{u}_{3}-{v}_{3}}{{u}_{2}-{u}_{3}}\cdot \mathbf{a}\left[z\partial \right]$  $\mathbf{c}\left[z\partial +1\right]\cdot \left(z\partial +{u}_{1}-{u}_{3}+1\right)=\left(z\partial +{u}_{1}-{v}_{3}\right)\cdot \mathbf{c}\left[z\partial \right]$  $\mathbf{a}\left[z\partial +1\right]\cdot \left(z\partial +{u}_{1}-{u}_{3}\right)+\left({u}_{2}-{u}_{3}\right)\cdot \mathbf{c}\left[z\partial +1\right]=\left(z\partial +{u}_{1}-{v}_{3}\right)\cdot \mathbf{a}\left[z\partial \right]$  $\mathbf{a}\left[z\partial \right]+\mathbf{b}\left[z\partial \right]\cdot \left(z\partial +{u}_{1}-{u}_{3}\right)-\left({u}_{2}-{u}_{3}\right)\cdot \mathbf{c}\left[z\partial \right]=\mathbf{a}\left[z\partial -1\right]+\left(z\partial +{u}_{1}-{v}_{3}\right)\cdot \mathbf{b}\left[z\partial -1\right]$  The solution of these equations has the form $\mathbf{a}\left[z\partial \right]=\frac{\Gamma \left(z\partial +{u}_{1}-{v}_{3}+1\right)}{\Gamma \left(z\partial +{u}_{1}-{u}_{3}+1\right)};\mathbf{b}\left[z\partial \right]=\frac{{u}_{3}-{v}_{3}}{{u}_{2}-{u}_{3}}\frac{\Gamma \left(z\partial +{u}_{1}-{v}_{3}+1\right)}{\Gamma \left(z\partial +{u}_{1}-{u}_{3}+1\right)}$  $\mathbf{c}\left[z\partial \right]=\frac{{u}_{3}-{v}_{3}}{{u}_{2}-{u}_{3}}\frac{\Gamma \left(z\partial +{u}_{1}-{v}_{3}\right)}{\Gamma \left(z\partial +{u}_{1}-{u}_{3}+1\right)}$  Collect everything together we obtain the expression for the operator ${\mathcal{ℛ}}_{3}$  from the Proposition. All calculations for the operator ${\mathcal{ℛ}}_{1}$  are very similar.
It remains to prove the equivalence of defining equation ( 4.2 ) to the system ( 4.6 ) and derive the explicit formula for the operator ${\mathcal{ℛ}}_{2}$  . First we show that the system ( 4.6 ) is the direct consequence of the eq. ( 4.2 ). Let us make the shift ${u}_{1}\to {u}_{1}+\lambda ,{u}_{2}\to {u}_{2}+\lambda ,{u}_{3}\to {u}_{3}+\mu ,{v}_{1}\to {v}_{1}+\nu ,{v}_{2}\to {v}_{2}+\lambda ,{v}_{3}\to {v}_{3}+\lambda$  in the defining equation ( 4.2 ) for the operator ${\mathcal{ℛ}}_{2}$  . The $\mathcal{ℛ}$  -operator is invariant under this shift and $L$  -operators transform as follows ${L}_{1}\to {L}_{1}+\lambda \cdot 1l+\left(\mu -\lambda \right)\left(\begin{array}{ccc}0& 0& 0\\ 0& 0& 0\\ -{z}_{1}+\frac{{\theta }_{1}{\overline{\theta }}_{1}}{2}& {\theta }_{1}& 1\end{array}\right);{L}_{2}\to {L}_{2}+\lambda \cdot 1l+\left(\nu -\lambda \right)\left(\begin{array}{ccc}1& 0& 0\\ -{\overline{\theta }}_{2}& \lambda & 0\\ {z}_{2}+\frac{{\theta }_{2}{\overline{\theta }}_{2}}{2}& 0& 0\end{array}\right)$  After all one obtains the equation which contains the arbitrary parameters $\lambda$  , $\mu$  and $\nu$  and as consequence we derive the system ( 4.6 ). Next we show that from the systems of equations ( 4.6 ) follows eq. ( 4.2 ). This will be almost evident if we rewrite these equations in equivalent form using the $s\ell \left(2|1\right)$  -invariance of the $L$  -operator and the commutativity of ${\mathcal{ℛ}}_{2}$  and ${z}_{1}-\frac{{\theta }_{1}{\overline{\theta }}_{1}}{2},{\theta }_{1},{z}_{2}+\frac{{\theta }_{2}{\overline{\theta }}_{2}}{2},{\overline{\theta }}_{2}$  .
First of all it is useful to make the transformation ${\mathcal{ℛ}}_{2}={\mathbb{S}}^{-1}\cdot \mathbf{R}\cdot \mathbb{S};\mathbb{S}={e}^{\frac{{\theta }_{1}{\overline{\theta }}_{1}}{2}{\partial }_{1}}\cdot {e}^{-\frac{{\theta }_{2}{\overline{\theta }}_{2}}{2}{\partial }_{2}}$  so that $\mathbf{R}$  commutes with ${z}_{1},{\theta }_{1},{z}_{2},{\overline{\theta }}_{2}$  now. The corresponding transformation for the $L$  -operators can be easily derived using factorized representation ( 3.2 ). Next it is possible to make the two similarity transformations of the defining equation ( 4.2 ) using simple matrices which commute with operator $\mathbf{F}$  . After all these transformations the defining equation ( 4.2 ) for the $\mathbf{F}$  -operator in factorized form looks as follows $\mathbf{R}\cdot {\mathbf{l}}_{1}\left({u}_{1},{u}_{2},{u}_{3}\right)\cdot \mathbf{M}\cdot {\mathbf{l}}_{2}\left({v}_{1},{v}_{2},{v}_{3}\right)={\mathbf{l}}_{1}\left({u}_{1},{v}_{2},{u}_{3}\right)\cdot \mathbf{M}\cdot {\mathbf{l}}_{2}\left({v}_{1},{u}_{2},{v}_{3}\right)\cdot \mathbf{R};\mathbf{M}\equiv \left(\begin{array}{ccc}1& 0& 0\\ -{\overline{\theta }}_{2}& 1& 0\\ -{z}_{12}-{\theta }_{1}{\overline{\theta }}_{2}& {\theta }_{1}& 1\end{array}\right)$  ${\mathbf{l}}_{1}\left({u}_{1},{u}_{2},{u}_{3}\right)\equiv \left(\begin{array}{ccc}1& 0& 0\\ -{\overline{\theta }}_{1}& 1& 0\\ 0& 0& 1\end{array}\right)\left(\begin{array}{ccc}{u}_{1}& -{\partial }_{{\overline{\theta }}_{1}}& -{\partial }_{1}\\ 0& {u}_{2}-1& {\partial }_{{\theta }_{1}}-{\overline{\theta }}_{1}{\partial }_{1}\\ 0& 0& {u}_{3}\end{array}\right)\left(\begin{array}{ccc}1& 0& 0\\ {\overline{\theta }}_{1}& 1& 0\\ 0& 0& 1\end{array}\right)$  ${\mathbf{l}}_{2}\left({v}_{1},{v}_{2},{v}_{3}\right)\equiv \left(\begin{array}{ccc}1& 0& 0\\ 0& 1& 0\\ 0& -{\theta }_{2}& 1\end{array}\right)\left(\begin{array}{ccc}{v}_{1}& -{\partial }_{{\overline{\theta }}_{2}}+{\theta }_{2}{\partial }_{2}& -{\partial }_{2}\\ 0& {v}_{2}-1& {\partial }_{{\theta }_{2}}\\ 0& 0& {v}_{3}\end{array}\right)\left(\begin{array}{ccc}1& 0& 0\\ 0& 1& 0\\ 0& {\theta }_{2}& 1\end{array}\right)$  Next step we rewrite the defining equation for the transformed operator $\mathbf{r}={e}^{{\theta }_{1}{\partial }_{{\theta }_{2}}}\cdot {e}^{{\overline{\theta }}_{2}{\partial }_{{\overline{\theta }}_{1}}}\cdot \mathbf{R}\cdot {e}^{-{\theta }_{1}{\partial }_{{\theta }_{2}}}\cdot {e}^{-{\overline{\theta }}_{2}{\partial }_{{\overline{\theta }}_{1}}}$  in the form
 $\begin{array}{c}\mathbf{r}\cdot {\mathbf{L}}_{1}\left({u}_{1},{u}_{2},{u}_{3}\right)\mathbf{m}{\mathbf{L}}_{2}\left({v}_{1},{v}_{2},{v}_{3}\right)={\mathbf{L}}_{1}\left({u}_{1},{v}_{2},{u}_{3}\right)\mathbf{m}{\mathbf{L}}_{2}\left({v}_{1},{u}_{2},{v}_{3}\right)\cdot \mathbf{r}\end{array}$ (4.18)
where ${\mathbf{L}}_{1}\left({u}_{1},{u}_{2},{u}_{3}\right)=\left(\begin{array}{ccc}{u}_{1}-1+{\overline{\theta }}_{1}{\partial }_{{\overline{\theta }}_{1}}& -{\partial }_{{\overline{\theta }}_{1}}& -{\partial }_{1}\\ \left({u}_{2}-{u}_{1}\right){\overline{\theta }}_{1}& {u}_{2}-1+{\overline{\theta }}_{1}{\partial }_{{\overline{\theta }}_{1}}& {\partial }_{{\theta }_{1}}-{\partial }_{{\theta }_{2}}-{\overline{\theta }}_{2}{\partial }_{1}\\ 0& 0& {u}_{3}\end{array}\right)$  ${\mathbf{L}}_{2}\left({v}_{1},{v}_{2},{v}_{3}\right)=\left(\begin{array}{ccc}{v}_{1}& -{\partial }_{{\overline{\theta }}_{2}}+{\partial }_{{\overline{\theta }}_{1}}+{\theta }_{1}{\partial }_{2}& -{\partial }_{2}\\ 0& {v}_{2}-{\theta }_{2}{\partial }_{{\theta }_{2}}& -{\partial }_{{\theta }_{2}}\\ 0& \left({v}_{3}-{v}_{2}\right){\theta }_{2}& {v}_{3}-{\theta }_{2}{\partial }_{{\theta }_{2}}\end{array}\right)$  To derive the system of equations which is equivalent to the system ( 4.6 ) written in terms of $\mathbf{r}$  we repeat the same trick with the shift of parameters and obtain
 $\begin{array}{c}\mathbf{r}\cdot \left[{\mathbf{L}}_{1}\left({u}_{1},{u}_{2},{u}_{3}\right)\mathbf{m}+\mathbf{m}{\mathbf{L}}_{2}\left({v}_{1},{v}_{2},{v}_{3}\right)\right]=\left[{\mathbf{L}}_{1}\left({u}_{1},{v}_{2},{u}_{3}\right)\mathbf{m}+\mathbf{m}{\mathbf{L}}_{2}\left({v}_{1},{u}_{2},{v}_{3}\right)\right]\cdot \mathbf{r}\end{array}$ (4.19)
This system results in a simple equations
 $\begin{array}{c}\mathbf{r}\left({\partial }_{{\theta }_{1}}-{\overline{\theta }}_{1}{\partial }_{1}\right)=\left({\partial }_{{\theta }_{1}}-{\overline{\theta }}_{1}{\partial }_{1}\right)\mathbf{r};\mathbf{r}\left({\partial }_{{\overline{\theta }}_{2}}-{\theta }_{1}{\partial }_{2}\right)=\left({\partial }_{{\overline{\theta }}_{2}}-{\theta }_{1}{\partial }_{2}\right)\mathbf{r};\mathbf{r}\left({\partial }_{1}+{\partial }_{2}\right)=\left({\partial }_{1}+{\partial }_{2}\right)\mathbf{r}\end{array}$ (4.20)
 $\begin{array}{c}\mathbf{r}\left({\theta }_{2}{\partial }_{{\theta }_{2}}+\left({z}_{12}+{\theta }_{1}{\overline{\theta }}_{2}\right){\partial }_{2}\right)=\left({\theta }_{2}{\partial }_{{\theta }_{2}}+\left({z}_{12}+{\theta }_{1}{\overline{\theta }}_{2}\right){\partial }_{2}\right)\mathbf{r}\end{array}$ (4.21)
 $\begin{array}{c}\mathbf{r}\left({\overline{\theta }}_{1}{\partial }_{{\overline{\theta }}_{1}}+\left({z}_{12}+{\theta }_{1}{\overline{\theta }}_{2}\right){\partial }_{1}\right)=\left({\overline{\theta }}_{1}{\partial }_{{\overline{\theta }}_{1}}+\left({z}_{12}+{\theta }_{1}{\overline{\theta }}_{2}\right){\partial }_{1}\right)\mathbf{r}\end{array}$ (4.22)
 $\begin{array}{c}\mathbf{r}\left({\overline{\theta }}_{1}{\partial }_{{\overline{\theta }}_{1}}-{\theta }_{2}{\partial }_{{\theta }_{2}}\right)=\left({\overline{\theta }}_{1}{\partial }_{{\overline{\theta }}_{1}}-{\theta }_{2}{\partial }_{{\theta }_{2}}\right)\mathbf{r}\end{array}$ (4.23)
 $\begin{array}{c}\mathbf{r}\left(\left({z}_{12}+{\theta }_{1}{\overline{\theta }}_{2}\right){\partial }_{{\theta }_{2}}+\left({u}_{2}-{u}_{1}\right){\overline{\theta }}_{1}\right)=\left(\left({z}_{12}+{\theta }_{1}{\overline{\theta }}_{2}\right){\partial }_{{\theta }_{2}}+\left({v}_{2}-{u}_{1}\right){\overline{\theta }}_{1}\right)\mathbf{r}\end{array}$ (4.24)
 $\begin{array}{c}\mathbf{r}\left(\left({z}_{12}+{\theta }_{1}{\overline{\theta }}_{2}\right){\partial }_{{\overline{\theta }}_{1}}+\left({v}_{3}-{v}_{2}\right){\theta }_{2}\right)=\left(\left({z}_{12}+{\theta }_{1}{\overline{\theta }}_{2}\right){\partial }_{{\overline{\theta }}_{1}}+\left({v}_{3}-{v}_{2}\right){\theta }_{2}\right)\mathbf{r}\end{array}$ (4.25)
Returning to the system ( 4.18 )(it is the system ( 4.2 ) written in terms of $\mathbf{r}$  ) we note that it is possible to factorize the matrix $diag\left({v}_{1};1;1\right)$  from the right and the matrix $diag\left(1;1;{u}_{3}\right)$  from the left so that ${v}_{1},{u}_{3}$  -dependence disappear from equation. The obtained system of equations is the consequence of the system ( 4.20 )-( 4.25 ). The proof is by direct calculation. Now we are going to the solution of the defining system of equations. The general solution of equations ( 4.20 )-( 4.23 ) has the form $\mathbf{r}=\mathbf{a}+\mathbf{b}\cdot {\overline{\theta }}_{1}{\partial }_{{\overline{\theta }}_{1}}+\mathbf{c}\cdot {\theta }_{2}{\partial }_{{\theta }_{2}}+\mathbf{d}\cdot \left({z}_{12}+{\theta }_{1}{\overline{\theta }}_{2}\right){\partial }_{{\overline{\theta }}_{1}}{\partial }_{{\theta }_{2}}+\mathbf{e}\cdot {\overline{\theta }}_{1}{\theta }_{2}\cdot {\partial }_{{\overline{\theta }}_{1}}{\partial }_{{\theta }_{2}}$  where $\mathbf{a},\mathbf{b},\mathbf{c},\mathbf{d},\mathbf{e}$  are some constants. The equations ( 4.24 ) and ( 4.25 ) fix these constants $\mathbf{a}=\frac{\left({u}_{2}-{u}_{1}\right)\left({v}_{2}-{v}_{3}\right)}{{v}_{2}-{u}_{2}}\cdot \mathbf{d};\mathbf{b}=\left({v}_{2}-{v}_{3}\right)\cdot \mathbf{d};\mathbf{c}=\left({u}_{1}-{u}_{2}\right)\cdot \mathbf{d};\mathbf{e}=\left({u}_{2}-{v}_{2}\right)\cdot \mathbf{d}$  Collect everything together we obtain the expression for the operator ${\mathcal{ℛ}}_{2}$  from the Proposition.

5 Conclusions

We have shown that the general R-matrix can be represented as the product of the simple ”building blocks” – $\mathcal{ℛ}$  -operators. In the first paper [1we have demonstrated how this factorization arises in the simplest situations of the symmetry algebra $s\ell \left(2\right)$  and $s\ell \left(3\right)$  . In the present paper we have showed that the same factorization take place for the R-matrix with supersymmetry algebra $s\ell \left(2|1\right)$  . It seems that this phenomenon is quite general and all results can be generalized to the symmetry algebra $s\ell \left(n\right)$  and to the supersymmetry algebra $s\ell \left(n|m\right)$  .

6 Acknowledgments

I would like to thank D.Karakhanyan, R.Kirschner, G.Korchemsky, P.Kulish and A.Manashov for the stimulating discussions and critical remarks on the different stages of this work. This work was supported by the grant 03-01-00837 of the Russian Foundation for Fundamental Research. Appendix In this Appendix we calculate the matrix elements of the $\mathcal{ℛ}$  -operators and as consequence obtain the matrix elements of R-matrix. It is additional check of the main results and after all we recover the formulae from the paper [8.
All lowest weights in the space ${V}_{{\ell }_{1},{b}_{1}}\otimes {V}_{{\ell }_{2},{b}_{2}}$  are divided on two sets. There are the even lowest weights ${\Phi }_{n}^{±}\equiv {\left({Z}_{12}±\frac{1}{2}{\theta }_{12}{\overline{\theta }}_{12}\right)}^{n};{D}_{1}^{±}{\Phi }_{n}^{±}=0,S{\Phi }_{n}^{±}=\left(n+{\ell }_{1}+{\ell }_{2}\right){\Phi }_{n}^{±},B{\Phi }_{n}^{±}=\left({b}_{1}+{b}_{2}\right){\Phi }_{n}^{±}$  and the odd lowest weights ${\Psi }_{n}^{-}\equiv {\theta }_{12}{Z}_{12}^{n};{\Psi }_{n}^{+}\equiv {\overline{\theta }}_{12}{Z}_{12}^{n};S{\Psi }_{n}^{±}=\left(n+{\ell }_{1}+{\ell }_{2}+\frac{1}{2}\right){\Psi }_{n}^{±},B{\Psi }_{n}^{±}=\left({b}_{1}+{b}_{2}±\frac{1}{2}\right){\Psi }_{n}^{±}$  In this section we shall calculate the action of $\mathcal{ℛ}$  -operators on these lowest weights.
Operator ${\mathcal{ℛ}}_{3}$  We have factorized representation for the operator ${\mathcal{ℛ}}_{3}$  ${\mathcal{ℛ}}_{3}\sim {\mathbb{S}}_{3}^{-1}\cdot {\mathbf{r}}_{3}\cdot {\mathbb{S}}_{3};{\mathbf{r}}_{3}\equiv \frac{\Gamma \left({z}_{1}{\partial }_{1}+{u}_{1}-{v}_{3}+1\right)}{\Gamma \left({z}_{1}{\partial }_{1}+{u}_{1}-{u}_{3}+1\right)}\left(\frac{{u}_{2}-{u}_{3}}{{u}_{3}-{v}_{3}}+{\theta }_{1}{\partial }_{{\theta }_{1}}\right)+\frac{\Gamma \left({z}_{1}{\partial }_{1}+{u}_{1}-{v}_{3}\right)}{\Gamma \left({z}_{1}{\partial }_{1}+{u}_{1}-{u}_{3}+1\right)}{z}_{1}{\partial }_{{\theta }_{1}}{\partial }_{{\overline{\theta }}_{1}}.$  In the explicit form the action of the operators ${\mathbb{S}}_{3}$  and ${\mathbb{S}}_{3}^{-1}$  is ${\mathbb{S}}_{3}\Phi \left({z}_{1},{\theta }_{1},{\overline{\theta }}_{1}|{z}_{2},{\theta }_{2},{\overline{\theta }}_{2}\right)=\Phi \left({z}_{1}+{z}_{2}+\frac{{\theta }_{2}{\overline{\theta }}_{1}}{2}-\frac{{\theta }_{1}{\overline{\theta }}_{2}}{2}-\frac{{\theta }_{1}{\overline{\theta }}_{1}}{2},{\theta }_{1}+{\theta }_{2},{\overline{\theta }}_{1}+{\overline{\theta }}_{2}|{z}_{2},{\theta }_{2},{\overline{\theta }}_{2}\right)$  ${\mathbb{S}}_{3}^{-1}\Phi \left({z}_{1},{\theta }_{1},{\overline{\theta }}_{1}|{z}_{2},{\theta }_{2},{\overline{\theta }}_{2}\right)=\Phi \left({z}_{1}-{z}_{2}+\frac{{\theta }_{1}{\overline{\theta }}_{1}}{2}+\frac{{\theta }_{2}{\overline{\theta }}_{2}}{2}-{\theta }_{2}{\overline{\theta }}_{1},{\theta }_{1}-{\theta }_{2},{\overline{\theta }}_{1}-{\overline{\theta }}_{2}|{z}_{2},{\theta }_{2},{\overline{\theta }}_{2}\right)$  First we calculate the action of ${\mathbb{S}}_{3}$  ${\mathbb{S}}_{3}{\Phi }_{n}^{+}={z}_{1}^{n};{\mathbb{S}}_{3}{\Phi }_{n}^{-}={\left({z}_{1}-{\theta }_{1}{\overline{\theta }}_{1}\right)}^{n}$  then the action of ${\mathbf{r}}_{3}$  ${\mathbf{r}}_{3}\cdot {z}_{1}^{n}=\frac{{u}_{2}-{u}_{3}}{{u}_{3}-{v}_{3}}\frac{\Gamma \left(n+{u}_{1}-{v}_{3}+1\right)}{\Gamma \left(n+{u}_{1}-{u}_{3}+1\right)}\cdot {z}_{1}^{n}$  ${\mathbf{r}}_{3}\cdot {\left({z}_{1}-{\theta }_{1}{\overline{\theta }}_{1}\right)}^{n}=\frac{{u}_{2}-{v}_{3}}{{u}_{3}-{v}_{3}}\frac{\Gamma \left(n+{u}_{1}-{v}_{3}\right)}{\Gamma \left(n+{u}_{1}-{u}_{3}\right)}\cdot {\left({z}_{1}-{\theta }_{1}{\overline{\theta }}_{1}\right)}^{n}+\left({u}_{2}-{u}_{1}\right)\cdot \frac{\Gamma \left(n+{u}_{1}-{v}_{3}\right)}{\Gamma \left(n+{u}_{1}-{u}_{3}+1\right)}\cdot {z}_{1}^{n}$  and finally action of ${\mathbb{S}}_{3}^{-1}$  ${\Phi }_{n}^{+}={\mathbb{S}}_{3}^{-1}\cdot {z}_{1}^{n};{\Phi }_{n}^{-}={\mathbb{S}}_{3}^{-1}{\left({z}_{1}-{\theta }_{1}{\overline{\theta }}_{1}\right)}^{n}$  so that one obtains ${\mathcal{ℛ}}_{3}\cdot {\Phi }_{n}^{+}=\frac{{u}_{2}-{u}_{3}}{{u}_{3}-{v}_{3}}\frac{\Gamma \left(n+{u}_{1}-{v}_{3}+1\right)}{\Gamma \left(n+{u}_{1}-{u}_{3}+1\right)}\cdot {\Phi }_{n}^{+}$  ${\mathcal{ℛ}}_{3}\cdot {\Phi }_{n}^{-}=\frac{{u}_{2}-{v}_{3}}{{u}_{3}-{v}_{3}}\frac{\Gamma \left(n+{u}_{1}-{v}_{3}\right)}{\Gamma \left(n+{u}_{1}-{u}_{3}\right)}\cdot {\Phi }_{n}^{-}+\left({u}_{2}-{u}_{1}\right)\cdot \frac{\Gamma \left(n+{u}_{1}-{v}_{3}\right)}{\Gamma \left(n+{u}_{1}-{u}_{3}+1\right)}\cdot {\Phi }_{n}^{+}$  For the odd lowest weights all is simpler ${\mathbb{S}}_{3}{\Psi }_{n}^{+}={\overline{\theta }}_{1}\cdot {z}_{1}^{n};{\mathbb{S}}_{3}{\Psi }_{n}^{-}={\theta }_{1}\cdot {z}_{1}^{n}$  ${\mathbf{r}}_{3}\cdot {\overline{\theta }}_{1}{z}_{1}^{n}=\frac{{u}_{2}-{u}_{3}}{{u}_{3}-{v}_{3}}\frac{\Gamma \left(n+{u}_{1}-{v}_{3}+1\right)}{\Gamma \left(n+{u}_{1}-{u}_{3}+1\right)}\cdot {\overline{\theta }}_{1}{z}_{1}^{n};{\mathbf{r}}_{3}\cdot {\theta }_{1}{z}_{1}^{n}=\frac{{u}_{2}-{v}_{3}}{{u}_{3}-{v}_{3}}\frac{\Gamma \left(n+{u}_{1}-{v}_{3}+1\right)}{\Gamma \left(n+{u}_{1}-{u}_{3}+1\right)}\cdot {\theta }_{1}{z}_{1}^{n}$  and finally we have ${\mathcal{ℛ}}_{3}\cdot {\Psi }_{n}^{+}=\frac{{u}_{2}-{u}_{3}}{{u}_{3}-{v}_{3}}\frac{\Gamma \left(n+{u}_{1}-{v}_{3}+1\right)}{\Gamma \left(n+{u}_{1}-{u}_{3}+1\right)}\cdot {\Psi }_{n}^{+};{\mathcal{ℛ}}_{3}\cdot {\Psi }_{n}^{-}=\frac{{u}_{2}-{v}_{3}}{{u}_{3}-{v}_{3}}\frac{\Gamma \left(n+{u}_{1}-{v}_{3}+1\right)}{\Gamma \left(n+{u}_{1}-{u}_{3}+1\right)}\cdot {\Psi }_{n}^{-}$  Operator ${\mathcal{ℛ}}_{2}$  We have factorized representation for the operator ${\mathcal{ℛ}}_{2}$  ${\mathcal{ℛ}}_{2}\sim {\mathbb{S}}^{-1}\cdot {\mathbf{r}}_{2}\cdot \mathbb{S}$  ${\mathbf{r}}_{2}\equiv \frac{\left({u}_{2}-{u}_{1}\right)\left({v}_{2}-{v}_{3}\right)}{{v}_{2}-{u}_{2}}+\left({u}_{2}-{u}_{1}\right){\theta }_{12}{\partial }_{{\theta }_{2}}+\left({v}_{2}-{v}_{3}\right){\overline{\theta }}_{12}{\partial }_{{\overline{\theta }}_{1}}+\left({z}_{12}+{\theta }_{1}{\overline{\theta }}_{2}\right){\partial }_{{\overline{\theta }}_{1}}{\partial }_{{\theta }_{2}}+\left({v}_{2}-{u}_{2}\right){\theta }_{12}{\overline{\theta }}_{12}{\partial }_{{\overline{\theta }}_{1}}{\partial }_{{\theta }_{2}}$  In the explicit form the action of the operators $\mathbb{S}$  and ${\mathbb{S}}^{-1}$  is $\mathbb{S}\Phi \left({z}_{1},{\theta }_{1},{\overline{\theta }}_{1}|{z}_{2},{\theta }_{2},{\overline{\theta }}_{2}\right)=\Phi \left({z}_{1}+\frac{{\theta }_{1}{\overline{\theta }}_{1}}{2},{\theta }_{1},{\overline{\theta }}_{1}|{z}_{2}-\frac{{\theta }_{2}{\overline{\theta }}_{2}}{2},{\theta }_{2},{\overline{\theta }}_{2}\right)$  ${\mathbb{S}}^{-1}\Phi \left({z}_{1},{\theta }_{1},{\overline{\theta }}_{1}|{z}_{2},{\theta }_{2},{\overline{\theta }}_{2}\right)=\Phi \left({z}_{1}-\frac{{\theta }_{1}{\overline{\theta }}_{1}}{2},{\theta }_{1},{\overline{\theta }}_{1}|{z}_{2}+\frac{{\theta }_{2}{\overline{\theta }}_{2}}{2},{\theta }_{2},{\overline{\theta }}_{2}\right)$  First we calculate the action of $\mathbb{S}$  $\mathbb{S}{\Phi }_{n}^{+}={\left({z}_{12}+{\theta }_{1}{\overline{\theta }}_{2}+{\theta }_{12}{\overline{\theta }}_{12}\right)}^{n};\mathbb{S}{\Phi }_{n}^{-}={\left({z}_{12}+{\theta }_{1}{\overline{\theta }}_{2}\right)}^{n}$  then the action of ${\mathbf{r}}_{2}$  ${\mathbf{r}}_{2}\cdot {\left({z}_{12}+{\theta }_{1}{\overline{\theta }}_{2}+{\theta }_{12}{\overline{\theta }}_{12}\right)}^{n}=\frac{\left({u}_{2}-{v}_{3}\right)\left({v}_{2}-{u}_{1}\right)}{{v}_{2}-{u}_{2}}\cdot {\left({z}_{12}+{\theta }_{1}{\overline{\theta }}_{2}+{\theta }_{12}{\overline{\theta }}_{12}\right)}^{n}-\left({u}_{1}-{v}_{3}+n\right)\cdot {\left({z}_{12}+{\theta }_{1}{\overline{\theta }}_{2}\right)}^{n}$  ${\mathbf{r}}_{2}\cdot {\left({z}_{12}+{\theta }_{1}{\overline{\theta }}_{2}\right)}^{n}=\frac{\left({u}_{2}-{u}_{1}\right)\left({v}_{2}-{v}_{3}\right)}{{v}_{2}-{u}_{2}}\cdot {\left({z}_{12}+{\theta }_{1}{\overline{\theta }}_{2}\right)}^{n}$  and finally one obtains ${\mathcal{ℛ}}_{2}\cdot {\Phi }_{n}^{+}=\frac{\left({u}_{2}-{v}_{3}\right)\left({v}_{2}-{u}_{1}\right)}{{v}_{2}-{u}_{2}}\cdot {\Phi }_{n}^{+}-\left({u}_{1}-{v}_{3}+n\right)\cdot {\Phi }_{n}^{-};{\mathcal{ℛ}}_{2}\cdot {\Phi }_{n}^{-}=\frac{\left({u}_{2}-{u}_{1}\right)\left({v}_{2}-{v}_{3}\right)}{{v}_{2}-{u}_{2}}\cdot {\Phi }_{n}^{-}$  For the odd lowest weights all is simpler $\mathbb{S}{\Psi }_{n}^{+}={\overline{\theta }}_{12}\cdot {\left({z}_{12}+{\theta }_{1}{\overline{\theta }}_{2}\right)}^{n};\mathbb{S}{\Psi }_{n}^{-}={\theta }_{12}\cdot {\left({z}_{12}+{\theta }_{1}{\overline{\theta }}_{2}\right)}^{n}$  ${\mathbf{r}}_{2}\cdot {\overline{\theta }}_{12}\cdot {\left({z}_{12}+{\theta }_{1}{\overline{\theta }}_{2}\right)}^{n}=\frac{\left({v}_{2}-{u}_{1}\right)\left({v}_{2}-{v}_{3}\right)}{{v}_{2}-{u}_{2}}\cdot {\overline{\theta }}_{12}\cdot {\left({z}_{12}+{\theta }_{1}{\overline{\theta }}_{2}\right)}^{n}$  ${\mathbf{r}}_{2}\cdot {\theta }_{12}\cdot {\left({z}_{12}+{\theta }_{1}{\overline{\theta }}_{2}\right)}^{n}=\frac{\left({u}_{2}-{u}_{1}\right)\left({u}_{2}-{v}_{3}\right)}{{v}_{2}-{u}_{2}}\cdot {\theta }_{12}\cdot {\left({z}_{12}+{\theta }_{1}{\overline{\theta }}_{2}\right)}^{n}$  and finally we have ${\mathcal{ℛ}}_{2}\cdot {\Psi }_{n}^{+}=\frac{\left({v}_{2}-{u}_{1}\right)\left({v}_{2}-{v}_{3}\right)}{{v}_{2}-{u}_{2}}\cdot {\Psi }_{n}^{+};{\mathcal{ℛ}}_{2}\cdot {\Psi }_{n}^{-}=\frac{\left({u}_{2}-{u}_{1}\right)\left({u}_{2}-{v}_{3}\right)}{{v}_{2}-{u}_{2}}\cdot {\Psi }_{n}^{-}.$  Operator ${\mathcal{ℛ}}_{1}$  We have factorized representation for the operator ${\mathcal{ℛ}}_{1}$  ${\mathcal{ℛ}}_{1}\sim {\mathbb{S}}_{1}^{-1}\cdot {\mathbf{r}}_{1}\cdot {\mathbb{S}}_{1};{\mathbf{r}}_{1}\equiv \frac{\Gamma \left({z}_{2}{\partial }_{2}+{u}_{1}-{v}_{3}+1\right)}{\Gamma \left({z}_{2}{\partial }_{2}+{v}_{1}-{v}_{3}+1\right)}\left(\frac{{v}_{1}-{v}_{2}}{{u}_{1}-{v}_{1}}+{\overline{\theta }}_{2}{\partial }_{{\overline{\theta }}_{2}}\right)-\frac{\Gamma \left({z}_{2}{\partial }_{2}+{u}_{1}-{v}_{3}\right)}{\Gamma \left({z}_{2}{\partial }_{2}+{v}_{1}-{v}_{3}+1\right)}{z}_{2}{\partial }_{{\theta }_{2}}{\partial }_{{\overline{\theta }}_{2}}$  and the action of the operators ${\mathbb{S}}_{1}$  and ${\mathbb{S}}_{1}^{-1}$  in explicit form is ${\mathbb{S}}_{1}\Phi \left({z}_{1},{\theta }_{1},{\overline{\theta }}_{1}|{z}_{2},{\theta }_{2},{\overline{\theta }}_{2}\right)=\Phi \left({z}_{1},{\theta }_{1},{\overline{\theta }}_{1}|{z}_{2}+{z}_{1}+\frac{{\theta }_{2}{\overline{\theta }}_{2}}{2}+\frac{{\theta }_{1}{\overline{\theta }}_{2}}{2}-\frac{{\theta }_{2}{\overline{\theta }}_{1}}{2},{\theta }_{2}+{\theta }_{1},{\overline{\theta }}_{2}+{\overline{\theta }}_{1}\right)$  ${\mathbb{S}}_{1}^{-1}\Phi \left({z}_{1},{\theta }_{1},{\overline{\theta }}_{1}|{z}_{2},{\theta }_{2},{\overline{\theta }}_{2}\right)=\Phi \left({z}_{1},{\theta }_{1},{\overline{\theta }}_{1}|{z}_{2}-{z}_{1}-\frac{{\theta }_{1}{\overline{\theta }}_{1}}{2}-\frac{{\theta }_{2}{\overline{\theta }}_{2}}{2}+{\theta }_{2}{\overline{\theta }}_{1},{\theta }_{2}-{\theta }_{1},{\overline{\theta }}_{2}-{\overline{\theta }}_{1}\right)$  First we calculate the action of ${\mathbb{S}}_{1}$  ${\mathbb{S}}_{1}{\Phi }_{n}^{+}=\left(-{z}_{2}{\right)}^{n};{\mathbb{S}}_{1}{\Phi }_{n}^{-}={\left(-{z}_{1}-{\theta }_{2}{\overline{\theta }}_{2}\right)}^{n}$  then the action of ${\mathbf{r}}_{1}$  ${\mathbf{r}}_{1}\cdot \left(-{z}_{2}{\right)}^{n}=\frac{{v}_{1}-{v}_{2}}{{u}_{1}-{v}_{1}}\frac{\Gamma \left(n+{u}_{1}-{v}_{3}+1\right)}{\Gamma \left(n+{v}_{1}-{v}_{3}+1\right)}\cdot \left(-{z}_{2}{\right)}^{n}$  ${\mathbf{r}}_{1}\cdot {\left(-{z}_{2}-{\theta }_{2}{\overline{\theta }}_{2}\right)}^{n}=\frac{{u}_{1}-{v}_{2}}{{u}_{1}-{v}_{1}}\frac{\Gamma \left(n+{u}_{1}-{v}_{3}\right)}{\Gamma \left(n+{v}_{1}-{v}_{3}\right)}\cdot {\left(-{z}_{2}-{\theta }_{2}{\overline{\theta }}_{2}\right)}^{n}+\left({v}_{3}-{v}_{2}\right)\cdot \frac{\Gamma \left(n+{u}_{1}-{v}_{3}\right)}{\Gamma \left(n+{v}_{1}-{v}_{3}+1\right)}\cdot \left(-{z}_{2}{\right)}^{n}$  and finally one obtains ${\mathcal{ℛ}}_{1}\cdot {\Phi }_{n}^{+}=\frac{{v}_{1}-{v}_{2}}{{u}_{1}-{v}_{1}}\frac{\Gamma \left(n+{u}_{1}-{v}_{3}+1\right)}{\Gamma \left(n+{v}_{1}-{v}_{3}+1\right)}\cdot {\Phi }_{n}^{+}$  ${\mathcal{ℛ}}_{1}\cdot {\Phi }_{n}^{-}=\frac{{u}_{1}-{v}_{2}}{{u}_{1}-{v}_{1}}\frac{\Gamma \left(n+{u}_{1}-{v}_{3}\right)}{\Gamma \left(n+{v}_{1}-{v}_{3}\right)}\cdot {\Phi }_{n}^{-}+\left({v}_{3}-{v}_{2}\right)\cdot \frac{\Gamma \left(n+{u}_{1}-{v}_{3}\right)}{\Gamma \left(n+{v}_{1}-{v}_{3}+1\right)}\cdot {\Phi }_{n}^{+}$  For the odd lowest weights we have ${\mathcal{ℛ}}_{1}\cdot {\Psi }_{n}^{+}=\frac{{u}_{1}-{v}_{2}}{{u}_{1}-{v}_{1}}\frac{\Gamma \left(n+{u}_{1}-{v}_{3}+1\right)}{\Gamma \left(n+{v}_{1}-{v}_{3}+1\right)}\cdot {\Psi }_{n}^{+};{\mathcal{ℛ}}_{1}\cdot {\Psi }_{n}^{-}=\frac{{v}_{1}-{v}_{2}}{{u}_{1}-{v}_{1}}\frac{\Gamma \left(n+{u}_{1}-{v}_{3}+1\right)}{\Gamma \left(n+{v}_{1}-{v}_{3}+1\right)}\cdot {\Psi }_{n}^{-}$  Operator R The matrix elements for the $\mathbb{R}$  -operator are obtained from the matrix elements of $\mathcal{ℛ}$  -operators by the formula $\mathbb{R}\left(u;v\right)\sim {\mathcal{ℛ}}_{1}\left({u}_{1};{v}_{1},{u}_{2},{u}_{3}\right){\mathcal{ℛ}}_{2}\left({u}_{1},{u}_{2};{v}_{2},{u}_{3}\right){\mathcal{ℛ}}_{3}\left({u}_{1},{u}_{2},{u}_{3};{v}_{3}\right)$  The result of calculations is the following $\mathbb{R}{\Phi }_{n}^{+}\sim R\cdot \left\{\left({u}_{2}-{u}_{3}\right)\left({v}_{2}-{v}_{1}\right)\cdot \frac{\Gamma \left(n+{u}_{1}-{v}_{3}+1\right)}{\Gamma \left(n+{v}_{1}-{u}_{3}+1\right)}\cdot {\Phi }_{n}^{+}+\left({u}_{2}-{v}_{2}\right)\cdot \frac{\Gamma \left(n+{u}_{1}-{v}_{3}+1\right)}{\Gamma \left(n+{v}_{1}-{u}_{3}\right)}\cdot {\Phi }_{n}^{+}\right\}$  $\mathbb{R}{\Phi }_{n}^{-}\sim R\cdot \left\{\left({u}_{2}-{u}_{1}\right)\left({v}_{2}-{v}_{3}\right)\cdot \frac{\Gamma \left(n+{u}_{1}-{v}_{3}\right)}{\Gamma \left(n+{v}_{1}-{u}_{3}\right)}\cdot {\Phi }_{n}^{-}+C\cdot \frac{\Gamma \left(n+{u}_{1}-{v}_{3}\right)}{\Gamma \left(n+{v}_{1}-{u}_{3}+1\right)}\cdot {\Phi }_{n}^{+}\right\}$  where $C=\left({u}_{2}-{v}_{3}\right)\left({v}_{2}-{u}_{3}\right)\left({u}_{1}-{v}_{1}\right)+\left({v}_{2}-{u}_{1}\right)\left({v}_{1}-{u}_{2}\right)\left({v}_{3}-{u}_{3}\right)+\left({u}_{1}-{v}_{1}\right)\left({u}_{2}-{v}_{2}\right)\left({u}_{3}-{v}_{3}\right)$  $\mathbb{R}{\Psi }_{n}^{+}\sim R\cdot \left({v}_{2}-{u}_{1}\right)\left({v}_{2}-{u}_{3}\right)\cdot \frac{\Gamma \left(n+{u}_{1}-{v}_{3}+1\right)}{\Gamma \left(n+{v}_{1}-{u}_{3}+1\right)}\cdot {\Psi }_{n}^{+}$  $\mathbb{R}{\Psi }_{n}^{-}\sim R\cdot \left({u}_{2}-{v}_{1}\right)\left({u}_{2}-{v}_{3}\right)\cdot \frac{\Gamma \left(n+{u}_{1}-{v}_{3}+1\right)}{\Gamma \left(n+{v}_{1}-{u}_{3}+1\right)}\cdot {\Psi }_{n}^{+}$  We extract the common normalization factor $R\equiv \frac{\left({u}_{2}-{u}_{1}\right)\left({u}_{2}-{u}_{3}\right)}{\left({u}_{1}-{v}_{1}\right)\left({u}_{2}-{v}_{2}\right)\left({u}_{3}-{v}_{3}\right)}.$  After substitution of parameters in explicit form ${u}_{1}=u+{b}_{1}+{\ell }_{1},{v}_{1}=v+{b}_{2}+{\ell }_{2};{u}_{2}=u+2{b}_{1},{v}_{2}=v+2{b}_{2};{u}_{3}=u+{b}_{1}-{\ell }_{1},{v}_{3}=v+{b}_{2}-{\ell }_{2}$  we recover the formulae for the matrix elements of R-operator from the paper [8.
References.

1. S.E. Derkachov ”Factorization of the R-matrix.I.” math.QA/0503396
2. N.Beisert ”The Dilatation Operator of N=4 Super Yang-Mills Theory and Integrability”, Phys.Rept. 405, (2005) 1, hep-th/0407277
3. L.Dolan, C.Nappi and E.Witten ”Yangian symmetry in D=4 superconformal Yang-Mills theory”, hep-th/0401243 L.Dolan and C.Nappi ”Spin models and superconformal Yang-Mills theory”, hep-th/0411020
4. R.Kirschner, Parton interaction in super Yang-Mills theory, JHEP 0407 (2004) 064
5. A.Belitsky, S.Derkachov, G.Korchemsky and A.Manashov, ”Dilatation operator in (super-)Yang-Mills theories on the light-cone”, Nucl.Phys.B 708, (2005) 115, hep-th/0409120 A.Belitsky, S.Derkachov, G.Korchemsky and A.Manashov, ”Quantum Integrability in the super Yang-Mills theory on the light-cone”, Phys.Lett.B 594, (2004) 385, hep-th/0403085
6. M.Scheunert, W.Nahm and V.Rittenberg , J.Math.Phys.18 (1977) 155
7. P.D.Jarvis, H.S.Green , J.Math.Phys.20 (1979) 2115
8. S.Derkachov, D.Karakhanyan, R.Kirschner ”Heisenberg spin chains based on $s\ell \left(2|1\right)$  symmetry”, Nucl.Phys. B 583, (2000) 691
9. M.Marcu, J.Math.Phys.21 (1980) 1277 , J.Math.Phys.21 (1980) 1284
10. D.Arnaudon, C.Chryssomalakos, L.Frappat, J.Math.Phys.36 (1995) 5262 L. Frappat, P. Sorba, A. Sciarrino, DICTIONARY ON LIE SUPERALGEBRAS hep-th/9607161
11. P.P. Kulish and E.K.Sklyanin , ”On the solutions of the Yang-Baxter equation” Zap.Nauchn.Sem. LOMI 95 (1980) 129
12. P.P. Kulish , Zap.Nauchn.Sem. LOMI 145 (1985) 140 , J.Soviet. Math. 35 (1986) 1111, ”Yang-Baxter equation and reflection equations in integrable models”, hep-th/9507070
13. P.P. Kulish and E.K. Sklyanin , ”Quantum spectral transform method. Recent developments”, Lect. Notes in Physics, v 151, (1982) , 61, L.D. Faddeev, ”How Algebraic Bethe Anstz works for integrable model”, Les-Houches lectures 1995, hep-th/9605187, E.K.Sklyanin,”Quantum Inverse Scattering Method.Selected Topics”, in ”Quantum Group and Quantum Integrable Systems” (Nankai Lectures in Mathematical Physics), ed. Mo-Lin Ge,Singapore:World Scientific,1992,pp.63-97; hep-th/9211111
14. H.Frahm, M.P.Pfannmüller and A.M.Tsvelik , Phys. Rev. Lett. 81, (1998) 2116
15. P.P. Kulish, N.Yu.Reshetikhin and E.K.Sklyanin, ”Yang-Baxter equation and representation theory”, Lett.Math.Phys. 5 (1981) 393-403