We propose, in the following, an alternative approach bypassing the diagonalization step. We first express E int in terms of the trace of √ Λ: by recognizing the trace invariance under any change of basis, the interaction energy is written as shown in Eq.( ) which, at this point is still exact. For simplicity, the MBD lower script is dropped from now on.