Dear Lammps Users,
I would like to couple Gayberne potential with DPD/tstat as a thermostat in my system, which helps preserve the hydrodynamic behavior of the system. However, It didn’t work. According to my understanding, since the first term in DPD/tstat is from a pairwise interaction,which is only dependent on distance. This being said, if I simply use LJ potential, dpd/tstat will definitely work; while if I am using Gayberne potential, with one more dependence arising from the orientation, which causes the DPD/tstat to fail in this system.
I am wondering is there any possibility at all that if I would like to preserve hydrodynamic interaction while I can still use Gayberne Potential?
Thanks very much in advance for any thoughts and ideas!
Sincerely,
LQ