[lammps-users] custom improper

I am attempting to use an united atom force field with a different form for the improper force. It uses the form


I have modified improper_harmonic.cpp but I am unsure if I did this correctly. I have attached the modified cpp file. Thanks

Ben Hanson

improper_PMMA.cpp (7.59 KB)