The computation of the third term in Equation
(
106
), which corresponds to the interaction between two quadrupoles,
, can be found in Ref.
[
21
]
.