Hello Aiswarya, Your program should look like this, in pseudo-code:
pseudo-code: do n_frames_of_your_trajectory g_mindist -f traj.trr -s .tpr -n (your_protein (first-group) , your_solvent (second-group)).ndx -o atm-pair.out -d your-distance-value -b (one frame of your trajectory) -e (the same frame) processing the atm-pair.out list and write them into an index file : 1.ndx, with each water in one group take just these atoms and residues, which are at your specific distance at your 'location' do measure the angle in relation to your protein for this frame for each specified water g_angle -f traj.trr -s .tpr -n 1.ndx -b (one frame of your trajectory) -e (the same frame) enddo read the angles from the output take just the one water with specific angle out of your list. (By using g_mindist -d , and g_angle you specify one specific point in relation to your protein) (By using -b and -e you specify one specific time) save the residue number of this water and your frame calculate the time of this frame enddo now you have an array of residue numbers at this particular location for each time analyse this array, then you know the maximum frame-time of one water-molecule at this specific location. You will not get any information about one particular time - just an average. This must be your residence time. Your application needs some little programming. >>> 24.08.11 18.14 Uhr >>> Emanuel, The residence time which am referring is the water molecule residing on the particular location at a specific time. Am not able to write a correct equation due to some character missing when typing from phone. AiswaryaSent from my BlackBerry® on Reliance Mobile, India's No. 1 Network. Go for it! From: aiswarya.pa...@gmail.com Date: Wed, 24 Aug 2011 16:13:05 +0000 To: Gromacs mail ReplyTo: aiswarya.pa...@gmail.com Subject: Re: Antw: Re: [gmx-users] Hydrogen existence time Emanuel, The residence time which am referring is the water molecule residing on the particular location at a specific time. Am not able to write a correct equation due to small character missing when typing from phone. AiswaryaSent from my BlackBerry® on Reliance Mobile, India's No. 1 Network. Go for it! From: "Emanuel Peter" Date: Wed, 24 Aug 2011 16:16:37 +0200 To: ; Subject: Antw: Re: [gmx-users] Hydrogen existence time As I told you, you can calculate the self-diffusion coefficient by using g_msd. Following: Large diffusion-coeff. --> means : The residence time is low. Small diffusion-coeff. --> means : The residence time is large !!!! Mean-square-displacement means: Displacement of molecules per time. No displacement of molecule at all : if(msd == 0) means: Infinite residence time. So I say from my idiotic point of view: msd ~ residence time ( I still do not know what you mean with this because this is usually related to chemical reactors in engineering and NOT in protein simulation ) In your index-file you can index those molecules, which you need. Your equation does not make sense. >>> 24.08.11 16.01 Uhr >>> Mark, The residence time equation goes like this- P(r)= function (N(t) - N(t+r))dt Where function(x-y) takes the value of 1 when x=y(x not equal to 0 and y not equal to 0) ie the site occupied and when x not equal to y or both x and y =0. N(t) index of water molecule found in the hydration site at time t. Aiswarya ------Original Message------ From: Mark Abraham To: Aiswarya Pawar Subject: Re: [gmx-users] Hydrogen existence time Sent: 24 Aug 2011 6:09 PM On 24/08/2011 10:31 PM, aiswarya.pa...@gmail.com wrote: > Any idea how can I calculate the residence time of water. First, seek define what you mean to calculate. If you can't write an equation down for it, you can't calculate it. Mark Sent from my BlackBerry® on Reliance Mobile, India's No. 1 Network. Go for it!
-- gmx-users mailing list gmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/Search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. Can't post? Read http://www.gromacs.org/Support/Mailing_Lists