To compute the vdW polarization density, MBD model as implemented in the libMBD library was employed. The MBD method maps atoms onto a model Hamiltonian of charged quantum harmonic (Drude) oscillators (QDO) centered at R A , with frequen-cies ω A , charges q A and masses m A . The MBD Hamiltonian for a system of N coupled oscillators reads