For calculating a local density matrix or in other words taking a partial trace has a neat trick (IMHO):
The Trick:
- Write the state given in density matrix form, if not already. Which means if $\vert \psi \rangle$ is the state, the density matrix is $\vert \psi \rangle \langle \psi \vert$.
- Turn all the outer-products into inner-products for the dimensions you don't want.
Example
- $$\vert \psi \rangle =\frac{|00\rangle + |01\rangle + |10\rangle}{\sqrt{3}}$$
$$\vert \psi \rangle \langle \psi \vert = \big(\frac{|00\rangle + |01\rangle + |10\rangle}{\sqrt{3}}\big)\big( \frac{\langle00| + \langle01| + \langle10\vert}{\sqrt{3}} \big) $$
$$=\big( \frac{\vert00\rangle\langle00\vert+\vert00\rangle\langle01\vert+\vert00\rangle\langle10\vert+\vert01\rangle\langle00\vert+\vert01\rangle\langle01\vert+\vert01\rangle\langle10\vert+\vert10\rangle\langle00\vert+\vert10\rangle\langle01\vert+\vert10\rangle\langle10\vert}{3}\big)$$
- Since you want the state of 2nd system, you would turn all the outer-products for 1st system to inner-prodcuts(This implies all the terms of form, $\langle0\vert1\rangle=\langle1\vert0\rangle=0$ and $\langle0\vert0\rangle=\langle1\vert1\rangle=1$,
$$\rho_2 = \big( \frac{1*\vert0\rangle\langle0\vert+1*\vert0\rangle\langle1\vert+0*\vert0\rangle\langle0\vert+1*\vert1\rangle\langle0\vert+1*\vert1\rangle\langle1\vert+0*\vert1\rangle\langle0\vert+0*\vert0\rangle\langle0\vert+0*\vert0\rangle\langle1\vert+1*\vert0\rangle\langle0\vert}{3}\big)$$
$$\rho_2 = \big( \frac{\vert0\rangle\langle0\vert+\vert0\rangle\langle1\vert+\vert1\rangle\langle0\vert+\vert1\rangle\langle1\vert+\vert0\rangle\langle0\vert}{3}\big)$$
$$\rho_2 = \big( \frac{2\vert0\rangle\langle0\vert+\vert0\rangle\langle1\vert+\vert1\rangle\langle0\vert+\vert1\rangle\langle1\vert}{3}\big)$$
Which is the state in the example you quote. Now, this will become second nature when you become little experienced and you wouldnt even be doing all this algebra!
Example on 3-qubit system:
- I deliberately choose an state with 2 terms just to show the proof of principle, you can work out a complicated state if you wish so!,
$$\vert \psi \rangle = \sqrt{\frac{2}{3}}\vert 011 \rangle + \sqrt{\frac{1}{3}}\vert100\rangle$$
$$ \vert \psi \rangle \langle \psi \vert =
\frac{2}{3}\vert011\rangle\langle011\vert
+\frac{1}{3}\vert100\rangle\langle100\vert
+\frac{\sqrt{2}}{3}\vert011\rangle\langle100\vert
+\frac{\sqrt{2}}{3}\vert100\rangle\langle011\vert$$
- Lets say you want the state of 2nd qubit, turn all the outer-products for 1 and 3 qubit to inner-products.
$$\rho_2=\frac{2}{3}\langle0\vert0\rangle\vert1\rangle\langle1\vert\langle1\vert1\rangle
+\frac{1}{3}\langle1\vert1\rangle\vert0\rangle\langle0\vert\langle0\vert0\rangle
+\frac{\sqrt{2}}{3}\langle0\vert1\rangle\vert0\rangle\langle0\vert\langle1\vert0\rangle
+\frac{\sqrt{2}}{3}\langle1\vert0\rangle\vert0\rangle\langle0\vert\langle0\vert1\rangle$$
$$\rho_2 = \frac{2}{3}\vert0\rangle\langle0\vert + \frac{1}{3}\vert1\rangle\langle1\vert$$
Explaination of trick:
Trace on a single qubit state(in computational basis), $\sigma$ is given by sum of diagonal elements,
$$Tr(\rho) = \sum_{i|i\in\{0,1\}}\langle i \vert\sigma\vert i\rangle.$$
Similarly for 2-qubit state $\rho$,
$$Tr(\rho) = \sum_{i,j |i,j\in\{0,1\} }\langle ij\vert\rho\vert ij\rangle.$$
Now, given this 2-qubit state if you want to know the local state of 1st qubit you would take a partial trace over 2nd qubit:
$$Tr_2(\rho) = \sum_{j|j\in\{0,1\}}\langle j \vert\rho\vert j\rangle.$$
In above you can see that if there was an element of form $\vert0\rangle\langle1\vert,\vert1\rangle\langle0\vert$ for 2nd qubit, they will go to zero while the terms of form $\vert0\rangle\langle0\vert,\vert1\rangle\langle1\vert$ goes to 1, and hence the trick!
and if you want to know the local state of 2nd qubit you would take a partial trace over 1st qubit:
$$Tr_1(\rho) = \sum_{i|i\in\{0,1\}}\langle i \vert\rho\vert i\rangle.$$