# Is a permutation of coordinates or labels really equivalent?

To construct a N-body anti-symmetric wave function some derivations start with the requirement that the N-body wave function should be anti-symmetric under a permutation of coordinates, other derivations start with the requirement that the total wave function should be anti-symmetric under a permutation of labels or states, for example this derivation . I can see the equivalence is trivial if you assume you can freely change the order of the wavefunctions. This is true in the case the wavefunctions are just scalars but what if you are dealing with multicomponent wavefunctions. For example relativistic wavefunctions which are 4x1 matrices.

A simple example of my problem is the following, suppose we want to write down a general 2 body wavefunction in terms of single particle wavefunctions, \begin{align} \Psi(x_1,x_2) \propto \varphi_{a}(x_1)\varphi_b(x_2) \end{align} Where $x_i$ denote the spatial,spin,... coordinates, and $a,b$ denote single particle eigenstates. If we permute the coordinates in order to derive a antisymmetric wavefunction we get, \begin{align} \frac{1}{\sqrt{2}} \left(\varphi_{a}(x_1)\varphi_b(x_2) - \varphi_{a}(x_2)\varphi_b(x_1) \right) \end{align} while a permutation in states gives, \begin{align} \frac{1}{\sqrt{2}} \left(\varphi_{a}(x_1)\varphi_b(x_2) - \varphi_{b}(x_1)\varphi_a(x_2) \right) \end{align} Obviously above expressions are identical for scalar wavefunctions. When I try to check the normalisation of the two particle wavefunction (assuming the single particle states are orthormal and general $N \times 1$ matrices) I get for the label permutation, \begin{align} \int | \Psi(x_1,x_2) |^{2} &= \frac{1}{2} \int \varphi_{b}^{\dagger}(x_2) \varphi_{a}^{\dagger}(x_1) \varphi_{a}(x_1)\varphi_b(x_2) - \frac{1}{2} \int \varphi_{b}^{\dagger}(x_2) \varphi_{a}^{\dagger}(x_1)\varphi_{b}(x_1)\varphi_a(x_2) \\ &- \frac{1}{2} \int \varphi_{a}^{\dagger}(x_2) \varphi_{b}^{\dagger}(x_1) \varphi_{a}(x_1)\varphi_b(x_2) + \frac{1}{2} \int \varphi_{a}^{\dagger}(x_2) \varphi_{b}^{\dagger}(x_1)\varphi_{b}(x_1)\varphi_a(x_2) \\ &= \frac{1}{2} - 0 - 0 + \frac{1}{2} = 1 \end{align} Which is the expected result (note that the integral is assumed to integrate over all continuous degrees of freedom, or sum over the discrete ones). For the coordinate permutation however I get, \begin{align} \int | \Psi(x_1,x_2) |^{2} &= \frac{1}{2} \int \varphi_{b}^{\dagger}(x_2) \varphi_{a}^{\dagger}(x_1) \varphi_{a}(x_1)\varphi_b(x_2) - \frac{1}{2}\int \varphi_{b}^{\dagger}(x_2) \varphi_{a}^{\dagger}(x_1)\varphi_{a}(x_2)\varphi_b(x_1) \\ &- \frac{1}{2}\int \varphi_{b}^{\dagger}(x_1) \varphi_{a}^{\dagger}(x_2) \varphi_{a}(x_1)\varphi_b(x_2) + \frac{1}{2}\int \varphi_{b}^{\dagger}(x_1) \varphi_{a}^{\dagger}(x_2)\varphi_{a}(x_2)\varphi_b(x_1) \\ \end{align} Where it is not clear to my why the cross terms, \begin{align} \propto \int \textrm{d}x_1 \textrm{d}x_2 \left[ \varphi_{b}^{\dagger}(x_1) \varphi_b(x_2) \right] \left[ \varphi_{a}^{\dagger}(x_2) \varphi_{a}(x_1) \right] \end{align} should vanish. (Even if they did my main question still stands namely if you can always freely change the order of the wavefunctions). Note that the choice of permutation in coordinates or labels is equivalent to the choice whether to expand to rows or columns in a Slater determinant.

