Suppose I have 2 fermions in a potential $V(x)$. Both particles are moving in one dimension: the $x$ axis. Then, neglecting the interaction between the particles, the spatial wave function of the system would be of the form $$\psi_{n_{1}}(x_{1})\psi_{n_{2}}(x_{2}) $$
Now, if I'm considering particles with spin 1/2, the notation $\alpha(1)$ indicates that the particle 1 has spin up, and $\beta(2)$ denotes the particle 2 having spin down.
Now, I want to write the complete wave function, a function of the form $$\psi_{n_{1}n_{2}s_{1}s_{2}}(x_{1},x_{2},s_{1},s_{2})=\psi_{n_{1}}(x_{1})\psi_{n_{2}}(x_{2})F(\alpha,\beta)$$ where $F(\alpha,\beta)$ is a function of the spin of the system.
To this end, I have that the only physically possible functions $F(\alpha,\beta)$ are:
Symmetric: $\chi_{\alpha}:=\alpha(1)\alpha(2),\quad\chi_{\beta}:=\beta(1)\beta(2),\quad\chi_{+}:=\frac{1}{\sqrt{2}}\left[\alpha(1)\beta(2)+\alpha(2)\beta(1)\right]$
Antisymmetric: $\chi_{-}:=\frac{1}{\sqrt{2}}\left[\alpha(1)\beta(2)-\alpha(2)\beta(1)\right]$
In order to write down the complete wave function with spin, I understand I have to consider the energy levels. For example the ground state: $\psi_{1}(x_{1})\psi_{1}(x_{2})$.
If $\psi_{1}(x_{1})\psi_{1}(x_{2})$ is symmetric (as I understand it is), then I must multiply this function times the antisymmetric function $\chi_{-}$ (in order to get an antisymmetric wave function, for two fermions).
If $\psi_{1}(x_{1})\psi_{1}(x_{2})$ is antisymmetric (and I understand this is impossible, since the ground state is not degenerate), then I'd have 3 wave functions, obtained by multiplying $\psi_{1}(x_{1})\psi_{1}(x_{2})$ times $\chi_{\alpha}$, $\chi_{\beta}$ and $\chi_{+}$.
Now, for the 1st excited level, say $\psi_{2}(x_{1})\psi_{1}(x_{2})$, my question is, what happens when this function is not symmetric neither antisymmetric?
I mean, I could built a symmetric
$$f_S=\frac{1}{\sqrt{2}}\left[\psi_{2}(x_{1})\psi_{1}(x_{2})+\psi_{2}(x_{2})\psi_{1}(x_{1})\right]$$
or an antisymmetric
$$f_A=\frac{1}{\sqrt{2}}\left[\psi_{2}(x_{1})\psi_{1}(x_{2})-\psi_{2}(x_{2})\psi_{1}(x_{1})\right]$$ wave function. But which one of these must I choose? Or, must I calculate the resultant complete wave functions with both? Then, when I count the states with the 1st excited energy, I'd have 4 instead of 1 or 3.