GROMACS version: 2020.4
GROMACS modification: No
I have a small box of chlorobenzene molecules and I’m using OPLS-AA FF and I got the parameters and charges from the LigParGen server. When I run the production MD in NPT, if I use Parrinello-Rahman barostat and Nose-Hoover thermostat the system will blow up. Conversely, If I use the Berendsen barostat and V-rescale thermostat everything is fine.
I tried more restrict energy minimization and equilibration on NVT and NPT with more time (200ps) and also a prior NPT simulation (200ps) without constraint in Berendsen barostat followed by a production NPT in P-R and it always blows up in few steps.
If there is something wrong with my simulation why is it all good in Berendsen? It gives reasonable density and dielectric constant too.
I have attached the mdp and topol files and em.gro (as em.dat) in case you want to take a look at them.