Mahnam a écrit :
In God We Trust
Hello Dear Dr David van der Spoel
Thank you for your quick response. I had used constraints before md simulation. I have done energy minimization and position restraint (20ps) (with constraints=all bonds )and then md simulation at 100 K, 200K and 300K for 100ps (with constraint=none) and then I extended my simulation at 300K for 700ps (e.g. 1000ps md simulation) with time step 1fs and everything is ok, only if I change time step to 2fs, it crashes.

In order to use 2 ps time steps, you need to constrain the h-bonded interactions, using "h-bonds" instead of none should solve you problem.

Consider also coupling counter-ions and solvent as mentionned by Justin A. Lemkul.

Cheers,
Stéphane

--
Stéphane Téletchéa, PhD.                  http://www.steletch.org
Unité Mathématique Informatique et Génome http://migale.jouy.inra.fr/mig
INRA, Domaine de Vilvert                  Tél : (33) 134 652 891
78352 Jouy-en-Josas cedex, France         Fax : (33) 134 652 901
_______________________________________________
gmx-users mailing list    gmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

Reply via email to