Esteban Vöhringer wrote:
  Hi,

I was looking into the code, how the velocities are generated in grompp, once the option genvel is used in the .mdp file. I don't understand why the number of degrees of freedom which are used to calculate the scaling factor for the random velocities is 3N and not 3N-Nc ( Nc= Number of constraitns; N=Number of atoms).
is there some special reason?
because the constraints are not taking into account during vel. generation. before the first MD step the constraints are applied and the temperature is approximately correct.

Tank you in advance for your help

Esteban


------------------------------------------------------------------------

_______________________________________________
gmx-users mailing list    [email protected]
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


--
David van der Spoel, Ph.D.
Molec. Biophys. group, Dept. of Cell & Molec. Biol., Uppsala University.
Box 596, 75124 Uppsala, Sweden. Phone:  +46184714205. Fax: +4618511755.
[EMAIL PROTECTED]       [EMAIL PROTECTED]   http://folding.bmc.uu.se
_______________________________________________
gmx-users mailing list    [email protected]
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