Jian Dai wrote:
Hi, dear all (especially for David):
When reading the paper "Thermodynamics of hydrogen bonding in hydrophilic and hydrophobic media" and its supplementary material and the code "gmx_hbond.c" for the analysis of hydrogen bond, I got questions: 1. When calculating the autocorrelation function of the interrupted lifetime using -ac option, is the error of the calculation given by "Chi^2" in the output?
Don't know.

2. As can be seen from the supplementary material of the paper, the forward and backward rate k and k' are dependent on the temperature, but from the code, I can not identify the link between "temp" and "k" since "temp" is only used in "calc_dg()", which means "temp" and "k" are not related. Am I right?

You have to do simulations at different Temperature.'

Do also check this paper which we are busy implementing in the program:
Reversible geminate recombination of hydrogen-bonded water  molecule pair
Omer Markovitcha and Noam Agmonb THE JOURNAL OF CHEMICAL PHYSICS 129, 084505 2008


Thank you.
Jian



--
David.
________________________________________________________________________
David van der Spoel, PhD, Professor of Biology
Dept. of Cell and Molecular Biology, Uppsala University.
Husargatan 3, Box 596,          75124 Uppsala, Sweden
phone:  46 18 471 4205          fax: 46 18 511 755
[email protected]    [email protected]   http://folding.bmc.uu.se
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
--
gmx-users mailing list    [email protected]
http://lists.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