Thank you for your answer Mark. I was not able to find this
information in manual. If it is present, could you please refer me to
the relevant section. If it is indeed not there, it would be very
helpful to include it in the manual.

Sincerely,
Ravi Bhadauria



On Thu, Nov 10, 2011 at 6:18 AM, Mark Abraham <mark.abra...@anu.edu.au> wrote:
> On 9/11/2011 9:44 PM, Ravi Bhadauria wrote:
>>
>> Hello users,
>>
>> I have few conflicting answers from the user mailing list about the
>> following question.
>>
>> Using leap-frog integrator, what is the velocity that is WRITTEN in
>> the trajectory file for a frame corresponding to r(t)? Is it v(t-dt/2)
>> or v(t). Has the protocol for writing velocity in trajectory file been
>> consistent from version 3.x to 4.5.x?
>
> With leap-frog, it is v(t-dt/2) in 4.5.x. v(t) is never even computed (see
> do_update_md() in src/mdlib/update.c). I don't know if Erik's statement was
> accurate in 2002, but it isn't accurate now.
>
> Mark
>
>>
>> Here are few posts hinting its v(t):
>> http://lists.gromacs.org/pipermail/gmx-users/2002-July/001969.html
>> http://lists.gromacs.org/pipermail/gmx-users/2003-September/006956.html
>>
>> Here are couple of posts hinting its v(t-dt/2):
>> http://lists.gromacs.org/pipermail/gmx-users/2006-June/022261.html
>> http://lists.gromacs.org/pipermail/gmx-users/2003-March/004827.html
>>
>> Please let me know. The problem is with certain class of simulations
>> of bounded flows, these two definitions provide contrasting results
>> for suitably averaged velocity.
>>
>> Sincerely
>> Ravi Bhadauria
>
> --
> 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
>
--
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

Reply via email to