Dear gromacs users and developers, I have a question about the methodology that gromacs h_bond module compute H_bond correlation function, C(t) from trajectory.
I want to compute hydrogen bond life time between a carbonyl oxygen of single Etoac and water. For example, I made a system which contain a single Etoac molecules and about 2000 water molecules. then I used h_bond module to get hbmap.xpm and hbond.ndx and I analyzed those two so that I can get existence function of each carbonyl oxygen - water pair. In detail, during 500 ps simulation, 23 water molecules formed hydrogen bond with carbonyl oxygen of Etoac at least one time. So, I computed 23 existence functions, each having 500 ps length. Then I use those 23 set of existence function to compute h_bond correlation function. I'll upload Matlab code (corr_031116.m )with which I computed correlation function and please tell me what the problem is.. So, What I want to know is that How Gromacs practically compute Correlation function, in detail.. Thank you for reading my email and Thank you for your reply in advance..
-- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-requ...@gromacs.org.