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>