The permutation operator acts on the basis as follows $$P_{12}{\mathbf{e}}_i\otimes {\mathbf{e}}_k=(-1{)}^{\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=(-\frac{1}{2},-\frac{1}{2});{\Lambda }_3=\Lambda =(\ell ,b)$  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}(u-v)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}(u-v)$$  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}}_{\pm },{\mathbf{v}}_{\pm },{\mathbf{w}}_{\pm }$  are $s\ell \left(2\right|1)$  -generators in the chiral representation ( 2.2 ) and $S,B,S_{\pm },V_{\pm },W_{\pm }$  are generators in the generic representation. The algebra $s\ell \left(2\right|1)$  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}(u-v){\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}(u-v)$$  The explicit expression for the second Lax-operator is the same but now $\mathbf{s},\mathbf{b},{\mathbf{s}}_{\pm },{\mathbf{v}}_{\pm },{\mathbf{w}}_{\pm }$  are $s\ell \left(2\right|1)$  -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\Rightarrow {\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}_{-}\Rightarrow {\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(u_1,u_2,u_3)$  in the functional representation $V_{\Lambda }$  is

$$\begin{array}{c}L(u_1,u_2,u_3)=\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 +(u_2-u_1)\overline{\theta };L_{32}=-\left(z+\frac{\theta \overline{\theta }}{2}\right){\partial }_{\overline{\theta }}-\frac{1}{2}\theta z\partial +(u_3-u_2)\theta$$  $$L_{31}=z^2\partial +z(\theta {\partial }_{\theta }+\overline{\theta }{\partial }_{\overline{\theta }})+(u_1-u_3)z+(u_1+u_3-2u_2)\frac{\theta \overline{\theta }}{2}$$  There exists the useful factorized representation for the Lax-operator

$$\begin{array}{c}L(u_1,u_2,u_3)\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^{\pm }$  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\right|1)$  -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}(u-v){\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}(u-v)$$  change the numeration of the representation spaces ${\Lambda }_2\to {\Lambda }_1=({\ell }_1,b_1);{\Lambda }_3\to {\Lambda }_2=({\ell }_2,b_2)$  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(u-v)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(u-v)$$  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=(z_k,{\theta }_k,{\overline{\theta }}_k)$  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}[Z_1,Z_2]=\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 [1] and [8] . The defining equation which is similar to [1, 8] is

