From: Chris Harrison (charris5_at_gmail.com)
Date: Thu Nov 04 2010 - 06:49:40 CDT
Dear Jun Zhang,
The specific nature of your system determines which case is best suited for
you. Chandler's paper examines liquid water, in which case it is a system
of homogeneous pairs, and <> over all pairs describes a statistical
microscopic property of the bulk system; in this case it makes sense to
average over all pairs. In the case of a protein system, which
contains heterogeneous pairs that are not a regular lattice, <> over all
h-bond pairs does not seem particularly meaningful. In such an example, you
should most likely restrict <> to a particular h-bonded pair over your
trajectory, unless you are examining for example the water-to-water
h-bonding properties of water solvating the protein.
Best,
Chris
On Thu, Nov 4, 2010 at 12:55 AM, Jun Zhang <coolrainbow_at_yahoo.cn> wrote:
> Hello everyone:
>
> I would like to calculate the HB correlation function according to:
>
> c(t)=<h(0)h(t)>/<h>
>
> in "Nature,379:55-57". However, I have been convoluted by the meaning of
> <>:
> Is it an average over all atom pair at a particular time or over the
> entire trajectory for a particular pair? I prefer the second definition, but
> I am not sure. Please help me. Thank you in advance.
>
> Jun Zhang
> Nankai University
> coolrainbow_at_yahoo.cn
>
>
>
>
>
>
-- Chris Harrison, Ph.D. Theoretical and Computational Biophysics Group NIH Resource for Macromolecular Modeling and Bioinformatics Beckman Institute for Advanced Science and Technology University of Illinois, 405 N. Mathews Ave., Urbana, IL 61801 char_at_ks.uiuc.edu Voice: 217-244-1733 http://www.ks.uiuc.edu/~char Fax: 217-244-6078
This archive was generated by hypermail 2.1.6 : Wed Feb 29 2012 - 15:56:19 CST