= Solution
The <canonical anticommutation relations> realize the internal states as a <fermionic Fock space> $\Lambda^\bullet\mathbb C^3$. Choosing a <Clifford vacuum> $|0\rangle$ annihilated by all $\psi_a$, the states are functions multiplying $|0\rangle$, $\psi_a^\dagger|0\rangle$, $\psi_a^\dagger\psi_b^\dagger|0\rangle$ with $a<b$, or $\psi_1^\dagger\psi_2^\dagger\psi_3^\dagger|0\rangle$. Thus the four <fermion-number sectors in supersymmetric quantum mechanics> are
$$
\boxed{\mathcal H_n=L^2(\mathbb R^3)\otimes\Lambda^n\mathbb C^3,\quad n=0,1,2,3,}
$$
with internal dimensions $1,3,3,1$. The <fermion number operator> $N=\sum_a\psi_a^\dagger\psi_a$ has value $n$ on $\mathcal H_n$. The sectors $n=0,2$ have even <fermion parity>, and $n=1,3$ have odd <fermion parity>.
Write $q=\sum_a\psi_a^\dagger A_a$, so that the <supercharge> is $Q=q+q^\dagger$. The <commutator> $[A_a,A_b]$ is zero because the two <mixed partial derivatives> of $\chi$ agree. The <canonical anticommutation relations> then imply $q^2=(q^\dagger)^2=0$ and
$$
H=qq^\dagger+q^\dagger q,\qquad [H,N]=[H,q]=[H,q^\dagger]=[H,Q]=0.
$$
The <Hamiltonian operator> preserves each sector, while $q$ raises and $q^\dagger$ lowers its degree by one. For a state in its operator domain, $\langle v,Hv\rangle=\|qv\|^2+\|q^\dagger v\|^2\geq0$, so the <spectrum> is nonnegative. If $Hv=Ev$ with $E>0$, then $Qv/\sqrt E$ is a state of the same energy and opposite <fermion parity>, and $Q(Qv/\sqrt E)=\sqrt E v$. This is an <odd involution pairing bosonic and fermionic states> and proves equal even and odd multiplicities at every positive <eigenvalue>.
There is also a more precise adjacent-sector pairing. On the positive-energy <eigenspace>,
$$
v=E^{-1}qq^\dagger v+E^{-1}q^\dagger qv.
$$
The two summands are <orthogonal> since $q^2=(q^\dagger)^2=0$. A state $u$ in the image of $q^\dagger$ obeys $q^\dagger u=0$, $q^\dagger qu=Eu$, so $q/\sqrt E$ maps it to a normalized partner in the next sector, with inverse $q^\dagger/\sqrt E$. Thus positive-energy doublets lie in $\mathcal H_0\oplus\mathcal H_1$, $\mathcal H_1\oplus\mathcal H_2$, or $\mathcal H_2\oplus\mathcal H_3$. In particular, every positive-energy state of $\mathcal H_0$ has a partner in $\mathcal H_1$, and every one of $\mathcal H_3$ has a partner in $\mathcal H_2$; the remaining middle-sector states pair with each other. There is no requirement that all four individual sector spectra coincide. Zero-energy states satisfy $qv=q^\dagger v=Qv=0$ and can be unpaired. No confinement assumption was given: a continuous <spectrum> is possible, with the same even/odd pairing on its positive spectral subspaces.
For real $\chi$, the <adjoint operator> is $A_a^\dagger=-\partial_a-\partial_a\chi$. In $\mathcal H_0$, every $\psi_a^\dagger\psi_b$ vanishes; in the one-dimensional internal space of $\mathcal H_3$, it equals $\delta_{ab}$. The two scalar <Hamiltonian operators> are consequently
$$
\boxed{H_0=\sum_a A_a^\dagger A_a=-\Delta+|\nabla\chi|^2+\Delta\chi,\qquad H_3=\sum_a A_aA_a^\dagger=-\Delta+|\nabla\chi|^2-\Delta\chi.}
$$
These are the extreme sectors of the <gradient superpotential Hamiltonian in supersymmetric quantum mechanics>.
For radial $\chi$ and a normalized state $v_0=f(r)|0\rangle$ with $E>0$, $q^\dagger v_0=0$ and $A_af=(x_a/r)(f'-\chi'f)$. The <radial superpartner in three-dimensional supersymmetric quantum mechanics> therefore has the simple form
$$
\boxed{v_1=\frac{f'(r)-\chi'(r)f(r)}{\sqrt E}\,\psi_r^\dagger|0\rangle,\qquad \psi_r^\dagger=\sum_a\frac{x_a}{r}\psi_a^\dagger.}
$$
It belongs to $\mathcal H_1$; $Hv_1=Ev_1$ follows from $[H,Q]=0$, and its <wavefunction normalization> follows from $\|Qv_0\|^2=\langle v_0,Hv_0\rangle=E$. Also $Qv_1=\sqrt E v_0$. At the origin the formula is interpreted by regular extension; its degree-one components carry the radial <unit vector> rather than being three independent radial scalars.
For a <zero mode in the fermion-vacuum sector>, the factorization of $H_0$ gives $0=\sum_a\|A_af\|^2$. Hence $\partial_a f=(\partial_a\chi)f$, or $\partial_a(e^{-\chi}f)=0$. The requested <wavefunction> is
$$
\boxed{f(r)=Ce^{\chi(r)},\qquad |C|^{-2}=4\pi\int_0^\infty r^2e^{2\chi(r)}\,dr.}
$$
A nonzero <normalizable wavefunction> exists only if this integral is finite. The positive sign in the exponential is fixed by the given convention $A_a=\partial_a-\partial_a\chi$.
Back to article page
