Yes. It would be worth going through this yourself, but here's my calculation:
The pure spin state for an anti-correlated pair of spin-1/2 particles is: [itex]\frac{1}{\sqrt{2}} (|U_z\rangle |D_z\rangle - |D_z\rangle |U_z\rangle)[/itex]. The corresponding density matrix is
[itex]\sum_{a i b j} \rho_{a i b j} |a\rangle |i \rangle \langle b| \langle j|[/itex]
where [itex]a, b[/itex] refers to the first particle, and [itex]i, j[/itex] refers to the second particle, and all the indices range over the possibilities [itex]U_z, D_z[/itex], and where the nonzero coefficients [itex]\rho_{a i b j}[/itex] are given by:
- [itex]\rho_{U_z D_z U_z D_z} = \rho_{D_z U_z D_z U_z} =\frac{1}{2}[/itex]
- [itex]\rho_{U_z D_z D_z U_z} = \rho_{D_z U_z U_z D_z} =-\frac{1}{2}[/itex]
The reduced matrix for particle 2 is:
[itex]\sum_{i j} \rho_{i j} |i\rangle \langle j|[/itex]
where [itex]\rho_{ij} \equiv \sum_{a} \rho_{a i a j}[/itex]. The nonzero values for [itex]\rho_{ij}[/itex] are:
- [itex]\rho_{U_z U_z} = \frac{1}{2}[/itex]
- [itex]\rho_{D_z D_z} = \frac{1}{2}[/itex]
So [itex]\rho = \frac{1}{2} |U_z\rangle \langle U_z| + \frac{1}{2} |D_z\rangle \langle D_z|[/itex]
Using the usual two-component spinors: [itex]|U_z\rangle = \left( \begin{array} \\ 1 \\ 0 \end{array} \right)[/itex] and [itex]|D_z\rangle = \left( \begin{array} \\ 0 \\ 1 \end{array} \right)[/itex], this means:
[itex]\rho = \ \left( \begin{array} \\ \frac{1}{2} & 0 \\ 0 & \frac{1}{2} \end{array} \right)[/itex].
Now, go through the whole thing again using the basis
[itex]|U_x\rangle = \frac{1}{\sqrt{2}} (|U_z\rangle + |D_z\rangle)[/itex]
[itex]|D_x\rangle = \frac{1}{\sqrt{2}} (|U_z\rangle - |D_z\rangle)[/itex]
In this basis, the pure state is [itex]\frac{1}{\sqrt{2}}(|U_x\rangle |D_x\rangle - |D_x\rangle |U_x\rangle)[/itex]. (You can work this out.) This is exactly the same as for the z-basis, except the names of the states have changed.
The composite density matrix is analogous to the case for the z-basis. Then when you perform the trace, you end up with the reduced matrix:
[itex]\rho = \frac{1}{2} |U_x\rangle \langle U_x| + \frac{1}{2} |D_x\rangle \langle D_x|[/itex]
Using the spinor representation: [itex]|U_x\rangle = \left(\begin{array} \\ \frac{1}{\sqrt{2}} \\ \frac{1}{\sqrt{2}} \end{array} \right)[/itex] and [itex]|D_x\rangle = \left(\begin{array} \\ \frac{1}{\sqrt{2}} \\ \frac{-1}{\sqrt{2}} \end{array} \right)[/itex], this means:
[itex]|U_x\rangle \langle U_x| = \left( \begin{array} \\ \frac{1}{\sqrt{2}} \\ \frac{1}{\sqrt{2}} \end{array} \right) \left( \begin{array} \\ \frac{1}{\sqrt{2}} & \frac{1}{\sqrt{2}} \end{array} \right) = \left( \begin{array} \\ \frac{1}{2} & \frac{1}{2} \\ \frac{1}{2} & \frac{1}{2} \end{array} \right)[/itex][itex]|D_x\rangle \langle D_x| = \left( \begin{array} \\ \frac{1}{\sqrt{2}} \\ \frac{-1}{\sqrt{2}} \end{array} \right) \left( \begin{array} \\ \frac{1}{\sqrt{2}} & \frac{-1}{\sqrt{2}} \end{array} \right) = \left( \begin{array} \\ \frac{1}{2} & \frac{-1}{2} \\ \frac{-1}{2} & \frac{1}{2} \end{array} \right)[/itex]
[itex]\rho = \frac{1}{2} |U_x\rangle \langle U_x| + \frac{1}{2} |D_x\rangle \langle D_x| = \left( \begin{array} \\ \frac{1}{2} & 0 \\ 0 & \frac{1}{2} \end{array} \right)[/itex]
You end up with the same reduced density matrix, regardless of which basis you work with.