schroedinger equation - Coupled many body quantum harmonic oscillator in 3 Dimension

Can coupled many body quantum harmonic oscillator be solved exactly ? $$H=\sum_{i}\dfrac{\vec{p}_{i}^{2}}{2m}+\dfrac{1}{2}\sum_{k,j}K(\vec{R_{k}}-\vec{R}_{j})^{2}$$. Is there any reference to solve this exactly.