2\]) to determine the matrices. However, this method is not feasible for high dimensional systems due to the computational complexity. To overcome the problem, we first seek diagonal matrices, which are generally easier to manipulate. To this end, we look for operators $M'_i$ such that $$\begin{aligned}
&&\rho^d_B=\sum\limits_\mu M'_\mu \rho^d_A {M'_\mu}^\dagger,\label{rhod}\\
&&\sum\limits_\mu {M'_\mu}^\dagger M'_\mu=I. \label{mm3}\end{aligned}$$ To solve the above equations, we obtain the matrix operators $M'_\mu $ ($\mu=0,1,2,...N-1$) with the entries $$\begin{aligned}
 (M'_\mu)_{ij}= \sqrt{p^B_i}\delta_{i,j-\mu}+\sqrt{p^B_i}\delta_{i,j-\mu+N}~~,
\label{m}\end{aligned}$$ where $i,j=1,2,...N$[@Tong2].

With the matrices $M'_\mu $, we can construct the Kraus operators $M_\mu $ satisfying (\[rhoa\]) and (\[mm2\]). To see this, we use the relations, $$\begin{aligned}
\rho^d_A=U_A^\dagger\rho_AU_A,\label{rhoab3}\\
\rho^d_B=U_B^\dagger\rho_BU_B. \label{rhoab4}\end{aligned}$$ Substituting (\[rhoab3\]) and (\[rhoab4\]) into (\[rhod\]), we have $$\begin{aligned}
U_B^\dagger\rho_BU_B=\sum\limits_\mu M'_\mu U_A^\dagger\rho_A
U_A{M'_\mu}^\dagger,\end{aligned}$$ that is, $$\begin{aligned}
\rho_B=\sum\limits_\mu U_B M'_\mu U_A^\dagger\rho_A (U_BM'_\mu
U_A^\dagger)^\dagger.\end{aligned}$$ Let $$\begin{aligned}
M_\mu=U_BM'_\mu U_A^\dagger, \label{mmu}\end{aligned}$$ it is obvious that $M_\mu
2\ ]) to determine the matrices. However, this method is not feasible for high dimensional system due to the computational complexity. To get the best the problem, we first seek aslant matrices, which are generally easier to fudge. To this end, we look for operator $ M'_i$ such that $ $ \begin{aligned } 
 & & \rho^d_B=\sum\limits_\mu M'_\mu \rho^d_A { M'_\mu}^\dagger,\label{rhod}\\ 
 & & \sum\limits_\mu { M'_\mu}^\dagger M'_\mu = I. \label{mm3}\end{aligned}$$ To solve the above equations, we prevail the matrix operators $ M'_\mu $ ($ \mu=0,1,2,... N-1 $) with the entries $ $ \begin{aligned } 
  (M'_\mu)_{ij}= \sqrt{p^B_i}\delta_{i, j-\mu}+\sqrt{p^B_i}\delta_{i, j-\mu+N}~~, 
 \label{m}\end{aligned}$$ where $ i, j=1,2,... N$[@Tong2 ]. 

 With the matrices $ M'_\mu $, we can reconstruct the Kraus hustler $ M_\mu $ satisfying (\[rhoa\ ]) and (\[mm2\ ]). To visualize this, we use the relations, $ $ \begin{aligned } 
 \rho^d_A = U_A^\dagger\rho_AU_A,\label{rhoab3}\\ 
 \rho^d_B = U_B^\dagger\rho_BU_B. \label{rhoab4}\end{aligned}$$ Substituting (\[rhoab3\ ]) and (\[rhoab4\ ]) into (\[rhod\ ]), we suffer $ $ \begin{aligned } 
 U_B^\dagger\rho_BU_B=\sum\limits_\mu M'_\mu U_A^\dagger\rho_A 
 U_A{M'_\mu}^\dagger,\end{aligned}$$ that is, $ $ \begin{aligned } 
 \rho_B=\sum\limits_\mu U_B M'_\mu U_A^\dagger\rho_A (U_BM'_\mu 
 U_A^\dagger)^\dagger.\end{aligned}$$ Let $ $ \begin{aligned } 
 M_\mu = U_BM'_\mu U_A^\dagger, \label{mmu}\end{aligned}$$ it is obvious that $ M_\mu
2\]) to determine the matrices. Mowever, this metkid is iot feaaible fof high dimensional systems dne ti the computational complexkty. To ovvrcome thw privlem, we ficat seek diagohwl mctcices, which are generally easier to manhpjlcte. To this end, we look for operatorf $M'_i$ suvh that $$\begin{alidned}
&&\ghj^d_B=\shm\limits_\mu M'_\mu \rho^d_A {M'_\mu}^\dagger,\labem{rhod}\\
&&\suk\limits_\mu {M'_\mu}^\cagger M'_\mu=I. \label{mm3}\end{aligjed}$$ Ho solve the above equations, qe ofrain the matfix operators $M'_\mu $ ($\mu=0,1,2,...N-1$) with the entries $$\begin{aligned}
 (O'_\mu)_{ij}= \sqrt{p^B_i}\deotq_{i,j-\lo}+\sqrt{p^B_i}\delte_{i,j-\mu+N}~~,
\jabel{m}\end{alinmed}$$ whare $i,j=1,2,...N$[@Yong2].

With the mstrmces $M'_\mu $, we can construcv the Kraus operatorf $M_\mu $ sadiafying (\[rhoa\]) and (\[mn2\]). Ro seg thiv, we yse tht rxlafions, $$\hegmn{aligned}
\rhk^d_A=U_A^\dagger\eho_AU_A,\label{rhoab3}\\
\rho^c_B=T_V^\dagger\rho_BU_B. \label{whjab4}\end{aligned}$$ Substituting (\[rhoab3\]) and (\[rhmab4\]) into (\[rhod\]), we have $$\begib{aligned}
U_B^\dagger\rho_BU_H=\sum\limitf_\mu M'_\mu U_A^\dagger\rho_A
U_A{M'_\mu}^\dagger,\end{aligned}$$ that is, $$\begii{auigucq}
\fyo_H=\sum\limits_\mu U_B M'_\mu U_A^\dagger\rho_A (U_BM'_\mu
U_A^\daggew)^\sabgvr.\end{aligned}$$ Let $$\negin{aligned}
M_\mu=U_BM'_\ku U_S^\qagger, \label{mou}\end{alighed}$$ it is obvious hhat $M_\mo
2\]) to determine the matrices. However, this not for high systems due to the we first seek matrices, which are easier to manipulate. To this end, look for operators $M'_i$ such that $$\begin{aligned} &&\rho^d_B=\sum\limits_\mu M'_\mu \rho^d_A {M'_\mu}^\dagger,\label{rhod}\\ &&\sum\limits_\mu {M'_\mu}^\dagger \label{mm3}\end{aligned}$$ To solve the above equations, we obtain the matrix operators $M'_\mu $ with entries (M'_\mu)_{ij}= \label{m}\end{aligned}$$ where $i,j=1,2,...N$[@Tong2]. With the matrices $M'_\mu $, we can construct the Kraus operators $M_\mu $ (\[rhoa\]) and (\[mm2\]). To see this, we use relations, $$\begin{aligned} \rho^d_A=U_A^\dagger\rho_AU_A,\label{rhoab3}\\ \rho^d_B=U_B^\dagger\rho_BU_B. Substituting (\[rhoab3\]) and (\[rhoab4\]) into we $$\begin{aligned} U_B^\dagger\rho_BU_B=\sum\limits_\mu U_A^\dagger\rho_A that $$\begin{aligned} \rho_B=\sum\limits_\mu U_B U_A^\dagger\rho_A (U_BM'_\mu U_A^\dagger)^\dagger.\end{aligned}$$ Let $$\begin{aligned} M_\mu=U_BM'_\mu U_A^\dagger, \label{mmu}\end{aligned}$$ it is obvious that $M_\mu
2\]) to determine the matrices. HowEver, this meThod iS noT feAsIble For hIgh dimensional SYsteMs due to the computationaL compLeXIty. TO OvErcomE the proBLeM, WE fiRsT sEek DiAGoNal maTriCes, whicH are generaLly EaSier to manipuLAtE. To this end, We lOok for operatOrs $m'_i$ such ThAt $$\bEGin{alIgnEd}
&&\rho^D_B=\sum\lIMits_\mu m'_\mu \rho^d_A {M'_\Mu}^\DAgger,\lABel{rhod}\\
&&\SUM\lImitS_\mu {M'_\mu}^\dagger M'_\mu=I. \lABeL{Mm3}\end{aligned}$$ To Solve tHe ABoVE EquAtiOns, we obtaiN tHe matRIx operaTOrS $m'_\MU $ ($\mu=0,1,2,...n-1$) With the entrieS $$\begin{alignED}
 (M'_\mU)_{ij}= \sqrT{p^b_i}\dELta_{i,j-\mU}+\sqrt{P^B_I}\DelTa_{i,j-\mu+N}~~,
\labeL{m}\enD{aligned}$$ wHere $i,j=1,2,...n$[@tong2].

WitH The matrIces $M'_\mU $, we Can ConsTRuCt The krAUs oPErAtoRS $M_\mU $ satisfyInG (\[rHoa\]) anD (\[mm2\]). TO SEE This, We uSe thE relaTions, $$\begin{aliGneD}
\rho^D_a=U_A^\DaggeR\rho_Au_A,\laBeL{rhoaB3}\\
\rho^d_B=u_B^\dagGeR\rho_BU_B. \label{rhoAb4}\enD{aligned}$$ SUbsTiTutInG (\[rhoaB3\]) And (\[rhoAb4\]) iNto (\[Rhod\]), we hAve $$\begiN{AliGnED}
u_b^\dAgger\rho_BU_B=\sum\limiTs_\MU m'_\mU U_A^\daggeR\rho_A
U_a{m'_\mU}^\dAGger,\end{aLiGneD}$$ thaT IS, $$\begiN{aliGNeD}
\rho_B=\sum\Limits_\MU U_b M'_\Mu U_A^\dagGeR\rho_A (U_bM'_\Mu
U_a^\daGger)^\dAGger.\End{aliGned}$$ Let $$\bEgin{aLIgned}
M_\mu=U_BM'_\mu U_a^\Dagger, \label{mmU}\EnD{ALiGNed}$$ iT is Obvious that $m_\mu
2\]) to determine the matr ices. Howe ver,thi s m et hodis n ot feasible fo r  hig h dimensional systemsdue t ot he c o mp utati onal co m pl e x ity .To  ov er c om e the  pr oblem,we first s eek  d iagonal matr i ce s, which a regenerally ea sie r to m an ipu l ate.Tothisend, w e  lookfor opera to r s $M'_ i $ sucht h at  $$\ begin{aligned}
&& \ rh o ^d_B=\sum\limi ts_\mu  M ' _\ m u  \r ho^ d_A {M'_\m u} ^\dag g er,\lab e l{ r h o d}\ \ 
&&\sum\limit s_\mu {M'_\ m u}^ \dagge rM'_ \ mu=I.\labe l{ m m3} \end{aligne d}$$  To solve  the a b ove equ a tions,we obt ain  th e ma t ri xope ra t ors $M '_\ m u $  ($\mu=0 ,1 ,2 ,...N -1$) w i t h th e e ntri es $$ \begin{aligne d} (M' _ \mu )_{ij }= \s qrt{ p^ B_i}\ delta_ {i,j- \m u}+\sqrt{p^B_i} \del ta_{i,j-\ mu+ N} ~~, 
\ label { m}\end {al ign ed}$$ w here $i , j=1 ,2 , . . .N $[@Tong2].

