\section{Introduction}

There exist just two quantum affine algebras of rank $2$: $\sltwo$ and $\advadva$. The algebra $\sltwo$ and its representation theory are both very well studied. The situation with $\advadva$ is somewhat different, although the algebra itself and its representation are interesting from the standpoints of both mathematics and mathematical physics.

The algebra $\advadva$ first appeared in works on physics. The fundamental representation of the $\Rc$-matrix of $\advadva$ was obtained in \cite{IK} as the $\Rc$-matrix of the quantum version of the Shabat-Mikhailov model, also known as the Izergin-Korepin model. The same representation was later obtained in \cite{J} among fundamental representations of $\Rc$-matrices of other nonexceptional affine Lie algebras. The Bethe ansatz technique was later extended onto $\advadva$ in \cite{T}. Finally, the ideas of the thermodynamic Bethe ansatz for $\advadva$ were developed in \cite{FRS}, and some finite dimensional representations were obtained there.

The algebra $\advadva$ was also investigated from an algebraic standpoint. In \cite{KT}, the Cartan-Weyl basis for $\advadva$ was established, and the universal $\Rc$-matrix was written in terms of infinite products of elements of the Cartan-Weyl basis. A classification of finite dimensional representations of algebra $\advadva$ was obtained in \cite{CP} by means of the Drinfeld polynomials. Later, an integral formula for the universal $\Rc$-matrix for $\advadva$ appeared in \cite{DK}, where $\advadva$ was treated as a topological Hopf algebra with the Drinfeld comultiplication. In the same work, the Serre relations in $\advadva$ were represented in terms of zeros and poles of products of the Drinfeld currents. Finally, integral representations for factors of the universal $\Rc$-matrix for $\advadva$ were derived in \cite{KS}, where $\advadva$ was endowed with the standard Hopf structure.

Quantum affine algebras allow three different realizations with different Hopf structures. The first one is the ``standard'' realization given by the Chevalley generators and relations determined by the corresponding Cartan matrix. The standard realization has a small number of generators, but is unfortunately very difficult to use in applications. The second realization, the Drinfeld ``new realization'', was first established in \cite{D} by means of generating functions (the Drinfeld ``currents'') and relations on them. The Drinfeld realization allows using methods of complex analysis. Moreover, it makes it possible to give the classification of the finite dimensional representations of the quantum affine algebra. Finally, the third one is the $RLL$ realization, based on the Faddeev-Reshetikhin-Takhtajan-Semenov-Tian-Shansky approach, where generators are combined into $L$-operators, satisfying the famous Yang-Baxter equation (see \cite{FRT}, \cite{RS}). The simplicity of comultiplication in the $RLL$ realization allows to construct new representations as tensor products of already known ones. The $RLL$ realization is therefore widely used in physical models.

It is universally acknowledged that the three realizations are isomorphic, although precise proofs hardly exist for any algebra other then of the $U_q(\widehat sl_n)$ type. For $U_q(\widehat sl_n)$, an isomorphism between the standard and the Drinfeld realizations was established in \cite{DF}, whereas links between the Drinfeld and the $RLL$ approaches were studied in \cite{DK2}.

In the case of quantum twisted affine algebras, the $RLL$ realization requires additional relations. Although it is believed that the three realizations are also isomorphic for the twisted algebras, there is no full understanding of what the exact isomorphism should look like. An isomorphism between the standard and the Drinfeld realizations for $\advadva$ was established in \cite{KS}, and an isomorphism between the $RLL$ and the Drinfeld realizations was partly derived in \cite{YZ}. A more complete bibliography can be found in \cite{H}.

In this work, we obtain a full description of the three realizations of algebra $\advadva$ (without grading element and with zero central charge), the isomorphisms between them, and the links between the three Hopf structures. We also pay special attention to the following fact: each realization has a ``minimal'' set of generators (or almost minimal in the twisted case) and an extended set of generators. These two sets should be linked by an analogue of the PBW-theorem if there exists one. In the standard realization, these sets are the Chevalley generators and the Cartan-Weyl basis. In the Drinfeld realization, they are the Drinfeld currents and the so-called ``composite currents'' (see \cite{DK3}). Finally, in the $RLL$-realization, the Gaussian coordinates immediately above or below the diagonal form the ``minimal'' set, and all the Gaussian coordinates form the extended set. We note that these extended sets of generators are crucial for calculations of the universal weight function. Here, we obtain a link between projections of the composite currents and the Gaussian coordinates for $\advadva$ as was done for $U_q(\widehat\glg_n)$ in \cite{KP}.

The paper is organized as follows. In the first section we introduce the three realizations. In Section $2$, we obtain the isomorphisms between them. In Section $3$, we represent elements of the extended set of Gaussian coordinates in terms of composite currents.

\section{Realizations of $\advadva$}

{\sl For simplicity of exposition, we consider the algebra $\advadva$ without grading element and with zero central charge. All the statements also hold for an algebra with grading element and arbitrary central charge, but the formulas become more involved.}

\subsection{Drinfeld realization}

Let $\AD$ be the associative algebra generated by elements
$$
e_n, f_n,\;\, n\in\Z,\quad a_n,\;\, n\in\Z\setminus\hc{0},\quad {\rm and} \quad k^{\pm 1},
$$
subject to certain commutation relations. The relations are given as formal power series identities for generating functions (currents)
$$
e(z)=\sum_{k\in\Z}e_kz^{-k},\quad f(z)=\sum_{k\in\Z}f_kz^{-k},
\quad K^\pm(z)=k^{\pm 1}\exp\left(\pm(q-q^{-1})\sum_{n>0}a_{\pm n}z^{\mp n}\right)
$$
as follows:
\pbea
(z-q^2w)(qz+w)e(z)e(w) &=& (q^2z-w)(z+qw)e(w)e(z),                           \\
(q^2z-w)(z+qw)f(z)f(w) &=& (z-q^2w)(qz+w)f(w)f(z),                           \\
K^+(z)e(w)K^+(z)^{-1} &=& \a(w/z)e(w),                                       \\
K^+(z)f(w)K^+(z)^{-1} &=& \a(w/z)^{-1}f(w),                                  \\
K^-(z)e(w)K^-(z)^{-1} &=& \a(z/w)^{-1}e(w),                                  \\
K^-(z)f(w)K^-(z)^{-1} &=& \a(z/w)f(w),                                       \\
K^\pm(z)K^\pm(w) &=& K^\pm(w)K^\pm(z),                                       \\
K^-(z)K^+(w) &=& K^+(w)K^-(z),                                               \\
e(z)f(w)-f(w)e(z) &=& \frac1{q-q^{-1}}\hr{\d(z/w)K^+(w)-\d(z/w)K^-(z)}.
\peea
where
$$
\a(x)=\dfrac{(q^2-x)(q^{-1}+x)}{(1-q^2x)(1+q^{-1}x)},
$$
and $\d\left(\dfrac zw\right)$ is a formal Laurent series, given by
$$
\d(z/w) = \sum\limits_{n\in\Z}(z/w)^n.
$$

The generating functions $e(z)$ and $f(z)$ also satisfy the cubic Serre relations (see \cite{D}):
\pbea
\Sym_{z_1,z_2,z_3} \hr{q^{-3}z_1 - (q^{-2}+q^{-1})z_2 + z_3} e(z_1)e(z_2)e(z_3) &=& 0, \\
\Sym_{z_1,z_2,z_3} \hr{q^{-3}z_1^{-1} - (q^{-2}+q^{-1})z_2^{-1} + z_3^{-1}} f(z_1)f(z_2)f(z_3) &=& 0, \\
\Sym_{z_1,z_2,z_3} \hr{q^3 z_1^{-1} - (q^{2}+q^{})z_2^{-1} + z_3^{-1}} e(z_1)e(z_2)e(z_3) &=& 0, \\
\Sym_{z_1,z_2,z_3} \hr{q^3 z_1 - (q^{2}+q^{})z_2 + z_3} f(z_1)f(z_2)f(z_3) &=& 0.
\peea

The Hopf algebra structure on $\AD$ can be defined as follows:
\pbea
\D_\Dc(e(z)) &=& e(z)\otimes 1+K^-(z)\otimes e(z), \\
\D_\Dc(f(z)) &=& 1\otimes f(z)+f(z)\otimes K^+(z), \\
\D_\Dc(K^\pm(z)) &=& K^\pm(z)\otimes K^\pm(z),     \\
S_\Dc(e(z)) &=& -(K^-(z))^{-1}e(z),                \\
S_\Dc(f(z)) &=& -f(z)(K^+(z))^{-1},                \\
S_\Dc(K^\pm(z)) &=& (K^\pm(z))^{-1},               \\
\e_\Dc(e(z)) &=& 0,                                \\
\e_\Dc(f(z)) &=& 0,                                \\
\e_\Dc(K^\pm(z)) &=& 1.
\peea
where $\D_\Dc$, $\e_\Dc$ and $S_\Dc$ are the comultiplication, the counit and the antipode maps, respectively. We call $\D_\Dc$ a \emph{Drinfeld comultiplication}. Here, we must note that $\AD$ is a topological bialgebra and $\D_\Dc$ is a map from $\AD$ to a topological extension of its tensor square (see \cite[Section 2]{EKP} for details).

\subsection{Chevalley realization}

Another realization of $\advadva$ is given in terms of Chevalley generators. Let the algebra $\ACh$ be the associative algebra generated by elements $e_{\pm\a}$, $e_{\pm(\d-2\a)}$, $k^{\pm1}_{\a}$, $k^{\pm1}_{\d-2\a}$, satisfying the following relations:
\pbea
k_{\a}e_{\pm\a}k^{-1}_{\a} = q^{\pm1}e_{\pm\a}, &\quad&
k_{\a}e_{\pm(\d-2\a)}k^{-1}_{\a} = q^{\mp 2}e_{\pm(\d-2\a)},
\\
k_{\d-2\a}e_{\pm{\a}}k^{-1}_{\d-2\a} = q^{\mp 2}e_{\pm{\a}}, &\quad&
k_{\d-2\a}e_{\pm{(\d-2\a)}}k^{-1}_{\d-2\a} = q^{\pm4}e_{\pm{(\d-2\a)}},
\\
k_{\a}^2 k_{\d-2\a} = 1, &\quad&
\hs{e_{\pm\a},e_{\mp(\d-2\a)}} = 0,
\\
\hs{e_\a,e_{-\a}} = \dfrac{k_\a-k_\a^{-1}}{q-q^{-1}}, &\quad&
\hs{e_{\d-2\a},e_{-(\d-2\a)}} = \dfrac{k_{\d-2\a}-k_{\d-2\a}^{-1}}{q-q^{-1}},
\\
(\ad_{q}e_{\pm\a})^5 e_{\pm(\d-2\a)} = 0, &\quad&
(\ad_{q}e_{\pm(\d-2\a)})^2 e_{\pm\a} = 0,
\peea
where
\pbea
(\ad_{q}e_{\pm\a})(x) &=& e_{\pm\a}x - k_\a^{\pm 1}xk_\a^{\mp 1}e_{\pm\a},
\\
(\ad_{q}e_{\pm(\d-2\a)})(x) &=& e_{\pm(\d-2\a)}x - k_{\d-2\a}^{\pm 1}xk_{\d-2\a}^{\mp 1}e_{\pm(\d-2\a)}.
\peea

The Hopf algebra structure associated with this realization can be defined by
\pbea
\D(e_{\a})=e_{\a}\otimes 1+k_{\a}\otimes e_{\a},
&\quad&
\D(e_{\d-2\a})=e_{\d-2\a}\otimes 1+k_{\d-2\a}\otimes e_{\d-2\a},
\\
\D(e_{-\a})=1\otimes e_{-\a}+e_{-\a}\otimes k^{-1}_{\a},
&\quad&
\D(e_{-(\d-2\a)})=1\otimes e_{-(\d-2\a)}+e_{-(\d-2\a)}\otimes k^{-1}_{\d-2\a},
\\
\D(k_{\a})\,=\,k_{\a}\otimes k_{\a},
&\quad&
\D(k_{\d-2\a})\,=\,k_{\d-2\a}\otimes k_{\d-2\a},
\\
\e(e_{\pm\a})=0,
&\quad&
\e(e_{\pm(\d-2\a)})=0,
\\
\e(k^{\pm1}_{\a})=1,
&\quad&
\e(k^{\pm1}_{\d-2\a})=1,
\\
S(e_{\a})=-k^{-1}_{\a}e_{\a},
&\quad&
S(e_{\d-2\a})=-k^{-1}_{\d-2\a}e_{\d-2\a},
\\
S(e_{-\a})=-e_{-\a}k_{\a},
&\quad&
S(e_{-(\d-2\a)})=-e_{-(\d-2\a)}k_{\d-2\a},
\\
S(k^{\pm1}_{\a}) =  k^{\mp1}_{\a},
&\quad&
S(k^{\pm1}_{\d-2\a}) =  k^{\mp1}_{\d-2\a},
\peea
where $\D$, $\e$ and $S$ are the comultiplication, the counit and the antipode maps, respectively. We call $\D$ the \emph{standard comultiplication}. We also consider the opposite comultiplication
$$
\D^{op}=\si\circ\D, \qquad\text{where}\qquad \si(u\otimes v) = v\otimes u.
$$

Now, let us recall the construction of the Cartan-Weyl basis, obtained for $\advadva$ in \cite{KT}\footnote{To adapt the results of \cite{KT} to our case, one should replace $q$-commutators with $q^{-1}$-commutators in the construction of the Cartan-Weyl basis.}. The twisted affine algebra $A_2^{(2)}$ has the root system $\D=\D_+\cup\D_-$, where
$$
\D_+ = \a\cup\hc{n\d,\;\pm\a+n\d,\;\pm2\a+(2n-1)\d\;|\;n\in\N}, \qquad \text{and} \qquad \D_- = -\D_+.
$$
The roots $\ga\in\D_+$ are called \emph{positive}, the roots $\ga\in\hc{n\d,\;n\in\Z}$ are called \emph{imaginary}, and the roots $\ga\in\D\setminus\hc{n\d,\;n\in\Z}$ are called \emph{real}. We define a scalar product on $\D$ by
$$
\hr{\a,\a} = 1, \qquad \hr{\a,\d} = \hr{\d,\d} = 0.
$$
Now, let
$$
q_\ga = q^{(\ga,\ga)}, \qquad (a)_q = \frac{q^a-1}{q-1}, \qquad (n)_q! = (1)_q\cdot\dots\cdot(n)_q,
$$
and
$$
\exp_q(x) = \sums{n\ge0}\frac{x^n}{(n)_q!} = 1+x+\frac{x^2}{(2)_q!}+\dots\,.
$$

