One can directly jump to the question and skip the first two paragraphs (motivation for the question)
Consider a system in a thermal gas of $N$ particles. If I want to study the reduced dynamics of the system alone, I need to find the partial trace (over the environment) of the total density matrix of the system and the $N$ environmental particles $\hat \rho_{tot}$.
If however the change in reduced density matrix for a system in a thermal gas by a single collision is $\Delta \hat \rho^{red}$ then for a very dilute gas where collisions can be looked at as "particle by particle" the total change in reduced density matrix should simply be $N \times \Delta \hat \rho^{red}$ and therefore, if one calculates the change in reduced density matrix for a system from system + just a "single colliding particle" density matrix (as opposed to the total density matrix), it should be enough.
I however want to see if I can reach $N \times \Delta \rho^{red}$ starting from writing the total density matrix for all $N$ particles.
So, I am asking the following question so that I can try to do that.
My question is:
If $$\hat \rho_{tot} = (\hat \rho)^{\otimes N}$$
then if I know $$\hat \rho^{red} = Tr_{E}(\hat \rho)$$
Is there a straight-forward way to find: $$\hat \rho_{tot}^{red} = Tr_{E}(\underbrace{\hat \rho \otimes \hat \rho ..... \otimes \hat \rho }_{N \ times} )$$
from the knowledge of $\hat \rho^{red}$?
Is simply $$\hat \rho_{tot}^{red} = (\hat \rho^{red})^{\otimes N}$$?