Hamiltonian is quadratic in the momentum operator, and the quadrupole moment is quadratic in position operator. From Heisenberg equations, it follows that the third derivative of the quadrupole moment vanishes (it was obtained by direct calculations for the LG packet in [1]}), thus, one can see that

Q,int(t) =Q,int(0)+Q(2)t2
