How Do Indices Determine Positions in a Quantum Density Matrix?

[doggydan42]

I am reading Leonard Susskind's Theoretical Minimum book on Quantum Mechanics. Excercise 7.4 is as follows:

Calculate the density matrix for $|\Psi\rangle = \alpha|u\rangle + \beta|d\rangle$.

Answer:
$$ \psi(u) = \alpha, \quad \psi^*(u) = \alpha^* \\
\psi(d) = \beta, \quad \psi^*(d) = \beta^*\\
\rho_{a'a} = \begin{pmatrix}
\alpha^*\alpha & \alpha^*\beta \\
\beta^*\alpha & \beta^*\beta
\end{pmatrix}$$

From my understanding, if we use $\rho_{a'a} = \psi(a')\psi^*(a)$, then I notice that the matrix is equivalent to:

$$\rho_{a'a} = \begin{pmatrix}
\psi^*(u)\psi(u) & \psi^*(u)\psi(d) \\
\psi^*(d)\psi(u) & \psi^*(d)\psi(d)
\end{pmatrix} = \begin{pmatrix}
\rho_{uu} & \rho_{du} \\
\rho_{ud} & \rho_{dd}
\end{pmatrix}
$$

I was wondering how the indices affect the position of each $\rho_{a'a}$ in the density matrix.

In other words, why is the matrix not:

$$
\rho_{a'a} = \begin{pmatrix}
\alpha^*\alpha & \beta^*\alpha \\
\alpha^*\beta & \beta^*\beta
\end{pmatrix}$$

 

[Orodruin]

Because it is not defined that way and because it is not symmetric (it is Hermitian).

 

[doggydan42]

Because it is not defined that way and because it is not symmetric (it is Hermitian).

But both matrices are Hermitian, would it be defined such that:

$$\rho_{aa'} = \begin{pmatrix}
\rho_{11} & \rho_{12} \\
\rho_{21} & \rho_{22}
\end{pmatrix}$$

And u would correlate to the 1 while d correlates to the 2, since $|u\rangle = \begin{pmatrix} 1 \\ 0 \end{pmatrix}$ while $|d\rangle$ is orthonormal?

 

[DrClaude]

The density operator for state $| \Psi \rangle$ is defined as $\hat{\rho} = | \Psi \rangle \langle \Psi|$. When you go to a density matrix representation, you can choose the order of the components, e.g., for $|\Phi \rangle = a | u \rangle + b | d \rangle$,
$$
\begin{pmatrix} a \\ b \end{pmatrix} \mbox{ or }\begin{pmatrix} b \\ a \end{pmatrix}.
$$
Once that order is chosen (the one on the left here), the position of the elements of the density matrix is fixed, so that the result of $\hat{\rho} |\Phi \rangle$ is the same as the corresponding matrix-vector multiplication.

 

[doggydan42]

The density operator for state $| \Psi \rangle$ is defined as $\hat{\rho} = | \Psi \rangle \langle \Psi|$. When you go to a density matrix representation, you can choose the order of the components, e.g., for $|\Phi \rangle = a | u \rangle + b | d \rangle$,
$$
\begin{pmatrix} a \\ b \end{pmatrix} \mbox{ or }\begin{pmatrix} b \\ a \end{pmatrix}.
$$
Once that order is chosen (the one on the left here), the position of the elements of the density matrix is fixed, so that the result of $\hat{\rho} |\Phi \rangle$ is the same as the corresponding matrix-vector multiplication.

That makes sense, since when you find the eigen values of of the density matrix with the state vector and make it depend on the order of components, you get 1 for both options.

Thank you.

 

[vanhees71]

Well, let's do the calculation in Dirac's way, which is much simpler, because you don't need to remember so much. The statistical operator of a pure state, given by one representative state ket is
$$\hat{\rho}=|\psi \rangle \langle \psi|=(\alpha |u \rangle + \beta |d \rangle) (\alpha^* \langle u| + \beta^* \langle d|).$$
Now by definition in a basis representation, it's written as
$$\hat{\rho} = \sum_{ij} |i \rangle \langle j | \rho_{ij}.$$
Now you multiply out the ket-bra products, ans you get
$$\rho_{uu}=\alpha \alpha^*, \quad \rho_{ud}=\alpha \beta^*, \quad \rho_{du}=\beta \alpha^*, \quad \rho_{dd}=\beta \beta^*.$$
There's no need to write this in matrix form, but usually in a basis representation you write the vectors as column vectors and the co-vectors as row vectors. In this sense we map
$$|\psi \rangle \rightarrow \psi=\begin{pmatrix}\alpha \\ \beta \end{pmatrix}.$$
The corresponding co-vector is then
$$\psi^{\dagger}=(\alpha^* \quad \beta^*).$$
Then the corresponding projector is according to the usual rules of matrix multiplication
$$\hat{\rho}=\psi \psi^{\dagger} = \begin{pmatrix} \rho_{uu} & \rho_{ud} \\ \rho_{du} & \rho_{dd} \end{pmatrix}=\begin{pmatrix}\alpha \\ \beta \end{pmatrix} (\alpha^* \quad \beta^*) = \begin{pmatrix} \alpha \alpha^* & \alpha \beta^* \\
\beta \alpha^* & \beta \beta^* \end{pmatrix}.$$