In Page 365, Operator-sum representation, Chapter 8, Quantum Computation and Quantum Information by Nielsen and Chuang, it is given that
Given a trace-preserving quantum operation expressed in the operator-sum representation, $\mathcal{E}(ρ)=\sum_k E_kρE_k^\dagger$, we can construct a physical model for it in the following way. From $(8.10)$, we want $U$ to satisfy $$ E_k=\langle e_k|U|e_0\rangle $$ where $U$ is some unitary operator, and $|e_k\rangle$ are orthonormal basis vectors for the environment system. Such a $U$ is conveniently represented as the block matrix $$ U=\begin{bmatrix} [E_1] & \cdots & \cdots \\\\ [E_2] & \cdots & \cdots \\\\ \vdots & \ddots & \vdots \\ \end{bmatrix} $$ in the basis $|e_k\rangle$. Note that the operation elements $E_k$ only determine the first block column of this matrix (unlike elsewhere, here it is convenient to have the first label of the states be the environment, and the second, the principal system). Determination of the rest of the matrix is left up to us; we simply choose the entries such that U is unitary. By the results of Chapter 4, $U$ can be implemented by a quantum circuit.
We have a principal system $Q$ and an environment $E$ and $U$ is a unitary operator acting on the combined system $QE$. \begin{align} \mathcal{E}(\rho)&=tr_E\bigg(U(\rho\otimes|e_0\rangle\langle e_0|)U^\dagger\bigg)\\ &=\sum_k(I\otimes\langle e_k|)(U(\rho\otimes|e_0\rangle\langle e_0|)U^\dagger)(I\otimes |e_k\rangle)\\ &=\sum_k\langle e_k|(U(\rho\otimes|e_0\rangle\langle e_0|)U^\dagger)|e_k\rangle\\ &=\sum_k(I\otimes\langle e_k|)U(I\otimes|e_0\rangle)\rho(I\otimes\langle e_0|)U^\dagger(I\otimes|e_k\rangle)\\ &=\sum_k\langle e_k|U|e_0\rangle\rho\langle e_0|U^\dagger|e_k\rangle\\ &=\sum_k E_k \rho E_k^\dagger \end{align} where $E_k=(I\otimes\langle e_k|)U(I\otimes|e_0\rangle)=\langle e_k|U|e_0\rangle$ is an operator on the state space of the principal system $Q$, $|e_k\rangle$ are the orthonormal basis vectors of the environment, $|e_0\rangle$ be the initial state of the environment.
Please check Prove that $tr\Big(\sum_k E_k^\dagger E_k\rho\Big)=1$ for all $\rho$ implies $\sum_k E_k^\dagger E_k=I$ for the derivation.
In the problem, $E_k=(I\otimes\langle e_k|)U(I\otimes|e_0\rangle)=\langle e_k|U|e_0\rangle$ where $|e_k\rangle$ are some orthonormal basis vectors and $|e_0\rangle$ be the initial state of the environment $E$.
How do we say that $U$ is represented as the block matrix $$ U=\begin{bmatrix} [E_1] & \cdots & \cdots\\ [E_2] & \cdots & \cdots\\ [E_3] & \cdots & \cdots\\ \vdots & \vdots & \vdots \end{bmatrix} $$ in the basis $|e_k\rangle$ ?
My observation:
Given an operator $A:V\to V$ $$ A=\begin{bmatrix}a_{00}&a_{01}\\a_{10}&a_{11}\end{bmatrix}=\begin{bmatrix}\langle0|A|0\rangle&\langle0|A|1\rangle\\\langle1|A|0\rangle&\langle1|A|1\rangle\end{bmatrix}=\sum_{i,j}a_{j}|i\rangle\langle j| $$ If we had an operator $A:V\to V$ and $\{|v_i\rangle\}$ constitute a basis of $V$ then
Let $\{|v_i\rangle\}$ be an orthonormal basis then the change of basis matrix is $P_v=V^{-1}=V^\dagger=\begin{bmatrix}\langle v_0|\\\langle v_1|\end{bmatrix}$
The matrix $A$ in the orthonormal basis $\{|v_i\rangle\}$ is, $$ B=[A]_v=P_vAP_v^{-1}=\begin{bmatrix}\langle v_0|\\\langle v_1|\end{bmatrix}A\begin{bmatrix} |v_0\rangle& |v_1\rangle\end{bmatrix}=\begin{bmatrix}\langle v_0|A|v_0\rangle&\langle v_0|A|v_1\rangle\\\langle v_1|A|v_0\rangle&\langle v_1|A|v_1\rangle\end{bmatrix} $$ $$ A=P_v^\dagger BP_v=\begin{bmatrix} |v_0\rangle& |v_1\rangle\end{bmatrix}\begin{bmatrix}\langle v_0|A|v_0\rangle&\langle v_0|A|v_1\rangle\\\langle v_1|A|v_0\rangle&\langle v_1|A|v_1\rangle\end{bmatrix}\begin{bmatrix}\langle v_0|\\\langle v_1|\end{bmatrix}\\ =\sum_{i,j}\langle v_i|A|v_j\rangle|v_i\rangle\langle v_j| $$ That means, $$ A=\sum_{i,j}a_{ij}|i\rangle\langle j|=\sum_{i,j}b_{ij}|v_i\rangle\langle v_j| $$ where $a_{ij}=\langle i|A|j\rangle$ and $b_{ij}=\langle v_i|A|v_j\rangle$
This is clear!
So the $(i,j)^{th}$ term of $A$ in the basis $|v_i\rangle$ is $b_{ij}=\langle v_i|A|v_j\rangle$ and $B=[A]_v=\sum_{i,j}b_{ij}|i\rangle\langle j|$
Unlike $\langle v_i|W|v_j\rangle$ which is a number and $|v_j\rangle$ is a vector, the term $E_k=(I\otimes\langle e_k|)U(I\otimes|e_0\rangle)=\langle e_k|U|e_0\rangle$ is an operator acting on the first system (principal system $Q$) and $I\otimes|e_0\rangle=|e_0\rangle,I\otimes|e_k\rangle=|e_k\rangle$ are matrices.
My Attempt
Thanks @glS for the hint.
\begin{align} E_k&=(I\otimes\langle e_k|)U(I\otimes|e_0\rangle)=\langle e_k|U|e_0\rangle\\ \end{align} If we interchange the order of systems, i.e., the first system is the environment $E$ and second is the principal system $Q$, then
Decomposing $U$ and $(\langle e_k|\otimes I)U(|e_0\rangle\otimes I)$ in terms of the orthonormal basis vectors of the respective systems as,
\begin{align} &(\langle e_k|\otimes I)U(|e_0\rangle\otimes I)=\\ &=(\langle e_k|\otimes I)(\sum_{m,n,m',n'}\mu_{m,n,m',n'}(|e_n\rangle\otimes|p_m\rangle)(\langle e_{n'}|\otimes\langle p_{m'}|))(|e_{0}\rangle\otimes I)\\ &=\langle e_k|\otimes I\bigg[\sum_{m,n,m',n'}\mu_{m,n,m',n'}|e_n\rangle\langle e_{n'}|\otimes|p_m\rangle\langle p_{m'}|\bigg]|e_{0}\rangle\otimes I\\ &=\sum_{m,n,m',n'}\mu_{m,n,m',n'}\langle e_k|e_n\rangle\langle e_{n'}|e_{0}\rangle\otimes|p_m\rangle\langle p_{m'}|\\ &=\sum_{m,n,m',n'}\mu_{m,n,m',n'}\delta_{k,n}\delta_{n',0}\otimes|p_m\rangle\langle p_{m'}|\\ &=\sum_{m,m'}\mu_{m,k,m',0}|p_m\rangle\langle p_{m'}|\\ \end{align}
\begin{align} U&=\sum_{m,n,m',n'}\mu_{m,n,m',n'}(|e_n\rangle\otimes|p_m\rangle)(\langle e_{n'}|\otimes\langle p_{m'}|)\\ &=\sum_{m,n,m',n'}\mu_{m,n,m',n'}|e_n\rangle\langle e_{n'}|\otimes|p_m\rangle\langle p_{m'}|\\ &=\sum_{n,n'}|e_n\rangle\langle e_{n'}|\otimes\sum_{m,m'}\mu_{m,n,m',n'}|p_m\rangle\langle p_{m'}|\\ \end{align}
If I define the change of basis matrix $P=P_E\otimes P_p$ then (extending the ideas in the previous section) the operators in the new basis are, $$ (\langle e_k|\otimes I)U(|e_0\rangle\otimes I)'=\sum_{m,m'}\mu_{m,k,m',0}|m\rangle\langle {m'}|\\ $$ \begin{align} U'&=\sum_{m,n,m',n'}\mu_{m,n,m',n'}|n\rangle\langle {n'}|\otimes|m\rangle\langle {m'}|\\ &=\sum_{n,n'}|n\rangle\langle {n'}|\otimes \sum_{m,m'}\mu_{m,n,m',n'}|m\rangle\langle {m'}|\\ &=\sum_{n,n'}|n\rangle\langle {n'}|\otimes\sigma_{n,n'} \end{align} where $\{|n\rangle\}$ and $\{|m\rangle\}$ are the standard basis such that $|n\rangle\langle {n'}|$ is the matrix with $1$ as its $(n,n')$ entry, and $\sigma_{n,n'}=\sum_{m,m'}\mu_{m,n,m',n'}|m\rangle\langle {m'}|$ is the $(n,n')^{th}$ block entry of $U'$.
Therefore, the $(k,0)^{th}$ block entry of $U'$ is
$\sigma_{k,0}=\sum_{m,m'}\mu_{m,k,m',0}|m\rangle\langle {m'}|=(\langle e_k|\otimes I)U(|e_0\rangle\otimes I)'$