I expect that in many-body problems of electrons, spin should cause same-spin-electrons to repel more strongly than opposite spin electrons because the Pauli exclusion principle is the observation that a slater determinant with repeated columns vanish identically. HF often capture the physics of many-body problems, and therefore we expect that there should be some repulsion due to spin in theories that HF works well. I am trying to pinpoint how this repulsion (1) enters mathematically in the exact solution (as it does not appear in the Hamiltonian) and (2) how the repulsion is expressed in $|\psi|^2$. The hamiltonian of an atom is an example that I have been using to understand this. The simplest is helium:
$$ H = -\frac{1}{2}\nabla^2_1 -\frac{1}{2}\nabla^2_2 - \frac{1}{r_1} - \frac{1}{r_2} + \frac{1}{r_{12}}. $$ As it is a spin-independent theory, any solution of the time-independent Schrödinger equation cannot depend on spin through the Schrödinger equation. I expect that any dependence of spin must enter by additional restrictions of $\psi$. Wikipedia claims in the third paragraph of this section that if $i,j = \uparrow,\downarrow$, a solution is a matrix of functions $\psi_{ij} : \mathbb{R}^3 \times \mathbb{R}^3 \to \mathbb{C}$.
This raises a couple of questions:
Does Wikipedia demand anything more than that the spatial part be antisymmetric? I suppose so, otherwise there would be no reason to label the different spin configurations individually. Wikipedia claims that the $\psi_{ij}$ are restricted by the pauli exclusion principle. How could that possibly provide a condition $\psi_{ij}$ here?
How can I see that these restrictions show up as a repulsion in $|\psi_{ii}|^2$?
It would be helpful if you could be as mathematically precise as possible.