Solution (source code)

= Solution

\b[Distinguishable noninteracting particles.] Interpret $V$ as the same external one-particle potential for each particle; the total potential is $V(\mathbf r_1)+V(\mathbf r_2)$. The <Time-dependent Schrödinger equation> is
$$
\boxed{i\hbar\frac{\partial\Psi}{\partial t}=
\left[-\frac{\hbar^2}{2m}(\nabla_1^2+\nabla_2^2)
+V(\mathbf r_1)+V(\mathbf r_2)\right]\Psi.}
$$
Let $\widehat h\phi_a=\epsilon_a\phi_a$ and $\widehat h\phi_b=\epsilon_b\phi_b$, where $\widehat h=-\hbar^2\nabla^2/(2m)+V$. Applying $\widehat h_1+\widehat h_2$ to the product multiplies it by $\epsilon_a+\epsilon_b$. The time <derivative> of $e^{-iEt/\hbar}$ multiplies it by $E$. Thus the proposed stationary product is a solution provided
$$
\boxed{E=\epsilon_a+\epsilon_b.}
$$
The <tensor product> separates the two independent one-particle dynamics.

\b[Identical spin-one-half particles.] Identical <fermions> require antisymmetry under exchange of the complete position-and-spin labels. First take $a,b$ to be distinct orthonormal spatial orbitals, and define
$$
\Phi_\pm(1,2)=\frac{\phi_a(\mathbf r_1)\phi_b(\mathbf r_2)
\pm\phi_b(\mathbf r_1)\phi_a(\mathbf r_2)}{\sqrt2}.
$$
The symmetric spatial state must multiply the antisymmetric <spin-one-half singlet state>; the antisymmetric spatial state must multiply one of the three symmetric <spin-one-half triplet states>. The four normalized possibilities, each with the common time factor $e^{-iEt/\hbar}$, are
$$
\begin{array}{c|c|c}
\text{spatial factor}&\text{spin factor}&(S,m_S)\\\hline
\Phi_+&\dfrac{|\uparrow\downarrow\rangle-|\downarrow\uparrow\rangle}{\sqrt2}&(0,0)\\[4pt]
\Phi_-&|\uparrow\uparrow\rangle&(1,1)\\
\Phi_-&\dfrac{|\uparrow\downarrow\rangle+|\downarrow\uparrow\rangle}{\sqrt2}&(1,0)\\[4pt]
\Phi_-&|\downarrow\downarrow\rangle&(1,-1)
\end{array}
$$
These are the required product states under <fermionic exchange symmetry>. If $a=b$, the antisymmetric spatial combination vanishes. Only the product $\phi_a(\mathbf r_1)\phi_a(\mathbf r_2)$ times the singlet remains, as required by the <Pauli exclusion principle>. A nonorthogonal pair would require different normalization factors; the displayed $1/\sqrt2$ assumes an orthonormal stationary-state basis.

\b[Weak repulsive <Coulomb interaction>.] For fixed orbitals $a,b$, first-order <perturbation theory> gives spatial <energy> shifts $J\pm K$. Here $J$ is the direct interaction <energy>, while the exchange integral for repulsion strength $g>0$ is
$$
K=g\iint\frac{\phi_a^*(\mathbf r)\phi_b(\mathbf r)
\phi_b^*(\mathbf r')\phi_a(\mathbf r')}{|\mathbf r-\mathbf r'|}\,d^3r\,d^3r'\ge0.
$$
The positivity follows from the positive <Fourier transform> of the Coulomb kernel. Consequently <Coulomb exchange splitting of two orbitals> gives
$$
\boxed{E_{\rm triplet}=\epsilon_a+\epsilon_b+J-K
\le E_{\rm singlet}=\epsilon_a+\epsilon_b+J+K.}
$$
All three triplet <spin> components remain degenerate for a spin-independent <Coulomb interaction>. Their antisymmetric spatial factor vanishes at coincidence and reduces close encounters. The inequality is strict when $K>0$; for negligible orbital overlap, the splitting may vanish. This comparison fixes one particle in each of the specified orbitals and does not rank singlets with a different orbital occupation.

\b[Spin-zero bosons.] For one particle in each of the same distinct orthonormal orbitals, the <spin> factor is symmetric and <bosonic exchange symmetry> requires $\Phi_+$. At coincident positions,
$$
|\Phi_+(\mathbf r,\mathbf r)|^2=2|\phi_a(\mathbf r)|^2|\phi_b(\mathbf r)|^2,
\qquad
|\Psi_D(\mathbf r,\mathbf r)|^2=|\phi_a(\mathbf r)|^2|\phi_b(\mathbf r)|^2.
$$
Thus the <bosonic coincidence enhancement in two orthogonal orbitals> is
$$
\boxed{\frac{P_B(\mathbf r,\mathbf r)}{P_D(\mathbf r,\mathbf r)}=2.}
$$
These are joint <probability densities>, or small-detector-volume limits; an exactly specified point has zero <probability> in a continuous distribution. If both particles occupy the same orbital, the properly normalized product is already symmetric and the corresponding enhancement relative to distinguishable particles in that same orbital is one. The paper's usual factor-two answer presumes distinct orthogonal input orbitals.