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