More generally, it should be noted that kQEq is not primarily intended as a stand-alone molecular dipole model. Since it is based on an energy functional, it can be used to model long-range electrostatic interactions in combination with local interatomic potentials. Indeed, it remains an open question whether dipole moments alone provide enough information for this purpose. Fortunately, the current approach can easily be expanded to also include higher moments, electrostatic potentials or reference partial charges. This will be explored in future work.