where the exact average is taken by considering the original
Hamiltonian(27). To evaluate the above expression, we have
approximated the average by applying the diagonalized harmonic
Hamiltonian. Once ϕi and Si
are uncoupled operators and ϕi has a
Gaussian distribution, we have