Suppose we have a qubit in state $| \Psi \rangle = \alpha | 0 \rangle + \beta | 1 \rangle$
Suppose we expose this to decoherence, which we will express as the state $| R \rangle$ such that $$| 0 \rangle| R \rangle \rightarrow | 0 \rangle| R_0 \rangle$$ $$| 1 \rangle| R \rangle \rightarrow | 1 \rangle| R_1 \rangle$$
Where $| R \rangle, | R_0 \rangle$ and $| R_1 \rangle$ are all normalised states.
I'm trying to work how the density operator of the qubit changes when we apply $R$. If we consider the density operator for $|\Psi \rangle$ $$\rho=|\Psi \rangle \langle\Psi |= \Big( \begin{matrix} \alpha^2 & \alpha \beta \\ \alpha \beta & \beta^2 \end{matrix} \Big)$$
Assuming alpha and beta are real.
Next, we apply $|R\rangle$ to out qubit, $$\rho=|\Psi \rangle \langle\Psi | \rightarrow |\Psi \rangle |R\rangle \langle\Psi | \langle R|$$ $$=\alpha^2 |0 \rangle |R_0\rangle \langle 0 | \langle R_0|+\alpha\beta( |0 \rangle |R_0\rangle \langle 1 | \langle R_1|+ |1 \rangle |R_1\rangle \langle 0 | \langle R_0| )+\beta^2 |1 \rangle |R_1\rangle \langle 1 | \langle R_1|$$
Next would we take the reduced density operator of $|\Psi \rangle$ to find it's density operator? This would produce a $2 \times 2$ density operator, which I'm hoping to express in terms of $\rho=|\Psi \rangle \langle\Psi |$
Is this the correct way to think about the relation between density operators and reduced density operators for entangled states? What about unentangled states?
Also would anyone be able to briefly explain what decoherence is and why we can describe it as such a operator $|R \rangle$? Any help would be greatly appreciated.