Dear Gigo: that's a good comprehensive testing and report. Please let us know what you find out from those authors. Their paper was short on methods (unless I missed it... I didn't check for any SI), so perhaps they did something non-standard and didn't report it.
I think at this point it is a good idea for you to file a redmine issue. It's not a gromacs error per se, but if this is true then pdb2gmx or grompp should give a warning or error for the combination of oplsaa and tip5p. Chris. gigo gigo at ibb.waw.pl Sun Sep 29 00:59:42 CEST 2013 Previous message: [gmx-users] OPLS/AA + TIP5P, anybody? Next message: [gmx-users] Restarting simulation -s. files? Messages sorted by: [ date ] [ thread ] [ subject ] [ author ] Dear Chris, I am really grateful for your help. This is what I did, with additional LJ terms on LP1 and LP2 of tip5p: - 5000 steps of steepest descent with positions restraints on protein and flexible water (flexibility like in tip4p), - 5000 steps of steep, no restraints, flexible water, - 5000 steps of cg, flexible water, - 100000 steps of MD with posres and constrained bond lengths, very weak temp coupling (v-rescale), starting without generating of initial velocities, so its heating very slowly to 300K, no pressure coupling - 100000 steps of MD with posres, no constraints, v-rescale, nvt - 100000 steps of MD, no posres, nvt, - 100000 steps of MD, v-rescale, Berendsen pressure coupling, - 100000 steps of MD, v-rescale, Parrinello-Rahman pressure coupling. Output of this chain, after removing LJ from LP, became input for 4 simulations just like the last one from the above chain, with or without posres and constraints turned on. Results: 1) Posres off, constraints off : "step 0" "WARNING: Listed nonbonded interaction between particles 1157 and 1163" <== Why does it say "particles" and not "atoms"? Nevermind, its lysine on the surface of the protein, one of atoms is dimensionless hydrogen. (...) "step 63: Water molecule starting at atom 7923 can not be settled." (...) "Segmentation fault" 2) Posres off, constraints on : Warning like above + LINCS warnings, all for charged aminoacids on surface of the protein Segfault at step 63 3) Posres on, constraints off, had to add refcoord_scaling = COM : "WARNING: Listed nonbonded interaction between particles 3075 and 3079" <== arginine on the surface, dimensionless hydrogen 4) Posres on, constraints on, refcoord_scaling = COM : Same warnings, several other positively charged dimensionless hydrogens listed, waters could not be settled, segfault. I tried to run few other peptides and proteins with tip5p, 100% of crashes. > Also, tip5p has been used successfully with the charmm foce field: > http://pubs.acs.org/doi/abs/10.1021/ct700053u Yes, there are no dimensionless charged hydrogens there except of these on ... tip3p, tip4p and tip5p water. > (...) > http://pubs.acs.org/doi/abs/10.1021/ct300180w > (...) > Are you using halogens with oplsaa-x? No, I use OPLSAA distributed with Gromacs. > Standard oplsaa (non-x) and tip5p seem to be a fine combination: > http://www.sciencedirect.com/science/article/pii/S0006349507713863 > > you might want to contact the authors of that paper to see if they > ever had such problems. Thank you, I will ask them. I am starting to be sure, though, that it is not possible to run any simulations with tip5p and proteins (containing arginine for example) without some tricks, strengthening/constraining bonds at least. It happened during the minimization, that the distance between HHxy and NHx in arginine grew up to 1.48A, while HHxy was landing on LP of the water. Repeatability is 100%. Regards, Grzegorz -- gmx-users mailing list gmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/Search before posting! * Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists