Harmonic Oscillator Wave Functions in Nuclear Spectroscopy

Abstract
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 σ1, σ2 and r1-r2. Since the coordinates appear only in the form r1-r2, it is convenient to write the two-body wave function in terms of the relative coordinate, r=r1-r2, and the center-of-mass coordinate, 2R=r1+r2. The eigenfunction in the new coordinates can be determined by noting that if the two particles are in the same oscillator level, then ψ(r1, r2), which is an eigenfunction of H1+H2, is also an eigenfunction of H1 alone. Transforming H1 to relative and center-of-mass coordinates implies that the operator p·P+mkr·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.