This leads to the separable Hamiltonian, with the relative motion part given by (rij=r)
ˆHr=−∇2r+14ω2r2+1r,plus a standard Harmonic Oscillator problem for the center-of-mass motion. This system has analytical solutions in two and three dimensions (M. Taut 1993 and 1994).