$$\begin{array}{c}{\mathbb{R}}_{{\Lambda }_1{\Lambda }_2}^{-1}(v-u)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}(v-u)\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}(-u)$  .

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}(-u)$  but for the more complicated algebras the action of this automorphism is nontrivial. To proceed in close analogy with [1, 8] we shall use the defining equation ( 3.4 ) so that we derive the expression for the operator ${\mathbb{R}}_{{\Lambda }_1{\Lambda }_2}^{-1}(v-u)$  .

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}(v-u)={\mathbb{P}}_{12}{\check{\mathbb{R}}}_{{\Lambda }_1{\Lambda }_2}(u;v)$  and solve the defining equation for the $\check{\mathbb{R}}$  -operator. The main defining equation for the $\check{\mathbb{R}}$  -operator is $$\check{\mathbb{R}}(u;v)L_1(u_1,u_2,u_3)L_2(v_1,v_2,v_3)=L_1(v_1,v_2,v_3)L_2(u_1,u_2,u_3)\check{\mathbb{R}}(u;v)$$  $$u_1=u+b_1+{\ell }_1,u_2=u+2b_1,u_3=u+b_1-{\ell }_1;v_1=v+b_2+{\ell }_2,v_2=v+2b_2,v_3=v+b_2-{\ell }_2$$  The operator $\check{\mathbb{R}}$  interchanges all parameters in the product of two $L$  -operators and similar to the $s\ell \left(3\right)$  -case $\check{\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(u_1,u_2,u_3)L_2(v_1,v_2,v_3)=L_1(v_1,u_2,u_3)L_2(u_1,v_2,v_3){\mathcal{ℛ}}_1\end{array}$$

(4.1)

$${\mathcal{ℛ}}_1={\mathcal{ℛ}}_1\left(u_1\right|v_1,v_2,v_3);{\mathcal{ℛ}}_1(u_1|v_1,v_2,v_3)={\mathcal{ℛ}}_1(u_1+\lambda |v_1+\lambda ,v_2+\lambda ,v_3+\lambda )$$  and these requirements fix the operator ${\mathcal{ℛ}}_1$  up to overall normalization constant $${\mathcal{ℛ}}_1\sim {\mathbb{S}}_1^{-1}\cdot \left[\frac{\Gamma (z_2{\partial }_2+u_1-v_3+1)}{\Gamma (z_2{\partial }_2+v_1-v_3+1)}\left(f_1+{\overline{\theta }}_2{\partial }_{{\overline{\theta }}_2}\right)-\frac{\Gamma (z_2{\partial }_2+u_1-v_3)}{\Gamma (z_2{\partial }_2+v_1-v_3+1)}{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(u_1,u_2,u_3)L_2(v_1,v_2,v_3)=L_1(u_1,v_2,u_3)L_2(v_1,u_2,v_3){\mathcal{ℛ}}_2\end{array}$$

(4.2)

$${\mathcal{ℛ}}_2={\mathcal{ℛ}}_2(u_1,u_2|v_2,v_3);{\mathcal{ℛ}}_2(u_1,u_2|v_2,v_3)={\mathcal{ℛ}}_1(u_1+\lambda ,u_2+\lambda |v_2+\lambda ,v_3+\lambda )$$  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}+(z_{12}+{\theta }_1{\overline{\theta }}_2){\partial }_{{\overline{\theta }}_1}{\partial }_{{\theta }_2}+(u_2-v_2){\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(u_1,u_2,u_3)L_2(v_1,v_2,v_3)=L_1(u_1,u_2,v_3)L_2(v_1,v_2,u_3){\mathcal{ℛ}}_3\end{array}$$

(4.3)

$${\mathcal{ℛ}}_3={\mathcal{ℛ}}_3(u_1,u_2,u_3|v_3);{\mathcal{ℛ}}_3(u_1,u_2,u_3\left|v_3\right)={\mathcal{ℛ}}_3(u_1+\lambda ,u_2+\lambda ,u_3+\lambda |v_3+\lambda )$$  and these requirements fix the operator ${\mathcal{ℛ}}_3$  up to overall normalization constant $${\mathcal{ℛ}}_3\sim {\mathbb{S}}_3^{-1}\cdot \left[\frac{\Gamma (z_1{\partial }_1+u_1-v_3+1)}{\Gamma (z_1{\partial }_1+u_1-u_3+1)}\left(f_3+{\theta }_1{\partial }_{{\theta }_1}\right)+\frac{\Gamma (z_1{\partial }_1+u_1-v_3)}{\Gamma (z_1{\partial }_1+u_1-u_3+1)}{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 $\check{\mathbb{R}}$  -operator can be factorized as follows $$\check{\mathbb{R}}(u;v)={\mathcal{ℛ}}_1(u_1;v_1,u_2,u_3){\mathcal{ℛ}}_2(u_1,u_2;v_2,u_3){\mathcal{ℛ}}_3(u_1,u_2,u_3;v_3)$$ 

There exist six equivalent ways to represent $\check{\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 $\check{\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\right|1)$  -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(u_1,u_2,u_3)+L_2(v_1,v_2,v_3)\right]=\left[L_1(v_1,u_2,u_3)+L_2(u_1,v_2,v_3)\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(u_1,u_2,u_3)+L_2(v_1,v_2,v_3)\right]=\left[L_1(u_1,v_2,u_3)+L_2(v_1,u_2,v_3)\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(u_1,u_2,u_3)+L_2(v_1,v_2,v_3)\right]=\left[L_1(u_1,u_2,v_3)+L_2(v_1,v_2,u_3)\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\right|1)$  -generators written in a compact form. In explicit notations we have for ${\Lambda }_1=({\ell }_1,b_1)$  and ${\Lambda }_2=({\ell }_2,b_2)$  $$\mathcal{ℛ}:V_{{\Lambda }_1}\otimes V_{{\Lambda }_2}\to V_{{\Lambda }_1^{\prime }}\otimes V_{{\Lambda }_2^{\prime }}$$  $${\mathcal{ℛ}}_1:{\Lambda }_1^{\prime }=({\ell }_1-{\xi }_1,b_1+{\xi }_1);{\Lambda }_2^{\prime }=({\ell }_2+{\xi }_1,b_2-{\xi }_1);{\xi }_1=\frac{u_1-v_1}{2}$$  $${\mathcal{ℛ}}_2:{\Lambda }_1^{\prime }=({\ell }_1,b_1-{\xi }_2);{\Lambda }_2^{\prime }=({\ell }_2,b_2+{\xi }_2);{\xi }_2=u_2-v_2$$  $${\mathcal{ℛ}}_3:{\Lambda }_1^{\prime }=({\ell }_1+{\xi }_3,b_1+{\xi }_3);{\Lambda }_2^{\prime }=({\ell }_2-{\xi }_3,b_2-{\xi }_3);{\xi }_3=\frac{u_3-v_3}{2}$$  The $s\ell \left(2\right|1)$  -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}(v-u)={\mathbb{P}}_{12}{\check{\mathbb{R}}}_{{\Lambda }_1{\Lambda }_2}(u;v)$  is automatically $s\ell \left(2\right|1)$  -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+(\mu -\lambda )\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)+(\nu -\lambda )\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\right|1)$  -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(u_1,u_2,u_3)\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(u_1,u_2,v_3)\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\right|1)$  -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(u_1,u_2,u_3)\mathbf{L}(v_1,v_2,v_3)=L_1(u_1,u_2,v_3)\mathbf{L}(v_1,v_2,u_3)\cdot \mathbf{R}\end{array}$$

(4.9)

where $$\mathbf{L}(v_1,v_2,v_3)\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(u_1,u_2,u_3)+\mathbf{L}(v_1,v_2,v_3)\right]=\left[L_1(u_1,u_2,v_3)+\mathbf{L}(v_1,v_2,u_3)\right]\cdot \mathbf{R}\end{array}$$

