According to the Copenhagen interpretation, the process of measurement is described by the collapse of the wavefunction, which is a non-unitary process. The Hermiticity of a Hamiltonian guarantees that the time dependent wavefunction undergoes a unitary transformation, so the Schrödinger equation can not be applied to the measurement process, and the Born rule should be used instead. The Many Worlds Interpretation avoids this issue by removing the collapse of the wavefunction and holds that the measurement process can be explained by an ordinary interaction between two quantum systems. This results in a final entangled state between the system to be measured and the measurement device. As I understand it, the total state evolves according to the evolution operator \begin{equation} U=\lim_{t\rightarrow\infty}\mathcal{T}e^{-i/\hbar\int_{t_0}^t Hdt} \end{equation} where $H$ is the interaction Hamiltonian describing the measurement process.
I have tried deriving the unitary operator $U$ which describes the measurement, but I did not obtain a unitary operator. Let the total wavefunction be written as \begin{equation} |\Psi\rangle=|\psi\rangle\otimes|\phi\rangle\; , \end{equation} where $|\psi\rangle=\sum_n|\psi_n\rangle$ is the state of the measured system and $|\psi_n\rangle$ are the eigenstates of the measured observable. In this notation, state $|\phi_n\rangle$ describes the state of the measurement device which has measured state $|\psi_n\rangle$ of the system. The measurement operator thus acts upon the states as \begin{equation} U|\psi_n\rangle\otimes|\phi\rangle=|\psi_n\rangle\otimes|\phi_n\rangle\; , \quad \forall |\phi\rangle\; . \end{equation}
Using the above relation, the matrix elements of $U$ are given by: \begin{equation} U_{n_1,n_2;m_1,m_2}=\langle\psi_{n_1}|\otimes\langle\phi_{n_2}|U|\psi_{m_1}\rangle\otimes|\phi_{m_2}\rangle=\langle\psi_{n_1}|\otimes\langle\phi_{n_2}|\psi_{m_1}\rangle\otimes|\phi_{m_1}\rangle=\delta_{n_1,m_1}\delta_{n_2,m_1}\; . \end{equation}
A unitary operator satisfies that $UU^\dagger=1$, that is to say $(UU^\dagger)_{n_1,n_2;m_1,m_2}=\delta_{n_1,m_1}\delta_{n_2,m_2}$. I have tried computing this property and obtained the following result: \begin{equation} (UU^\dagger)_{n_1,n_2;m_1,m_2}=\sum_{a,b}U_{n_1,n_2;a,b}(U_{m_1,m_2;a,b})^*=\sum_{a,b}\delta_{n_1,a}\delta_{n_2,a}\delta_{m_1,a}\delta_{m_2,a}\\ =\sum_{b}\delta_{n_1,m_1}\delta_{n_2,m_2}\delta_{n_1,n_2}\neq\delta_{n_1,m_1}\delta_{n_2,m_2}\; . \end{equation}
What is the reason for this inconsistency? I do not know if I made any wrong assumption about the MWI or maybe I did not define $U$ correctly. All the explanations I have read about the MWI do not delve into the mathematics behind the measurement, is there any book that covers this topic mathematically?
Edit: After posting the question I have realised that this derivation is equivalent to the quantum cloning problem, which is forbidden by the no cloning theorem. Therefore, there must be something wrong with the assumptions made. Maybe there are multiple degenerate states $|\phi_{n,g}\rangle$ corresponding to the measurement of $|\psi_n\rangle$?