Hello. We are trying to model polymer chains in which there are both bonded and non-bonded potentials. We want the bonded atoms to interact via the bond potential and the non-bonded atoms to interact via LJ. We have not been able to figure out an efficient way to do this. Any suggestions?