(4.10)

$$\begin{array}{c}\mathbf{R}\cdot L_1(u_1,u_2,u_3)\left(\begin{array}{cc}1& 0\\ 0& 1\\ 0& 0\end{array}\right)=L_1(u_1,u_2,v_3)\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(u_1,u_2,u_3)+\mathbf{L}(v_1,v_2,v_3)\right]=\left[L_1(u_1,u_2,v_3)+\mathbf{L}(v_1,v_2,u_3)\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(v_1;v_2-1;1)$  from the right $$\mathbf{R}\cdot L_1(u_1,u_2,u_3)\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(u_1,u_2,v_3)\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(u_1,u_2,u_3)\left(\begin{array}{c}{\partial }_1\\ -V_1^{-}\\ v_3\end{array}\right)=L_1(u_1,u_2,v_3)\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 $[\mathbf{R},D_2^{\pm }]=[\mathbf{R},{\partial }_2]=[\mathbf{R},W_1^{-}]=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 $[\mathbf{R},z_2]=[\mathbf{R},{\theta }_2]=[\mathbf{R},{\overline{\theta }}_2]=[\mathbf{R},D_2^{\pm }]=[\mathbf{R},{\partial }_2]=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(u_1-u_3)+\theta \overline{\theta }(u_3-u_2)\right)=$$ 

