Тепловое дипольное взаимодействие осцилляторов

Найдите среднюю потенциальную энергию дальнодействующего электростатического взаимодействия двух одинаковых изотропных трёхмерных гармонических осцилляторов при произвольной температуре $T$. Дипольный момент каждого равен $\mathbf d=q\mathbf s$, где $\mathbf s$ — отклонение от положения равновесия.