I've encounter two different approaches to the Unruh effect and I feel like they are not consistent with one another.
Bogoliubov Transformation
In this approach the basic statement is that the vacuum state for a Minkowski observer $|0\rangle_M$ is written as a superposition of excited states in the basis of the Rindler observer
\begin{equation} |0\rangle_M=\sum_{n_1=0}^\infty \dots \sum_{n_r=0}^\infty f(n_1,\dots,n_r) |n_1,\dots,n_r\rangle_R \end{equation}
Schematically, $n_i$ represents the number of excitations in the $i$ state. $r$ represents the last state but it could be infinite too. In this approach, it is straightforward to calculate the expectation value of the number of particles in the Rindler state $k$ $\hat{n}_k^{\text{(Rindler)}}$ in the Minkowski vacuum
\begin{equation} \langle0|_M \hat{n}_k^{\text{(Rindler)}}|0\rangle_M=\sum_{n_1=0}^\infty \dots \sum_{n_k=0}^\infty \dots \sum_{n_r=0}^\infty f^2 (n_1,\dots,n_k,\dots,n_r) n_k \end{equation}
given the right function $f$ it is possible for this expectation value to give a thermal number of particles, that is:
\begin{equation} \langle0|_M \hat{n}_k^{\text{(Rindler)}}|0\rangle_M=\frac{1}{e^{\frac{k}{2\pi a}}-1} \end{equation}
where $a$ is the acceleration of the Rindler obserber. Note that I don't really care about the function $f$. I'm just checking that it is possible to have a thermal expectation value for the number operator of one basis if we calculate it respect to a state from a different basis. Also, note that the fact that I'm writting the Minkowski vacuum as a sum of particle states in the Rindler basis means (I think) that the Rindler states form a complete basis or in other words, that the two Hilbert spaces have the same size.
Tracing the left wedge in Rindler space
This method is a bit more abstract and much more general. One starts with the Minkowski vacuum and traces over everything not covered by the Rindler coordinates, leaving us with the reduced density matrix on that region of spacetime, which is
\begin{equation} \rho_{\text{Rindler}}=e^{-2\pi K} \end{equation}
where $K$ is the charge associated with boosts. Now, if I remember quantum mechanics correctly, if a state is pure, it is pure for every choice of basis so the statement $|0\rangle_M$ becomes a thermal (i.e. mixed) state in the Rindler basis can't be true. It can only be true if we are also tracing over the patch of spacetime not covered by Rindler coordinates.
So the question is: How are these two methods giving the same thermal result? where is the "tracing out" in the first one? Is it implicitely using that the Rindler basis is smaller than the Minkowski basis?
Note: Even if the Rindler basis is indeed smaller than the Minkowski basis, it doesn't look like it is necessary in the derivation I provided. The question remains: If we calculate the expectation value of the number of particles operator of one basis respect to a pure state on a different basis for the same Hilbert space, can we obtain a thermal distribution?