  A flux qubit can interact strongly when it is capacitively coupled to other
circuit elements. This interaction can be separated in two parts, one acting on
the qubit subspaces and one in which excited states mediate the interaction.
The first term dominates the interaction between the flux qubit and an
LC-resonator, leading to ultrastrong couplings of the form
$\sigma^y(a+a^\dagger),$ which complement the inductive
$\sigma^xi(a^\dagger-a)$ coupling. However, when coupling two flux qubits
capacitively, all terms need to be taken into account, leading to complex
non-stoquastic ultrastrong interaction of the $\sigma^y\sigma^y$,
$\sigma^z\sigma^z$ and $\sigma^x\sigma^x$ type. Our theory explains all these
interactions, describing them in terms of general circuit properties---coupling
capacitances, qubit gaps, inductive, Josephson and capactive energies---, that
apply to a wide variety of circuits and flux qubit designs.
