Hi all,
I want to simulate the viscosity evolutions of water with the increase of surfactant concentrations using Lammps. Since I need to keep simulating box orthogonal, I adopt Muller - Plathe method and “fix viscosity” in Lammps. In the manual, this fix should not use a thermostat, so I chose “NVE” after a 5 ns NPT equilibration, but temperature are still fluctuating dur this process, and finally increase 10 degC.
(1) Is it necessary to control temperature during the swaping process?
(2) I want to know the viscosity changes of water, so I only need to swap the momentum of water, rather than the whole system, right?
Regards,
Jun