# [gmx-users] Formula for cosine content in essential dynamics

Berk Hess gmx3 at hotmail.com
Mon Aug 27 15:11:11 CEST 2007

```

>From: "Steve Spronk" <spronk at umich.edu>
>Reply-To: Discussion list for GROMACS users <gmx-users at gromacs.org>
>To: <gmx-users at gromacs.org>
>Subject: [gmx-users] Formula for cosine content in essential dynamics
>Date: Fri, 24 Aug 2007 17:04:04 -0400
>
>Hi,
>
>
>
>I am working on essential dynamics analysis of my protein, and I have a
>basic question on the calculation of the cosine content of the principal
>components.
>
>
>
>The way that I understand it, the first few principal components may
>resemble cosine functions, with the period equal to half the principal
>component index i.  However, in the formula for cosine component listed in
>Hess's papers, the Gromacs manual, and all other sources I've checked, the
>cosine function is
>
>
>
>cos(i*pi*t)
>
>
>
>which is not actually a cosine that has a period of half the index.  I must
>be misunderstanding something, because if you want to compare to a cos
>function that passes through half a period over the whole time of the
>simulation, shouldn't the formula include this instead:
>
>
>
>cos(i*pi*t/T)
>
>
>
>where T is the length of the simulation?  This way, it also appears to me
>that the units work out correctly.
>
>
>
>I'm confused.  Can you help me understand the formula better?
>
>
>
>Thanks,
>
>Steve
>

I probably forgot to add a factor of 1/T in the manual.
I'll fix it.

Thanks,

Berk.

>_______________________________________________
>gmx-users mailing list    gmx-users at 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 gmx-users-request at gromacs.org.