EDIT: upon request a more elaborate version. Suppose we want to write down a general wave function for 2 identical fermions $\Psi$. This wavefunction has to be antisymmetric. Suppose the single particle states are given by $$\varphi_{\alpha_i}(r_i,s_i) \,\,\,\,\,\ i \in \{1,2\}$$ Where $\alpha_i$ denotes the quantum numbers of a state, while $r_i,s_i$ denote the spatial and spin coordinates. A way to antisymmeterize $N$-body wavefunction in terms of single particle wavefunctions is often introduced by the concept of a Slater determinant. In the case of two particles a Slater determinant looks like \begin{align} \Psi_{\alpha_1,\alpha_2}(r_1,s_1,r_2,s_2) = \frac{1}{\sqrt{2}} \left| \begin{array}{c c} \varphi_{\alpha_1}(r_1,s_1) & \varphi_{\alpha_2}(r_1,s_1) \\ \varphi_{\alpha_1}(r_2,s_2) & \varphi_{\alpha_2}(r_2,s_2) \end{array} \right|. \end{align} Now we can calculate this determinant by expanding along a row or along a column. Which is equivalent to the choice whether to permute the quantum numbers of the single particle states or the the coordinates of the single particle states. More explicitly, expanding the matrix along the first row (or equivalently, permuting the label of the single particle states) we get, $$\frac{1}{\sqrt{2}} \left(\varphi_{\alpha_1}(r_1,s_1)\varphi_{\alpha_2}(r_2,s_2) - \varphi_{\alpha_2}(r_1,s_1)\varphi_{\alpha_1}(r_2,s_2) \right).$$ While an expansion along the first column gives (equivalent to a coordinate permutation), $$\frac{1}{\sqrt{2}} \left(\varphi_{\alpha_1}(r_1,s_1)\varphi_{\alpha_2}(r_2,s_2) - \varphi_{\alpha_1}(r_2,s_2)\varphi_{\alpha_2}(r_1,s_1) \right).$$ The two above expression certainly look to be equivalent. However if we consider the wave functions to be $N$-component vectors I get some strange results. Suppose we want to calculate the normalization of the two particle wave function. We would then write the following $$\sum_{s_1} \sum_{s_2} \int \textrm{d}r_1 \textrm{d}r_2 \Psi^{\dagger}(r_1,s_1,r_2,s_2) \Psi(r_1,s_1,r_2,s_2)$$ Note that I have introduced the $\dagger$ instead of the complex conjugate $*$ as the normalization should give a scalar. If we now replace $\Psi$ with the expression for the "label" permutation case we get, \begin{align*} \sum_{s_1} \sum_{s_2} & \int \textrm{d}r_1 \textrm{d}r_2 \Psi^{\dagger}(r_1,s_1,r_2,s_2) \Psi(r_1,s_1,r_2,s_2) \\ &= \frac{1}{2} \sum_{s_1} \sum_{s_2} \int \textrm{d}r_1 \textrm{d}r_2 \left(\varphi_{\alpha_1}(r_1,s_1)\varphi_{\alpha_2}(r_2,s_2) - \varphi_{\alpha_2}(r_1,s_1)\varphi_{\alpha_1}(r_2,s_2) \right)^{\dagger} \\ & \times \left(\varphi_{\alpha_1}(r_1,s_1)\varphi_{\alpha_2}(r_2,s_2) - \varphi_{\alpha_2}(r_1,s_1)\varphi_{\alpha_1}(r_2,s_2) \right) \\ &= \frac{1}{2} \sum_{s_1} \sum_{s_2} \int \textrm{d}r_1 \textrm{d}r_2 \left( \varphi^{\dagger}_{\alpha_2}(r_2,s_2) \varphi^{\dagger}_{\alpha_1}(r_1,s_1) \varphi_{\alpha_1}(r_1,s_1)\varphi_{\alpha_2}(r_2,s_2) \right. \\ & \,\, - \varphi^{\dagger}_{\alpha_2}(r_2,s_2) \varphi^{\dagger}_{\alpha_1}(r_1,s_1) \varphi_{\alpha_2}(r_1,s_1)\varphi_{\alpha_1}(r_2,s_2) \\ & \,\, - \varphi^{\dagger}_{\alpha_1}(r_2,s_2) \varphi^{\dagger}_{\alpha_2}(r_1,s_1) \varphi_{\alpha_1}(r_1,s_1)\varphi_{\alpha_2}(r_2,s_2) \\ & \,\, \left. + \varphi^{\dagger}_{\alpha_1}(r_2,s_2) \varphi^{\dagger}_{\alpha_2}(r_1,s_1) \varphi_{\alpha_2}(r_1,s_1)\varphi_{\alpha_1}(r_2,s_2) \right) \\ &= \frac{1}{2} \sum_{s_1} \int \textrm{d}r_1 \varphi^{\dagger}_{\alpha_1}(r_1,s_1) \varphi_{\alpha_1}(r_1,s_1) \sum_{s_2} \int \textrm{d}r_2 \varphi^{\dagger}_{\alpha_2}(r_2,s_2) \varphi_{\alpha_2}(r_2,s_2) \\ &-\frac{1}{2}\sum_{s_1} \int \textrm{d}r_1 \varphi^{\dagger}_{\alpha_1}(r_1,s_1) \varphi_{\alpha_2}(r_1,s_1) \sum_{s_2} \int \textrm{d}r_2 \varphi^{\dagger}_{\alpha_2}(r_2,s_2) \varphi_{\alpha_1}(r_2,s_2) \\ &-\frac{1}{2}\sum_{s_1} \int \textrm{d}r_1 \varphi^{\dagger}_{\alpha_2}(r_1,s_1) \varphi_{\alpha_1}(r_1,s_1) \sum_{s_2} \int \textrm{d}r_2 \varphi^{\dagger}_{\alpha_1}(r_2,s_2) \varphi_{\alpha_2}(r_2,s_2) \\ &+\frac{1}{2}\sum_{s_1} \int \textrm{d}r_1 \varphi^{\dagger}_{\alpha_2}(r_1,s_1) \varphi_{\alpha_2}(r_1,s_1) \sum_{s_2} \int \textrm{d}r_2 \varphi^{\dagger}_{\alpha_1}(r_2,s_2) \varphi_{\alpha_1}(r_2,s_2) \end{align*} Note that we have used the general expression $(AB)^{\dagger} = B^{\dagger} A^{\dagger}$ as well as the fact that the product $\varphi^{\dagger} \varphi$ contracts to a scalar and can be safely dragged trough the expressions. Now if one assumes the single particle states are orthogonal we have the following relation, \begin{align} \sum_{s_i}\int \textrm{d} r_i \varphi^{\dagger}_{\alpha_m}(r_i,s_i) \varphi_{\alpha_n}(r_i,s_i) = \delta_{m,n} \end{align} So we get, \begin{align} \sum_{s_1} \sum_{s_2} \int \textrm{d}r_1 \textrm{d}r_2 &\Psi^{\dagger}(r_1,s_1,r_2,s_2) \Psi(r_1,s_1,r_2,s_2) \\ & =\frac{1}{2} \delta_{\alpha_1,\alpha_1} \delta_{\alpha_2,\alpha_2} - \frac{1}{2} \delta_{\alpha_1,\alpha_2} \delta_{\alpha_2,\alpha_1} - \frac{1}{2} \delta_{\alpha_2,\alpha_1} \delta_{\alpha_1,\alpha_2} + \frac{1}{2} \delta_{\alpha_2,\alpha_2} \delta_{\alpha_1,\alpha_1} \\ &=1 \end{align} Which is to be expected. If we start from the second expression for $\Psi$ however we get the following, \begin{align*} \sum_{s_1} \sum_{s_2} & \int \textrm{d}r_1 \textrm{d}r_2 \Psi^{\dagger}(r_1,s_1,r_2,s_2) \Psi(r_1,s_1,r_2,s_2) \\ &= \frac{1}{2} \sum_{s_1} \sum_{s_2} \int \textrm{d}r_1 \textrm{d}r_2 \left(\varphi_{\alpha_1}(r_1,s_1)\varphi_{\alpha_2}(r_2,s_2) - \varphi_{\alpha_1}(r_2,s_2)\varphi_{\alpha_2}(r_1,s_1) \right)^{\dagger} \\ & \times \left(\varphi_{\alpha_1}(r_1,s_1)\varphi_{\alpha_2}(r_2,s_2) - \varphi_{\alpha_1}(r_2,s_2)\varphi_{\alpha_2}(r_1,s_1) \right) \\ &= \frac{1}{2} \sum_{s_1} \sum_{s_2} \int \textrm{d}r_1 \textrm{d}r_2 \left( \varphi^{\dagger}_{\alpha_2}(r_2,s_2) \varphi^{\dagger}_{\alpha_1}(r_1,s_1) \varphi_{\alpha_1}(r_1,s_1)\varphi_{\alpha_2}(r_2,s_2) \right. \\ & \,\, - \varphi^{\dagger}_{\alpha_2}(r_2,s_2) \varphi^{\dagger}_{\alpha_1}(r_1,s_1) \varphi_{\alpha_1}(r_2,s_2)\varphi_{\alpha_2}(r_1,s_1) \\ & \,\, - \varphi^{\dagger}_{\alpha_2}(r_1,s_1) \varphi^{\dagger}_{\alpha_1}(r_2,s_2) \varphi_{\alpha_1}(r_1,s_1)\varphi_{\alpha_2}(r_2,s_2) \\ & \,\, \left. + \varphi^{\dagger}_{\alpha_2}(r_1,s_1) \varphi^{\dagger}_{\alpha_1}(r_2,s_2) \varphi_{\alpha_1}(r_2,s_2)\varphi_{\alpha_2}(r_1,s_1) \right) \\ \end{align*} It is clear that the first and the fourth term will give $1/2 + 1/2 = 1$ as they are exactly the same as in the previous expression. However I do not see why the cross term should be necessarily zero. So it would seem that starting from the different antisymmeterization expression one gets different results. Hence my phrase "order matters" as the following appears to be true $$\frac{1}{\sqrt{2}} \left(\varphi_{\alpha_1}(r_1,s_1)\varphi_{\alpha_2}(r_2,s_2) - \varphi_{\alpha_2}(r_1,s_1)\varphi_{\alpha_1}(r_2,s_2) \right) \neq \frac{1}{\sqrt{2}} \left(\varphi_{\alpha_1}(r_1,s_1)\varphi_{\alpha_2}(r_2,s_2) - \varphi_{\alpha_1}(r_2,s_2)\varphi_{\alpha_2}(r_1,s_1) \right)$$ A final note about the Dirac notation. If one would start from the Dirac formalism and introduce $| \Psi \rangle$ as $$| \Psi \rangle = \frac{1}{\sqrt{2}} \left( | \alpha_1 \rangle | \alpha_2 \rangle - | \alpha_2 \rangle | \alpha_1 \rangle \right)$$ A projection into coordinate space gives the "label" permutation expression, $$\langle r_1,s_1 ; r_2,s_2 | \Psi \rangle = \frac{1}{\sqrt{2}} \left( \langle r_1,s_1| \alpha_1 \rangle \langle r_2,s_2 | \alpha_2\rangle - \langle r_1,s_1| \alpha_2 \rangle \langle r_2,s_2 | \alpha_1\rangle \right) \\ = \frac{1}{\sqrt{2}} \left(\varphi_{\alpha_1}(r_1,s_1)\varphi_{\alpha_2}(r_2,s_2) - \varphi_{\alpha_2}(r_1,s_1)\varphi_{\alpha_1}(r_2,s_2) \right)$$ The expression for the normalization starting from the Dirac formalism gives, $$\langle \Psi| \Psi \rangle = \frac{1}{2} \left( \langle \alpha_1 | \alpha_1 \rangle \langle \alpha_2 | \alpha_2 \rangle - \langle \alpha_1 | \alpha_2 \rangle \langle \alpha_2 | \alpha_1 \rangle - \langle \alpha_2 | \alpha_1 \rangle \langle \alpha_1 | \alpha_2 \rangle + \langle \alpha_2 | \alpha_2 \rangle \langle \alpha_1 | \alpha_1 \rangle \right)$$ Using the unit identities: $1 = \sum_s | s \rangle \langle s |$ and, $1 = \int \textrm{d} r | r \rangle \langle r |$ we get, $$\langle \Psi| \Psi \rangle = \frac{1}{2} \sum_{s_1} \sum_{s_2} \int \textrm{d} r_1 \int \textrm{d} r_2 \\ \left( \langle \alpha_1 | r_1,s_1 \rangle \langle r_1,s_1 | \alpha_1 \rangle \langle \alpha_2 | r_2,s_2 \rangle \langle r_2,s_2| \alpha_2 \rangle - \langle \alpha_1 | r_1,s_1 \rangle \langle r_1,s_1 | \alpha_2 \rangle \langle \alpha_2 | r_2,s_2 \rangle \langle r_2,s_2| \alpha_1 \rangle - \langle \alpha_2 | r_1,s_1 \rangle \langle r_1,s_1| \alpha_1 \rangle \langle \alpha_1 | r_2,s_2 \rangle \langle r_2,s_2| \alpha_2 \rangle + \langle \alpha_2 | r_1,s_1 \rangle \langle r_1,s_1| \alpha_2 \rangle \langle \alpha_1 | r_2,s_2 \rangle \langle r_2,s_2| \alpha_1 \rangle \right)$$ Which is when using $\langle r_i, s_i | \alpha_i \rangle = \varphi_{\alpha_i}(r_i,s_i)$ exactly the same as the normalization expression I get for the label permutation case, which works out nicely to $1$ as expected.

