Hello everyone,

I want to use different functional forms for bond stretching and bond angles in my forcefield document for a DPD simulation. There is a paper (doi: 10.3389/fmolb.2021.619381) which replaced harmonic potential with leonard jones potential for bond stretching in a coarse-grained MD run. There is also another paper (DOI https://doi.org/10.1039/D1SM01194D) which modelled the bond angle differently for a DPD run (equation 6c in paper). I was thinking that I can use tabulated values of energy and bond distances/bond angles using these new expressions in the forcefield document. Is there any simpler method to implement it? 

Thanks,

Ricky