Added the calculation of electrostatic interaction between subsystems.