We fix the normal ordering on $\D$
$$
\a,\; 2\a+\d,\; \a+\d, \; 2\a+3\d, \; \a+2\d, \; \dots,
\; \d, \; 2\d, \; 3\d, \; \dots, \;
2\d-\a, \; 3\d-2\a, \; \d-\a, \; \d-2\a.
$$
In accordance with the procedure of the construction of the Cartan-Weyl basis described in \cite{KT}, we set
\beq
\begin{split}
&e_{\d-\a}=a[e_{\a},e_{\d-2\a}]_{q^{-1}}, \qquad
e'_{\d}=b[e_{\a},e_{\d-\a}]_{q^{-1}}, \qquad
e_{\a+n\d}=b[e_{\a+(n-1)\d},e'_{\d}]_{q^{-1}}, \\
&e_{\d-\a+n\d}=b[e'_{\d},e_{\d-\a+(n-1)\d}]_{q^{-1}}, \qquad
e'_{n\d}=b[e_{\a+(n-1)\d},e_{\d-\a}]_{q^{-1}}, \\
&e_{2\a+(2n+1)\d}=a[e_{\a+n\d},e_{\a+(n+1)\d}]_{q^{-1}}, \qquad
e_{\d-2\a+2n\d}=a[e_{\d-\a+n\d}, e_{\d-\a+(n-1)\d}]_{q^{-1}},
\end{split}
\eeq
where
$$ \label{ab}
a=\frac{1}{\sqrt{q+q^{-1}}}, \qquad b=\frac{1}{\sqrt{q+1+q^{-1}}},
$$
and
$$
[e_\b,e_\ga]_{q^{-1}} = e_\b e_\ga-q^{-\hr{\b,\ga}}e_\ga e_\b.
$$
holds for all $\b, \ga \in \D$. Finally, we define \emph{imaginary} roots $e_{n\d}$ by means of the Schur polynomials:
\beq
e'_{n\d}=\sum_{p_{1}+2p_{2}+\dots +np_{n}=n}
\frac{\hr{(q-q^{-1})b^{-1}}^{\sum p_{i}-1}}{p_{1}!\ldots p_{n}!}
e_{\d}^{p_1}e_{2\d}^{p_2}\dots e_{n\d}^{p_n}.
\eeq
Monomials in the elements $e_{n\d}$, $e_{n\d\pm\a}$ and $k_\a^{\pm1}$ form a linear basis of $\ACh$.

\subsection{$RLL$ realization}

Let
\beq \label{R}
R(x) =
\begin{pmatrix}
1 & 0 & 0 & 0 & 0 & 0 & 0 & 0 & 0 \\
0 & a & 0 & d & 0 & 0 & 0 & 0 & 0 \\
0 & 0 & b & 0 & e & 0 & p & 0 & 0 \\
0 & f & 0 & a & 0 & 0 & 0 & 0 & 0 \\
0 & 0 & g & 0 & c & 0 & e & 0 & 0 \\
0 & 0 & 0 & 0 & 0 & a & 0 & d & 0 \\
0 & 0 & q & 0 & g & 0 & b & 0 & 0 \\
0 & 0 & 0 & 0 & 0 & f & 0 & a & 0 \\
0 & 0 & 0 & 0 & 0 & 0 & 0 & 0 & 1
\end{pmatrix},
\eeq
where
\begin{align*}
a &= \frac{q(x-1)}{q^2x-1},         & b &= \frac{q^2(qx+1)(x-1)}{(q^2x-1)(q^3x+1)},       &
c &= \frac{q(x-1)}{q^2x-1} + \frac{(q^2-1)(q^3+1)x}{(q^2x-1)(q^3x+1)},                    \\
d &= \frac{q^2-1}{(q^2x-1)},        & e &= \frac{(q^2-1)(x-1)q^{1/2}}{(q^2x-1)(q^3x+1)},  &
p &= \frac{(q^2-1)(q^3x+qx-q+1)}{(q^2x-1)(q^3x+1)},                                       \\
f &= \frac{(q^2-1)x}{(q^2x-1)},     & g &= \frac{(1-q^2)x(x-1)q^{5/2}}{(q^2x-1)(q^3x+1)}, &
q &= \frac{(q^2-1)x(q^3x-q^2x+q^2+1)}{(q^2x-1)(q^3x+1)}.
\end{align*}
We consider the generating functions
\pbeq
\begin{split}
l_{ij}^\pm(z) = \suml{n=0}{\infty}l_{ij}^\pm[\pm n]z^{\mp n}, &\qquad 1 \le i,j \le 3 \\
l_{ij}^+[0]=l_{ji}^-[0] = 0, &\qquad 1\le i<j \le 3
\end{split}
\peeq
and the shifted generating functions
$$
\tilde l_{ij}^\pm(z) = q^{\frac12(i-j)} l_{4-j,4-i}^\pm(-q^{-3}z).
$$
We define the \emph{$L$-operators} $L^\pm(z)$ and their shifts $\tilde L^\pm(z)$ as
$$
L^\pm(z) = \hr{l_{ij}^\pm(z)}_{i,j=1}^3, \qquad \tilde L^\pm(z) = \hr{\tilde l_{ij}^\pm(z)}_{i,j=1}^3.
$$
Finally, we recall the notion of $q$-determinant:
$$
\det\nolimits_q(L^\pm(z)) =
\sums{\tau\in \Sg_3}(-q)^{\sgn(\tau)}l^\pm_{1\tau(1)}(z)l^\pm_{2\tau(2)}(q^2z)l^\pm_{3\tau(3)}(q^4z),
$$
where $\Sg_n$ is the permutation group on $n$ elements.

Define $\AR$ as the associative algebra with the generators $l_{ij}^\pm[\pm n],\,n\in\N,\,1\le i,j \le 3$ and $l_{ij}^+[0]$, $l_{ji}^-[0]$ for $1\le j\le i\le 3$, subject to the following relations:
\beq \label{q-det}
\det\nolimits_q(L^\pm(z)) = 1,
\eeq
\beq \label{symmetry}
L^\pm(z)\tilde L^\pm(z) = I_3,
\eeq
\beq \label{pm_identity}
l_{ii}^+[0]l_{ii}^-[0] = l_{ii}^-[0]l_{ii}^+[0] = 1,
\eeq
\beq
\begin{split} \label{RLL}
R\hr{\frac zw}L_1^\pm(z)L_2^\pm(w) &= L_2^\pm(w)L_1^\pm(z)R\hr{\frac zw}, \\
R\hr{\frac zw}L_1^+(z)L_2^-(w) &= L_2^-(w)L_1^+(z)R\hr{\frac zw}.
\end{split}
\eeq
where $L_1(z) = L(z)\otimes I_3$, $L_2(z) = I_3 \otimes L(z)$, and $I_n$ denotes the $n\times n$ identity matrix.

The Hopf algebra structure on $\AR$ is given by
$$
\D_R(l_{ij}^\pm) = \suml{k=1}{3}l_{ik}^\pm \otimes l_{kj}^\pm, \qquad
S_R(L^\pm) = (L^\pm)^{-1}, \qquad
\e_R(L^\pm) = I_3,
$$
where $\D_R$, $S_R$ and $\e_R$ are the comultiplication, the antipode and the counit maps, respectively.

The $L$-operators admit the Gaussian decomposition
\beq \label{Gauss}
L^\pm(z) =
\begin{pmatrix}
1 & f^\pm_1(z) & f^\pm_{13}(z) \\
0 & 1          & f^\pm_2(z)    \\
0 & 0          & 1
\end{pmatrix}
\begin{pmatrix}
k^\pm_1(z) & 0          & 0          \\
0          & k^\pm_2(z) & 0          \\
0          & 0          & k^\pm_3(z)
\end{pmatrix}
\begin{pmatrix}
1             & 0          & 0 \\
e^\pm_1(z)    & 1          & 0 \\
e^\pm_{31}(z) & e^\pm_2(z) & 1
\end{pmatrix},
\eeq
where
\begin{align*}
e^\pm_i(z) &= \sums{n\ge0}e^\pm_i[n]z^{\mp n},\quad i=1,2,   &e^\pm_{31}(z) &= \sums{n\ge0}e^\pm_{31}[n]z^{\mp n}, \\
f^\pm_i(z) &= \sums{n\ge0}f^\pm_i[n]z^{\mp n},\quad i=1,2,   &f^\pm_{13}(z) &= \sums{n\ge0}f^\pm_{13}[n]z^{\mp n}, \\
k^\pm_i(z) &= \sums{n\ge0}k_i^\pm[n]z^{\mp n},\quad i=1,2,3.
\end{align*}
and
$$
e^-_1[0] = e^-_2[0] = e^-_{31}[0] = 0, \qquad f^+_1[0] = f^+_2[0] = f^+_{13}[0] = 0.
$$
We have thus defined the Gaussian generators
$$
f^\pm_1[n],\;f^\pm_2[n],\;f^\pm_{13}[n], \qquad
e^\pm_1[n],\;e^\pm_2[n],\;e^\pm_{31}[n], \qquad
k^\pm_1[n],\;k^\pm_2[n],\;k^\pm_3[n], \quad n\ge0,
$$
subject to relations \rf{q-det} -- \rf{RLL}. Here $f^\pm_1[n],\;f^\pm_2[n],\;e^\pm_1[n],\;e^\pm_2[n]$ form the ``minimal'' set of generators, while $f^\pm_{13}[n]$ and $e^\pm_{31}[n]$ form the extended one.

\section{Links between realizations}

\subsection{Universal $\Rc$-matrix}

Recall that the universal $\Rc$-matrix of a quasitriangular Hopf algebra $\Ac$ is an element $\Rc\in\Ac\otimes\Ac$, subject to the following relations:
\beq \label{intertwining}
\D^{op}(x)=\Rc\D(x)\Rc^{-1}
\eeq
for any $x\in\Ac$, and
\beq \label{quasi}
(\D\otimes\id)\Rc=\Rc_{13}\Rc_{23}, \qquad (\id\otimes\Delta)\Rc=\Rc_{13}\Rc_{12},
\eeq
where
$$
\Rc_{12}=\Rc\otimes1, \qquad \Rc_{23}=1\otimes \Rc, \qquad \Rc_{13} = (\si\otimes\,\id)(\Rc_{23})
$$
are elements of $\Ac\otimes\Ac\otimes\Ac$.

Now we fix the representation $\pi_z$ defined on Chevalley generators as
\begin{align*}
\pi_{z}\hr{e_\a} &=
\begin{pmatrix}
0&q^\frac14&0 \\
0&0&-q^{-\frac14} \\
0&0&0
\end{pmatrix},
&
\pi_{z}\hr{e_{\d-2\a}} &=
\begin{pmatrix}
0&0&0 \\
0&0&0 \\
\sqrt{q+q^{-1}}&0&0
\end{pmatrix}z, \\
\pi_{z}\hr{e_{-\a}} &=
\begin{pmatrix}
0&0&0 \\
q^{-\frac14}&0&0 \\
0&-q^\frac14&0
\end{pmatrix},
&
\pi_{z}\hr{e_{-\d+2\a}} &=
\begin{pmatrix}
0&0&\sqrt{q+q^{-1}} \\
0&0&0 \\
0&0&0
\end{pmatrix}z^{-1}, \\
\pi_{z}\hr{k_\a} &=
\begin{pmatrix}
q&0&0 \\
0&1&0 \\
0&0&q^{-1}
\end{pmatrix},
&
\pi_{z}\hr{k_{\d-2\a}} &=
\begin{pmatrix}
q^{-2}&0&0 \\
0&1&0 \\
0&0&q^2
\end{pmatrix}.
\end{align*}
After the representation is chosen, we obtain the $R$-matrix
\footnote{
The scaling of the basis of the $\advadva$-module is chosen in such a way that $R$-matrix $R(x)$ becomes symmetric with respect to the secondary diagonal. In that basis, the $R$-matrix coincides with the one obtained in \cite{J}.
}
$$
R\hr{\frac zw} = \hr{\pi_z\otimes\pi_w}\Rc
$$
directly from the definition of $\Rc$ by applying $\pi_z\otimes\pi_w$ to relation \rf{intertwining} and substituting the Chevalley generators for $x$. Condition \rf{intertwining} defines $R(x)$ up to a scalar.

\subsection{Chevalley and Drinfeld realizations}

Following \cite{KT}, we represent the universal $\Rc$-matrix of $\ACh$ as
\footnote{The formula for the $\Rc$-matrix differs from the formula in \cite{KT} because $q$-commutators are replaced with $q^{-1}$-commutators. Moreover, our formula differs from the one in \cite{KS}; the $\Rc$-matrix in \cite{KS} is inverse to the $\Rc$-matrix here.}:
\beq \label{R_matrix}
\begin{split}
\Rc = \Rc_-\Kc q^{-h \otimes h}\Rc_+, &\qquad
\Rc_- = \hr{\prodl{\ga<\d}\ra \exp_{q_\ga}
\hr{(q^{-1}-q)e_\ga \otimes e_{-\ga}}}, \\
\Kc = \exp\hr{-\sums{n>0}{S_{n}}}, &\qquad
\Rc_+ = \hr{\prodl{\ga>\d}\ot\exp_{q_\ga}
\hr{(q^{-1}-q)e_\ga k_\a^{-(\a,\ga)}\otimes k_\a^{(\a,\ga)} e_{-\ga}}},
\end{split}
\eeq
where $\ga$ ranges over real positive roots, $h$ is defined by
$$
q^{\pm h} = k^{\pm 1},
$$
and $S_{n}$ is given by the formula
$$
S_{n}=\frac{n\hr{q-q^{-1}}^2\hr{q+1+q^{-1}}(e_{n\d}\otimes e_{-n\d})}
{\hr{q^n-q^{-n}}\hr{q^{n}+(-1)^{n+1}+q^{-n}}}.
$$

Now, we can find the scalar prefactor for $R(x)$. Applying $\pi_z \otimes \pi_w$ to $\Kc \cdot q^{-h \otimes h}$, we obtain the scalar
\beq \label{prefactor}
q^{-1}\exp\hr{-\sums{n>0}\frac{\hr{q^n-q^{-n}}x^n}{\hr{q^n+(-1)^{n+1}+q^{-n}}n}}.
\eeq
Therefore, $R(x)$ is equal to matrix \rf{R} times scalar prefactor \rf{prefactor}.

\begin{rem}
Scalar \rf{prefactor} expands into the product
$$
q^{-1}\prodl{k=0}{\infty}\hr{\frac{(1+q^{3k+2}x)(1-q^{3k+3}x)}{(1+q^{3k}x)(1-q^{3k+1}x)}}^{(-1)^k},
$$
in the domain $|q|<1$ and into the product
$$
q^{-1}\prodl{k=0}{\infty}\hr{\frac{(1+q^{-3k}x)(1-q^{-(3k+1)}x)}{(1+q^{-(3k+2)}x)(1-q^{-(3k+3)}x)}}^{(-1)^k},
$$
in the domain $|q|>1$.
\end{rem}

\begin{theor} \label{theor Ch-Dr}
The isomorphism between associative algebras $\ACh$ and $\AD$ is established by the following map:
\beq \label{Ch-Dr}
\begin{array}{rclrcl}
k_{\d-2\a}  &\mapsto& k^{-2},                             &\qquad k_{\a} &\mapsto& k,     \\
e_{\d-2\a}  &\mapsto& a(qf_1f_0-f_0f_1)k^{-2},            &\qquad e_{\a} &\mapsto& e_{0},   \\
e_{-\d+2\a} &\mapsto& ak^2(q^{-1}e_0e_{-1}-e_{-1}e_0),    &\qquad e_{-\a} &\mapsto& f_{0},
\end{array}
\eeq
where $a$ is defined in \rf{ab}.
\end{theor}
\begin{proof}
A straightforward verification and the inductive construction of the Cartan-Weyl basis show that the map \rf{Ch-Dr} above is a surjective homomorphism. More precisely, every Drinfeld generator can be obtained as
\beq \label{CW-Dr}
\begin{array}{rcllrcll}
e_{\a+n\d}      &\mapsto& e_n,    &n\ge0,    & \qquad -k_\a^{-1}e_{\a-n\d} &\mapsto& e_{-n}, &n>0,                 \\
e_{-\a-n\d}     &\mapsto& f_{-n}, &n\ge0,    & \qquad -e_{-\a+n\d}k_\a     &\mapsto& f_n,    &n>0,                 \\
k_\a            &\mapsto& k,      &          & \qquad -b^{-1}e_{n\d}       &\mapsto& a_n,    &n\in\Z\setminus\{0\}.
\end{array}
\eeq
Reversing all the arrows in map \rf{CW-Dr}, we obtain a map inverse to \rf{Ch-Dr}. Finally, map \rf{CW-Dr} with reversed arrows is also an epimorphism, and the statement of the theorem follows.
\end{proof}
\begin{rem} \label{rem} The above map is an isomorphism of associative algebras. It preserves the counit but does not respect comultiplication maps $\D$ and $\D_\Dc$.
\end{rem}
As for the link between the Drinfeld and the standard comultiplications, we have the following
\begin{prop}
The tensor $\RR_-$ is a cocycle for $\D_\Dc$ so that for any $x\in\ACh$
$$
\D(x)=\hr{\Rc_-}^{-1}\D_\Dc(x)\Rc_-.
$$
\end{prop}
The proof is given in \cite{KS}. Another expression for $\Rc_-$ in terms of the Drinfeld generators can also be found there.

\subsection{From Chevalley to $RLL$ realization}

According to \cite{D}, the universal $R$-matrix of any quasitriangular Hopf algebra satisfies the quantum Yang-Baxter equation
\beq \label{qybe}
\Rc_{12}\Rc_{13}\Rc_{23}=\Rc_{23}\Rc_{13}\Rc_{12}.
\eeq
\begin{theor} \label{theor RLL-Ch}
An isomorphism between the associative algebras $\AR$ and $\ACh$ can is given by the maps
\beq \label{Loper}
\begin{split}
L^-(z) &\mapsto (\pi_z\otimes\id)\Rc, \\
L^+(z) &\mapsto (\pi_z\otimes\id)\Rc_{21}^{-1}.
\end{split}
\eeq
\end{theor}

\begin{proof}
Applying $\;\pi_z\otimes\pi_w\otimes\id\;$, $\;\pi_z\otimes\id\otimes\pi_w\;$ and $\;\id\otimes\pi_z\otimes\pi_w\;$ to relation \rf{qybe}, we derive commutation relations \rf{pm_identity} and \rf{RLL}. Using formulas \rf{R_matrix} and \rf{Loper}, we obtain the Gaussian decomposition \rf{Gauss}.

The matrix $R(-q^3)$ has a one-dimensional kernel, which is spanned by vector
$$
v = q^{-\frac12} v_1\otimes v_3 + v_2\otimes v_2 + q^\frac12 v_3\otimes v_1,
$$
where $\hc{v_1,v_2,v_3}$ is the basis of the representation $\pi_z$. It follows that $v$ is an eigenvector of operator $L_1(z)L_2(-q^{-3}z)$. Let $\la(z)$ be an eigenvalue of $v$. One can check that $\hr{\pi_z\otimes\pi_{-q^{-3}z}}\D(x)$ vanishes on $v$ for every $x\in\ACh$ (it is sufficient to varify the vanishing condition for the Chevalley generators only). Applying $\hr{\pi_z\otimes\pi_{-q^{-3}z}\otimes\id}$ to the first equation in \rf{quasi}, we obtain $\la(z)=1$. Hence, $v$ is stable under $L_1(z)L_2(-q^{-3}z)$, which turns out to be equivalent to the relation \rf{symmetry}:
$$
\begin{pmatrix}
l_{11}(z) & l_{12}(z) & l_{13}(z) \\
l_{21}(z) & l_{22}(z) & l_{23}(z) \\
l_{31}(z) & l_{32}(z) & l_{33}(z)
\end{pmatrix}
\begin{pmatrix}
          l_{33}(-q^{-3}z)  &  q^{-\frac12} l_{23}(-q^{-3}z)  &  q^{-1}       l_{13}(-q^{-3}z) \\
q^\frac12 l_{32}(-q^{-3}z)  &               l_{22}(-q^{-3}z)  &  q^{-\frac12} l_{12}(-q^{-3}z) \\
q         l_{31}(-q^{-3}z)  &  q^\frac12    l_{21}(-q^{-3}z)  &               l_{11}(-q^{-3}z)
\end{pmatrix}
=I_3.
$$
Next, it is possible to check that
$$
\det\nolimits_q(L^\pm(z)) = k_1^\pm(z)k_2^\pm(q^{-2}z)k_3^\pm(q^{-4}z).
$$
Following the procedure of the construction of the Cartan-Weyl basis we get
$$
\pi_z\hr{e_{\pm n\d}} = b\frac{q^n-q^{-n}}{n(q-q^{-1})}
\begin{pmatrix}
-z^{\pm n} & 0 & 0 \\
0 & (q^2z)^{\pm n} - (-qz)^{\pm n} & 0 \\
0 & 0 & (-q^3z)^{\pm n}
\end{pmatrix}, \quad n>0.
$$
Applying $\pi_z$ to either first or second tensor component of $\Kc$, we arrive at
\beq \label{add2}
\begin{split}
k_1^\pm(z) &= k_3^\pm(-q^{-3}z)^{-1}, \\
k_2^\pm(z) &= k_3^\pm(-q^{-1}z)k_3^\pm(q^{-2}z)^{-1},
\end{split}
\eeq
which implies the desired relation \rf{q-det}.

Since all the Chevalley generators are in the image of map \rf{Loper}, we obtain an epimorphism from $\AR$ to $\ACh$ (and, therefore, to $\AD$). In Theorem \rfs{theor Dr-RLL} below, we establish an epimorphism from $\AD$ to $\AR$, which is inverse to \rf{Loper}. Hence, the theorem is proved modulo the result of Theorem \rfs{theor Dr-RLL}.
\end{proof}

\begin{prop}
Let $\phi(x)\in\ACh$ be the image of $x\in\AR$ under isomorphism \rf{Loper}. Then
$$
\phi\hr{\D_R(x)} = \si\circ\D(\phi(x)).
$$
\end{prop}

\begin{proof}
The proof can be obtained by applying $\pi_z\otimes\id\otimes\id$ to the relation
$$
(\id\otimes\D^{op})\Rc = \Rc_{12}\Rc_{13},
$$
which is equivalent to \rf{quasi}.
\end{proof}

\subsection{From $RLL$ to Drinfeld realization}

First of all, we express all generators of the $RLL$ realization in terms of $e_1^\pm(z)$, $f_1^\pm(z)$, $k_1^\pm(z)$. The expressions
\beq \label{k2k3}
\begin{split}
k_2^\pm(z) &= k_1^\pm(-qz)k_1^\pm(q^2z)^{-1}, \\
k_3^\pm(z) &= k_1^\pm(-q^3z)^{-1}
\end{split}
\eeq
follow directly from \rf{add2}. Relations \rf{RLL} imply
\begin{multline*}
(q^2-1)z(q^3z+w)L_{23}^\pm(z)L_{33}^\pm(w) + q(z-w)(q^3z+w)L_{33}^\pm(z)L_{23}^\pm(w) = \\ = (q^2z-w)(q^3z+w)L_{23}^\pm(w)L_{33}^\pm(z),
\end{multline*}
which results in
$$
q(z-w)k_3^\pm(z)f_2^\pm(w) = (q^2z-w)f_2^\pm(w)k_3^\pm(z) - (q^2-1)zf_2^\pm(z)k_3^\pm(z).
$$
Setting $w=q^2z$, we get
$$
f_2^\pm(q^2z) = q^{-1}k_3^\pm(z)^{-1}f_2^\pm(z)k_3^\pm(z).
$$
On the other hand, relation \rf{symmetry} implies
$$
f_1^\pm(z) = -q^{-\frac12}k_3^\pm(-q^{-3}z)^{-1}f_2^\pm(-q^{-3}z)k_3^\pm(-q^{-3}z).
$$
We hence obtain
\beq \label{e2f2}
\begin{split}
f_2^\pm(z) &= -q^{-\frac12}f_1^\pm(-qz), \\
e_2^\pm(z) &= -q^{\frac12}e_1^\pm(-qz),
\end{split}
\eeq
where the expression for $e_2^\pm(z)$ is derived in the similar way.
The expressions for $e_{31}^\pm(z)$ and $f_{13}^\pm(z)$ are considered in Section $3$.

\begin{rem}
It is also possible to use the explicit formula for $L^\pm(z)$ to derive relations \rf{e2f2}, as we did for relations \rf{add2}.
\end{rem}

Now let
\pbeq
\begin{split}
E_i(z) &= e_i^+(z) - e_i^-(z), \\
F_i(z) &= f_i^+(z) - f_i^-(z).
\end{split}
\peeq
Then the equalities
\beq \label{E2F2}
\begin{split}
E_2(z) &= -q^\frac 12 E_1(-qz), \\
F_2(z) &= -q^{-\frac 12} F_1(-qz).
\end{split}
\eeq
hold. Moreover, relations \rf{q-det}--\rf{RLL} imply
\beq \label{kk}
\begin{split}
k_1^\pm(z)k_1^\pm(w) &= k_1^\pm(w)k_1^\pm(z), \\
k_1^-(z)k_1^+(w) &= k_1^+(w)k_1^-(z),
\end{split}
\eeq

\beq \label{kEkkFk}
\begin{split}
k_1^\pm(z)E_1(w)k_1^\pm(z)^{-1} &= \frac{q(z-w)}{q^2z-w}E_1(w), \\
k_1^\pm(z)^{-1}F_1(w)k_1^\pm(z) &= \frac{z-q^2w}{q(z-w)}F_1(w),
\end{split}
\eeq

\beq \label{EEFF}
\begin{split}
(qz+w)(z-q^2w)E_1(z)E_1(w) &= (z+qw)(q^2z-w)E_1(w)E_1(z), \\
(z+qw)(q^2z-w)F_1(z)F_1(w) &= (qz+w)(z-q^2w)F_1(w)F_1(z),
\end{split}
\eeq

\beq \label{EF}
\hs{E_1(z),F_1(w)} = (q-q^{-1})\d\hr{\frac zw}\hr{k_1^+(w)k_2^+(w)^{-1} - k_1^-(z)k_2^-(z)^{-1}}.
\eeq

\begin{theor} \label{theor Dr-RLL}
An isomorphism\footnote{As in remark \rf{rem}, the mapping above is an isomorphism of associative algebras. It preserves the counit but does not respect comultiplication maps $\D$ and $\D_\Dc$.} between the algebras $\AD$ and $\AR$ is given by the following maps:
\beq \label{Dr-RLL}
\begin{array}{lcl}
q^{-\frac 14} \hr{q-q^{-1}}e(qz)   &\mapsto &E_1(z),                       \\
q^{\frac 14}\hr{q-q^{-1}}f(qz)     &\mapsto &F_1(z),                       \\
K^\pm(qz)                          &\mapsto &k_1^\pm(z)k_2^\pm(z)^{-1}.
\end{array}
\eeq
\end{theor}

\begin{proof}
Formulas \rf{kk}--\rf{EF} imply that the map \rf{Dr-RLL} is a homomorphism. For all $n\ge0$ Gaussian coordinates $k_1^\pm[n]$ can be obtained step by step from the product $k_1^\pm(z)k_2^\pm(z)^{-1}$, which together with relations \rf{k2k3} and \rf{E2F2} imply that map \rf{Dr-RLL} is an epimorphism. Finally, one can check that the composition of maps \rf{Loper}, \rf{Ch-Dr} and \rf{Dr-RLL} is the identity map.
\end{proof}

\section{Link with composite Drinfeld currents}
Here we revise some definitions, given in \cite{KS} and then express the generating functions $f_{13}^\pm(z)$ and $e_{31}^\pm(z)$ in terms of the Drinfeld currents.

Recall that in any quantum affine algebra there exist two types of Borel subalgebras. Borel subalgebras of the first type come from the Drinfeld realization. Let $U_F$ denote the subalgebra of $\advadva$ generated by $k^{\pm 1},f_n,\,n\in\Z;\;a_n,\,n>0$, and let $U_E$ denote the subalgebra of $\advadva$ generated by $k^{\pm 1},e_n,\,n\in\Z;\;a_n,\,n<0$. The ``current'' Borel subalgebra $U_F$ contains subalgebra $U_f$ generated by $f_n,\,n\in\Z$. The ``current'' Borel subalgebra $U_E$ contains subalgebra $U_e$ generated by $e_n,\,n\in\Z$.

Borel subalgebras of the second type are obtained via the Chevalley realization. Let $U_q(\bgt_+)$ and $U_q(\bgt_-)$ denote a pair of subalgebras of $\advadva$ generated by
\pbe
e_{\a},\;\,e_{\d-2\a},\;\,k^{\pm1}_{\a} \qquad\text{and}\qquad
e_{-\a},\;\,e_{-(\d-2\a)},\;\,k^{\pm1}_{\a}
\pee
respectively. In terms of the Drinfeld realization, these subalgebras are generated by
\pbe
k^{\pm 1},\;\,e_0,\;\,qf_1f_0-f_0f_1 \qquad\text{and}\qquad k^{\pm 1},\;\,f_0,\;\,q^{-1}e_0e_{-1}-e_{-1}e_0.
\pee
respectively.

Let $U_F^+$, $U_f^-$, $U_e^+$ and $U_E^-$ denote the following intersections
of Borel subalgebras:

\pbea
U_f^-=U_F\cap U_q(\bgt_-), &\qquad U_F^+=U_F\cap U_q(\bgt_+), \\
U_e^+=U_E\cap U_q(\bgt_+), &\qquad U_E^-=U_E\cap U_q(\bgt_-).
\peea
The superscript sign indicates the Borel subalgebra $U_q(\bgt_\pm)$ containing the given algebra, and the subscript letter indicates ``current'' Borel subalgebra $U_F$ or $U_E$ that it is intersected with. These letters are capitals if the subalgebra contains imaginary root generators $a_n$ and are lower case otherwise.

Let $P^+$ and $P^-$ denote projection operators such that for any $f_+\in U_F^+$ and for any $f_-\in U_f^-$
\be\label{Pdef}
P^+(f_-f_+)=\e(f_-)f_+, \qquad P^-(f_-f_+)=f_-\e(f_+).
\ee
Also let $P^{*+}$ and $P^{*-}$ denote projecton operators such that for any $e_+\in U_e^+$ and for any $e_-\in U_E^-$
\be\label{Pdefa}
P^{*+}(e_+e_-)=e_+\e(e_-), \qquad P^{*-}(e_+e_-)=\e(e_+)e_-.
\ee

Finally, let us introduce composite currents $s(z)$ and $r(z)$ as follows:
\begin{align*}
s(z) &= \res\limits_{w=-q^{-1} z}f(z)f(w)\frac{dw}w, \\
r(z) &= \res\limits_{w=-qz}e(w)e(z)\frac{dw}w.
\end{align*}

\begin{rem}
The currents $s(z)$ and $\tilde s(z)$ were defined in \cite{KS}. The definition of $s(z)$ coincides with the one given above, while $\tilde s(z) = -s(-qz)$ or equivalently
$$
\tilde s(z) = \res\limits_{w=-qz}f(w)f(z)\frac{dw}w.
$$
An involution $\iota$ of the algebra $\advadva$ was also defined in \cite{KS}:
$$
\iota(e_n) = f_{-n}, \qquad \iota(f_n) = e_{-n}, \qquad \iota(a_n) = a_{-n}, \qquad \iota(K_0) = K_0^{-1}.
$$
Therefore, the current $r(z)$ satisfies the equality
$$
r(z) = \iota\hr{-\tilde s\hr{z^{-1}}}.
$$
\end{rem}

\begin{prop} \label{gauss_currents}
Isomorphism \rf{Dr-RLL} maps the Gaussian coordinates $e_{31}^\pm(z)$ and $f_{13}^\pm(z)$ to the expressions
\beq \label{e31f13}
\begin{split}
e_{31}^\pm(z) &\mapsto
\mp(1-q) \hr{ q^{-1} \hs{P^{*\pm}(e(-q^2z)),e_0}_{q^{-1}} + \hs{-q^2ze_{-1},P^{*\pm}(e(-q^2z))\mp e_0}_{q^{-1}} }, \\
f_{13}^\pm(z) &\mapsto
\pm(1-q) \hr{ q\hs{P^\pm(f(-q^2z)),f_0}_{q^{-1}} + \hs{(-q^2z)^{-1}f_1,P^\pm(f(-q^2z))\pm f_0}_{q^{-1}} }.
\end{split}
\eeq
\end{prop}

\begin{proof}
We consider only the case of $f_{13}^+(z)$ because the other cases are similar.
Relations \rf{RLL} imply
\begin{multline*}
(q^2-1)z(q^3z+w)l_{13}^+(z)l_{22}^-(w) + q(z-w)(q^3z+w)l_{23}^+(z)l_{12}^-(w) = \\ =
(q^2-1)w(q^3z+w)l_{13}^-(w)l_{22}^+(z) + q(z-w)(q^3z+w)l_{12}^-(w)l_{23}^+(z).
\end{multline*}
Assuming $w=0$, we get
$$
(q-q^{-1})l_{13}^+(z) = \hs{l_{12}^-(0),l_{23}^+(z)},
$$
or, equivalently,
$$
(q-q^{-1})f_{13}^+(z)k_3^+(z) = f_1^-[0]f_2^+(z)k_3^+(z) - f_2^+(z)k_3^+(z)f_1^-[0].
$$
Using \rf{RLL} once again, we obtain
$$
k_3^+(z)f_1^-[0] = qf_1^-[0]k_3^+(z) + (q-q^{-1})q^{\frac12}f_2^+(z)k_3^+(z),
$$
and therefore derive
$$
f_{13}^+(z) = \frac1{q-q^{-1}}\hr{q^\frac12 f_1^+(-qz)f_1^-[0] - q^{-\frac12}f_1^-[0]f_1^+(-qz)} - q^{-\frac12}f_1^+(-qz)^2.
$$

Now, let us recall a result from \cite{KS}: on one hand, we have
\begin{multline*}
P(f(z_1)f(z_2))= f^+(z_1)\hr{f^+(z_2) - \dfrac{q^2-1}{q^2-z_2/z_1}f^+(z_1)} - \\ -
\dfrac{(1+q^3)(1-z_2/z_1)}{(1+q)(q^2-z_2/z_1)(1+qz_2/z_1)}P(s(z_1)),
\end{multline*}
and, on the other hand,
$$
P(f(z_1)f(z_2)) = \frac{(q^2-z_1/z_2)(q^{-1}+z_1/z_2)}{(1-q^2z_1/z_2)(1+q^{-1}z_1/z_2)}P(f(z_2)f(z_1)).
$$
Equating the right-hand sides of these expressions, multiplying them by $z_2$, and letting $z_2$ tend to infinity, we obtain
$$
(q-q^{-1)}P(f(z))^2 = \hr{f_1P(f(z))-q^{-1}P(f(z))f_1}z^{-1} + \frac{1+q^3}{q(1+q)}P(s(z)) + \hr{f_1f_0-q^{-1}f_0f_1}z^{-1}.
$$
Finally, since isomorphism \rf{Dr-RLL} maps $f_1^\pm(-qz)$ to $\pm q^{\frac14}(q-q^{-1})P^\pm(f(-q^2z))$, we derive the desired formula \rf{e31f13}.
\end{proof}

\begin{theor} \label{theor}
Isomorphism \rf{Dr-RLL} links the Gaussian coordinates $f_{13}^+(z)$, $e_{31}^+(z)$ and the composite Drinfeld currents $s(z)$, $r(z)$ by the following expressions:
\beq
\begin{split}
P^+(s(-q^2z))    &= \frac{1}{(q-q^{-1})(q-1+q^{-1})}f_{13}^+(z), \\
P^-(s(-q^2z))    &= \frac{1}{(q-q^{-1})(q-1+q^{-1})}\hr{f_1^-(z)f_2^-(z) - f_{13}^-(z)}, \\
P^{*+}(r(-q^2z)) &= \frac{1}{(q-q^{-1})(q-1+q^{-1})}e_{31}^+(z), \\
P^{*-}(r(-q^2z)) &= \frac{1}{(q-q^{-1})(q-1+q^{-1})}\hr{e_2^-(z)e_1^-(z) - e_{31}^-(z)}.
\end{split}
\eeq
\end{theor}

\begin{proof}
The statement for $s(z)$ follows directly from Proposition \rfs{gauss_currents} and Theorem~1 in \cite{KS}. For $r(z)$ the involution $\iota$ should also be used.
\end{proof}

\begin{rem}
Since $f_1^-(z)f_2^-(z) - f_{13}^-(z)$ and $e_2^-(z)e_1^-(z) - e_{31}^-(z)$ are the corner elements in the Gaussian decomposition of the inverse $L$-operator $\hr{L^-(z)}^{-1}$, it might make sense to use the elements of the decomposition of the inverse $L$-operator $\hr{L^-(z)}^{-1}$ as the Gaussian coordinates instead of the elements of the decomposition of usual $L$-operator $L^-(z)$.
\end{rem}

\section*{Acknowledgements}
The author thanks Sergey Khoroshkin for the invaluable help and his constant interest to the work. The author is also grateful to Stanislav Pakuliak for many fruitful discussions. The author was supported by RFBR grant 08-01-00667, joint CNRS-RFBR grant 09-01-93106-NCNIL, and by Federal Agency for Science and Innovations of Russian Federation under contract 14.740.11.0081.


\section{Introduction}
\label{sec_one}

There is a rich literature on $ ^6 He $  based on an effective
$ \alpha-n-n $ 3-body problem \cite{rep}. Besides pair interactions also
ad hoc 3-body forces are used.
 The Pauli principle is approximately incorporated by projecting out
"Pauli forbidden states" for the neutrons inside the $\alpha$-particle
wave function. While that approach catches presumably the
halo structure of the two loosely bound neutrons clearly because of
its strongly restricted ansatz
 it is not suited  to  probe modern nucleon-nucleon and
three-nucleon forces, like the ones derived recently through effective
field theory and based on chiral symmetry \cite{epelbaum}.

Nevertheless some approaches already exist which directly attack the
6-nucleon problem and beyond
 with realistic nuclear forces, namely in the method of
no-core shell model (NCSM) \cite{ncsm} and the Greens function
\cite{gfmc} Monte Carlo treatment.
 In \cite{ncsm} chiral two-nucleon and three-nucleon forces were used
 and applied to
 $^7Li$ with some under binding. In \cite{gfmc} the AV18
 nucleon-nucleon interaction  and a Urbana three-nucleon
   force was used
again leading to  some under binding, now for $^6He$.
 There is also the stochastic variational
   Monte Carlo method \cite{varga} which, however, still applied
   simplified forces.

The exploration of chiral
 forces goes on, also including explicitely the $\Delta$-degree
of freedom \cite{evgeny}, which calls for an increased effort to
establish rigorous approaches beyond $A=4$.
 The achievments in \cite{ncsm} and \cite{gfmc}
 demonstrate that a direct treatment of 6 nucleons is feasible on
present day computers and therefore we felt that  another approach,
the exact formulation within the Yakubovsky equations, is timely.
About 20 years ago an analogous step turned out to be very fruitful,
namely the exact formulation of the $\alpha$-particle within
the Yakubovsky scheme \cite{kamada1,kamada2}. This pioneering
study opened the way to a nowadays
standard treatment \cite{kamada3,kamada4,nogga,nogga1}
     and allows the inclusion of the most modern two- and three-nucleon
     forces and even first estimates of four-nucleon forces \cite{nogga2}.

In section \ref{sec_two} we apply the
Yakubovsky equations \cite{yakubovsky}  to
the six-body problem using the basic notation for sub clusters \cite{mybook}.
        In section \ref{sec_three} we add
the identity of the nucleons which leads
to a set of 5 coupled Yakubovsky equations related to 5 different
sequential sub clusterings of 6 particles.
 In view of the expectation for the dominant  structure of $^6He$,
namely an $\alpha$-core and two loosely bound neutrons, we stop the sequential
sub clustering with 3 fragments, though the additional step with
two fragments could be easily performed.

In section \ref{sec_four} and the Appendices we provide technicalities which we
consider  useful for a numerical performance.
Finally we summarize in section \ref{sec_five}.

\section{The Yakubovsky approach to 6 particles}
\label{sec_two}

We use the standard notation $ a_n $ to denote the various members of
n fragments for a total of N particles. Here N=6 and $ a_2, a_3, a_4,
a_5 $ denote two-, three-,  up to 5-body fragmentations. 5-body
fragmentations necessarily have one pair left and thus $ a_5 $ can be
used to point to a specific pair.

$ a_3 \subset a_2 $ means that the three-body fragments $ a_3 $
consist of  sub clusters out of the two  fragments in
$ a_2 $ or $ a_3  \supset a_4 $ means that the sub clusters $ a_3 $
when broken up lead to the sub clusters  $ a_4 $.

The bound state $ \Psi$ obeys the homogeneous equation
\begin{eqnarray}
\Psi = G_0 \sum_{ a_5} V_{a_5} \Psi
\label{tot}
\end{eqnarray}
where $ G_0 $ is the 6-particle free Greens operator and $ V_{a_5} $
is a pair force.  The first step is the summation of each pair force to
infinite order. Defining
  \begin{eqnarray}
  \psi_{a_5} \equiv G_0 V_{a_5} \Psi\label{fad}
  \end{eqnarray}
one obtains like for 3 particles
\begin{eqnarray}
\psi_{a_5}  =   G_0 t_{a_5} \sum_{b_5} \overline{\delta_{ a_5 b_5}} \psi_{b_5}
\end{eqnarray}
where $ t_{a_5}$ is a two-body t-operator obeying the Lippmann
Schwinger equation
\begin{eqnarray}
t_{a_5} = V_{a_5} + V_{a_5} G_0 t_{a_5}
\end{eqnarray}
and $ \overline{\delta_{ a_5 b_5}} \equiv 1 - \delta_{ a_5 b_5} $.

Next one defines  new components
\begin{eqnarray}
\psi_{a_5 a_4}  \equiv   G_0 t_{a_5} \sum_{ b_5 \subset a_4}
\overline{\delta_{ a_5 b_5}} \psi_{b_5}\label{5}
\end{eqnarray}
where all pairs $ a_5, b_5$ are sub clusters of the fragments in $ a_4
$. Clearly
\begin{eqnarray}
\psi_{a_5}  =   \sum_{ a_4 \supset a_5} \psi_{ a_5 a_4}\label{6}
\end{eqnarray}

That relation (\ref{6})  is used to obtain a closed set of equations
for $ \psi_{a_5 a_4} $:
\begin{eqnarray}
\psi_{a_5 a_4}  =   G_0 t_{a_5} \sum_{ b_5 \subset a_4}
\overline{\delta_{ a_5 b_5}} \sum_{ b_4 \supset b_5} \psi_{b_5 b_4}
\end{eqnarray}

One separates now the components $  \psi_{a_5 a_4} $ for a given $ a_4
$ from the rest :
  \begin{eqnarray}
& & \psi_{a_5 a_4} - G_0 t_{a_5} \sum_{ b_5 \subset a_4}
\overline{\delta_{ a_5 b_5}} \psi_{b_5 a_4}
 =   G_0 t_{a_5} \sum_{ b_5 \subset a_4} \overline{\delta_{ a_5
    b_5}}
\sum_{ b_4 \supset b_5} \overline{\delta_{ a_4 b_4}} \psi_{b_5 b_4}
\label{7}
\end{eqnarray}

Let us define for a fixed $ a_4 $ the column vectors $ \psi^{a_4} $
and $\psi^{(a_4)} $ with the components
\begin{eqnarray}
 (\psi^{a_4})_{a_5} \equiv  \psi_{a_5 a_4}\label{comp}
 \end{eqnarray}
and
\begin{eqnarray}
 (\psi^{(a_4)})_{b_5} \equiv  \sum_{ b_4 \supset b_5}
\overline{\delta_{ a_4 b_4}} \psi_{b_5 b_4}\label{comp1}
 \end{eqnarray}

Then  introducing the matrix $ C^{ a_4} $ with the elements
$ C^{a_4}_{ a_5 b_5} \equiv t_{a_5} \overline{\delta_{ a_5 b_5}} $
Eq.(\ref{7}) reads
  \begin{eqnarray}
( 1- G_0 C^{a_4}) \psi^{a_4} = G_0 C^{a_4} \psi^{(a_4)}
\end{eqnarray}
 or
  \begin{eqnarray}
 \psi^{a_4} & = & ( 1- G_0 C^{a_4})^{-1} G_0 C^{a_4}  \psi^{(a_4)}
 \equiv G_0 T^{a_4} \psi^{(a_4)}\label{9}
 \end{eqnarray}
Apparently  $ T^{a_4} $ obeys
\begin{eqnarray}
 T^{a_4}   =   C^{a_4} + C^{a_4} G_0 T^{a_4}\label{10}
\end{eqnarray}

In explicite notation (\ref{9}) and (\ref{10}) read
\begin{eqnarray}
\psi_{a_5 a_4}  =   G_0 \sum_{ b_5 \subset a_4} T^{a_4}_{a_5 b_5}
\psi^{(a_4)}_{b_5}
=   G_0 \sum_{ b_5 \subset a_4} T^{a_4}_{a_5 b_5}
 \sum_{ b_4 \supset b_5}
\overline{\delta_{ a_4 b_4}} \psi_{b_5 b_4}
\label{11}
\end{eqnarray}
\begin{eqnarray}
T^{a_4}_{a_5 b_5}  =   t_{a_5} \overline{\delta_{ a_5 b_5}} +\sum_{
  c_5 \subset a_4}
t_{a_5} \overline{\delta_{ a_5 c_5}} G_0 T^{a_4}_{c_5 b_5}\label{15}
\end{eqnarray}

Note, there are two types of T-matrices. For $ a_4$ of the type
$ 123,4,5,6$ $ T_{ a_5 b_5}^{a_4}$ is a $ 3 \times 3 $ matrix and for
$ a_4$ of the type $ 12,34;5,6$  $ T_{ a_5 b_5}^{a_4}$ is a $ 2 \times 2 $ matrix.

Next we further decompose the right hand side of (\ref{11})  according
to  3-body  fragments $a_3$:
\begin{eqnarray}
 \psi_{a_5 a_4}^{ a_3} \equiv  \sum_{ b_5 \subset a_4} G_0 T_{a_5
   b_5}^{a_4}
\sum_{ b_4 \supset b_5, b_4 \subset a_3} \overline{\delta_{ a_4 b_4}}
\psi_{b_5 b_4} \label{16}
 \end{eqnarray}
and again
\begin{eqnarray}
 \psi_{a_5 a_4}  =\sum_{a_3 \supset a_4} \psi_{a_5 a_4}^{ a_3}
 \label{14}\end{eqnarray}
is an obvious consequence.

Using again (\ref{14}) Eq. (\ref{16}) can be rewritten as
\begin{eqnarray}
 \psi_{a_5 a_4}^{ a_3}  = \sum_{ b_5 \subset a_4} G_0 T_{a_5
   b_5}^{a_4}
\sum_{ b_4 \supset b_5, b_4 \subset a_3} \overline{\delta_{ a_4 b_4}}
  \sum_{b_3 \supset b_4} \psi_{b_5 b_4}^{ b_3}
\end{eqnarray}

Analogous to (\ref{7}) one separates $ b_3 = a_3 $ from $ b_3 \ne a_3 $ and gets
\begin{eqnarray}
 & & \psi_{a_5 a_4}^{ a_3}  - G_0 \sum_{ b_5 \subset a_4} T_{a_5
    b_5}^{a_4}
\sum_{ b_4 \supset b_5, b_4 \subset a_3} \overline{\delta_{ a_4 b_4}}
 \psi_{b_5 b_4}^{ a_3}
  =  G_0 \sum_{ b_5 \subset a_4} T_{a_5 b_5}^{a_4} \sum_{ b_4
   \supset b_5, b_4
\subset a_3} \overline{\delta_{ a_4 b_4}}
 \psi_{b_5 b_4}^{ (a_3)}
 \label{19}\end{eqnarray}
where we  defined
   \begin{eqnarray}
   \psi_{b_5 b_4}^{ (a_3)} \equiv \sum_{ b_3 \supset b_4}
\overline{\delta_{ a_3 b_3}}\psi_{b_5 b_4}^{ b_3}
   \end{eqnarray}

Using a matrix notation it is easily seen that (\ref{19}) leads to
\begin{eqnarray}
\psi^{a_3}_{a_5 a_4}  =  G_0 \sum_{ b_4  \subset a_3}
\sum_{ b_5 \subset b_4}D^{a_3}_{a_5,a_4;b_5 b_4}\psi^{(a_3)}_{b_5 b_4}\label{125}
  \end{eqnarray}
where $ D^{a_3}_{a_5,a_4;b_5 b_4} $ obeys
\begin{eqnarray}
   D^{a_3}_{a_5,a_4;b_5 b_4}  =    T^{a_4}_{a_5 b_5}
   \overline{\delta_{ a_4 b_4}}
+ \sum_{ c_4  \subset a_3} \sum_{ c_5 \subset a_4}
    T^{a_4}_{a_5 c_5} \overline{\delta_{ a_4 c_4}} G_0 D^{a_3}_{c_5,c_4;b_5 b_4}
   \end{eqnarray}

Note for $ a_3 = 1234,5,6 $ there are 18 pairs of $ a_5,a_4$,
for $ a_3 = 123,45,6$ there are 9 pairs of $ a_5,a_4$ and for $ a_3 = 12,34,56$
there are 6 pairs of $ a_5,a_4$. This defines the dimensions of
the different D-matrices. In the following the single nucleons in
$a_n$ will no longer be displayed.

\section{ Implementation of the identity of the nucleons}
\label{sec_three}

We start from (\ref{11}) choosing the case $ a_5 =12$ and $ a_4 =12,34$
 and obtain
\begin{eqnarray}
\psi_{12; 12,34} = G_0  T_{12, 12}^{12,34} \psi^{(12,34)}_{12} +
G_0  T_{12, 34}^{12,34} \psi^{(12,34)}_{34}\label{26}
\end{eqnarray}
which according to (\ref{comp1}) is
\begin{eqnarray}
\psi_{12; 12,34} & = &  G_0  T_{12, 12}^{12,34}  ( \psi_{12,123} +
\psi_{12,124}
+ \psi_{12,125} + \psi_{12,126}\cr
& + &  \psi_{12; 12,35} + \psi_{12; 12,36} + \psi_{12; 12,45}+
\psi_{12; 12,46}
+ \psi_{12; 12,56})\cr
& + & G_0  T_{12, 34}^{12,34} ( \psi_{34,134} + \psi_{34,234} +
\psi_{34,345}
+ \psi_{34,346}\cr
& + &  \psi_{34; 34,15} + \psi_{34; 34,16} + \psi_{34; 34,25}+
\psi_{34; 34,26}
+ \psi_{34; 34,56})\label{27}
\end{eqnarray}

It is easily seen, going back to the definitions (\ref{fad}) and
(\ref{5}) together with the antisymmetry requirement for the total
state $ \Psi$ that
\begin{eqnarray}
  & &    \psi_{34,134} + \psi_{34,234} + \psi_{34,345} + \psi_{34,346}
 +   \psi_{34; 34,15} + \psi_{34; 34,16} + \psi_{34; 34,25}+
  \psi_{34; 34,26}
+ \psi_{34; 34,56}\cr
& =&   P_{13}P_{24} (  \psi_{12,123} + \psi_{12,124} + \psi_{12,125}
+ \psi_{12,126}\cr
& + &  \psi_{12; 12,35} + \psi_{12; 12,36} + \psi_{12; 12,45}+
\psi_{12; 12,46}
+ \psi_{12; 12,56})\label{28}
\end{eqnarray}

Therefore (\ref{27}) turns into
\begin{eqnarray}
& & \psi_{12; 12,34} = G_0 (  T_{12, 12}^{12,34} +  T_{12, 34}^{12,34}
P_{13}P_{24})  ( \psi_{12,123} + \psi_{12,124} + \psi_{12,125} + \psi_{12,126}\cr
& + &  \psi_{12; 12,35} + \psi_{12; 12,36} + \psi_{12; 12,45}+
\psi_{12; 12,46}
+ \psi_{12; 12,56})\label{29}
\end{eqnarray}

The coupled equations (\ref{15}) written out for $ a_4 = 12,34$
and $ a_5 $ or $ b_5 $ equal to $12$ or $34$  yield when acting with
$ \tilde P \equiv P_{13}P_{24}$ from both sides
\begin{eqnarray}
  \tilde P T_{12, 12}^{12,34} \tilde P & = &  T_{34,34}^{12,34}\cr
  \tilde P  T_{12, 34}^{12,34} \tilde P & = &   T_{34,12}^{12,34}
  \end{eqnarray}

Then defining
\begin{eqnarray}
T^{12,34} \equiv  T_{12, 12}^{12,34} + T_{12, 34}^{12,34} \tilde P\label{3.6}
\end{eqnarray}
it follows that $  T^{12,34} $ obeys
\begin{eqnarray}
T^{12,34}   =   t_{12}  \tilde P  + t_{12} G_0 \tilde P T^{12,34}\label{141}
\end{eqnarray}

Therefore (\ref{29}) simplifies to
\begin{eqnarray}
& & \psi_{12; 12,34} = G_0   T^{12,34}  ( \psi_{12,123} +
  \psi_{12,124}
+ \psi_{12,125} + \psi_{12,126}\cr
& + &  \psi_{12; 12,35} + \psi_{12; 12,36} + \psi_{12; 12,45}+
\psi_{12; 12,46}
+ \psi_{12; 12,56})\label{33}
\end{eqnarray}

Starting again from (\ref{11}) but now for $ a_5 =12$ and $ a_4 = 123$
one obtains
\begin{eqnarray}
& & \psi_{12,123} = G_0 T^{123}_{12,12}  \psi^{(123)}_{12} + G_0
  T^{123}_{12,23}
\psi^{(123)}_{23}
+ G_0 T^{123}_{12,31}  \psi^{(123)}_{31}\cr
& =&    G_0  T^{123}_{12,12}  ( \psi_{ 12,124} + \psi_{12,125} + \psi_{12,126}\cr
  &+&   \psi_{ 12; 12,34} + \psi_{ 12; 12,35} + \psi_{ 12; 12,36}
 +   \psi_{ 12; 12,45} + \psi_{ 12; 12,46} + \psi_{ 12; 12,56})\cr
& + & G_0 T^{123}_{12,23} ( \psi_{ 23,234} + \psi_{23,235} + \psi_{23,236}\cr
&  + &  \psi_{ 23; 23,14} + \psi_{ 23; 23,15} + \psi_{ 23; 23,16}
 +   \psi_{ 23; 23,45} + \psi_{ 23; 23,46} + \psi_{ 23; 23,56})\cr
& + & G_0 T^{123}_{12,31} ( \psi_{ 31,314} + \psi_{31,315} + \psi_{31,316}\cr
&  + &  \psi_{ 31; 31,24} + \psi_{ 31; 31,25} + \psi_{ 31; 31,26}
 +   \psi_{ 31; 31,45} + \psi_{ 31; 31,46} + \psi_{ 31; 31,56})\cr
& =&    G_0 ( T^{123}_{12,12} +  T^{123}_{12,23} P_{12}P_{23} +
T^{123}_{12,31}
P_{13}P_{23})
   ( \psi_{ 12,124} + \psi_{12,125} + \psi_{12,126}\cr
&  + &  \psi_{ 12; 12,34} + \psi_{ 12; 12,35} + \psi_{ 12; 12,36}
 +   \psi_{ 12; 12,45} + \psi_{ 12; 12,46} + \psi_{ 12; 12,56})\label{34}
\end{eqnarray}
where we used permutation properties similar as in (\ref{28}).

The corresponding coupled sets (\ref{15}) for $a_4=123$ and using
relations like $P_{13}P_{23}t_{23}P_{23}P_{13}=t_{12}$ reveals that
\begin{eqnarray}
T^{123} \equiv T^{123}_{12,12} +  T^{123}_{12,23} P_{12}P_{23}
+ T^{123}_{12,31} P_{13}P_{23}\label{145}
\end{eqnarray}
obeys the equation
\begin{eqnarray}
T^{123} = t_{12} P +  t_{12} P G_0 T^{123}\label{36}
\end{eqnarray}
where $ P \equiv P_{12}P_{23} + P_{13} P_{23} $.

Then (\ref{34}) simplifies to
\begin{eqnarray}
& & \psi_{12,123}  =    G_0  T^{123}  ( \psi_{ 12,124} + \psi_{12,125}
+ \psi_{12,126}\cr
&  + &  \psi_{ 12; 12,34} + \psi_{ 12; 12,35} + \psi_{ 12; 12,36}
 +   \psi_{ 12; 12,45} + \psi_{ 12; 12,46} + \psi_{ 12; 12,56})\label{3.12}
\end{eqnarray}

The next step is to  decompose $ \psi_{12,123}$  according to
(\ref{16}). For $a_5 =12,a_4 = 123$ the possible  $a_3$'s  are:
  $ 1234 - 1235 - 1236 - 123,45 - 123,46 - 123,56$.

Lets  begin with
\begin{eqnarray}
  & & \psi_{12,123}^{ 1234} \equiv G_0 T_{12, 12}^{123} (
  \psi_{12,124}
+ \psi_{12; 12,34})\cr
  & + &  G_0 T_{12, 23}^{123} ( \psi_{23,234} + \psi_{23; 23,14})
   +  G_0 T_{12, 31}^{123} ( \psi_{31,134} + \psi_{31; 31,24})
\label{3.13}
  \end{eqnarray}

Since
  \begin{eqnarray}
  \psi_{23,234} + \psi_{23; 23,14} = P_{12}P_{23} ( \psi_{12,124}
+ \psi_{12; 12,34})\cr
  \psi_{31,134} + \psi_{31; 31,24} = P_{13}P_{23} ( \psi_{12,124}
+ \psi_{12; 12,34})
  \end{eqnarray}
(\ref{3.13}) simplifies according to (\ref{36}) to
   \begin{eqnarray}
   \psi_{12,123}^{ 1234} = G_0 T^{123} ( \psi_{12,124} + \psi_{12; 12,34})\label{40}
  \end{eqnarray}

Similarily
  \begin{eqnarray}
   \psi_{12,123}^{ 1235} & = &  G_0 T^{123} ( \psi_{12,125} + \psi_{12; 12,35})\cr
   \psi_{12,123}^{ 1236} & = &  G_0 T^{123} ( \psi_{12,126} + \psi_{12; 12,36})
  \end{eqnarray}

Again  using symmetry properties one gets
\begin{eqnarray}
  & & \psi_{12,123}^{ 123,45} = G_0 T^{123}  \psi_{12;12,45}\label{42}
  \end{eqnarray}
  \begin{eqnarray}
  & & \psi_{12,123}^{ 123,46} = G_0 T^{123}  \psi_{12;12,46}
\end{eqnarray}
  \begin{eqnarray}
  & & \psi_{12,123}^{ 123,56} = G_0 T^{123}  \psi_{12;12,56}
  \end{eqnarray}
All summed up
\begin{eqnarray}
& & \psi_{12,123} = \psi_{12,123}^{ 1234} + \psi_{12,123}^{ 1235}
+ \psi_{12,123}^{ 1236} + \psi_{12,123}^{ 123,45}
 +   \psi_{12,123}^{ 123,46} + \psi_{12,123}^{ 123,56}\label{44}
\end{eqnarray}
agrees with (\ref{3.12}).

Next we decompose (\ref{11}) according to (\ref{16}) for $a_5 =12,a_4=12,34$.
The possible $a_3$'s are:
    $1234 - 125,34 - 126,34 - 12,345 - 12,346 - 12,34,56$,
which are now regarded in turn.
\begin{eqnarray}
    & & \psi_{12;12,34}^{1234} = G_0  T_{12, 12}^{12,34} (
  \psi_{12,123}
+ \psi_{12,124})
      +   G_0  T_{12, 34}^{12,34} ( \psi_{34,234} + \psi_{34,134})
    \end{eqnarray}

Since
    \begin{eqnarray}
     P_{13}P_{24} ( \psi_{12,123} + \psi_{12,124}) = \psi_{34,234} + \psi_{34,134}
    \end{eqnarray}
one can use (\ref{3.6}) and gets
  \begin{eqnarray}
    \psi_{12;12,34}^{1234} & = &  G_0 ( T_{12, 12}^{12,34} +  T_{12,34}^{12,34}
P_{13}P_{24})  ( \psi_{12,123} + \psi_{12,124})\cr
    & = & G_0  T^{12,34}  ( \psi_{12,123} + \psi_{12,124})\label{47}
    \end{eqnarray}

Next
    \begin{eqnarray}
    \psi_{12;12,34}^{125,34} & = &  G_0  T_{12,
      12}^{12,34}\psi_{12,125}
+ G_0   T_{12, 34}^{12,34} ( \psi_{34; 15,34} + \psi_{34;25,34})\cr
    \psi_{12;12,34}^{ 345,12} & = &  G_0  T_{12, 12}^{12,34}(\psi_{12;
      12,35}
+ \psi_{12;12,45})+ G_0   T_{12, 34}^{12,34}  \psi_{34,345}
    \end{eqnarray}

The two amplitudes $\psi_{12;12,34}^{125,34}$ and $
\psi_{12;12,34}^{345,12}$ can not be related by permutations, but
their sum can be used
\begin{eqnarray}
     \psi_{12;12,34}^{125,34} +  \psi_{12;12,34}^{ 345,12} &=&
G_0  T_{12, 12}^{12,34}( \psi_{12,125} + \psi_{12; 12,35} + \psi_{12;12,45})\cr
    & + & G_0   T_{12, 34}^{12,34} ( \psi_{34; 15,34} +
\psi_{34;25,34}+ \psi_{34,345})
    \end{eqnarray}
in the sense
\begin{eqnarray}
   P_{13}P_{24}     ( \psi_{12,125} + \psi_{12; 12,35} +
   \psi_{12;12,45})
= \psi_{34; 15,34} + \psi_{34;25,34}+ \psi_{34,345}
\end{eqnarray}

This leads to
\begin{eqnarray}
   \psi_{12;12,34}^{125,34} +  \psi_{12;12,34}^{ 345,12} = G_0
   T^{12,34}( \psi_{12,125}  + \psi_{12; 12,35} + \psi_{12;12,45})\label{51}
    \end{eqnarray}

Similarily
\begin{eqnarray}
   \psi_{12;12,34}^{126,34} + \psi_{12;12,34}^{346,12}
=  G_0  T^{12,34}( \psi_{12,126}  + \psi_{12; 12,36} + \psi_{12;12,46})\label{52}
    \end{eqnarray}
and finally
\begin{eqnarray}
    \psi_{12;12,34}^{12,34,56} = G_0  T_{12, 12}^{12,34} \psi_{12;
      12,56}
+ G_0  T_{12, 34}^{12,34} \psi_{34; 34,56} = G_0  T^{12,34}
\psi_{12; 12,56}\label{53}
    \end{eqnarray}

Thus, Eqs.(\ref{47}), (\ref{51})-(\ref{53}),    summarizes   to
 \begin{eqnarray}
    & & \psi_{12;12,34} = \psi_{12;12,34}^{1234} +
   \psi_{12;12,34}^{125,34}
+ \psi_{12;12,34}^{345,12}
     +    \psi_{12;12,34}^{126,34} + \psi_{12;12,34}^{346,12}+
\psi_{12;12,34}^{12,34,56}\label{54}
\end{eqnarray}
which when written out agrees with (\ref{33}).

The two amplitudes $ \psi_{12,123}^{1234}$ and $\psi_{12;12,34}^{1234}$
expressed in (\ref{40}) and (\ref{47}) are connected to each other as shown now.
 The expression (\ref{44}) can easily be converted to $ \psi_{12,124}$
and using in addition (\ref{54}) one finds
\begin{eqnarray}
& & \psi_{12,123}^{ 1234} - G_0 T^{123} ( \psi_{12,124}^{ 1234}
+ \psi_{12;12,34}^{ 1234})\cr
 & = &  G_0 T^{123} ( \psi_{12,124}^{ 1245} + \psi_{12,124}^{ 1246}
+ \psi_{12,124}^{ 124,35}
 +  \psi_{12,124}^{ 124,36 } + \psi_{12,124}^{ 124,56 }
+  \psi_{12;12,34}^{125,34 }\cr
 & + &  \psi_{12;12,34}^{345,12 } +  \psi_{12;12,34}^{126,34 }
+ \psi_{12;12,34}^{346,12 }
   +   \psi_{12;12,34}^{ 12,34, 56 })\label{55}
\end{eqnarray}

Correspondingly  (\ref{47})  yields
\begin{eqnarray}
& & \psi_{12;12,34}^{ 1234} - G_0  T^{12,34} (  \psi_{12,123}^{ 1234}
+\psi_{12,124}^{ 1234}) \cr
 & = &  G_0  T^{12,34}  ( \psi_{12,123}^{ 1235} +  \psi_{12,123}^{
  123,45 }
+ \psi_{12,123}^{ 1236}
  +   \psi_{12,123}^{123,46} +   \psi_{12,123}^{123,56}
+   \psi_{12,124}^{ 1245}\cr
 &  + &   \psi_{12,124}^{ 124,35 } + \psi_{12,124}^{ 1246}
+   \psi_{12,124}^{124,36} + \psi_{12,124}^{124,56})
\label{56}
\end{eqnarray}

With $ \psi_{12,124}^{ 1234}  = - P_{34} \psi_{12,123}^{ 1234}$  we
can put (\ref{55}) and (\ref{56}) into a matrix form:
\begin{eqnarray}
    & & \left(
   \begin{array}{rcl}
   \psi_{12,123}^{ 1234}\\
   \psi_{12;12,34}^{ 1234}\\
   \end{array}
   \right) - G_0 \left(
   \begin{array}{rcl}
   -T^{123} P_{34}  & T^{123} \\
    T^{12,34} ( 1-P_{34}) & 0 \\
   \end{array}
   \right)
   \left(
   \begin{array}{rcl}
   \psi_{12,123}^{ 1234}\\
   \psi_{12;12,34}^{ 1234}\\
   \end{array}
   \right)
   \cr
   & = &    G_0  \left(
   \begin{array}{rcl}
    T^{123} &(&  \psi_{12,124}^{ 1245} + \psi_{12,124}^{ 1246}
+ \psi_{12,124}^{ 124,35}
 +  \psi_{12,124}^{ 124,36 } + \psi_{12,124}^{ 124,56 }
+   \psi_{12;12,34}^{125,34 }\cr
 & + &  \psi_{12;12,34}^{345,12 } +  \psi_{12;12,34}^{126,34 }
+ \psi_{12;12,34}^{346,12 }
  +    \psi_{12;12,34} ^{ 12,34, 56 })\\
    T^{12,34} &(&  \psi_{12,123}^{ 1235} +  \psi_{12,123}^{ 123,45 }
+ \psi_{12,123}^{ 1236} + \psi_{12,123}^{123,46}
  +   \psi_{12,123}^{123,56}  +   \psi_{12,124}^{ 1245} \cr
&+&  \psi_{12,124}^{ 124,35 } + \psi_{12,124}^{ 1246}
  +   \psi_{12,124}^{124,36} + \psi_{12,124}^{124,56})\\
   \end{array}
   \right)
\label{57}
\end{eqnarray}

Since
\begin{eqnarray}
   \psi_{12,124} ^{1245} & = &  - P_{34} \psi_{12,123}^{1235}\cr
   \psi_{12,124} ^{124, 35} & = &  - P_{34} \psi_{12,123}^{123, 45}
   \end{eqnarray}
the right hand side of (\ref{57}) can be  factored and (\ref{57})
achieves the form
\begin{eqnarray}
 & &    \left(
   \begin{array}{rcl}
   \psi_{12,123}^{ 1234}\\
   \psi_{12;12,34}^{ 1234}\\
   \end{array}
   \right) - G_0 \left(
   \begin{array}{rcl}
   -T^{123} P_{34}  & T^{123} \\
    T^{12,34} ( 1-P_{34}) & 0 \\
   \end{array}
   \right)
   \left(
   \begin{array}{rcl}
   \psi_{12,123}^{ 1234}\\
   \psi_{12;12,34}^{ 1234}\\
   \end{array}
   \right)\cr
  & = &  G_0 \left(
   \begin{array}{rcl}
   T^{123}(-P_{34})  & T^{123} \\
    T^{12,34} ( 1 - P_{34})  & 0\\
   \end{array}
   \right)
   \left(
   \begin{array}{rcl}
    & & \psi_{12,123}^{ 1235} + \psi_{12,123}^{ 123,45 }+ \psi_{12,123}^{ 1236}\cr
      &+&   \psi_{12,123}^{ 123,46 } + \psi_{12,123}^{ 123,56 }  \\
    & &  \psi_{12;12,34}^{125,34 } + \psi_{12;12,34}^{345,12 }
+   \psi_{12;12,34}^{126,34 }\cr
     &+&   \psi_{12;12,34}^{346,12 } +  \psi_{12;12,34}^{12,34,56 }\\
   \end{array}
   \right)
\end{eqnarray}

The right hand side can be reduced applying permutations and one obtains
\begin{eqnarray}
 & &    \left(
   \begin{array}{rcl}
   \psi_{12,123}^{ 1234}\\
   \psi_{12;12,34}^{ 1234}\\
   \end{array}
   \right) - G_0 \left(
   \begin{array}{rcl}
   -T^{123} P_{34}  & T^{123} \\
    T^{12,34} ( 1-P_{34}) & 0 \\
   \end{array}
   \right)
   \left(
   \begin{array}{rcl}
   \psi_{12,123}^{ 1234}\\
   \psi_{12;12,34}^{ 1234}\\
   \end{array}
   \right)\cr
  & = &  G_0 \left(
   \begin{array}{rcl}
   T^{123}(-P_{34})  & T^{123} \\
    T^{12,34} ( 1 - P_{34})  & 0\\
   \end{array}
   \right)
         \left(
   \begin{array}{rcl}
    - ( P_{45} + P_{46}) \psi_{ 12,123}^{1234}
+ ( 1 - P_{56} - P_{46}) \psi_{12,123}^{ 123,45 }\\
     (1 - P_{56}) (\psi_{12;12,34}^{125,34 } + \psi_{12;12,34}^{345,12
}) +  \psi_{12;12,34}^{12,34,56 }\\
   \end{array}
   \right)
   \label{60}
   \end{eqnarray}
which can be put into the form     $\left(
   \begin{array}{rcl}
   \psi_{12,123}^{ 1234}\\
   \psi_{12;12,34}^{ 1234}\\
   \end{array}
   \right) $ with the result
\begin{eqnarray}
   \left(
   \begin{array}{rcl}
   \psi_{12,123}^{ 1234}\\
   \psi_{12;12,34}^{ 1234}\\
   \end{array}
   \right)
   \equiv G_0 \left(
   \begin{array}{rcl}
   D_{11} & D_{12} \\
   D_{21} & D_{22}\\
   \end{array}
   \right)
  \left(
   \begin{array}{rcl}
    - ( P_{45} + P_{46}) \psi_{ 12,123}^{1234} + ( 1 - P_{56}
- P_{46}) \psi_{12,123}^{ 123,45 }\\
     (1 - P_{56}) (\psi_{12;12,34}^{125,34 } +
 \psi_{12;12,34}^{345,12 }) +
\psi_{12;12,34}^{12,34,56 }\\
   \end{array}
   \right)
   \label{61}
\end{eqnarray}
where  the matrix $ D$ obeys
\begin{eqnarray}
  & &    \left(
   \begin{array}{rcl}
   D_{11} & D_{12} \\
   D_{21} & D_{22}\\
   \end{array}
   \right) =   \left(
   \begin{array}{rcl}
   -T^{123} P_{34}  & T^{123} \\
    T^{12,34} ( 1-P_{34}) & 0 \\
   \end{array}
   \right)\cr
   &+&  \left(
   \begin{array}{rcl}
   -T^{123} P_{34}  & T^{123} \\
    T^{12,34} ( 1-P_{34}) & 0 \\
   \end{array}
   \right)
   G_0
      \left(
   \begin{array}{rcl}
   D_{11} & D_{12} \\
   D_{21} & D_{22}\\
   \end{array}
   \right)
   \label{62}
\end{eqnarray}

For the numerical treatment, however,  we consider the structure (\ref{60})
to be more advantageous.

The right hand side  of (\ref{60}) contains new amplitudes.
After adequate permutation of (\ref{54}) one obtains from (\ref{42})
\begin{eqnarray}
\psi_{12,123}^{ 123,45}  & = & G_0 T^{123}  (\psi_{12;12,45}^{ 1245}
+ \psi_{12;12,45}^{ 123,45} + \psi_{12;12,45}^{ 345,12}
  +   \psi_{12;12,45}^{ 456,12} + \psi_{12;12,45}^{ 126,45}
+ \psi_{12;12,45}^{ 12,45,36})
\end{eqnarray}
or
\begin{eqnarray}
& & \psi_{12,123}^{ 123,45} - G_0 T^{123} (-( 1- P_{36}) P_{53} (
\psi_{12;12,34}^{ 125,34} + \psi_{12;12,34}^{ 345,12}))\cr
 &=&  G_0 T^{123} (  - P_{35} \psi_{12;12,34}^{ 1234}
+ \psi_{12;12,45}^{ 12,45,36})
\end{eqnarray}

Further  (\ref{51}) yields inserting the decomposition of the right
hand side related to (\ref{44}) and (\ref{54})  yields
\begin{eqnarray}
& & \psi_{12;12,34}^{ 125,34} + \psi_{12,12,34}^{ 345,12}
= G_0  T^{12,34}   (\psi_{12,125}^{1235} + \psi_{12,125}^{1245}
+ \psi_{12,125}^{125,34}\cr
 & + &  \psi_{12,125}^{1256} + \psi_{12,125}^{125,36}+  \psi_{12,125}^{125,46}
  +  \psi_{12;12,35}^{ 1235} + \psi_{12;12,35}^{ 124,35}
+ \psi_{12;12,35}^{ 345,12}\cr
&  +  & \psi_{12;12,35}^{ 356,12} + \psi_{12;12,35}^{ 126,35}
+ \psi_{12;12,35}^{ 12,56,34}
 +  \psi_{12;12,45}^{ 1245} + \psi_{12;12,45}^{ 123,45}
+ \psi_{12;12,45}^{ 345,12}\cr
&  +  & \psi_{12;12,45}^{ 456,12} + \psi_{12;12,45}^{ 126,45}
+ \psi_{12;12,45}^{ 12,34,56}
\end{eqnarray}

Here quite a few amplitudes can be related to previous ones by permutations
leading to:
\begin{eqnarray}
& & \psi_{12;12,34}^{ 125,34} + \psi_{12,12,34}^{ 345,12}  -  G_0  T^{12,34}
 (  P_{35} P_{56} - P_{35} - P_{46} - P_{45}) ( \psi_{12;12,34}^{
    125,34} + \psi_{12;12,34}^{ 345,12})\cr
 &  = &  G_0  T^{12,34} ( ( 1- P_{36} - P_{34} - P_{46})P_{34} P_{35}
\psi_{12,123}^{1234} + ( 1 - P_{34}) P_{46}P_{35}\psi_{12,123}^{123,45}\cr
&  -& ( P_{45} + P_{35}) \psi_{12;12,34}^{ 1234} - P_{35}
\psi_{12;12,34}^{ 12,34,56})
\end{eqnarray}

Finally we regard $ \psi_{12;12,34}^{ 12,34,56}$, which after
(\ref{53}) and  suitable permutations of (\ref{54}) turns into
\begin{eqnarray}
& & \psi_{12;12,34}^{ 12,34,56}  - G_0T^{12,34}P_{35} P_{46}
\psi_{12;12,34}^{ 12,34, 56}\cr
&   = & G_0   T^{12,34} ( P_{35} P_{46} \psi_{12;12,34}^{ 1234}
+ ( 1-P_{34})P_{35}P_{46}( \psi_{12;12,34}^{ 125,34} +  \psi_{12;12,34}^{ 345,12}))
\end{eqnarray}
where we used
\begin{eqnarray}
\psi_{12;12,56}^{ 12,56,34} & =  & G_0 T^{12,56} \psi_{12;12,34}
 =    P_{35} P_{46} G_0 T^{12,34}  \psi_{12;12,56} =
P_{35} P_{46} \psi_{12;12,34}^{ 12,34,56}
\end{eqnarray}

Thus we end up with 5 independent amplitudes:
$ \psi_{12,123}^{ 1234}$, $\psi_{12;12,34}^{ 1234}$, $\psi_{12,123}^{
  123,45}$, $\psi_{12;12,34}^{ 125,34} + \psi_{12,12,34}^{ 345,12}$,
and $\psi_{12;12,34}^{ 12,34,56}$, coupled in the equations
\begin{eqnarray}
 & &    \left(
   \begin{array}{rcl}
   \psi_{12,123}^{ 1234}\\
   \psi_{12;12,34}^{ 1234}\\
   \end{array}
   \right) - G_0 \left(
   \begin{array}{rcl}
   -T^{123} P_{34}  & T^{123} \\
    T^{12,34} ( 1-P_{34}) & 0 \\
   \end{array}
   \right)
   \left(
   \begin{array}{rcl}
   \psi_{12,123}^{ 1234}\\
   \psi_{12;12,34}^{ 1234}\\
   \end{array}
   \right)\cr
  & = &  G_0 \left(
   \begin{array}{rcl}
   T^{123}(-P_{34})  & T^{123} \\
    T^{12,34} ( 1 - P_{34})  & 0\\
   \end{array}
   \right)
         \left(
   \begin{array}{rcl}
    - ( P_{45} + P_{46}) \psi_{ 12,123}^{1234} + ( 1 - P_{56} -
    P_{46}) \psi_{12,123}^{ 123,45 }\\
     (1 - P_{56}) (\psi_{12;12,34}^{125,34 } + \psi_{12;12,34}^{345,12
    }) +  \psi_{12;12,34}^{12,34,56 }\\
   \end{array}
   \right)
   \label{70}
   \end{eqnarray}

\begin{eqnarray}
& & \psi_{12,123}^{ 123,45} - G_0 T^{123} (-( 1- P_{36}) P_{53} (
\psi_{12;12,34}^{ 125,34} + \psi_{12;12,34}^{ 345,12}))\cr
& = &  - G_0 T^{123} P_{35}( \psi_{12;12,34}^{ 1234}
+ \psi_{12;12,34}^{ 12,34,56})\label{71}
\end{eqnarray}
\begin{eqnarray}
& & \psi_{12;12,34}^{ 125,34} + \psi_{12,12,34}^{ 345,12}  -  G_0  T^{12,34}
 (  P_{35} P_{56} - P_{35} - P_{46} - P_{45}) ( \psi_{12;12,34}^{
    125,34} + \psi_{12;12,34}^{ 345,12})\cr
 &  = &  G_0  T^{12,34} ( ( 1- P_{36} - P_{34} - P_{46})P_{34} P_{35}
\psi_{12,123}^{1234} + ( 1 - P_{34}) P_{46}P_{35}\psi_{12,123}^{123,45}\cr
&  -& ( P_{45} + P_{35}) \psi_{12;12,34}^{ 1234} -P_{35}
\psi_{12;12,34}^{ 12,34,56})\label{72}
\end{eqnarray}

\begin{eqnarray}
& & \psi_{12;12,34}^{ 12,34,56}  - G_0T^{12,34}P_{35} P_{46}
\psi_{12;12,34}^{ 12,34, 56}
 =  G_0   T^{12,34} ( P_{35} P_{46} \psi_{12;12,34}^{ 1234}\cr
&  + &  ( 1-P_{34})P_{35}P_{46}( \psi_{12;12,34}^{ 125,34} +
\psi_{12;12,34}^{ 345,12}))\label{73}
\end{eqnarray}

It remains to establish the expression for the total state $ \Psi $.
We use (\ref{tot}), (\ref{5}), (\ref{44}), (\ref{54})
 together with permutations and obtain
\begin{eqnarray}
\Psi & = &  [  1 - P_{23} - P_{24} - P_{13} - P_{14}  + P_{13} P_{24}]
[ ( 1 - P_{34}) \psi_{12,123}^{1234} +  \psi_{12;12,34}^{1234}]\cr
& - &  [  1 - P_{23} - P_{24} - P_{13} - P_{14}  + P_{13} P_{24}]\cr
&& [ ( ( 1 - P_{34}) ( P_{45} + P_{46} + P_{35} +  P_{36})
- ( P_{35}P_{46} + P_{36} P_{45})) \psi_{12,123}^{1234}\cr
& + & ( P_{45} + P_{46} + P_{35} +  P_{36} -  P_{35} P_{46})
\psi_{12;12,34}^{1234}]\cr
& - & [  P_{25} + P_{26} + P_{15} + P_{16} -    P_{13} P_{25}
-  P_{13} P_{26} -  P_{14} P_{25} -  P_{14} P_{26} -  P_{15} P_{26}]\cr
& & [ ( 1 - P_{34}) \psi_{12,123}^{1234} +  \psi_{12;12,34}^{1234}\cr
& - & ( (  1 -P_{34}) ( P_{45} + P_{46} + P_{35} +  P_{36})- ( P_{35}
      P_{46}
+ P_{36} P_{45})) \psi_{12,123}^{1234}\cr
& - & ( P_{45} + P_{46} + P_{35} +  P_{36} -  P_{35} P_{46})
\psi_{12;12,34}^{1234}]\cr
& + & [  1 - P_{23} - P_{24} - P_{13} - P_{14}  + P_{13} P_{24}]
 [( 1 - P_{34}  -   P_{35} - P_{36}) (   1 -  P_{56} - P_{46})
\psi_{12,123}^{123,45}\cr
& + & (1 -  P_{45} - P_{46} - P_{35} -  P_{36}  + P_{35} P_{46})
 ( ( 1 - P_{56}) ( \psi_{12;12,34}^{125,34} +
\psi_{12;12,34}^{345,12})
+ \psi_{12;12,34}^{12,34,56})]\cr
&   - & [  P_{25} + P_{26} + P_{15} + P_{16}   -  P_{13} P_{25}
-  P_{13} P_{26} -  P_{14} P_{25} -  P_{14} P_{26} -  P_{15} P_{26}]\cr
& & [  ( 1 - P_{34}  -   P_{35} - P_{36}) (   1 -  P_{56} - P_{46})
\psi_{12,123}^{123,45}
+ (1 -  P_{45} - P_{46} - P_{35} -  P_{36}  + P_{35} P_{46})\cr
& & ( ( 1 - P_{56}) ( \psi_{12;12,34}^{125,34} +
\psi_{12;12,34}^{345,12})
+ \psi_{12;12,34}^{12,34,56})]
\label{TOT}
\end{eqnarray}

The first piece
\begin{eqnarray}
\Psi_4 & \equiv  &  [  1 - P_{23} - P_{24} - P_{13} - P_{14}  + P_{13} P_{24}]
[ ( 1 - P_{34}) \psi_{12,123}^{1234} +  \psi_{12;12,34}^{1234}]
\end{eqnarray}
has exactly the form for a 4-nucleon bound state
\cite{kamada1,kamada3}
based on the two Yakubovsky components
 $ \psi_{12,123}^{1234}$ and $  \psi_{12;12,34}^{1234} $. Now for
6 nucleons these two amplitudes depend not only
 on the momenta of nucleons $1$ to $4$ but  also on the  momenta
of nucleons 5 and 6. In the spirit of the effective 3-body
  model $ \alpha-n-n $ $\Psi_4$ factorizes into a product of the
$\alpha$-state and the wave function for the two neutrons.
   The next two pieces, still related to the same two Yakubovsky
   components,
antisymmetrize the nucleons 5 and 6
    in relation to the nucleons 1 to 4. The remaining two pieces, going
with the additional three Yakubovsky components
 $\psi_{12,123}^{123,45}$, $ \psi_{12;12,34}^{125,34} +
\psi_{12;12,34}^{345,12}$
and $ \psi_{12;12,34}^{12,34,56} $,
   go clearly beyond the effective 3-body model
and allow for additional sub clusterings of different  types.
  In addition one has to keep in mind that the most general form (\ref{TOT})
  allows many more distributions of neutrons and protons
  in the sub clusters than restricted in the effective $\alpha-n-n$ picture.

\section{ Technicalities for a numerical implementation}
\label{sec_four}

The 5 independent amplitudes require 5 different Jacobi momenta.
For $ \psi_{12,34}^{ 123} $
we choose
\begin{eqnarray}
    \vec a_1 & = &  \frac{1}{2} ( \vec k_1 - \vec k_2)\cr
 \vec a_2  & = &  \frac{1}{3} ( 2 \vec k_3 - \vec k_1 - \vec k_2)\cr
   \vec a_3  & = &  \frac{1}{4} ( 3 \vec k_4 - \vec k_1
- \vec k_2 -\vec k_3 )\cr
  \vec a_4  & = &  \frac{1}{2} ( \vec k_5 - \vec k_6)\cr
    \vec a_5  & = &  \frac{1}{3} ( 2 ( \vec k_5 + \vec k_6)
- \vec k_1 - \vec k_2 -\vec k_3 - \vec k_4 )\label{74}
\end{eqnarray}

Then  the individual momenta (under the condition $ \sum_i \vec k_i
=0$) in terms of those Jacobi momenta are
\begin{eqnarray}
  \vec k_1 & = &  \vec a_1 -   \frac{1}{2} \vec a_2 - \frac{1}{3}\vec
  a_3
- \frac{1}{4}\vec a_5\cr
  \vec k_2  & = &  - \vec a_1 -   \frac{1}{2} \vec a_2 -
  \frac{1}{3}\vec a_3
- \frac{1}{4}\vec a_5\cr
  \vec k_3  & = &  \vec a_2 -  \frac{1}{3}\vec a_3 - \frac{1}{4}\vec a_5\cr
  \vec k_4  & = &  \vec a_3 - \frac{1}{4}\vec a_5\cr
  \vec k_5  & = &  \vec a_4 + \frac{1}{2}\vec a_5
  \end{eqnarray}

The kinetic energy is
\begin{eqnarray}
   \sum_{i=1}^{6}\frac{k_i^2}{2m} = \frac{1}{m} ( a_1^2 + \frac{3}{4}
   a_2^2
+ \frac{2}{3} a_3^2 + a_4^2 + \frac{3}{8} a_5^2)
  \end{eqnarray}

The remaining choices of Jacobi momenta are given in Appendix \ref{appendix_a}.

In a partial wave decomposition the basis states suitable
for $ \psi_{12,34}^{ 123} $ are
\begin{eqnarray}
  & &  |a_1 a_2 a_3 a_4 a_5 \alpha_1>\cr
   &\equiv& | a_1 a_2 a_3 a_4 a_5, ( l_1 s_{12}) j_1
(l_2 \frac{1}{2})j_2 ( j_1 j_2) I_{3} (l_3\frac{1}{2}) j_4 ( I_3 j_4)
I_4  (l_5 s_{56})j_5 (l_4 j_5)I_{5}
    ( I_{4} I_{5}) J M> \cr
   & & | (t_{12}\frac{1}{2})t_3 ( t_3 \frac{1}{2}) t_4 ( t_4 t_{56}) T
M_T>
\label{77}
  \end{eqnarray}
Here the orbital angular momenta $l_i$ go with the $ \vec a_i$,
$s_{ij} $ are
two-body spins for particles $ij$, $j_i$ are total 1- and  2-body  angular
   momenta coupled out of orbital and spin angular momenta, $I_{3} $
and $ I_4$ are  total 3- and 4-body  angular momenta, $I_5$ the
     total angular momentum of particles 5 and 6 relative
to particles 1-4 and finally $ I_4$ and $I_5$ are coupled to $J$, the
conserved total 6-body angular momentum.  The second state in (\ref{77})
   refers to isospin in an obvious manner.

Using (\ref{77}) the amplitude $ \psi^{1234}_{123}$ has the representation
\begin{eqnarray}
  |\psi^{1234}_{123}> = \sum_{\alpha_1}\ \int \Pi_{i=1}^5 da_i a_i^2
| a_1 a_2 a_3 a_4 a_5 \alpha_1> < a_1 a_2 a_3 a_4 a_5 \alpha_1| \psi^{1234}_{123}>
\label{78}
\end{eqnarray}
where the discrete set of quantum  numbers $ \alpha_1$ runs   over
all values for a given $J$ and $T$.

Thereby antisymmetry requires that $ l_1 + s_{12} + t_{12}$
and $l_4 + s_{56} + t_{56}$ have to be odd.

Analogous basis states for the other 4 amplitudes and related Jacobi
momenta should be obvious. One other example is given below.

If one projects the coupled equations (\ref{70}) - (\ref{73}) from the
left onto the adequate basis states and expands the 5 amplitudes on the right
       hand side like in (\ref{78}) one faces the task to evaluate
the various kernels. As an example out of (\ref{70}) we take
\begin{eqnarray}
   & & < a_1 a_2 a_3 a_4 a_5 \alpha_1|\psi^{1234}_{12,123}>\cr
    &=& - G_0  < a_1 a_2 a_3 a_4 a_5 \alpha_1|  T^{123} P_{34}
    \sum_{\alpha_1'} \int \Pi_{i=1}^5 da_i' a_i^{'2} | a_1' a_2' a_3'
    a_4' a_5'
\alpha_1'> \cr
    & & < a_1' a_2' a_3' a_4' a_5' \alpha_1'| \psi^{1234}_{12,123}>\cr
    &+&  G_0  < a_1 a_2 a_3 a_4 a_5 \alpha_1|  T^{123}|
\sum_{\alpha_2'}
\int \Pi_{i=1}^5 db_i' b_i^{'2} | b_1' b_2' b_3' b_4' b_5' \alpha_2'>\cr
   & &  < b_1' b_2' b_3' b_4' b_5' \alpha_2'| \psi^{1234}_{12;12,34}> +
    \cdots
\label{79}
\end{eqnarray}
where the $ \cdots$ refer to the remaining operators and amplitudes
and the b-states are related to $ \psi_{12;12,34}^{1234}$.

Using techniques like the ones presented in
\cite{mybook,kamada1,nogga1,witala} it is straightforward, though tedious,
 to generate the kernels like
$ < a_1 a_2 a_3 a_4 a_5 \alpha_1|  T^{123} P_{34}
 | a_1' a_2' a_3' a_4' a_5' \alpha_1'> $
       or $ < b_1  b_2 b_3 b_4 b_5 \alpha_2| T^{12,34} ( 1-P_{34}) (
       P_{45}
+ P_{46}) | a_1' a_2' a_3' a_4' a_5' \alpha_1'>$.

In the course of the required recoupling among the different Jacobi
momenta the variables for the 5 amplitudes on the right hand sides are
in general linear combinations of intermediate integration variables
which include angles besides momentum magnitudes. An example
illustrates that situation:
\begin{eqnarray}
  H \equiv P_{45}  \psi_{12,123}^{1234}
  \end{eqnarray}

We project onto the basis states for Jacobi momenta of type b given in
Appendix \ref{appendix_a}:
\begin{eqnarray}
   | b > & \equiv &  | b_1 b_2 b_3 b_4 b_5; ( l_1 s_{12}) j_1 ( l_2
   s_{34}) j_2
( j_1 j_2 ) S ( L S ) I ( l_4 s_{56}) j_5 ( l_5 j_5) I_5 ( I I_5 ) JM>\cr
   & & | (t_{12} t_{34}) t_{1-4} ( t_{1-4} t_{56} ) T M_T>\label{81}
\end{eqnarray}
where the orbital angular momenta $l_i, i\ne 3$ go with $\vec b_i$
and
$L$ goes with $\vec b_3$. The 2-body spins are $s_{12}, s_{34}$ and
$s_{56}$,
$I$ is the total angular momentum for particles 1-4 and $I_5$ the
total
angular momentum of the pair 56 against the subsystem 1-4. They are
coupled
to the total angular momentum $J$, which is conserved. The isospin
coupling
should be obvious.

For the sake of simplicity we choose  s-waves which simplifies (\ref{81}) to
\begin{eqnarray}
    | b > & \equiv &   \delta_{S I} \delta_{ s_{12} j_1}\delta_{
      s_{34} j_2}
\delta_{ j_5 I_5} \delta_{ s_{56} j_5}\cr
  & & | b_1 b_2 b_3 b_4 b_5>
   | ( s_{12} s_{34} ) S ( S s_{56}) JM> |(t_{12} t_{34}) t_{1-4}
( t_{1-4} t_{56} ) T M_T>
  \end{eqnarray}

The basis states (\ref{77}) related to the Jacobi momenta of type a
and restricted to s-waves are
\begin{eqnarray}
| a > & = &  \delta_{ s_{12} j_1} \delta_{ j_2 \frac{1}{2}}
\delta_{j_4 \frac{1}{2}} \delta_{ j_5 s_{56}} \delta_{ j_5 I_5}
  | a_1 a_2 a_3 a_4 a_5> \cr
&& | ( s_{12} \frac{1}{2}) I_3
( I_3 \frac{1}{2} ) I_4 ( I_4 s_{56}) JM>
 | (t_{12} \frac{1}{2}) t_3 ( t_3 \frac{1}{2} )t_4 ( t_4 t_{56})  T
M_T>
\label{88}
\end{eqnarray}

Then $ \psi_{12,123}^{1234} $ has the representation
\begin{eqnarray}
   |\psi_{12,123}^{1234}> \equiv \sum_{\alpha_1} \delta \cdots
\int \Pi_{i=1}^5 da_i a_i^2 | a_1 a_2 a_3 a_4 a_5 > | spin>_a | isospin>_a
    < a|\psi_{12,123}^{1234}>\label{89}
  \end{eqnarray}
and H projected from the left is
\begin{eqnarray}
   < b| H> & = &  \delta_{S I} \delta_{ s_{12} j_1}\delta_{ s_{34}
     j_2}
\delta_{ j_5 I_5} \delta_{ s_{56} j_5}\cr
   & & < (t_{12} t_{34}) t_{1-4} ( t_{1-4} t_{56} ) T M_T|  <( s_{12}
s_{34} )
S ( S s_{56}) J M | <  b_1 b_2 b_3 b_4 b_5| P_{45}\cr
   & & \sum_{\alpha_1'} \delta \cdots \int \Pi_{i=1}^5 da_i a_i^2
| a_1 a_2 a_3 a_4 a_5 >
   | ( s_{12}' \frac{1}{2}) I_3' ( I_3' \frac{1}{2} ) I_4' ( I_4'
   s_{56}') JM>
\cr
& & | (t_{12}' \frac{1}{2}) t_3' ( t_3' \frac{1}{2} )t_4' ( t_4'
t_{56}')  T M_T>
< a|\psi_{12,123}^{1234}>\cr
& \equiv &  \delta \cdots  \sum_{\alpha_1} \delta \cdots  \int
\Pi_{i=1}^5
da_i a_i^2  {_b} < isospin| P_{45}^{isospin}| isospin>_a\cr
  & & _b < spin| P_{45}^{spin}| spin>_a < b_1 b_2 b_3 b_4 b_5|
P_{45}^{mom}
| a_1 a_2 a_3 a_4 a_5 >\cr
  & &  < a|\psi_{12,123}^{1234}>
\end{eqnarray}
The permutation is separated into the 3 spaces: isospin, spin,
and momentum, and  the $ \delta's\cdots$ are the strings of
Kronecker symbols.

The isospin- and spin-matrix elements can be calculated in a standard manner:
\begin{eqnarray}
& &   {_b} < isospin| P_{45}^{isospin}| isospin>_a
 =  \delta_{T T'} \delta_{M_T m_T'} \delta_{ t_{12} t_{12}'}
( - ) ^{ t_{12} + t_{1-4} + 1} \cr
&& \sqrt{ \hat t_3' \hat t_{34} \hat
    t_{1-4}
\hat t_{56} \hat t_4' \hat t_{56}'}
  \left\{{\begin{array}{*{20}c}
   t_{12} & \frac{1}{2} & t_3'  \\
   \frac{1}{2}  & t_{1-4}  & t_{34}  \\
\end{array}}\right\}
 \left\{{\begin{array}{*{20}c}
   t_3' & {\frac{1}{2}} & t_{1-4}  \\
 {\frac{1}{2}}    & {\frac{1}{2}} & t_{56}  \\
  t_4' & t_{56}' & T  \\
\end{array}}\right\}
\end{eqnarray}
\begin{eqnarray}
& &   {_b} < spin| P_{45}^{spin}| spin>_a
 =   \delta_{ s_{12} s_{12}'}
( - ) ^{ s_{12} + S + 1} \sqrt{ \hat I_3' \hat s_{34} \hat S \hat
    s_{56}
\hat I_4' \hat s_{56}'}\cr
 & & \left\{{\begin{array}{*{20}c}
   s_{12} & \frac{1}{2} & I_3'  \\
   \frac{1}{2}  & S  & s_{34}  \\
\end{array}}\right\}
 \left\{{\begin{array}{*{20}c}
   I_3' & {\frac{1}{2}} & S  \\
 {\frac{1}{2}}    & {\frac{1}{2}} & s_{56}  \\
  I_4' & s_{56}' & J  \\
\end{array}}\right\}
\end{eqnarray}

For the momentum space part we insert a complete basis and obtain
\begin{eqnarray}
  < b_1 b_2 b_3 b_4 b_5| P_{45}^{mom} & = &  \int d^3 b_1'
\cdots d^3 b_5'< b_1 b_2 b_3 b_4 b_5| \vec b_1' \vec b_2'\vec b_3'
\vec b_4'\vec b_5'>
   < \vec b_1' \vec b_2'\vec b_3'\vec b_4'\vec b_5'| P_{45}\cr
   &  = &  (\frac{1}{\sqrt{4 \pi}})^5 \int d \hat b_1' \cdots d
\hat  b_5'
   <  b_1 \hat b_1'   b_2 \hat  b_2' b_3 \hat  b_3' b_4 \hat  b_4' b_5
\hat  b_5'| P_{45}^{mom}
   \end{eqnarray}

The transposition $ P_{45}^{mom}  $ acting to the left leads to a
state with the same quantum  numbers but different meaning.
Particles 4 and 5 are interchanged. That state can be reexpressed
again in terms of the old b-state using the relation  between the  Jacobi
   momenta $ \vec b_1,\cdots \vec b_5 $ and the Jacobi momenta with
particles 4 and 5 interchanged. It results
\begin{eqnarray}
    < \vec b_1 \vec b_2 \vec b_3 \vec b_4 \vec b_5| P_{45}
    &=&   < \vec b_1, \frac{1}{2} ( \vec b_2 - \frac{1}{2} \vec b_3
-   \vec b_4 - \frac{3}{4} \vec b_5),
     \frac{1}{2} ( -\vec b_2 + \frac{3}{2} \vec b_3 -  \vec b_4
- \frac{3}{4} \vec b_5),\cr
   & & \frac{1}{2} ( -\vec b_2 - \frac{1}{2} \vec b_3 +  \vec b_4
- \frac{3}{4} \vec b_5),
    -   ( \vec b_2 + \frac{1}{2} \vec b_3 +  \vec b_4  -
\frac{1}{4}
\vec b_5)|
   \end{eqnarray}

Consequently
\begin{eqnarray}
& & < b_1 b_2 b_3 b_4 b_5|  P_{45}^{mom} =  (\frac{1}{\sqrt{4 \pi}})^5
\int d \hat b_1' \cdots d \hat b_5'\cr
& &   < b_1 \hat b_1', \frac{1}{2} ( b_2 \hat b_2' - \frac{1}{2} b_3
\hat b_3' -  b_4 \hat b_4' - \frac{3}{4}   b_5 \hat b_5'),
\frac{1}{2} ( - b_2 \hat b_2' + \frac{3}{2}  b_3 \hat b_3' - b_4 \hat
b_4'
- \frac{3}{4}  b_5 \hat b_5'),\cr
  & &   \frac{1}{2} ( - b_2 \hat b_2' - \frac{1}{2}  b_3 \hat b_3' +
b_4
\hat b_4'   - \frac{3}{4}  b_5 \hat b_5'),
-   (  b_2 \hat b_2' + \frac{1}{2} b_3 \hat b_3' +  b_4 \hat b_4'
-  \frac{1}{4} b_5 \hat b_5')|\label{87}
\end{eqnarray}

The state $ | \vec b_1 \cdots \vec b_5>$  can be reexpressed in terms
of the state $ | \vec a_1 \cdots \vec a_5>$  as
\begin{eqnarray}
  | \vec b_1 \vec b_2 \vec b_3 \vec b_4  \vec b_5> = |  \vec b_1 ,
  \frac{2}{3}
( \vec b_2 - \vec b_3), - \frac{1}{2} ( 2 \vec b_2 + \vec b_3),
   \vec b_4, \vec b_5>\label{91}
   \end{eqnarray}

Consequently using (\ref{87}) and ( \ref{91}) one obtains
\begin{eqnarray}
& & < b_1 b_2 b_3 b_4 b_5|  P_{45}^{mom} \int \Pi_{i=1}^5 da_i a_i^2
| a_1 a_2 a_3 a_4 a_5 > <  a_1 a_2 a_3 a_4 a_5 | \psi_{12,123}^{1234}>\cr
  & = &   \int d \hat  b_1' \cdots d \hat b_5'  <  b_1,  \frac{1}{3}
| 2 b_2 \hat b_2' - 2  b_3 \hat b_3'|,
   \frac{1}{2} |  \frac{1}{2} b_2 \hat b_2' + \frac{1}{4}  b_3 \hat
   b_3'
-  \frac{3}{4}  b_4  \hat b_4' - \frac{9}{8}   b_5 \hat b_5')|,\cr
  & &   \frac{1}{2} | - \frac{1}{2} b_2  \hat b_2' - \frac{1}{2}  b_3
\hat b_3' +  b_4 \hat b_4' - \frac{3}{4} b_5 \hat b_5')|,\cr
  & &    | b-2 \hat b_2' + \frac{1}{2}  b_3 \hat b_3' +  b_4 \hat b_4'
- \frac{1}{4} b_5 \hat b_5' | | \psi_{12,123}^{1234}>( \frac{1}{\sqrt{4 \pi}})^5
\end{eqnarray}

The angular integration over $ \hat b_1'$ yields directly $ 4 \pi $.
Further one can put $ \hat b_2' $ into the z-direction
  and $\hat b_3'$  into the x-z plane. Therefore
\begin{eqnarray}
& & < b_1 b_2 b_3 b_4 b_5|  P_{45}^{mom} \int \Pi_{i=1}^5 da_i a_i^2
| a_1 a_2 a_3 a_4 a_5 > <  a_1 a_2 a_3 a_4 a_5 | \psi_{12,123}^{1234}>\cr
 & = &  \sqrt{\pi}  \int  d cos \theta_3' d \hat  b_4' d \hat b_5'
    <  b_1,  \frac{1}{3} | 2 b_2 \hat b_2' - 2  b_3 \hat b_3'|,
     \frac{1}{2} |  \frac{1}{2} b_2 \hat b_2' + \frac{1}{4} b_3
\hat b_3'
-  \frac{3}{4}   b_4 \hat b_4' - \frac{9}{8}   b_5 \hat b_5')|,\cr
  & &   \frac{1}{2} | - \frac{1}{2} b_2  \hat b_2' - \frac{1}{2} b_3
\hat b_3'
+  b_4 \hat b_4' - \frac{3}{4} b_5 \hat b_5')|,
      | b_2 \hat b_2' + \frac{1}{2} b_3 \hat b_3' +  b_4 \hat b_4'
- \frac{1}{4} b_5 \hat b_5'| | \psi_{12,123}^{1234}>
\end{eqnarray}

This is an example where a 4-dimensional interpolation appears necessary.
Therefore if the momentum magnitudes are discretized choosing for each
one a certain grid, interpolations are inevitable.

In the 3- and 4-nucleon problems cubic Hermitean spline interpolation
turned out to be very efficient \cite{hueber}.
Thus a 5-dimensional interpolation, for instance,  has the form
\begin{eqnarray}
&&f_{ijklm} ( a_1 a_2 a_3 a_4 a_5)  =  \cr
&& \sum_{r=0}^3 \sum_{s = 0}^3 \sum_{t=0}^3 \sum_{u=0}^3
\sum_{v=0}^3 S_r( a_1 )S_s( a_2 )S_t( a_3 )S_u( a_4 )S_v( a_5 )
     f( a_{1r} a_{2s}a_{3t}a_{4u}a_{5v})
\label{80}
\end{eqnarray}
where $ ijklm $ denotes the 5-dimensional cubus around the point
$ a_1 a_2 a_3 a_4 a_5 $ and $ a_{ik} $ denotes the grid points for
the variable $a_i$.

For $ a_1 $,  for instance,
   the 4 grid points related to $a_i$ are   $ a_{10} < a_{11} \leq a_1
\leq  a_{12} < a_{13} $
 and similar for the other variables. Here for the sake of a simpler
notation we renumbered the grid points in that context.

Further $ f $ is the function to be interpolated
   and $ f_{ijklm} $ the interpolating one.  The spline functions
are given in Appendix \ref{appendix_b} and the conditions underlying that form
(\ref{80}) can be found in \cite{hueber}.

The coupled set (\ref{70})-(\ref{73}) in a matrix notation has  the
schematic structure
\begin{eqnarray}
        \eta ( E) \psi = K(E) \psi
\end{eqnarray}
where $E$ is the searched for energy eigenvalue at which the auxiliary
kernel eigenvalue $\eta(E)=1$. For 3- and 4-nucleon bound states
   a Lanczos type algorithm turned out to be very
efficient  \cite{stadler,nogga1}. Starting from an arbitrary
initial $\psi = \psi_0 $
 one generates by consecutive applications of $K$ a sequence of
 amplitudes $\psi_n$, which after orthogonalisation form a basis into which
    $\psi$ is expanded. It turn out that a reasonably small number
of K-applications (of the order of 10-20) is sufficient, which leads
to an algebraic eigenvalue problem
 of rather low dimension. Then the  energy is varied such that one
reaches $ \eta(E)=1 $.

If one regards the sub clustering underlying the 5 amplitudes only
two of them, $ \psi^{1234}_{12,123} $ and $ \psi^{1234}_{12;12,34} $,
  are related to the very approximative effective 3-body model of an
inert $ \alpha $-core and two neutrons. The total isospin
 quantum numbers of the $^6He $ ground state are $T=1$ and $M_T =-1$
(if we define the magnetic isospin quantum number of the neutron
as $-\frac{1}{2}$).

Now even these two amplitudes depending on the Jacobi momenta of type
a and b,  (\ref{74}) and (\ref{106}), respectively, are not restricted to
  $ t_4 =0 $ and $ t_{56}=1$ in case of $ \psi^{1234}_{12,123}$ but
also allow $ t_4 =1 $ and $ t_{56}=1 $, $ t_4 =2 $ and $ t_{56}=1 $
 and $ t_4 =1 $ and $t_{56}=0 $. Similarily the "deuteron-deuteron"
 like substructures of $  \psi^{1234}_{12;12,34}$ are  not restricted to
  $ t_{12}=t_{34}=0$ and $ t_{56}=1$ but also other two-body isospins
 are allowed which do not built up a $ t=0 $ $\alpha$-core. The amplitudes
    $ \psi^{123,45}_{12,123}$ going with Jacobi momenta of type c,
 (\ref{6.4}), refer to a 3-body together with a two-body sub
 clustering, which is
 also  not present in the effective 3-body model. The linear
 combinations  $ \psi^{125,34}_{12;12,34} + \psi^{345,12}_{12;12,34}$ refer
  again to a 3-body together  with a 2-body sub clustering, where
the underlying fragmentation related to 2-body fragments  differs from
  $ \psi^{123,45}_{12,123}$. And  again this is beyond the effective
3-body model. Finally  $ \psi^{12,34,56}_{12;12,34}$ allows for several
 additional 2-body sub clusters which are not contained in the
effective 3-body model, either.

If one would add another step in the Yakubovsky scheme, namely  to
2-body fragmentations $ a_2 $,  the resulting  amplitudes
$ \psi^{a_2}_{a_5,a_4,a_3}$
 would point to 5-body subclusters together with a single  nucleon or
to two 3-body subclusters (like $^3H-^3H$) in addition to
4-body and  2-body subclusters. Now all those additional structures
are of course also generated by the coupled system
(\ref{70})-(\ref{73})  ending with $a_3$ fragmentation, which we presented.
 In other words that Yakubovsky scheme is complete and delivers an
exact description of the 6-nucleon problem.

\section{Summary}
\label{sec_five}

The Yakubovsky equations have been derived long time ago by
O. A. Yakubovsky \cite{yakubovsky}. Therefore the application
presented here could have been given also long time ago. It is,
however, only now after the experiences with 3- and 4-nucleon
problems in the Faddeev-Yakubovsky schemes that the technical
 expertise has been developed in the last decades and the very strong
increase  of computer  power just recently achieved allows to attack
the 6-body problem in that exact formulation. Therefore we felt it is
timely to work out that scheme for that system. Another argument is
the development of nuclear forces in a systematic manner in the realm
of effective field theory and based on chiral symmetry. Two-, three-,
and four-nucleon forces have been derived
  consistently to each other and they are waiting to be applied in
light nuclear systems and checked against nature. Several tests in
that spirit already appeared \cite{ncsm,gfmc,nogga1,gitter} but for
the purpose of benchmarking the exact approach in the Yakubovsky
scheme is strongly recommended.

Here we restricted the formulation to two-nucleon forces only but
the inclusion of three-nucleon forces can easily be done like
pioneered in \cite{kamada2}.

In section \ref{sec_two} we used the general
basic and standard formulation,
where $a_n$ points to  n-body fragmentations for a system of $ N > n$ particles.

We worked out that scheme ending with $a_3$ in the spirit of the
usually applied though approximate effective $ \alpha -n -n $ 3-body model.

Though the step to $ a_2$ could be easily done, ending with $ a_3$
still includes exactly the whole dynamics.

The Pauli principle is then exactly incorporated in section
\ref{sec_three} leading to a set of 5 coupled equations for 5
 independent Yakubovsky components, which built up the total state.

The technical performance in a partial wave decomposition is only
touched in section \ref{sec_four}. The 5 different Jacobi momenta
as well as a necessary
 multi-dimensional interpolation scheme, like modified cubic
 Hermitean splines, are given.  For solving the high dimensional
 energy eigenvalue problem of the 5 coupled equations  we point
to  the Lanczos type algorithm, which turned out to be very efficient
in the 3- and 4-nucleon problem.  It remains to work out the partial
wave projected kernels, which is straightforward and can be carried
through along the lines cited above. An example for that is presented.

We expect that on the most modern supercomputers with parallel
architecture this formulation can be numerically mastered.

\begin{acknowledgments}
This work was partially supported
 by the Polish 2008-2011 science
funds as the research project No. N N202 077435
and by the Helmholtz
Association through funds provided to the virtual institute ``Spin
and strong QCD''(VH-VI-231).
\end{acknowledgments}