$$\begin{array}{c}=\left(z^2\partial +z\left(\theta {\partial }_{\theta }+\overline{\theta }{\partial }_{\overline{\theta }}\right)+z(u_1-v_3)+\theta \overline{\theta }(v_3-u_2)\right)\mathbf{r}\end{array}$$

(4.15)

$$\begin{array}{c}\mathbf{r}\left(z{\partial }_{\overline{\theta }}+(u_2-u_3)\theta \right)=\left(z{\partial }_{\overline{\theta }}+(u_2-v_3)\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 }+(u_2-u_1)\overline{\theta }\right)=\left(-z\left({\partial }_{\theta }+\overline{\theta }\partial \right)+\theta \overline{\theta }{\partial }_{\theta }+(u_2-u_1)\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 $[{\mathbf{S}}^{+},{\mathbf{W}}^{-}]={\mathbf{W}}^{+}$  . The general solution of the equations ( 4.14 ) is $$\mathbf{r}=\mathbf{a}[z\partial ]+\mathbf{b}[z\partial ]\cdot \theta {\partial }_{\theta }+\mathbf{c}[z\partial ]\cdot z{\partial }_{\theta }{\partial }_{\overline{\theta }}$$  The equations ( 4.15 ) and ( 4.16 ) results in recurrence relations $$\mathbf{a}[z\partial ]-\mathbf{a}[z\partial -1]=(u_2-u_3)\cdot \mathbf{c}[z\partial ];\mathbf{b}[z\partial ]=\frac{u_3-v_3}{u_2-u_3}\cdot \mathbf{a}[z\partial ]$$  $$\mathbf{c}[z\partial +1]\cdot (z\partial +u_1-u_3+1)=(z\partial +u_1-v_3)\cdot \mathbf{c}[z\partial ]$$  $$\mathbf{a}[z\partial +1]\cdot (z\partial +u_1-u_3)+(u_2-u_3)\cdot \mathbf{c}[z\partial +1]=(z\partial +u_1-v_3)\cdot \mathbf{a}[z\partial ]$$  $$\mathbf{a}[z\partial ]+\mathbf{b}[z\partial ]\cdot (z\partial +u_1-u_3)-(u_2-u_3)\cdot \mathbf{c}[z\partial ]=\mathbf{a}[z\partial -1]+(z\partial +u_1-v_3)\cdot \mathbf{b}[z\partial -1]$$  The solution of these equations has the form $$\mathbf{a}[z\partial ]=\frac{\Gamma (z\partial +u_1-v_3+1)}{\Gamma (z\partial +u_1-u_3+1)};\mathbf{b}[z\partial ]=\frac{u_3-v_3}{u_2-u_3}\frac{\Gamma (z\partial +u_1-v_3+1)}{\Gamma (z\partial +u_1-u_3+1)}$$  $$\mathbf{c}[z\partial ]=\frac{u_3-v_3}{u_2-u_3}\frac{\Gamma (z\partial +u_1-v_3)}{\Gamma (z\partial +u_1-u_3+1)}$$  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+(\mu -\lambda )\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+(\nu -\lambda )\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\right|1)$  -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(u_1,u_2,u_3)\cdot \mathbf{M}\cdot {\mathbf{l}}_2(v_1,v_2,v_3)={\mathbf{l}}_1(u_1,v_2,u_3)\cdot \mathbf{M}\cdot {\mathbf{l}}_2(v_1,u_2,v_3)\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(u_1,u_2,u_3)\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(v_1,v_2,v_3)\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(u_1,u_2,u_3)\mathbf{m}{\mathbf{L}}_2(v_1,v_2,v_3)={\mathbf{L}}_1(u_1,v_2,u_3)\mathbf{m}{\mathbf{L}}_2(v_1,u_2,v_3)\cdot \mathbf{r}\end{array}$$

(4.18)

