def yield_pair_energies(self, index1, index2):
strength = self.strengths[index1, index2]<EOL>distance = self.distances[index1, index2]<EOL>yield strength*distance**(-<NUM_LIT:6>), <NUM_LIT:1><EOL>
Yields pairs ((s(r_ij), v(bar{r}_ij))