Re: [gmx-users] Re: rotational correlation function

2010-09-24 Thread David van der Spoel
On 2010-09-23 21.10, Paymon Pirzadeh wrote: OK, To summurize: 1. I calculate the Ctot(t) with g_rotacf using N-H index files. 2. Use trjconv to get a trajectory of internal motions. => How I get the Cint(t)??? 3. Crot = Ctot/Cint But, is there any way to calculate the Crot(t) directly?

Re: [gmx-users] Re: rotational correlation function

2010-09-24 Thread Paymon Pirzadeh
OK, To summurize: 1. I calculate the Ctot(t) with g_rotacf using N-H index files. 2. Use trjconv to get a trajectory of internal motions. => How I get the Cint(t)??? 3. Crot = Ctot/Cint But, is there any way to calculate the Crot(t) directly? Is trjconv the solution again? regards, Paymo

Re: [gmx-users] Re: rotational correlation function

2010-09-23 Thread David van der Spoel
On 2010-09-23 20.38, Paymon Pirzadeh wrote: Great tips and advices. I appreciate your attention. But, what would the g_rms -fit rot+trans do? compute the rmsd. something different. I said trjconv. Paymon On Thu, 2010-09-23 at 20:07 +0200, David van der Spoel wrote: On 2010-09-23 18.14, P

Re: [gmx-users] Re: rotational correlation function

2010-09-23 Thread Paymon Pirzadeh
Great tips and advices. I appreciate your attention. But, what would the g_rms -fit rot+trans do? Paymon On Thu, 2010-09-23 at 20:07 +0200, David van der Spoel wrote: > On 2010-09-23 18.14, Paymon Pirzadeh wrote: > > Dear Dr. van der Spoel, > > Thanks a lot for your clarification. But some qu

Re: [gmx-users] Re: rotational correlation function

2010-09-23 Thread David van der Spoel
On 2010-09-23 18.14, Paymon Pirzadeh wrote: Dear Dr. van der Spoel, Thanks a lot for your clarification. But some questions yet: 1. If I do not provide the "-aver" option, will the program calculate the C(t) for individual N-H bond? Yes. 2. In a paper that I just read, JCP vol. 131 155103 (2

Re: [gmx-users] Re: rotational correlation function

2010-09-23 Thread Paymon Pirzadeh
Dear Dr. van der Spoel, Thanks a lot for your clarification. But some questions yet: 1. If I do not provide the "-aver" option, will the program calculate the C(t) for individual N-H bond? 2. In a paper that I just read, JCP vol. 131 155103 (2009), they claim they have calculated the Ctot(t) whi

Re: [gmx-users] Re: rotational correlation function

2010-09-22 Thread David van der Spoel
On 2010-09-23 01.27, Paymon Pirzadeh wrote: Dear Dr. Chaban, I meant the N-H bond vectors of protein backbone for calculation of rotational time correlation function to calculate the rotational diffusion constant of my protein. I need a protocol which walks me step by step through the procedure.

[gmx-users] Re: rotational correlation function

2010-09-22 Thread Vitaly Chaban
On Wed, Sep 22, 2010 at 7:47 PM, Paymon Pirzadeh wrote: > Well, I found a paper: J. Chem Phys. Vol. 131, 155103 (2009) in which > they have used GROMACS and its utilities to calculate the correlation > functions I need, however, it is not clear to me what utilities were > used and what steps were

Re: [gmx-users] Re: rotational correlation function

2010-09-22 Thread Vitaly Chaban
Hmm... Rotational time..? If you need diffusion constant, just use g_rotacf and then integrate the resulting VACF. If I recall correctly, you should make an index file with three numbers to define two vectors for the calculation. On Wed, Sep 22, 2010 at 7:27 PM, Paymon Pirzadeh wrote: > Dear Dr.

Re: [gmx-users] Re: rotational correlation function

2010-09-22 Thread Paymon Pirzadeh
Dear Dr. Chaban, I meant the N-H bond vectors of protein backbone for calculation of rotational time correlation function to calculate the rotational diffusion constant of my protein. I need a protocol which walks me step by step through the procedure. Thanks for your Attention. Paymon On Wed, 2

[gmx-users] Re: rotational correlation function

2010-09-22 Thread Vitaly Chaban
> Message: 2 > Date: Wed, 22 Sep 2010 12:07:27 -0600 > From: Paymon Pirzadeh > Subject: [gmx-users] rotational correlation function > To: gmx-users@gromacs.org > Message-ID: <1285178847.11669.66.ca...@paymon-desktop> > Content-Type: text/plain > > Hello, > Imagine if the total correlation function