• There might be some issue when exchanging Grassmanian coordinates ... Commented Jun 20, 2013 at 22:32
• What do you mean exactly with "some issue"? What would be the correct way to antisymmetrise the wavefunction? Looking at my simple example above I would say the permutation in labels seems to be more correct. But why isn't this equivalent to a coordinate permutation? Commented Jun 21, 2013 at 7:38
• "For example relativistic wavefunctions which are 4x1 matrices." .No. First, a 4x1 matrix is a vector... Second, a wafefunction is not a vector, and a wafefunction is not a matrix. Wavefunctions are complex quantities, so your 2 expressions are obviously equal. Commented Jun 21, 2013 at 8:47
• @camelthemammel I didn't mean to write everything in the question. I thought that writing down all the details by yourself will make everything clear for you. You are already considering "N-component vectors" by writing down spin indices. $\langle r_i, s_i | \alpha_i \rangle = \varphi_{\alpha_i}(r_i,s_i)$ -- is a number. It commutes. Commented Jul 5, 2013 at 9:27
• @camelthemammel Honestly, I'm starting to suspect that understanding is not your actual goal. Let me guess -- you are arguing with someone and don't want to end up loosing. That is why you are making more and more complex conglomerations in support of your claim. Commented Jul 5, 2013 at 9:32