where $${\mathbf{L}}_1(u_1,u_2,u_3)=\left(\begin{array}{ccc}{u}_1-1+{\overline{\theta }}_1{\partial }_{{\overline{\theta }}_1}& -{\partial }_{{\overline{\theta }}_1}& -{\partial }_1\\ (u_2-u_1){\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(v_1,v_2,v_3)=\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& (v_3-v_2){\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(u_1,u_2,u_3)\mathbf{m}+\mathbf{m}{\mathbf{L}}_2(v_1,v_2,v_3)\right]=\left[{\mathbf{L}}_1(u_1,v_2,u_3)\mathbf{m}+\mathbf{m}{\mathbf{L}}_2(v_1,u_2,v_3)\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}({\partial }_1+{\partial }_2)=({\partial }_1+{\partial }_2)\mathbf{r}\end{array}$$

(4.20)

$$\begin{array}{c}\mathbf{r}\left({\theta }_2{\partial }_{{\theta }_2}+(z_{12}+{\theta }_1{\overline{\theta }}_2){\partial }_2\right)=\left({\theta }_2{\partial }_{{\theta }_2}+(z_{12}+{\theta }_1{\overline{\theta }}_2){\partial }_2\right)\mathbf{r}\end{array}$$

(4.21)

$$\begin{array}{c}\mathbf{r}\left({\overline{\theta }}_1{\partial }_{{\overline{\theta }}_1}+(z_{12}+{\theta }_1{\overline{\theta }}_2){\partial }_1\right)=\left({\overline{\theta }}_1{\partial }_{{\overline{\theta }}_1}+(z_{12}+{\theta }_1{\overline{\theta }}_2){\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((z_{12}+{\theta }_1{\overline{\theta }}_2){\partial }_{{\theta }_2}+(u_2-u_1){\overline{\theta }}_1\right)=\left((z_{12}+{\theta }_1{\overline{\theta }}_2){\partial }_{{\theta }_2}+(v_2-u_1){\overline{\theta }}_1\right)\mathbf{r}\end{array}$$

(4.24)

$$\begin{array}{c}\mathbf{r}\left((z_{12}+{\theta }_1{\overline{\theta }}_2){\partial }_{{\overline{\theta }}_1}+(v_3-v_2){\theta }_2\right)=\left((z_{12}+{\theta }_1{\overline{\theta }}_2){\partial }_{{\overline{\theta }}_1}+(v_3-v_2){\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(v_1;1;1)$  from the right and the matrix $diag(1;1;u_3)$  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 (z_{12}+{\theta }_1{\overline{\theta }}_2){\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{(u_2-u_1)(v_2-v_3)}{v_2-u_2}\cdot \mathbf{d};\mathbf{b}=(v_2-v_3)\cdot \mathbf{d};\mathbf{c}=(u_1-u_2)\cdot \mathbf{d};\mathbf{e}=(u_2-v_2)\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 [1] we 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\right|1)$  . 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\right|m)$  .

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^{\pm }\equiv {\left(Z_{12}\pm \frac{1}{2}{\theta }_{12}{\overline{\theta }}_{12}\right)}^n;D_1^{\pm }{\Phi }_n^{\pm }=0,S{\Phi }_n^{\pm }=(n+{\ell }_1+{\ell }_2){\Phi }_n^{\pm },B{\Phi }_n^{\pm }=(b_1+b_2){\Phi }_n^{\pm }$$  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^{\pm }=\left(n+{\ell }_1+{\ell }_2+\frac{1}{2}\right){\Psi }_n^{\pm },B{\Psi }_n^{\pm }=\left(b_1+b_2\pm \frac{1}{2}\right){\Psi }_n^{\pm }$$  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 (z_1{\partial }_1+u_1-v_3+1)}{\Gamma (z_1{\partial }_1+u_1-u_3+1)}\left(\frac{u_2-u_3}{u_3-v_3}+{\theta }_1{\partial }_{{\theta }_1}\right)+\frac{\Gamma (z_1{\partial }_1+u_1-v_3)}{\Gamma (z_1{\partial }_1+u_1-u_3+1)}{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 (n+u_1-v_3+1)}{\Gamma (n+u_1-u_3+1)}\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 (n+u_1-v_3)}{\Gamma (n+u_1-u_3)}\cdot {\left(z_1-{\theta }_1{\overline{\theta }}_1\right)}^n+(u_2-u_1)\cdot \frac{\Gamma (n+u_1-v_3)}{\Gamma (n+u_1-u_3+1)}\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 (n+u_1-v_3+1)}{\Gamma (n+u_1-u_3+1)}\cdot {\Phi }_n^{+}$$  $${\mathcal{ℛ}}_3\cdot {\Phi }_n^{-}=\frac{u_2-v_3}{u_3-v_3}\frac{\Gamma (n+u_1-v_3)}{\Gamma (n+u_1-u_3)}\cdot {\Phi }_n^{-}+(u_2-u_1)\cdot \frac{\Gamma (n+u_1-v_3)}{\Gamma (n+u_1-u_3+1)}\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 (n+u_1-v_3+1)}{\Gamma (n+u_1-u_3+1)}\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 (n+u_1-v_3+1)}{\Gamma (n+u_1-u_3+1)}\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 (n+u_1-v_3+1)}{\Gamma (n+u_1-u_3+1)}\cdot {\Psi }_n^{+};{\mathcal{ℛ}}_3\cdot {\Psi }_n^{-}=\frac{u_2-v_3}{u_3-v_3}\frac{\Gamma (n+u_1-v_3+1)}{\Gamma (n+u_1-u_3+1)}\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{(u_2-u_1)(v_2-v_3)}{v_2-u_2}+(u_2-u_1){\theta }_{12}{\partial }_{{\theta }_2}+(v_2-v_3){\overline{\theta }}_{12}{\partial }_{{\overline{\theta }}_1}+(z_{12}+{\theta }_1{\overline{\theta }}_2){\partial }_{{\overline{\theta }}_1}{\partial }_{{\theta }_2}+(v_2-u_2){\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{(u_2-v_3)(v_2-u_1)}{v_2-u_2}\cdot {\left(z_{12}+{\theta }_1{\overline{\theta }}_2+{\theta }_{12}{\overline{\theta }}_{12}\right)}^n-(u_1-v_3+n)\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{(u_2-u_1)(v_2-v_3)}{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{(u_2-v_3)(v_2-u_1)}{v_2-u_2}\cdot {\Phi }_n^{+}-(u_1-v_3+n)\cdot {\Phi }_n^{-};{\mathcal{ℛ}}_2\cdot {\Phi }_n^{-}=\frac{(u_2-u_1)(v_2-v_3)}{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{(v_2-u_1)(v_2-v_3)}{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{(u_2-u_1)(u_2-v_3)}{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{(v_2-u_1)(v_2-v_3)}{v_2-u_2}\cdot {\Psi }_n^{+};{\mathcal{ℛ}}_2\cdot {\Psi }_n^{-}=\frac{(u_2-u_1)(u_2-v_3)}{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 (z_2{\partial }_2+u_1-v_3+1)}{\Gamma (z_2{\partial }_2+v_1-v_3+1)}\left(\frac{v_1-v_2}{u_1-v_1}+{\overline{\theta }}_2{\partial }_{{\overline{\theta }}_2}\right)-\frac{\Gamma (z_2{\partial }_2+u_1-v_3)}{\Gamma (z_2{\partial }_2+v_1-v_3+1)}{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^{+}=(-z_2{)}^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 (-z_2{)}^n=\frac{v_1-v_2}{u_1-v_1}\frac{\Gamma (n+u_1-v_3+1)}{\Gamma (n+v_1-v_3+1)}\cdot (-z_2{)}^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 (n+u_1-v_3)}{\Gamma (n+v_1-v_3)}\cdot {\left(-z_2-{\theta }_2{\overline{\theta }}_2\right)}^n+(v_3-v_2)\cdot \frac{\Gamma (n+u_1-v_3)}{\Gamma (n+v_1-v_3+1)}\cdot (-z_2{)}^n$$  and finally one obtains $${\mathcal{ℛ}}_1\cdot {\Phi }_n^{+}=\frac{v_1-v_2}{u_1-v_1}\frac{\Gamma (n+u_1-v_3+1)}{\Gamma (n+v_1-v_3+1)}\cdot {\Phi }_n^{+}$$  $${\mathcal{ℛ}}_1\cdot {\Phi }_n^{-}=\frac{u_1-v_2}{u_1-v_1}\frac{\Gamma (n+u_1-v_3)}{\Gamma (n+v_1-v_3)}\cdot {\Phi }_n^{-}+(v_3-v_2)\cdot \frac{\Gamma (n+u_1-v_3)}{\Gamma (n+v_1-v_3+1)}\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 (n+u_1-v_3+1)}{\Gamma (n+v_1-v_3+1)}\cdot {\Psi }_n^{+};{\mathcal{ℛ}}_1\cdot {\Psi }_n^{-}=\frac{v_1-v_2}{u_1-v_1}\frac{\Gamma (n+u_1-v_3+1)}{\Gamma (n+v_1-v_3+1)}\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}(u;v)\sim {\mathcal{ℛ}}_1(u_1;v_1,u_2,u_3){\mathcal{ℛ}}_2(u_1,u_2;v_2,u_3){\mathcal{ℛ}}_3(u_1,u_2,u_3;v_3)$$  The result of calculations is the following $$\mathbb{R}{\Phi }_n^{+}\sim R\cdot \left\{(u_2-u_3)(v_2-v_1)\cdot \frac{\Gamma (n+u_1-v_3+1)}{\Gamma (n+v_1-u_3+1)}\cdot {\Phi }_n^{+}+(u_2-v_2)\cdot \frac{\Gamma (n+u_1-v_3+1)}{\Gamma (n+v_1-u_3)}\cdot {\Phi }_n^{+}\right\}$$  $$\mathbb{R}{\Phi }_n^{-}\sim R\cdot \left\{(u_2-u_1)(v_2-v_3)\cdot \frac{\Gamma (n+u_1-v_3)}{\Gamma (n+v_1-u_3)}\cdot {\Phi }_n^{-}+C\cdot \frac{\Gamma (n+u_1-v_3)}{\Gamma (n+v_1-u_3+1)}\cdot {\Phi }_n^{+}\right\}$$  where $$C=(u_2-v_3)(v_2-u_3)(u_1-v_1)+(v_2-u_1)(v_1-u_2)(v_3-u_3)+(u_1-v_1)(u_2-v_2)(u_3-v_3)$$  $$\mathbb{R}{\Psi }_n^{+}\sim R\cdot (v_2-u_1)(v_2-u_3)\cdot \frac{\Gamma (n+u_1-v_3+1)}{\Gamma (n+v_1-u_3+1)}\cdot {\Psi }_n^{+}$$  $$\mathbb{R}{\Psi }_n^{-}\sim R\cdot (u_2-v_1)(u_2-v_3)\cdot \frac{\Gamma (n+u_1-v_3+1)}{\Gamma (n+v_1-u_3+1)}\cdot {\Psi }_n^{+}$$  We extract the common normalization factor $$R\equiv \frac{(u_2-u_1)(u_2-u_3)}{(u_1-v_1)(u_2-v_2)(u_3-v_3)}.$$  After substitution of parameters in explicit form $$u_1=u+b_1+{\ell }_1,v_1=v+b_2+{\ell }_2;u_2=u+2b_1,v_2=v+2b_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\right|1)$  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