Pairing energy calculations are generally carried out assuming the nucleon wave functions are those of an harmonic oscillator. The two-body interaction is assumed to be some function of ${\ensuremath{\sigma}}_{1}$, ${\ensuremath{\sigma}}_{2}$ and ${\mathbf{r}}_{1}$-${\mathbf{r}}_{2}$. Since the coordinates appear only in the form ${\mathbf{r}}_{1}$-${\mathbf{r}}_{2}$, it is convenient to write the two-body wave function in terms of the relative coordinate, r=${\mathbf{r}}_{1}$-${\mathbf{r}}_{2}$, and the center-of-mass coordinate, 2R=${\mathbf{r}}_{1}$+${\mathbf{r}}_{2}$. The eigenfunction in the new coordinates can be determined by noting that if the two particles are in the same oscillator level, then $\ensuremath{\psi}({\mathrm{r}}_{1}, {\mathrm{r}}_{2})$, which is an eigenfunction of ${H}_{1}+{H}_{2}$, is also an eigenfunction of ${H}_{1}$ alone. Transforming ${H}_{1}$ to relative and center-of-mass coordinates implies that the operator $\mathrm{p}\ifmmode\cdot\else\textperiodcentered\fi{}\mathrm{P}+mk\mathrm{r}\ifmmode\cdot\else\textperiodcentered\fi{}\mathrm{R}$ (where $m$ is the mass of the particle and $k$ is the spring constant of the oscillator) must give zero when operating on the wave function. This condition plus certain requirements arising from the radial form of the oscillator eigenfunctions is sufficient to determine the wave function in the new coordinate system.
Read full abstract