Dear gromacs users,<br><br>I did a g_hbond -ac analysis on a 10 ns trajectory of a polycarboxylate in water. Following is the output I obtained for the polycarboxylate H-bondign with water<br><br>ACF 36008/36008<br>Normalization for c(t) = 0.00709164 for gh(t) = 7.09124e-07<br>
Hydrogen bond thermodynamics at T = 298.15 K<br>Fitting parameters chi^2 = 0.000661899<br>Q = 0<br>--------------------------------------------------<br>Type Rate (1/ps) Time (ps) DG (kJ/mol) Chi^2<br>Forward 0.300 3.334 7.513 0.000661899<br>
Backward 1.608 0.622 3.350<br>One-way 0.089 11.284 10.536<br>Integral 0.039 25.374 12.544<br>Relaxation 0.062 16.080 11.413<br><br>I have read David paper (2006) on thermodynamics of H-bonding. In the paper two delta G values are mentioned. one is the Gibbs energy of activation (dG#) and another is the free energy of H-bonding. In the output given from the -ac analysis, which dG values do these correspond to? I assumed dG given for the forward rate correspond to the energy of activation. Is it correct? What is the value for free energy of H-bonding.<br>
<br>Another question I have is, for a polycarboxylic acid system in water (COOH groups as side chains) for a 10ns run, I got the following as the output where the life time of the H-bond (between the polyacid and water) is negative. Is this correct?<br>
<br>ACF 22057/22057<br>Normalization for c(t) = 0.0217714 for gh(t) = 4.35398e-06<br>Hydrogen bond thermodynamics at T = 298.15 K<br>Fitting parameters chi^2 = 0.0146697<br>Q = 0<br>--------------------------------------------------<br>
Type Rate (1/ps) Time (ps) DG (kJ/mol) Chi^2<br>Forward -0.271 -3.687 -666.000 0.0146697<br>Backward -2.291 -0.437 -666.000<br>One-way 0.101 9.882 10.207<br>Integral 0.034 29.376 12.907<br>
Relaxation 0.063 15.803 11.370<br><br><br>Thankyou for any help,<br><br>Dr. M. S. Sulatha<br>Dept. of Chemical Engineering<br>IIT-Madras<br>India<br>