When I do langevin dynamic simulation with fix langevin command, I found the large deviation of temperature from target temperature, about 300/± 10. My system include a fixed rodlike polyelectrolyte , two mobile semiflexible polyelectrolytes and 600 ions.
I use fix nve and langevin command.
ANY suggestion? Thanks in advance.