I suspect that the root of this (quite standard) confusion is that one thinks that wavefunctions are the basic notion of quantum mechanics. While you have to think about Hilbert spaces, vectors, scalar and tensor products.

If, say, you'd use Dirac notation : $$\mid\!\Psi\rangle = \frac{1}{\sqrt2}\left(\mid\!a\rangle\otimes\mid\!b\rangle\,-\mid\!b\rangle\otimes\mid\!a\rangle\right)$$ Then there is no "label/argument ambiguity" in the first place.

So your question doesn't make much sense if you think about basics of QM.

Finally, talking about normalization (I hope that you'll vaguely recognize yours in it): $$\langle\Psi\mid\!\Psi\rangle=\frac12\left(\langle a\mid\!a\rangle\langle b\mid\!b\rangle-\langle a\mid\!b\rangle\langle b\mid\!a\rangle-\langle b\mid\!a\rangle\langle a\mid\!b\rangle+\langle b\mid\!b\rangle\langle a\mid\!a\rangle\right)$$ Depends on mutual properties of $\mid\!a\rangle$ and $\mid\!b\rangle$. Usually they form a complete set (like in your case, since you called them "eigenstates"): $$\langle i\mid\!j\rangle = \delta_{ij}$$ Which leads to a unit normalization.

• Ok, I agree. But suppose we project your first equation into coordinate space, we then get something like $\langle x_1, x_2 | \Psi \rangle = 1/\sqrt{2} ( \langle x_1 | a \rangle \langle x_2 | b \rangle - \langle x_1 | b \rangle \langle x_2 | a \rangle)$. Which seems very familiar to my "label" permutation wave function. Hence is it ok to state that the "label" permutation is the (more) correct one? Commented Jun 21, 2013 at 13:22
• @camelthemammel , nice! It is not just "familiar" -- it is exactly what you are writing: $\phi_a(x_1) = \langle x_1\mid a\rangle$. You can get your expressions by inserting unit operator $1 = \int\,dx\, \mid x\rangle\langle x\mid$ here and there. Try it! Commented Jun 21, 2013 at 13:55
• Yes it is quite obvious that the expression are the same. Just to give my brain some closure; would you agree then that the "label permutation" is the way to go when antisymmetrising a wave function and the coordinate permutation equivalence only follows if the wave functions are scalars? Commented Jun 21, 2013 at 15:26
• Could you comment on my last remark as it is my main question ("equivalence" of label/coordinate permutation)? I'm not accepting your response as the answer yet because of this. Commented Jun 25, 2013 at 9:53
• @camelthemammel No, do not agree. You can swap coordinates as long as you swap spin indices too. $\frac{1}{\sqrt{2}} \left(\varphi_{a}(x_1,J_1)\varphi_b(x_2,J_2) - \varphi_{a}(x_2,J_2)\varphi_b(x_1,J_1) \right)$ You could have gotten that by yourself if you would have done the unnit operation insertion I've suggested. Commented Jun 25, 2013 at 16:53

You had the following cross term for the case of coordinate permutation:

\begin{align} \propto \int \textrm{d}x_1 \textrm{d}x_2 \left[ \varphi_{b}^{\dagger}(x_1) \varphi_b(x_2) \right] \left[ \varphi_{a}^{\dagger}(x_2) \varphi_{a}(x_1) \right]. \end{align}

However, this is not how you take the inner product of $\varphi_{a}(x_{1})\varphi_{b}(x_{2})$ and $\varphi_{a}(x_{2})\varphi_{b}(x_{1})$.

$\varphi_{a}(x_{1})$ and $\varphi_{b}(x_{1})$ belong to the Hilbert space of particle 1, which I denote $\mathcal{H}_{1}$, and similarly, $\varphi_{a}(x_{2}), \varphi_{b}(x_{2}) \in \mathcal{H}_{2}$. The objects $\varphi_{a}(x_1) \varphi_{b}(x_2)$ and $\varphi_{a}(x_2) \varphi_{b}(x_1)$ are members of $\mathcal{H}_{1}\otimes\mathcal{H}_{2}$, which is the tensor product of $\mathcal{H}_{1}$ and $\mathcal{H}_{2}$. Furthermore, as $\varphi_{a}(x_{1})$ and $\varphi_{b}(x_{2})$ belong to different Hilbert spaces, $\varphi_{a}(x_1) \varphi_{b} (x_2)$ and $\varphi_{b} (x_2)\varphi_{a}(x_1)$ really mean the same thing: there is no ordering issue even if $\varphi$ is a $N$-component vector.

When taking the inner product of the two objects $\varphi_{a}(x_{1})\varphi_{b}(x_{2})$ and $\varphi_{a}(x_{2})\varphi_{b}(x_{1})$ in $\mathcal{H}_{1}\otimes\mathcal{H}_{2}$, you should contract the parts in $\mathcal{H}_{1}$ together, and the parts in $\mathcal{H}_{2}$ together. Then, $$\int dx_{1} dx_{2} \big[\varphi_{a}(x_{1})\varphi_{b}(x_{2})\big]^{\dagger} \varphi_{a}(x_{2})\varphi_{b}(x_{1}) = \left[\int dx_{1} \varphi_{a}(x_{1})^{\dagger} \varphi_{b}(x_{1})\right]\left[\int dx_{2} \varphi_{b}(x_{2})^{\dagger} \varphi_{a}(x_{2})\right] = 0.$$