Abstract
Charge transfer is one of the mechanisms involved in non-covalent interactions. In molecular dimers, its contribution to the pairwise interaction energies has been studied extensively using a variety of interaction energy decomposition schemes. In polar interactions such as hydrogen bonds, it can contribute with ten or few tens of percents of the interaction energy. Less is known about its importance in higher-order interactions in many-body systems, mainly because of the lack of methods applicable to this problem. In this work, we extend our method for quantifying the charge transfer energy based on constrained DFT to many-body cases and apply it to model trimers extracted from molecular crystals. Our calculations show that charge transfer can account for a large fraction of the total three-body interaction energy. This also has implications for DFT calculations of many-body interactions in general as it is known that many DFT functionals struggle to correctly describe charge transfer effects.