I have a system containing two aromatic molecules. One of them contains a functional group with a lone pair. I want to calculate the interaction energy in details. I used PSI4 to perform SAPT2+ calculation and found the interaction energy classified into electrostatic, induction, exchange and dispersion.
Now I would like to calculate how much of the interaction energy is originated from dipole-pi interaction and how much is from pi-pi stacking. Is there any way to calculate this information?