With t he m at rices $M '_\mu$ ,we can cons tr uct  the K rausoper a to rs $M_\m u $ sa t is fy ing (\[ rh oa\])an d ( \[m m2\]) .  Tosee th is, we u se th e  relations, $$ \ begin{aligned } 
\ r h o^ d _A=U _A^ \dagger\rho _AU_ A ,\la bel{ r ho ab3 } \\
\r ho^d_ B= U _B ^ \dagger\rho_BU_B. \ la bel{rh oab4} \end{aligned} $$ Substit u t i ng (\[rh oab3 \ ]) and (\[rhoab4\ ]) in to (\[rhod \ ]), we h ave $ $\begin{ aligned}U _ B^\dagge r\r ho_ BU_ B=\ s u m\ limits_\mu M' _ \ mu U _A ^\dagge r\r ho_A
U_ A{M '_\ mu} ^\d ag ger,\end{ aligned} $$  t ha tis,  $$\b e gin{alig ne d}\r ho_ B=\su m \limit s_\mu  U_B  M '_ \ muU_A^\da g ge r \ rho_ A(U _BM' _\m uU_A^\ dagg e r)^ \dagger .\end{ali gne d }$$Le t$$\begi n{aligned}
M_ \m u=U_BM'_\m uU_A ^\dagg e r , \label {mmu}\end{aligned}$$ it is obvi ous  that  $M_ \mu
2\]) to_determine the_matrices. However, this method_is not_feasible_for high_dimensional_systems due to_the computational complexity._To overcome the problem,_we first seek_diagonal_matrices, which are generally easier to manipulate. To this end, we look for operators_$M'_i$_such that_$$\begin{aligned}
&&\rho^d_B=\sum\limits_\mu_M'_\mu_\rho^d_A {M'_\mu}^\dagger,\label{rhod}\\
&&\sum\limits_\mu {M'_\mu}^\dagger M'_\mu=I. \label{mm3}\end{aligned}$$_To solve the above equations,_we obtain_the matrix operators $M'_\mu $ ($\mu=0,1,2,...N-1$) with the_entries_$$\begin{aligned}
 (M'_\mu)_{ij}= \sqrt{p^B_i}\delta_{i,j-\mu}+\sqrt{p^B_i}\delta_{i,j-\mu+N}~~,
\label{m}\end{aligned}$$_where $i,j=1,2,...N$[@Tong2].

With the matrices $M'_\mu $, we can construct_the Kraus operators $M_\mu $ satisfying_(\[rhoa\]) and (\[mm2\])._To_see_this, we use the_relations, $$\begin{aligned}
\rho^d_A=U_A^\dagger\rho_AU_A,\label{rhoab3}\\
\rho^d_B=U_B^\dagger\rho_BU_B. \label{rhoab4}\end{aligned}$$ Substituting (\[rhoab3\]) and_(\[rhoab4\]) into (\[rhod\]), we have $$\begin{aligned}
U_B^\dagger\rho_BU_B=\sum\limits_\mu_M'_\mu U_A^\dagger\rho_A
U_A{M'_\mu}^\dagger,\end{aligned}$$ that is, $$\begin{aligned}
\rho_B=\sum\limits_\mu U_B M'_\mu_U_A^\dagger\rho_A (U_BM'_\mu
U_A^\dagger)^\dagger.\end{aligned}$$ Let $$\begin{aligned}
M_\mu=U_BM'_\mu U_A^\dagger, \label{mmu}\end{aligned}$$_it is obvious that $M_\mu