TODO:
1. create profile of each donor (#right)
2. right function to compute donation of all neighbors (array of ids to
atoms)
d. calculate number of neighbors