[gmx-users] mailing list migration
Hi, Today we are migrating the mailinglist server to a new machine, so there might be delays in the delivery of messages. Cheers, Rossen -- gmx-users mailing listgmx-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
Re: [gmx-users] Re: Bilayer COM removal issue: Large VCM
Hi Rajat, If you remove comm on the bilayer, there may be relative comm between leaflets. If that relative motion is significant and you switch to removing comm per leaflet, the program suddenly finds itself resetting the com over a large distance. About equilibration, you equilibrated with comm_grps = SOL DMPC, the system is not equilibrated for another scheme. You can solve this issue by regenerating velocities, or by running short cycles with the time step increasing from very small to normal. Hope it helps, Tsjerk On Wed, Nov 13, 2013 at 8:06 AM, rajat desikan rajatdesi...@gmail.comwrote: Hi All, Any suggestions? Thanks, On Mon, Nov 11, 2013 at 12:38 AM, rajat desikan rajatdesi...@gmail.com wrote: Hi All, I am experiencing a few problems in membrane simulations wrt COM removal. I downloaded a 400 ns pre-equilibrated Slipid-DMPC membrane with all the accompanying files. I then carried out the following steps: 1) energy minimization 2) NVT Eq - 100 ps 3) NPT Eq - 250 ps (Berendsen temp, Pres coupling) Then I used g_select to select the upper and lower DMPC leaflets. The then carried out a 250 ps NPT eq again. The only change was: comm-grps= SOL DMPC == comm-grps= SOL upper lower On every step in log file, I get the following message: *Step Time Lambda 124000 248.00.0 Large VCM(group lower): -0.00051, -0.00515, -0.00652, Temp-cm: 8.11828e+29 Energies (kJ/mol)U-BProper Dih. Improper Dih. LJ-14 Coulomb-147.23818e+044.19778e+046.46641e+024.54801e+03 -1.45245e+05 LJ (SR)LJ (LR) Disper. corr. Coulomb (SR) Coul. recip.2.79689e+04 -3.78407e+03 -2.10679e+03 -5.84134e+05 -8.87497e+04 PotentialKinetic En. Total EnergyTemperature Pres. DC (bar)-6.76497e+051.76468e+05 -5.00029e+05 3.10424e+02 -1.05704e+02 Pressure (bar) Constr. rmsd -1.85927e+02 6.42934e-06* *Large VCM(group lower): -0.00187, -0.00369, 0.00032, Temp-cm: 2.02076e+29 Large VCM(group lower): -0.00725, -0.00278, -0.00549, Temp-cm: 1.05988e+30Large VCM(group lower): 0.00020, 0.00308, -0.00176, Temp-cm: 1.48126e+29Large VCM(group lower): -0.00541, 0.00546, -0.00166, Temp-cm: 7.24656e+29 Large VCM(group lower): -0.00220, 0.00362, -0.00741, Temp-cm: 8.53812e+29Large VCM(group lower): 0.00140, -0.00160, 0.00029, Temp-cm: 5.39679e+28Large VCM(group lower): -0.00056, -0.00293, -0.00364, Temp-cm: 2.59422e+29 Large VCM(group lower): -0.00172, -0.00260, 0.00494, Temp-cm: 3.99945e+29Large VCM(group lower): 0.00252, 0.00594, 0.00068, Temp-cm: 4.93342e+29* *DD step 124999 vol min/aver 0.702 load imb.: force 1.3% pme mesh/force 0.636* I do not know what to make of it. There are no issues when I remove COM for the entire system. I have seen this issue come up a few times in the archives too, but I didn't find a satisfactory solution since the bilayer was very well equilibrated. I would appreciate any suggestions. Thank you. -- Rajat Desikan (Ph.D Scholar) Prof. K. Ganapathy Ayappa's Lab (no 13), Dept. of Chemical Engineering, Indian Institute of Science, Bangalore -- Rajat Desikan (Ph.D Scholar) Prof. K. Ganapathy Ayappa's Lab (no 13), Dept. of Chemical Engineering, Indian Institute of Science, Bangalore -- gmx-users mailing listgmx-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 -- Tsjerk A. Wassenaar, Ph.D. -- gmx-users mailing listgmx-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
Re: [gmx-users] help about ibi
Hi, Something went wrong earlier in your workflow. Check your log files, etc. Mark On Nov 13, 2013 3:57 AM, guozhicheng222 guozhicheng...@126.com wrote: Hi: When I am running the ibi procedure, I get the following error message: A coordinate in file conf.gro does not contain a '.' Additionally, I check the coordinate file of confout.gro in step_001. It showed that 'nan' symbol appeared in confout.gro. What is wrong with this? How can I fix it? I am very appreciating for anyone's help. Best Wishes! Zhicheng Guo -- gmx-users mailing listgmx-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 -- gmx-users mailing listgmx-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
Re: [gmx-users] Re: Bilayer COM removal issue: Large VCM
Hi Tsjerk, That was very sage advice! Thank you. I will try regenerating velocities and see if the motion goes away... On Wed, Nov 13, 2013 at 2:00 PM, Tsjerk Wassenaar tsje...@gmail.com wrote: Hi Rajat, If you remove comm on the bilayer, there may be relative comm between leaflets. If that relative motion is significant and you switch to removing comm per leaflet, the program suddenly finds itself resetting the com over a large distance. About equilibration, you equilibrated with comm_grps = SOL DMPC, the system is not equilibrated for another scheme. You can solve this issue by regenerating velocities, or by running short cycles with the time step increasing from very small to normal. Hope it helps, Tsjerk On Wed, Nov 13, 2013 at 8:06 AM, rajat desikan rajatdesi...@gmail.com wrote: Hi All, Any suggestions? Thanks, On Mon, Nov 11, 2013 at 12:38 AM, rajat desikan rajatdesi...@gmail.com wrote: Hi All, I am experiencing a few problems in membrane simulations wrt COM removal. I downloaded a 400 ns pre-equilibrated Slipid-DMPC membrane with all the accompanying files. I then carried out the following steps: 1) energy minimization 2) NVT Eq - 100 ps 3) NPT Eq - 250 ps (Berendsen temp, Pres coupling) Then I used g_select to select the upper and lower DMPC leaflets. The then carried out a 250 ps NPT eq again. The only change was: comm-grps= SOL DMPC == comm-grps= SOL upper lower On every step in log file, I get the following message: *Step Time Lambda 124000 248.00.0 Large VCM(group lower): -0.00051, -0.00515, -0.00652, Temp-cm: 8.11828e+29 Energies (kJ/mol)U-BProper Dih. Improper Dih. LJ-14 Coulomb-147.23818e+044.19778e+046.46641e+024.54801e+03 -1.45245e+05 LJ (SR)LJ (LR) Disper. corr. Coulomb (SR) Coul. recip.2.79689e+04 -3.78407e+03 -2.10679e+03 -5.84134e+05 -8.87497e+04 PotentialKinetic En. Total Energy Temperature Pres. DC (bar)-6.76497e+051.76468e+05 -5.00029e+05 3.10424e+02 -1.05704e+02 Pressure (bar) Constr. rmsd -1.85927e+02 6.42934e-06* *Large VCM(group lower): -0.00187, -0.00369, 0.00032, Temp-cm: 2.02076e+29 Large VCM(group lower): -0.00725, -0.00278, -0.00549, Temp-cm: 1.05988e+30Large VCM(group lower): 0.00020, 0.00308, -0.00176, Temp-cm: 1.48126e+29Large VCM(group lower): -0.00541, 0.00546, -0.00166, Temp-cm: 7.24656e+29 Large VCM(group lower): -0.00220, 0.00362, -0.00741, Temp-cm: 8.53812e+29Large VCM(group lower): 0.00140, -0.00160, 0.00029, Temp-cm: 5.39679e+28Large VCM(group lower): -0.00056, -0.00293, -0.00364, Temp-cm: 2.59422e+29 Large VCM(group lower): -0.00172, -0.00260, 0.00494, Temp-cm: 3.99945e+29Large VCM(group lower): 0.00252, 0.00594, 0.00068, Temp-cm: 4.93342e+29* *DD step 124999 vol min/aver 0.702 load imb.: force 1.3% pme mesh/force 0.636* I do not know what to make of it. There are no issues when I remove COM for the entire system. I have seen this issue come up a few times in the archives too, but I didn't find a satisfactory solution since the bilayer was very well equilibrated. I would appreciate any suggestions. Thank you. -- Rajat Desikan (Ph.D Scholar) Prof. K. Ganapathy Ayappa's Lab (no 13), Dept. of Chemical Engineering, Indian Institute of Science, Bangalore -- Rajat Desikan (Ph.D Scholar) Prof. K. Ganapathy Ayappa's Lab (no 13), Dept. of Chemical Engineering, Indian Institute of Science, Bangalore -- gmx-users mailing listgmx-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 -- Tsjerk A. Wassenaar, Ph.D. -- gmx-users mailing listgmx-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 -- Rajat Desikan (Ph.D Scholar) Prof. K. Ganapathy Ayappa's Lab (no 13), Dept. of Chemical Engineering, Indian Institute of Science, Bangalore -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users * Please
Re: [gmx-users] Re: Bilayer COM removal issue: Large VCM
An update to anyone interested: Regenerating velocities by itself did not solve the problem. I had to regenerate velocities and couple the upper and lower leaflets separately to the thermostat to equilibrate the system. To smoothen the equilibration process further, I used a 0.5 fs timestep instead of 2 fs (though this is probably unnecessary). Thank you once more, Tsjerk. Old .mdp: comm-grps= SOL DMPC tcoupl = v-rescale; Thermostat tc-grps = DMPC SOL ; Couple lipids and SOL separately tau-t= 0.1 0.1 ; Time constant for temperature coupling ref-t= 310 310 ; Desired temperature (K) New .mdp: comm-grps= SOL upper lower tcoupl = v-rescale; Thermostat, v-rescale is also fine tc-grps = upper lower SOL ; Couple lipid leaflets and SOL separately tau-t= 0.1 0.1 0.1 ; Time constant for temperature coupling ref-t= 310 310 310 ; Desired temperature (K) On Wed, Nov 13, 2013 at 4:07 PM, rajat desikan rajatdesi...@gmail.comwrote: Hi Tsjerk, That was very sage advice! Thank you. I will try regenerating velocities and see if the motion goes away... On Wed, Nov 13, 2013 at 2:00 PM, Tsjerk Wassenaar tsje...@gmail.comwrote: Hi Rajat, If you remove comm on the bilayer, there may be relative comm between leaflets. If that relative motion is significant and you switch to removing comm per leaflet, the program suddenly finds itself resetting the com over a large distance. About equilibration, you equilibrated with comm_grps = SOL DMPC, the system is not equilibrated for another scheme. You can solve this issue by regenerating velocities, or by running short cycles with the time step increasing from very small to normal. Hope it helps, Tsjerk On Wed, Nov 13, 2013 at 8:06 AM, rajat desikan rajatdesi...@gmail.com wrote: Hi All, Any suggestions? Thanks, On Mon, Nov 11, 2013 at 12:38 AM, rajat desikan rajatdesi...@gmail.com wrote: Hi All, I am experiencing a few problems in membrane simulations wrt COM removal. I downloaded a 400 ns pre-equilibrated Slipid-DMPC membrane with all the accompanying files. I then carried out the following steps: 1) energy minimization 2) NVT Eq - 100 ps 3) NPT Eq - 250 ps (Berendsen temp, Pres coupling) Then I used g_select to select the upper and lower DMPC leaflets. The then carried out a 250 ps NPT eq again. The only change was: comm-grps= SOL DMPC == comm-grps= SOL upper lower On every step in log file, I get the following message: *Step Time Lambda 124000 248.00.0 Large VCM(group lower): -0.00051, -0.00515, -0.00652, Temp-cm: 8.11828e+29 Energies (kJ/mol)U-BProper Dih. Improper Dih. LJ-14 Coulomb-147.23818e+044.19778e+046.46641e+024.54801e+03 -1.45245e+05 LJ (SR)LJ (LR) Disper. corr. Coulomb (SR) Coul. recip.2.79689e+04 -3.78407e+03 -2.10679e+03 -5.84134e+05 -8.87497e+04 PotentialKinetic En. Total Energy Temperature Pres. DC (bar)-6.76497e+051.76468e+05 -5.00029e+05 3.10424e+02 -1.05704e+02 Pressure (bar) Constr. rmsd -1.85927e+02 6.42934e-06* *Large VCM(group lower): -0.00187, -0.00369, 0.00032, Temp-cm: 2.02076e+29 Large VCM(group lower): -0.00725, -0.00278, -0.00549, Temp-cm: 1.05988e+30Large VCM(group lower): 0.00020, 0.00308, -0.00176, Temp-cm: 1.48126e+29Large VCM(group lower): -0.00541, 0.00546, -0.00166, Temp-cm: 7.24656e+29 Large VCM(group lower): -0.00220, 0.00362, -0.00741, Temp-cm: 8.53812e+29Large VCM(group lower): 0.00140, -0.00160, 0.00029, Temp-cm: 5.39679e+28Large VCM(group lower): -0.00056, -0.00293, -0.00364, Temp-cm: 2.59422e+29 Large VCM(group lower): -0.00172, -0.00260, 0.00494, Temp-cm: 3.99945e+29Large VCM(group lower): 0.00252, 0.00594, 0.00068, Temp-cm: 4.93342e+29* *DD step 124999 vol min/aver 0.702 load imb.: force 1.3% pme mesh/force 0.636* I do not know what to make of it. There are no issues when I remove COM for the entire system. I have seen this issue come up a few times in the archives too, but I didn't find a satisfactory solution since the bilayer was very well equilibrated. I would appreciate any suggestions. Thank you. -- Rajat Desikan (Ph.D Scholar) Prof. K. Ganapathy Ayappa's Lab (no 13), Dept. of Chemical Engineering, Indian Institute of Science, Bangalore -- Rajat Desikan (Ph.D Scholar)
[gmx-users] Invalid order for directive defaults
Dear Justin Very thanks for your reply. I created a new topol.top file as below: 1) I used once default directive. 2) I put cnt.itp file in working directory. 3) I copied pr.top and renamed it to topol.top. I added #include cnt.itp in the end of topol.top file. I modified [ molecules ] directive. -- begining of topol.top file is as follows: ; Include forcefield parameters #include charmm27.ff/forcefield.itp [ moleculetype ] ; Namenrexcl Protein_chain_A 3 [ atoms ] . . . . end of com.top file is as follows: ; Include Position restraint file #ifdef POSRES #include posre.itp #endif #include cnt.itp ; Include water topology #include charmm27.ff/tip3p.itp #ifdef POSRES_WATER ; Position restraint for each water oxygen [ position_restraints ] ; i funct fcxfcyfcz 11 1000 1000 1000 #endif ; Include topology for ions #include charmm27.ff/ions.itp [ system ] ; Name Protein in water [ molecules ] ; Compound#mols Protein_chain_A 1 CNT 1 SOL 1388 --- Previous error (Invalid order for directive defaults) was solved, but When I used grompp -f minim.mdp -c system.gro -p topol.top -o minim.tpr, I encountered with this error: ERROR1 [file cnt.itp, line 2861]: No default Angle types . . . . . . ERROR 1218 [file cnt.itp, line 4078]: No default Angle types Fatal error: There were 1218 errors in input file(s). Lines 2861-4078 are related to [ angles ] directive in cnt.itp file. How to solve this issue? Any help will highly appreciated. -- gmx-users mailing listgmx-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
[gmx-users] error while running pdb2gmx
Hello GROMACS users, I have phosphorylated Serine residue in my protein (140 residues) of interest, now when I run pdb2gmx I get this following error Atom OXT in residue ALA 140 was not found in rtp entry ALA with 6 atoms while sorting atoms. I checked aminoacid.rtp, there is no separate entry for OXT there.When I did the simulation for the same protein prior phosphorylation I did not get this error. What is the reason for this and how should I rectify this error? Please help me with this regard Regards, Hasthi -- gmx-users mailing listgmx-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
[gmx-users] Invalid order for directive defaults
Dear Justin My cnt is infinite. I obtained cnt.top by g_x2top and then modified cnt.top to cnt.itp. For obtaining cnt.top, I used following files: --- ffcnt.atp: CA 12.01100 ; aromatic C --- ffcnt.n2t: CCA0.0012.011 3C 0.141 C 0.141 C 0.141 CCA0.0012.011 2C 0.141 C 0.141 --- ffcntbon.itp: [ bondtypes ] ; i j funcb0 kb CA CA 3 0.1418 47890.0 21.867 [ angletypes ] ; i j k functh0 cth ub0 cub CA CA CA 2 120.00 562.20 [ dihedraltypes ] ; i j k l funcphi0cp mult CA CA CA CA 5 0.00 25.12 0.00 0.00 --- ffcntnonbon.itp: [ atomtypes ] ;name at.num masscharge ptype sigma epsi CA 6 12.011000.00A 0.385 0.4396 --- In cnt.itp file, angle section of file is as follows: [ angles ] ; aiajak functc0c1c2c3 2 1 8 1 2 1 287 1 8 1 287 1 1 2 3 1 1 210 1 3 210 1 2 3 289 1 2 3 406 1 289 3 406 1 5 417 1 5 4 320 1 17 4 320 1 4 5 6 1 . . . . . . I saw system.gro file by VMD, there are all angles defined above in [angle] directive. -- gmx-users mailing listgmx-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
[gmx-users] about my single point calculation
Hello Mark, I don't get any informing of your reply by e-mail, but get your reply searched by google. Anyway thanks very much for your reply! Yeah, I used two totally different mdp files for the single point calculation because I thought that gromcs would report the potential energy of my system if I used the option -rerun, no matter what mdp files I used. Gromacs-4.5.5 was used in the calculations. Some time later, I also tried as you mentioned in the e-mail below. Problem was the same. I used this mdp (attached minim.mdp) file for both 0-step minimization and single point calculation with rerun: The only difference is the integrator, steep for 0-step minimization while md for single point calculation. And the following lines are the output I got: Single point calculation with rerun : Step Time Lambda 00.00.0 Energies (kJ/mol) G96Bond G96Angle Improper Dih. LJ-14 Coulomb-14 3.97878e+051.44370e+041.06677e+043.93431e+01 9.36593e+01 LJ (SR)Coulomb (SR) Potential Kinetic En. Total Energy 2.77574e+015.17380e+024.23661e+050.0e+00 4.23661e+05 Temperature Pressure (bar) 0.0e+000.0e+00 0-step minimization: Step Time Lambda 00.00.0 Energies (kJ/mol) G96Bond G96AngleImproper Dih. LJ-14 Coulomb-14 5.75700e+011.78703e+018.25973e-02 -9.22001e+00 5.93500e+01 LJ (SR)Coulomb (SR) Potential Pressure (bar) -1.96671e+013.75963e+024.81949e+020.0e+00 Later, I also used the mdp file which was pasted on the forum (attached sp.mdp) to do single point calculation with rerun, The following is what I got: Step Time Lambda 00.00.0 Energies (kJ/mol) G96BondG96Angle Improper Dih. LJ-14 Coulomb-14 3.97878e+051.44370e+041.06677e+043.93431e+01 9.36593e+01 LJ (SR) Coulomb (SR) Potential Kinetic En. Total Energy 2.77575e+015.17379e+024.23661e+056.54617e+01 4.23726e+05 Conserved En.Temperature Pressure (bar) 4.23726e+051.45800e+020.0e+00 I found that those energies are pretty much the same as the one mentioned above. The differences between corresponding energies are huge, I still don't understand the difference. My system only contains 37 atom, the energies generated from the 0-step minimization seem more reasonable than those from single point calculation with rerun. Do you know the possible reasons which could result in the huge difference? Or I made some mistake? Thanks very much! All the best, Qinghua On Wed, Nov 6, 2013 at 4:07 PM, fantasticqhl fantastic...@gmail.com wrote: Dear Justin, I am sorry for the late reply. I still can't figure it out. It isn't rocket science - your two .mdp files describe totally different model physics. To compare things, change as few things as necessary to generate the comparison. So use the same input .mdp file for the MD vs EM single-point comparison, just changing the integrator line, and maybe unconstrained-start (I forget the details). And be aware of http://www.gromacs.org/Documentation/How-tos/Single-Point_Energy Mark Could you please send me the mdp file which was used for your single point calculations. I want to do some comparison and then solve the problem. Thanks very much! All the best, Qinghua -- View this message in context: http://gromacs.5086.x6.nabble.com/single-point-calculation-with-gromacs-tp5012084p5012295.html Sent from the GROMACS Users Forum mailing list archive at Nabble.com. -- gmx-users mailing listgmx-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? Readhttp://www.gromacs.org/Support/Mailing_Lists -- gmx-users mailing listgmx-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? Readhttp://www.gromacs.org/Support/Mailing_Lists define = ;-DPOSRES integrator = md ; molecular dynamics algorithm tinit = 0.0 ; start time and timestep in ps dt = 0.002; time step in ps nsteps = 2; number of steps for 1000ns run emtol = 100; convergence criterion emstep
[gmx-users] How to construct mixed lipid bilayer
Dear All, can anyone tell me how to construct mixed lipid bilayer in gromacs id possible also provide me the command to construct the mixed bilayer Thanks in advance Nikhil -- -- gmx-users mailing listgmx-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
[gmx-users] Recompile Gromacs 4.6.3
Hello, all I intend to make some modification on minimize.c in mdlib. Do I need to do cmake make make install all over again? Or is there a quick way for recompiling? Thanks for any tips. JhengWei Li Institute of Atomic and Molecular Sciences, Academia Sinica, Taipei 106, Taiwan -- gmx-users mailing listgmx-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
Re: [gmx-users] Calculating diffusion coefficient in three dimension
On 11/13/13 12:20 AM, Venkat Reddy wrote: Dear Justin and Piggot, Thanks for the suggestions. Actually, I have constructed a CG lipid vesicle by placing lipids in random conformation in a simulation box. My lipid system is heterogeneous, i.e., it has different types of lipids (POPC,CHOLESTEROL,PPC...etc). Some of them occupy the outer layer of vesicle (POPC), and some stay in the intermediate region (CHOL). So, I want to calculate the diffusion rates of these lipids. Since POPC forms the surface (polar heads interacting with water and their tails points to the core), I suppose we have to calculate 2D diffusion for POPC. For the lipids in the core, they can diffuse in 3-dimension. So, it requires a 3D diffusion coefficient for these core lipids. How to calculate 2D and 3D diffusion coeff.? Hope I am clear. 2D diffusion coefficients are what the -lateral option does. I really don't understand why you want a 2D value for anything with spherical symmetry. If there is an outer layer of a vesicle, that's as much a sphere as anything inside it. -Justin On Wed, Nov 13, 2013 at 12:58 AM, Thomas Piggot t.pig...@soton.ac.ukwrote: Hi Venkat, Can you make it a bit clearer what you actually want? If it is the diffusion of the lipids along the curved surface of the vesicle, rather than simply the overall 3D diffusion, this is not trivial to calculate as I don't believe g_msd will do this for you. This property has been studied before though, so I suggest you search the literature for papers simulating vesicles to see how the lipid diffusion was calculated. Cheers Tom On 11/12/2013 06:35 PM, Justin Lemkul wrote: On 11/12/13 1:33 PM, Venkat Reddy wrote: Thanks Justin. So, I have to calculate diffusion coefficient three times (x,y,z) and finally add-up together to get in 3D??? If you just want the overall diffusion constant, that's what g_msd does without any additional options. -Justin On Tue, Nov 12, 2013 at 11:25 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/12/13 12:30 PM, Venkat Reddy wrote: Dear Sir, Thanks for the quick reply. So, I have to declare -type no flag. Isn't it?? The options for the -type flag are x, y, or z. You said you wanted the diffusion coefficient in each spatial dimension. That is precisely what this option will do. and I have recently gone through Justin's membrane protein tutorial, where You mean my tutorial :) he has calculated diffusion coefficient for lipids in a membrane by creating an index group for a particular atom. So, here also shall I do the same thing? Moreover, mine is a coarse-grained system. Yes, a representative atom is usually what is passed to g_msd. -Justin On Tue, Nov 12, 2013 at 9:57 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/12/13 11:25 AM, Venkat Reddy wrote: Then, how to mention the direction for spherical particles Sir? Read g_msd -h again, paying specific attention to the -type flag. -Justin On Tue, Nov 12, 2013 at 7:28 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/12/13 8:55 AM, Venkat Reddy wrote: Thank you sir for the prompt reply. *g_msd -f traj.xtc -s topol.tpr -n msd.ndx -lateral z -o msd.xvg -tu ns* Here I am giving -lateral z (like for membrane simulations). Is it fine for spherical systems also? No. The system is a sphere, so what use is it to calculate motion perpendicular to z when you have lipids moving in all three spatial dimensions? A vesicle is very different from a membrane, in which the lipids move in a plane, thus making -lateral z useful. -Justin On Tue, Nov 12, 2013 at 5:32 PM, Dr. Vitaly Chaban vvcha...@gmail.com wrote: MSD is 3D by default. Dr. Vitaly V. Chaban On Tue, Nov 12, 2013 at 6:01 AM, Venkat Reddy venkat...@gmail.com wrote: Dear all, I am simulating a spherical lipid vesicle. I want to calculate the diffusion coefficient for each lipid component in 3D. How to calculate it using g_msd (or any other tool like g_velacc)? Thank you for your concern -- With Best Wishes Venkat Reddy Chirasani -- gmx-users mailing listgmx-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 -- gmx-users mailing listgmx-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 -- == Justin A. Lemkul,
Re: [gmx-users] Invalid order for directive defaults
On 11/13/13 5:51 AM, Atila Petrosian wrote: Dear Justin Very thanks for your reply. I created a new topol.top file as below: 1) I used once default directive. 2) I put cnt.itp file in working directory. 3) I copied pr.top and renamed it to topol.top. I added #include cnt.itp in the end of topol.top file. I modified [ molecules ] directive. -- begining of topol.top file is as follows: ; Include forcefield parameters #include charmm27.ff/forcefield.itp [ moleculetype ] ; Namenrexcl Protein_chain_A 3 [ atoms ] . . . . end of com.top file is as follows: ; Include Position restraint file #ifdef POSRES #include posre.itp #endif #include cnt.itp ; Include water topology #include charmm27.ff/tip3p.itp #ifdef POSRES_WATER ; Position restraint for each water oxygen [ position_restraints ] ; i funct fcxfcyfcz 11 1000 1000 1000 #endif ; Include topology for ions #include charmm27.ff/ions.itp [ system ] ; Name Protein in water [ molecules ] ; Compound#mols Protein_chain_A 1 CNT 1 SOL 1388 --- Previous error (Invalid order for directive defaults) was solved, but When I used grompp -f minim.mdp -c system.gro -p topol.top -o minim.tpr, I encountered with this error: ERROR1 [file cnt.itp, line 2861]: No default Angle types . . . . . . ERROR 1218 [file cnt.itp, line 4078]: No default Angle types Fatal error: There were 1218 errors in input file(s). Lines 2861-4078 are related to [ angles ] directive in cnt.itp file. How to solve this issue? In your previous setup, you were effectively trying to use CHARMM27 + some other force field related to the CNT. You can't do that. What you can do (and need to do) is incorporate the nonbonded and bonded parameters related to the CNT into the parent force field. You may be able to simply #include the ffnonbonded.itp and ffbonded.itp files in the topology. Your current approach has simply deleted necessary information. -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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
Re: [gmx-users] error while running pdb2gmx
On 11/13/13 6:02 AM, hasthi wrote: Hello GROMACS users, I have phosphorylated Serine residue in my protein (140 residues) of interest, now when I run pdb2gmx I get this following error Atom OXT in residue ALA 140 was not found in rtp entry ALA with 6 atoms while sorting atoms. I checked aminoacid.rtp, there is no separate entry for OXT there.When I did the simulation for the same protein prior phosphorylation I did not get this error. What is the reason for this and how should I rectify this error? Please help me with this regard Presumably you have modified the force field to include the phosphorylated residue, correct? Have you followed every one of the steps shown at http://www.gromacs.org/Documentation/How-tos/Adding_a_Residue_to_a_Force_Field#Adding_a_new_residue? If you need further help, we will need more information, which will include (at minimum): 1. Snippet of the PDB file containing the problematic residue 2. Your exact pdb2gmx command 3. The screen output of pdb2gmx (all of it, not just the error message) -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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
Re: [gmx-users] How to construct mixed lipid bilayer
Start with packmol if you want to start from scratch. Or else get a pretty equilibrated mixed lipid bilayer if available somewhere on web. On Nov 13, 2013 6:45 PM, Nikhil Agrawal nikhil.08...@gmail.com wrote: Dear All, can anyone tell me how to construct mixed lipid bilayer in gromacs id possible also provide me the command to construct the mixed bilayer Thanks in advance Nikhil -- -- gmx-users mailing listgmx-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 -- gmx-users mailing listgmx-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
Re: [gmx-users] How to construct mixed lipid bilayer
Hi Nikhil, The first step would be to determine what forcefield you are going to use for the lipids. If you are going to use Charmm or Slipids, you can use charmmgui (just google it). If you are planning to use the Gromos forcefields, you can check Prof. Tieleman's website or lipidbook for pure bilayers and then build your own from scratch using packmol... Hope this helps... On Wed, Nov 13, 2013 at 6:44 PM, Nikhil Agrawal nikhil.08...@gmail.comwrote: Dear All, can anyone tell me how to construct mixed lipid bilayer in gromacs id possible also provide me the command to construct the mixed bilayer Thanks in advance Nikhil -- -- gmx-users mailing listgmx-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 -- Rajat Desikan (Ph.D Scholar) Prof. K. Ganapathy Ayappa's Lab (no 13), Dept. of Chemical Engineering, Indian Institute of Science, Bangalore -- gmx-users mailing listgmx-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
Re: [gmx-users] How to construct mixed lipid bilayer
Let me correct myself :). Its pre-equilibrated. Not pretty equilibrated. :) On Nov 13, 2013 7:23 PM, Arun kumar V arun.tar...@gmail.com wrote: Start with packmol if you want to start from scratch. Or else get a pretty equilibrated mixed lipid bilayer if available somewhere on web. On Nov 13, 2013 6:45 PM, Nikhil Agrawal nikhil.08...@gmail.com wrote: Dear All, can anyone tell me how to construct mixed lipid bilayer in gromacs id possible also provide me the command to construct the mixed bilayer Thanks in advance Nikhil -- -- gmx-users mailing listgmx-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 -- gmx-users mailing listgmx-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
Re: [gmx-users] Recompile Gromacs 4.6.3
For just modifying a file, just doing make is sufficient. I would recommend not installing the modified version (since you can run the build/src/kernel/mdrun directly), or if you must install, to use the suffixing options available in the ccmake advanced mode. Mark On Nov 13, 2013 2:48 PM, Jheng Wei Li lijheng...@gmail.com wrote: Hello, all I intend to make some modification on minimize.c in mdlib. Do I need to do cmake make make install all over again? Or is there a quick way for recompiling? Thanks for any tips. JhengWei Li Institute of Atomic and Molecular Sciences, Academia Sinica, Taipei 106, Taiwan -- gmx-users mailing listgmx-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 -- gmx-users mailing listgmx-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
Re: [gmx-users] error while running pdb2gmx
Probably the default behaviour of pdb2gmx for termini is not appropriate for your input. Use pdb2gmx -ter and choose wisely Mark On Nov 13, 2013 12:03 PM, hasthi durgs7kr...@gmail.com wrote: Hello GROMACS users, I have phosphorylated Serine residue in my protein (140 residues) of interest, now when I run pdb2gmx I get this following error Atom OXT in residue ALA 140 was not found in rtp entry ALA with 6 atoms while sorting atoms. I checked aminoacid.rtp, there is no separate entry for OXT there.When I did the simulation for the same protein prior phosphorylation I did not get this error. What is the reason for this and how should I rectify this error? Please help me with this regard Regards, Hasthi -- gmx-users mailing listgmx-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 -- gmx-users mailing listgmx-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
Re: [gmx-users] Calculating diffusion coefficient in three dimension
Dear Justin, I have referred to an article (Vuorela T, Catte A, Niemela PS, Hall A, Hyvonen MT, et al. (2010) Role of Lipids in Spheroidal High Density Lipoproteins. PLoS Comput Biol 6(10): e1000964. doi:10.1371/journal.pcbi.1000964), where the authors have clearly described the fitting of 2D diffusion coefficient to the surface lipids (diffusion along the lipid-water interface ) and 3D diffusion coefficient to the core lipids. On Wed, Nov 13, 2013 at 7:19 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/13/13 12:20 AM, Venkat Reddy wrote: Dear Justin and Piggot, Thanks for the suggestions. Actually, I have constructed a CG lipid vesicle by placing lipids in random conformation in a simulation box. My lipid system is heterogeneous, i.e., it has different types of lipids (POPC,CHOLESTEROL,PPC...etc). Some of them occupy the outer layer of vesicle (POPC), and some stay in the intermediate region (CHOL). So, I want to calculate the diffusion rates of these lipids. Since POPC forms the surface (polar heads interacting with water and their tails points to the core), I suppose we have to calculate 2D diffusion for POPC. For the lipids in the core, they can diffuse in 3-dimension. So, it requires a 3D diffusion coefficient for these core lipids. How to calculate 2D and 3D diffusion coeff.? Hope I am clear. 2D diffusion coefficients are what the -lateral option does. I really don't understand why you want a 2D value for anything with spherical symmetry. If there is an outer layer of a vesicle, that's as much a sphere as anything inside it. -Justin On Wed, Nov 13, 2013 at 12:58 AM, Thomas Piggot t.pig...@soton.ac.uk wrote: Hi Venkat, Can you make it a bit clearer what you actually want? If it is the diffusion of the lipids along the curved surface of the vesicle, rather than simply the overall 3D diffusion, this is not trivial to calculate as I don't believe g_msd will do this for you. This property has been studied before though, so I suggest you search the literature for papers simulating vesicles to see how the lipid diffusion was calculated. Cheers Tom On 11/12/2013 06:35 PM, Justin Lemkul wrote: On 11/12/13 1:33 PM, Venkat Reddy wrote: Thanks Justin. So, I have to calculate diffusion coefficient three times (x,y,z) and finally add-up together to get in 3D??? If you just want the overall diffusion constant, that's what g_msd does without any additional options. -Justin On Tue, Nov 12, 2013 at 11:25 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/12/13 12:30 PM, Venkat Reddy wrote: Dear Sir, Thanks for the quick reply. So, I have to declare -type no flag. Isn't it?? The options for the -type flag are x, y, or z. You said you wanted the diffusion coefficient in each spatial dimension. That is precisely what this option will do. and I have recently gone through Justin's membrane protein tutorial, where You mean my tutorial :) he has calculated diffusion coefficient for lipids in a membrane by creating an index group for a particular atom. So, here also shall I do the same thing? Moreover, mine is a coarse-grained system. Yes, a representative atom is usually what is passed to g_msd. -Justin On Tue, Nov 12, 2013 at 9:57 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/12/13 11:25 AM, Venkat Reddy wrote: Then, how to mention the direction for spherical particles Sir? Read g_msd -h again, paying specific attention to the -type flag. -Justin On Tue, Nov 12, 2013 at 7:28 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/12/13 8:55 AM, Venkat Reddy wrote: Thank you sir for the prompt reply. *g_msd -f traj.xtc -s topol.tpr -n msd.ndx -lateral z -o msd.xvg -tu ns* Here I am giving -lateral z (like for membrane simulations). Is it fine for spherical systems also? No. The system is a sphere, so what use is it to calculate motion perpendicular to z when you have lipids moving in all three spatial dimensions? A vesicle is very different from a membrane, in which the lipids move in a plane, thus making -lateral z useful. -Justin On Tue, Nov 12, 2013 at 5:32 PM, Dr. Vitaly Chaban vvcha...@gmail.com wrote: MSD is 3D by default. Dr. Vitaly V. Chaban On Tue, Nov 12, 2013 at 6:01 AM, Venkat Reddy venkat...@gmail.com wrote: Dear all, I am simulating a spherical lipid vesicle. I want to calculate the diffusion coefficient for each lipid component in 3D. How to calculate it using g_msd (or any other tool like g_velacc)? Thank you for your concern -- With Best Wishes Venkat Reddy Chirasani -- gmx-users mailing listgmx-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
Re: [gmx-users] Calculating diffusion coefficient in three dimension
On 11/13/13 9:41 AM, Venkat Reddy wrote: Dear Justin, I have referred to an article (Vuorela T, Catte A, Niemela PS, Hall A, Hyvonen MT, et al. (2010) Role of Lipids in Spheroidal High Density Lipoproteins. PLoS Comput Biol 6(10): e1000964. doi:10.1371/journal.pcbi.1000964), where the authors have clearly described the fitting of 2D diffusion coefficient to the surface lipids (diffusion along the lipid-water interface Diffusion along a lipid-water interface is one thing. Trying to use g_msd do to it is another, because I don't think it will. It looks for a plane in the configuration and calculates relative to it. I suspect you will need to modify the code to implement a custom algorithm. -Justin ) and 3D diffusion coefficient to the core lipids. On Wed, Nov 13, 2013 at 7:19 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/13/13 12:20 AM, Venkat Reddy wrote: Dear Justin and Piggot, Thanks for the suggestions. Actually, I have constructed a CG lipid vesicle by placing lipids in random conformation in a simulation box. My lipid system is heterogeneous, i.e., it has different types of lipids (POPC,CHOLESTEROL,PPC...etc). Some of them occupy the outer layer of vesicle (POPC), and some stay in the intermediate region (CHOL). So, I want to calculate the diffusion rates of these lipids. Since POPC forms the surface (polar heads interacting with water and their tails points to the core), I suppose we have to calculate 2D diffusion for POPC. For the lipids in the core, they can diffuse in 3-dimension. So, it requires a 3D diffusion coefficient for these core lipids. How to calculate 2D and 3D diffusion coeff.? Hope I am clear. 2D diffusion coefficients are what the -lateral option does. I really don't understand why you want a 2D value for anything with spherical symmetry. If there is an outer layer of a vesicle, that's as much a sphere as anything inside it. -Justin On Wed, Nov 13, 2013 at 12:58 AM, Thomas Piggot t.pig...@soton.ac.uk wrote: Hi Venkat, Can you make it a bit clearer what you actually want? If it is the diffusion of the lipids along the curved surface of the vesicle, rather than simply the overall 3D diffusion, this is not trivial to calculate as I don't believe g_msd will do this for you. This property has been studied before though, so I suggest you search the literature for papers simulating vesicles to see how the lipid diffusion was calculated. Cheers Tom On 11/12/2013 06:35 PM, Justin Lemkul wrote: On 11/12/13 1:33 PM, Venkat Reddy wrote: Thanks Justin. So, I have to calculate diffusion coefficient three times (x,y,z) and finally add-up together to get in 3D??? If you just want the overall diffusion constant, that's what g_msd does without any additional options. -Justin On Tue, Nov 12, 2013 at 11:25 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/12/13 12:30 PM, Venkat Reddy wrote: Dear Sir, Thanks for the quick reply. So, I have to declare -type no flag. Isn't it?? The options for the -type flag are x, y, or z. You said you wanted the diffusion coefficient in each spatial dimension. That is precisely what this option will do. and I have recently gone through Justin's membrane protein tutorial, where You mean my tutorial :) he has calculated diffusion coefficient for lipids in a membrane by creating an index group for a particular atom. So, here also shall I do the same thing? Moreover, mine is a coarse-grained system. Yes, a representative atom is usually what is passed to g_msd. -Justin On Tue, Nov 12, 2013 at 9:57 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/12/13 11:25 AM, Venkat Reddy wrote: Then, how to mention the direction for spherical particles Sir? Read g_msd -h again, paying specific attention to the -type flag. -Justin On Tue, Nov 12, 2013 at 7:28 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/12/13 8:55 AM, Venkat Reddy wrote: Thank you sir for the prompt reply. *g_msd -f traj.xtc -s topol.tpr -n msd.ndx -lateral z -o msd.xvg -tu ns* Here I am giving -lateral z (like for membrane simulations). Is it fine for spherical systems also? No. The system is a sphere, so what use is it to calculate motion perpendicular to z when you have lipids moving in all three spatial dimensions? A vesicle is very different from a membrane, in which the lipids move in a plane, thus making -lateral z useful. -Justin On Tue, Nov 12, 2013 at 5:32 PM, Dr. Vitaly Chaban vvcha...@gmail.com wrote: MSD is 3D by default. Dr. Vitaly V. Chaban On Tue, Nov 12, 2013 at 6:01 AM, Venkat Reddy venkat...@gmail.com wrote: Dear all, I am simulating a spherical lipid vesicle. I want to calculate the diffusion coefficient for each lipid component in 3D. How to calculate it using g_msd (or any other tool like g_velacc)? Thank you for your concern -- With Best Wishes Venkat Reddy
[gmx-users] Invalid order for directive defaults
Dear Justin Thanks for your reply. In your previous setup, you were effectively trying to use CHARMM27 + some other force field related to the CNT. You can't do that. Thus, Gromacs is not appropriate for systems containing cnt. Is my deduction true? In my case, peptid + cnt + water molecules, what is your suggestion? Please guide me and explain more. How to do MD simulation of my system by gromacs? Any help will highly appreciated. -- gmx-users mailing listgmx-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
Re: [gmx-users] Invalid order for directive defaults
On 11/13/13 10:39 AM, Atila Petrosian wrote: Dear Justin Thanks for your reply. In your previous setup, you were effectively trying to use CHARMM27 + some other force field related to the CNT. You can't do that. Thus, Gromacs is not appropriate for systems containing cnt. Is my deduction true? Of course not. People simulate CNTs with Gromacs all the time. You just didn't construct the force field properly. In my case, peptid + cnt + water molecules, what is your suggestion? Please guide me and explain more. How to do MD simulation of my system by gromacs? You have missing parameters in the .top/.itp file. You have those parameters already in ffbonded.itp for the CNT. As long as those parameters are compatible with the peptide force field (CHARMM27), then you just need to add those parameters. Again, that may be as simple as adding #include cntffbonded.itp of whatever it is to the .top file after the #include statement for the parent force field. Your only problem was #including another force field that re-defined a [defaults] directive. That is syntactically illegal. Nothing else was inherently problematic, unless you're mixing incompatible force fields, but I haven't seen any evidence of that. -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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
[gmx-users] Invalid order for directive defaults
Dear Justin Thanks for your quick reply. I was confused. If I add #include ffcntbon.itp after #include cnt.itp in .top file, my problem was solved and error was solved? -- gmx-users mailing listgmx-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
Re: [gmx-users] Invalid order for directive defaults
On 11/13/13 11:53 AM, Atila Petrosian wrote: Dear Justin Thanks for your quick reply. I was confused. If I add #include ffcntbon.itp after #include cnt.itp in .top file, my problem was solved and error was solved? No. The parameters are at the force field level and thus have to be #included before any [moleculetype] is introduced (see Chapter 5 of the manual for required order of directives). If you do: #include charmm27.ff/forcefield.itp #include ffcntbon.itp you should be fine. -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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
[gmx-users] GROMACS 4.6.4 is released
Hi GROMACS users, GROMACS 4.6.4 is officially released. It contains numerous bug fixes, and some noteworthy simulation performance enhancements (particularly with GPUs!). We encourage all users to upgrade their installations from earlier 4.6-era releases. You can find the code, manual, release notes, installation instructions and test suite at the links below. Note that some tests have been added, and the manual has changed only in chapter 7 and appendix D. ftp://ftp.gromacs.org/pub/gromacs/gromacs-4.6.4.tar.gz ftp://ftp.gromacs.org/pub/manual/manual-4.6.4.pdf http://www.gromacs.org/About_Gromacs/Release_Notes/Versions_4.6.x#Release_notes_for_4.6.4 http://www.gromacs.org/Documentation/Installation_Instructions http://gromacs.googlecode.com/files/regressiontests-4.6.4.tar.gz Happy simulating! The GROMACS team -- gmx-users mailing listgmx-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
Re: [gmx-users] GROMACS 4.6.4 is released
Will a simulation from 4.6.1 continue running fine if I upgrade to 4.6.4? Hi GROMACS users, GROMACS 4.6.4 is officially released. It contains numerous bug fixes, and some noteworthy simulation performance enhancements (particularly with GPUs!). We encourage all users to upgrade their installations from earlier 4.6-era releases. You can find the code, manual, release notes, installation instructions and test suite at the links below. Note that some tests have been added, and the manual has changed only in chapter 7 and appendix D. ftp://ftp.gromacs.org/pub/gromacs/gromacs-4.6.4.tar.gz ftp://ftp.gromacs.org/pub/manual/manual-4.6.4.pdf http://www.gromacs.org/About_Gromacs/Release_Notes/Versions_4.6.x#Release_notes_for_4.6.4 http://www.gromacs.org/Documentation/Installation_Instructions http://gromacs.googlecode.com/files/regressiontests-4.6.4.tar.gz Happy simulating! The GROMACS team -- gmx-users mailing listgmx-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 -- gmx-users mailing listgmx-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
[gmx-users] GROMACS-4.6.3 CUDA version on multiple nodes each having 2 GPUs
Hello, I am trying to run MPI, OpenMP and CUDA enable GROMACS 4.6.3 on nodes having 12 cores (2 CPUs) and 2 GPUs (Tesla M2090) each. The problem is when I launch job GROMCAS is using only GPUs on first node come across and failing to use GPUs on other nodes. The command I used for two gpu enable nodes was, mpirun -np 2 mdrun -v -deffnm $configfile I tried with many other options but none of them worked. The one thing needs to remember here is that on all the nodes, GPUs got id 0 and 1 so -gpu_id option also didn't work. This old thread gave me some idea but I didn't understand it completely. http://lists.gromacs.org/pipermail/gmx-users/2013-March/079802.html Please suggests me the possible solutions for this issue. Thank you --Jignesh -- gmx-users mailing listgmx-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
[gmx-users] (no subject)
sir, I have a basic doubt about remd simulation. In remd is it possible to run 16 replicas in 8 processors? -- gmx-users mailing listgmx-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
[gmx-users] remd
sir, I have a basic doubt about remd simulation. In remd is it possible to run 16 replicas in 8 processors? -- gmx-users mailing listgmx-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
[gmx-users] shear viscosity
Dear users, g_energy -f *.edr -vis I have two questions about the results of eviscoi.xvg ( derivative of Einstein relation): 1.) I dont understand the unit of y-axis. It is kg.m^(-1).s^(-1).10^(-3) in B.Hess 2002 In eviscoi.xvg @yaxis label (kg m\S-1\N s\S-1\N ps) That is The unit of y-axis:kg.m^(-1).s^(-1).ps What is that? kg.m^(-1).s^(-1).ps equals to kg.m^(-1).s^(-1).10^(-3)? 2.) There are 5 columns in eviscoi.xvg. 1th is time. What are the rest? By the way, there is 216 water molecules in spc216.gro. But I want to calculate the shear viscosity of 512 water molecules. How can I get/derive 512 water molecules from spc216.gro? Can anyone give me some hint of this? Thanks in advance -- Ahmet Yıldırım -- gmx-users mailing listgmx-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
Re: [gmx-users] Re: g_analyze
Sorry, I attached the wrong file . Here's the average file generate from one of the files I sent in my last mail. I used the command g_analyze -f hbond_115-water.xvg -av hbond_115-water-avg.xvg. Here's the file obtained from this command :- https://www.dropbox.com/s/sovzk40cudznfjw/hbond_115-water-avg.xvg Now, if you see, the graph (in previous mail) and average file, both correlates well. I have a doubt about interpreting the result from g_analyze. The value 7.150740e+00 implies that on average 7 hydrogen bonds are formed during the simulation time of 5ns to 10ns. What does then the average file or its graph tells ?? On Mon, Nov 11, 2013 at 9:58 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/11/13 4:06 AM, bharat gupta wrote: In addition to my previous question, I have another question about g_analyze. When I used the hbond.xvg file to get the average and plotted the average.xvg file I found that the average value is round 4 to 5 according to the graph. But g_analyze in its final calculation gives 7.150 as the average values... Here's the link for the graph and result of average value calculated by g_analyze :- std. dev.relative deviation of standard - cumulants from those of set average deviation sqrt(n-1) a Gaussian distribition cum. 3 cum. 4 SS1 * 7.150740e+00 * 8.803173e-01 1.760635e-02 0.0620.163 SS2 1.490604e+00 1.164761e+00 2.329523e-02 0.4950.153 https://www.dropbox.com/s/1vqixenyerha7qq/115-water.png Here's the link hbond.xvg file and its averaged file https://www.dropbox.com/s/4n0m47o3mrjn3o8/hbond_115-water.xvg https://www.dropbox.com/s/4n0m47o3mrjn3o8/hbond_115-water.xvg Neither of these files produce output that corresponds to the PNG image above. Both files have values in 6-9 H-bond range and thus agree with the g_analyze output, which I can reproduce. I suspect you're somehow getting your files mixed up. -Justin On Mon, Nov 11, 2013 at 3:30 PM, bharat gupta bharat.85.m...@gmail.com wrote: thank you informing about g_rdf... Is it possible to dump the structure with those average water molecules interacting with the residues. I generated the hbond.log file which gives the details but I need to generate a figure for this ?? On Mon, Nov 11, 2013 at 10:40 AM, Justin Lemkul jalem...@vt.edu wrote: On 11/10/13 8:38 PM, bharat gupta wrote: But trjorder can be used to calculate the hydration layer or shell around residues ... Right ?? Yes, but I also tend to think that integrating an RDF is also a more straightforward way of doing that. With trjorder, you set some arbitrary cutoff that may or may not be an informed decision - with an RDF it is clear where the hydration layers are. -Justin On Mon, Nov 11, 2013 at 10:35 AM, Justin Lemkul jalem...@vt.edu wrote: On 11/10/13 8:30 PM, bharat gupta wrote: Thanks for your reply. I was missing the scientific notation part. Now everything is fine. Regarding trjorder, it doesn't measure h-bonds but gives the water nearest to protein. I wouldn't try to draw any sort of comparison between the output of trjorder and g_hbond. If you want to measure H-bonds, there's only one tool for that. -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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 -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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 --
Re: [gmx-users] Re: g_analyze
Hi, I tried g_select to dump the structure with the interacting water molecules, but I don't know know how to do that. I searched for some threads in the discussion but wasn't able to find anything related to my need. Can you explain how can I do that ? On Tue, Nov 12, 2013 at 7:39 AM, bharat gupta bharat.85.m...@gmail.comwrote: Sorry, I attached the wrong file . Here's the average file generate from one of the files I sent in my last mail. I used the command g_analyze -f hbond_115-water.xvg -av hbond_115-water-avg.xvg. Here's the file obtained from this command :- https://www.dropbox.com/s/sovzk40cudznfjw/hbond_115-water-avg.xvg Now, if you see, the graph (in previous mail) and average file, both correlates well. I have a doubt about interpreting the result from g_analyze. The value 7.150740e+00 implies that on average 7 hydrogen bonds are formed during the simulation time of 5ns to 10ns. What does then the average file or its graph tells ?? On Mon, Nov 11, 2013 at 9:58 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/11/13 4:06 AM, bharat gupta wrote: In addition to my previous question, I have another question about g_analyze. When I used the hbond.xvg file to get the average and plotted the average.xvg file I found that the average value is round 4 to 5 according to the graph. But g_analyze in its final calculation gives 7.150 as the average values... Here's the link for the graph and result of average value calculated by g_analyze :- std. dev.relative deviation of standard - cumulants from those of set average deviation sqrt(n-1) a Gaussian distribition cum. 3 cum. 4 SS1 * 7.150740e+00 * 8.803173e-01 1.760635e-02 0.0620.163 SS2 1.490604e+00 1.164761e+00 2.329523e-02 0.4950.153 https://www.dropbox.com/s/1vqixenyerha7qq/115-water.png Here's the link hbond.xvg file and its averaged file https://www.dropbox.com/s/4n0m47o3mrjn3o8/hbond_115-water.xvg https://www.dropbox.com/s/4n0m47o3mrjn3o8/hbond_115-water.xvg Neither of these files produce output that corresponds to the PNG image above. Both files have values in 6-9 H-bond range and thus agree with the g_analyze output, which I can reproduce. I suspect you're somehow getting your files mixed up. -Justin On Mon, Nov 11, 2013 at 3:30 PM, bharat gupta bharat.85.m...@gmail.com wrote: thank you informing about g_rdf... Is it possible to dump the structure with those average water molecules interacting with the residues. I generated the hbond.log file which gives the details but I need to generate a figure for this ?? On Mon, Nov 11, 2013 at 10:40 AM, Justin Lemkul jalem...@vt.edu wrote: On 11/10/13 8:38 PM, bharat gupta wrote: But trjorder can be used to calculate the hydration layer or shell around residues ... Right ?? Yes, but I also tend to think that integrating an RDF is also a more straightforward way of doing that. With trjorder, you set some arbitrary cutoff that may or may not be an informed decision - with an RDF it is clear where the hydration layers are. -Justin On Mon, Nov 11, 2013 at 10:35 AM, Justin Lemkul jalem...@vt.edu wrote: On 11/10/13 8:30 PM, bharat gupta wrote: Thanks for your reply. I was missing the scientific notation part. Now everything is fine. Regarding trjorder, it doesn't measure h-bonds but gives the water nearest to protein. I wouldn't try to draw any sort of comparison between the output of trjorder and g_hbond. If you want to measure H-bonds, there's only one tool for that. -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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 -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-users@gromacs.org
[gmx-users] installation error under openSuse 12.2
Date: Mon, 11 Nov 2013 08:27:05 -0800 (PST) From: kolnkempff kolnkem...@gmail.com Subject: [gmx-users] installation error under openSuse 12.2 To: gmx-users@gromacs.org Message-ID: 1384187225465-5012430.p...@n6.nabble.com Content-Type: text/plain; charset=us-ascii Dear gmx-users, I am trying to install gromacs-4.6.3 on an older Dell that is running openSuse 12.2 Using DGMX_BUILD_OWN_FFTW=ON failed for me so to get through cmake I had to compile fftw from scratch and I followed the recommendation of going directly to their website. Now I am at the make stage and get the following message: [ 67%] Built target gmxfftw make[2]: *** No rule to make target `//home/koln/bin/gromacs-4.6.3/build/src/contrib/fftw/gmxfftw-prefix/lib/libfftw3f.a', needed by `src/gmxlib/libgmx.so.8'. Stop. make[1]: *** [src/gmxlib/CMakeFiles/gmx.dir/all] Error 2 make: *** [all] Error 2 When I check for libfftw, I get: koln@linux-5bim:~/bin/gromacs-4.6.3/build2 rpm -qf /usr/lib/libfftw3f.* libfftw3-3-32bit-3.3.2-1.1.2.x86_64 Any suggestions for how to get past this obstacle would be greatly appreciated! This is an known issue in OpenSuse! We added a workaround a while ago: https://gerrit.gromacs.org/#/c/2540/ Until Gromacs release 4.6.4. you will have to change line 41 of src/contrib/fftw/CMakeLists.txt yourself or just install fftw through OpenSuse or by hand. Christoph Koln -- View this message in context: http://gromacs.5086.x6.nabble.com/installation-error-under-openSuse-12-2-tp5012430.html Sent from the GROMACS Users Forum mailing list archive at Nabble.com. -- Christoph Junghans Web: http://www.compphys.de -- gmx-users mailing listgmx-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
[gmx-users] Re: installation error under openSuse 12.2
Thank you so much Justin. On the one hand, I feel dumb because I could have sworn that I was using a clean build directory. On the other hand, I obviously lost track of what I was doing because your suggestion worked like a charm! Koln -- View this message in context: http://gromacs.5086.x6.nabble.com/installation-error-under-openSuse-12-2-tp5012430p5012436.html Sent from the GROMACS Users Forum mailing list archive at Nabble.com. -- gmx-users mailing listgmx-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
[gmx-users] Calculating diffusion coefficient in three dimension
Dear all, I am simulating a spherical lipid vesicle. I want to calculate the diffusion coefficient for each lipid component in 3D. How to calculate it using g_msd (or any other tool like g_velacc)? Thank you for your concern -- With Best Wishes Venkat Reddy Chirasani -- gmx-users mailing listgmx-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
Re: [gmx-users] hydrogen bond calculation problem
On 11/11/13 12:36 PM, Sushma Yadav wrote: Dear gromacs users, I have NaCl-water system of finite concentration .I want to calculate water-water hydrogen bond in the first solvation shell around ion..So I used g_hbond -shell option.For single ion in water,its easy to calculate..But say,for water-water hydrogen bonds around two Na+ ions,it is saying to Select one atom for shell. Although I have more than 1 Na+ ions and I want to see water-water hydrogen bond around all Na+ ions simultaneously.Is it possible to take solvation shell of each ion at a same time? According to the help description for the -shell option, no. Only one particle can be considered at a time. You'll have to analyze each Na+ separately, which is easily scripted. -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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
[gmx-users] Restarting a simulation after replacing an empty md.trr file
Dear Gromacs users, I am running a 50ns simulation of a protein having nearly 700 residues on 60 threads (Gromacs 4.6.3). At one point i got a disk space problem, so i have deleted the md.trr file and created an empty md.trr file. when i tried to restart the simulation from check point file on 100 threads, [ mdrun -v -deffnm md -cpi md.cpt -nt 100 ] i am getting a note and an error as fallows Reading checkpoint file md.cpt generated: #PME-nodes mismatch, current program: 100 checkpoint file: 60 Gromacs binary or parallel settings not identical to previous run. Continuation is exact, but is not guaranteed to be binary identical. ... Source code file: checkpoint.c, line: 1767 Fatal error: Can't read 1048576 bytes of 'md.trr' to compute checksum. The file has been replaced or its contents has been modified. please help me in overcoming this problem. Thanking you. -- Arun Kumar Somavarapu Project-JRF Dr. Pawan Gupta's lab Protein Science and Engineering Dept, Institute of Microbial Tecnology, Sec 39-A, Chandigarh - 160036. -- gmx-users mailing listgmx-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
Re: [gmx-users] Re: Reaction field zero and ions
On 11/11/13 12:08 PM, Williams Ernesto Miranda Delgado wrote: Hello If I did the MD simulation using PME and neutralized with ions, and I want to rerun this time with reaction field zero, is there any problem if I keep the ions? This is for LIE calculation. I am using AMBER99SB. Why do you think it necessary to delete them? -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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
Re: [gmx-users] Calculating diffusion coefficient in three dimension
MSD is 3D by default. Dr. Vitaly V. Chaban On Tue, Nov 12, 2013 at 6:01 AM, Venkat Reddy venkat...@gmail.com wrote: Dear all, I am simulating a spherical lipid vesicle. I want to calculate the diffusion coefficient for each lipid component in 3D. How to calculate it using g_msd (or any other tool like g_velacc)? Thank you for your concern -- With Best Wishes Venkat Reddy Chirasani -- gmx-users mailing listgmx-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 -- gmx-users mailing listgmx-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
Re: [gmx-users] Re: Reaction field zero and ions
There are no problems to have ions while using Reaction-Field treatment. Dr. Vitaly V. Chaban On Mon, Nov 11, 2013 at 7:06 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/11/13 12:08 PM, Williams Ernesto Miranda Delgado wrote: Hello If I did the MD simulation using PME and neutralized with ions, and I want to rerun this time with reaction field zero, is there any problem if I keep the ions? This is for LIE calculation. I am using AMBER99SB. Why do you think it necessary to delete them? -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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 -- gmx-users mailing listgmx-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
Re: [gmx-users] Re: g_analyze
On 11/11/13 5:39 PM, bharat gupta wrote: Sorry, I attached the wrong file . Here's the average file generate from one of the files I sent in my last mail. I used the command g_analyze -f hbond_115-water.xvg -av hbond_115-water-avg.xvg. Here's the file obtained from this command :- https://www.dropbox.com/s/sovzk40cudznfjw/hbond_115-water-avg.xvg Now, if you see, the graph (in previous mail) and average file, both correlates well. I have a doubt about interpreting the result from g_analyze. The value 7.150740e+00 implies that on average 7 hydrogen bonds are formed during the simulation time of 5ns to 10ns. What does then the average file or its graph tells ?? It's an average over sets. It is not equivalent to the output printed to the screen, nor is it supposed to. The value printed to the screen is the actual average of the data set of interest, as is intuitive from your values. An average of 4 is impossible if all the data points are in the range of 6-9. -Justin On Mon, Nov 11, 2013 at 9:58 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/11/13 4:06 AM, bharat gupta wrote: In addition to my previous question, I have another question about g_analyze. When I used the hbond.xvg file to get the average and plotted the average.xvg file I found that the average value is round 4 to 5 according to the graph. But g_analyze in its final calculation gives 7.150 as the average values... Here's the link for the graph and result of average value calculated by g_analyze :- std. dev.relative deviation of standard - cumulants from those of set average deviation sqrt(n-1) a Gaussian distribition cum. 3 cum. 4 SS1 * 7.150740e+00 * 8.803173e-01 1.760635e-02 0.0620.163 SS2 1.490604e+00 1.164761e+00 2.329523e-02 0.4950.153 https://www.dropbox.com/s/1vqixenyerha7qq/115-water.png Here's the link hbond.xvg file and its averaged file https://www.dropbox.com/s/4n0m47o3mrjn3o8/hbond_115-water.xvg https://www.dropbox.com/s/4n0m47o3mrjn3o8/hbond_115-water.xvg Neither of these files produce output that corresponds to the PNG image above. Both files have values in 6-9 H-bond range and thus agree with the g_analyze output, which I can reproduce. I suspect you're somehow getting your files mixed up. -Justin On Mon, Nov 11, 2013 at 3:30 PM, bharat gupta bharat.85.m...@gmail.com wrote: thank you informing about g_rdf... Is it possible to dump the structure with those average water molecules interacting with the residues. I generated the hbond.log file which gives the details but I need to generate a figure for this ?? On Mon, Nov 11, 2013 at 10:40 AM, Justin Lemkul jalem...@vt.edu wrote: On 11/10/13 8:38 PM, bharat gupta wrote: But trjorder can be used to calculate the hydration layer or shell around residues ... Right ?? Yes, but I also tend to think that integrating an RDF is also a more straightforward way of doing that. With trjorder, you set some arbitrary cutoff that may or may not be an informed decision - with an RDF it is clear where the hydration layers are. -Justin On Mon, Nov 11, 2013 at 10:35 AM, Justin Lemkul jalem...@vt.edu wrote: On 11/10/13 8:30 PM, bharat gupta wrote: Thanks for your reply. I was missing the scientific notation part. Now everything is fine. Regarding trjorder, it doesn't measure h-bonds but gives the water nearest to protein. I wouldn't try to draw any sort of comparison between the output of trjorder and g_hbond. If you want to measure H-bonds, there's only one tool for that. -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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 -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users * Please search the archive at
Re: [gmx-users] Restarting a simulation after replacing an empty md.trr file
On 11/12/13 1:15 AM, arun kumar wrote: Dear Gromacs users, I am running a 50ns simulation of a protein having nearly 700 residues on 60 threads (Gromacs 4.6.3). At one point i got a disk space problem, so i have deleted the md.trr file and created an empty md.trr file. when i tried to restart the simulation from check point file on 100 threads, [ mdrun -v -deffnm md -cpi md.cpt -nt 100 ] i am getting a note and an error as fallows Reading checkpoint file md.cpt generated: #PME-nodes mismatch, current program: 100 checkpoint file: 60 Gromacs binary or parallel settings not identical to previous run. Continuation is exact, but is not guaranteed to be binary identical. ... Source code file: checkpoint.c, line: 1767 Fatal error: Can't read 1048576 bytes of 'md.trr' to compute checksum. The file has been replaced or its contents has been modified. please help me in overcoming this problem. Use -noappend. If you mess with any of the files that mdrun produces, the file appending procedure (which is the default behavior) will fail. -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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
Re: [gmx-users] Re: g_analyze
On 11/11/13 6:56 PM, bharat gupta wrote: Hi, I tried g_select to dump the structure with the interacting water molecules, but I don't know know how to do that. I searched for some threads in the discussion but wasn't able to find anything related to my need. Can you explain how can I do that ? Start with g_select -select 'help all' and see what you can determine. Such selections are rather straightforward and have been explained several times on the list. If you need help, show us what you're doing and describe why it isn't what you want. It will ultimately save a lot of time. -Justin On Tue, Nov 12, 2013 at 7:39 AM, bharat gupta bharat.85.m...@gmail.comwrote: Sorry, I attached the wrong file . Here's the average file generate from one of the files I sent in my last mail. I used the command g_analyze -f hbond_115-water.xvg -av hbond_115-water-avg.xvg. Here's the file obtained from this command :- https://www.dropbox.com/s/sovzk40cudznfjw/hbond_115-water-avg.xvg Now, if you see, the graph (in previous mail) and average file, both correlates well. I have a doubt about interpreting the result from g_analyze. The value 7.150740e+00 implies that on average 7 hydrogen bonds are formed during the simulation time of 5ns to 10ns. What does then the average file or its graph tells ?? On Mon, Nov 11, 2013 at 9:58 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/11/13 4:06 AM, bharat gupta wrote: In addition to my previous question, I have another question about g_analyze. When I used the hbond.xvg file to get the average and plotted the average.xvg file I found that the average value is round 4 to 5 according to the graph. But g_analyze in its final calculation gives 7.150 as the average values... Here's the link for the graph and result of average value calculated by g_analyze :- std. dev.relative deviation of standard - cumulants from those of set average deviation sqrt(n-1) a Gaussian distribition cum. 3 cum. 4 SS1 * 7.150740e+00 * 8.803173e-01 1.760635e-02 0.0620.163 SS2 1.490604e+00 1.164761e+00 2.329523e-02 0.4950.153 https://www.dropbox.com/s/1vqixenyerha7qq/115-water.png Here's the link hbond.xvg file and its averaged file https://www.dropbox.com/s/4n0m47o3mrjn3o8/hbond_115-water.xvg https://www.dropbox.com/s/4n0m47o3mrjn3o8/hbond_115-water.xvg Neither of these files produce output that corresponds to the PNG image above. Both files have values in 6-9 H-bond range and thus agree with the g_analyze output, which I can reproduce. I suspect you're somehow getting your files mixed up. -Justin On Mon, Nov 11, 2013 at 3:30 PM, bharat gupta bharat.85.m...@gmail.com wrote: thank you informing about g_rdf... Is it possible to dump the structure with those average water molecules interacting with the residues. I generated the hbond.log file which gives the details but I need to generate a figure for this ?? On Mon, Nov 11, 2013 at 10:40 AM, Justin Lemkul jalem...@vt.edu wrote: On 11/10/13 8:38 PM, bharat gupta wrote: But trjorder can be used to calculate the hydration layer or shell around residues ... Right ?? Yes, but I also tend to think that integrating an RDF is also a more straightforward way of doing that. With trjorder, you set some arbitrary cutoff that may or may not be an informed decision - with an RDF it is clear where the hydration layers are. -Justin On Mon, Nov 11, 2013 at 10:35 AM, Justin Lemkul jalem...@vt.edu wrote: On 11/10/13 8:30 PM, bharat gupta wrote: Thanks for your reply. I was missing the scientific notation part. Now everything is fine. Regarding trjorder, it doesn't measure h-bonds but gives the water nearest to protein. I wouldn't try to draw any sort of comparison between the output of trjorder and g_hbond. If you want to measure H-bonds, there's only one tool for that. -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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 -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room
Re: [gmx-users] Re: g_analyze
Thanks justin for your replies. I understood the g_analyze related data. I tired g_analyze to dump the structures as you said. But, I didn't find any switch that can be used to dump the structure in pdb format. On Tue, Nov 12, 2013 at 10:15 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/11/13 6:56 PM, bharat gupta wrote: Hi, I tried g_select to dump the structure with the interacting water molecules, but I don't know know how to do that. I searched for some threads in the discussion but wasn't able to find anything related to my need. Can you explain how can I do that ? Start with g_select -select 'help all' and see what you can determine. Such selections are rather straightforward and have been explained several times on the list. If you need help, show us what you're doing and describe why it isn't what you want. It will ultimately save a lot of time. -Justin On Tue, Nov 12, 2013 at 7:39 AM, bharat gupta bharat.85.m...@gmail.com wrote: Sorry, I attached the wrong file . Here's the average file generate from one of the files I sent in my last mail. I used the command g_analyze -f hbond_115-water.xvg -av hbond_115-water-avg.xvg. Here's the file obtained from this command :- https://www.dropbox.com/s/sovzk40cudznfjw/hbond_115-water-avg.xvg Now, if you see, the graph (in previous mail) and average file, both correlates well. I have a doubt about interpreting the result from g_analyze. The value 7.150740e+00 implies that on average 7 hydrogen bonds are formed during the simulation time of 5ns to 10ns. What does then the average file or its graph tells ?? On Mon, Nov 11, 2013 at 9:58 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/11/13 4:06 AM, bharat gupta wrote: In addition to my previous question, I have another question about g_analyze. When I used the hbond.xvg file to get the average and plotted the average.xvg file I found that the average value is round 4 to 5 according to the graph. But g_analyze in its final calculation gives 7.150 as the average values... Here's the link for the graph and result of average value calculated by g_analyze :- std. dev.relative deviation of standard - cumulants from those of set average deviation sqrt(n-1) a Gaussian distribition cum. 3 cum. 4 SS1 * 7.150740e+00 * 8.803173e-01 1.760635e-02 0.0620.163 SS2 1.490604e+00 1.164761e+00 2.329523e-02 0.4950.153 https://www.dropbox.com/s/1vqixenyerha7qq/115-water.png Here's the link hbond.xvg file and its averaged file https://www.dropbox.com/s/4n0m47o3mrjn3o8/hbond_115-water.xvg https://www.dropbox.com/s/4n0m47o3mrjn3o8/hbond_115-water.xvg Neither of these files produce output that corresponds to the PNG image above. Both files have values in 6-9 H-bond range and thus agree with the g_analyze output, which I can reproduce. I suspect you're somehow getting your files mixed up. -Justin On Mon, Nov 11, 2013 at 3:30 PM, bharat gupta bharat.85.m...@gmail.com wrote: thank you informing about g_rdf... Is it possible to dump the structure with those average water molecules interacting with the residues. I generated the hbond.log file which gives the details but I need to generate a figure for this ?? On Mon, Nov 11, 2013 at 10:40 AM, Justin Lemkul jalem...@vt.edu wrote: On 11/10/13 8:38 PM, bharat gupta wrote: But trjorder can be used to calculate the hydration layer or shell around residues ... Right ?? Yes, but I also tend to think that integrating an RDF is also a more straightforward way of doing that. With trjorder, you set some arbitrary cutoff that may or may not be an informed decision - with an RDF it is clear where the hydration layers are. -Justin On Mon, Nov 11, 2013 at 10:35 AM, Justin Lemkul jalem...@vt.edu wrote: On 11/10/13 8:30 PM, bharat gupta wrote: Thanks for your reply. I was missing the scientific notation part. Now everything is fine. Regarding trjorder, it doesn't measure h-bonds but gives the water nearest to protein. I wouldn't try to draw any sort of comparison between the output of trjorder and g_hbond. If you want to measure H-bonds, there's only one tool for that. -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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
Re: [gmx-users] Re: g_analyze
On 11/12/13 8:33 AM, bharat gupta wrote: Thanks justin for your replies. I understood the g_analyze related data. I tired g_analyze to dump the structures as you said. But, I didn't find any switch that can be used to dump the structure in pdb format. Because that's not the function of g_analyze. Use trjconv -dump with a suitable index file (from g_select). -Justin On Tue, Nov 12, 2013 at 10:15 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/11/13 6:56 PM, bharat gupta wrote: Hi, I tried g_select to dump the structure with the interacting water molecules, but I don't know know how to do that. I searched for some threads in the discussion but wasn't able to find anything related to my need. Can you explain how can I do that ? Start with g_select -select 'help all' and see what you can determine. Such selections are rather straightforward and have been explained several times on the list. If you need help, show us what you're doing and describe why it isn't what you want. It will ultimately save a lot of time. -Justin On Tue, Nov 12, 2013 at 7:39 AM, bharat gupta bharat.85.m...@gmail.com wrote: Sorry, I attached the wrong file . Here's the average file generate from one of the files I sent in my last mail. I used the command g_analyze -f hbond_115-water.xvg -av hbond_115-water-avg.xvg. Here's the file obtained from this command :- https://www.dropbox.com/s/sovzk40cudznfjw/hbond_115-water-avg.xvg Now, if you see, the graph (in previous mail) and average file, both correlates well. I have a doubt about interpreting the result from g_analyze. The value 7.150740e+00 implies that on average 7 hydrogen bonds are formed during the simulation time of 5ns to 10ns. What does then the average file or its graph tells ?? On Mon, Nov 11, 2013 at 9:58 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/11/13 4:06 AM, bharat gupta wrote: In addition to my previous question, I have another question about g_analyze. When I used the hbond.xvg file to get the average and plotted the average.xvg file I found that the average value is round 4 to 5 according to the graph. But g_analyze in its final calculation gives 7.150 as the average values... Here's the link for the graph and result of average value calculated by g_analyze :- std. dev.relative deviation of standard - cumulants from those of set average deviation sqrt(n-1) a Gaussian distribition cum. 3 cum. 4 SS1 * 7.150740e+00 * 8.803173e-01 1.760635e-02 0.0620.163 SS2 1.490604e+00 1.164761e+00 2.329523e-02 0.4950.153 https://www.dropbox.com/s/1vqixenyerha7qq/115-water.png Here's the link hbond.xvg file and its averaged file https://www.dropbox.com/s/4n0m47o3mrjn3o8/hbond_115-water.xvg https://www.dropbox.com/s/4n0m47o3mrjn3o8/hbond_115-water.xvg Neither of these files produce output that corresponds to the PNG image above. Both files have values in 6-9 H-bond range and thus agree with the g_analyze output, which I can reproduce. I suspect you're somehow getting your files mixed up. -Justin On Mon, Nov 11, 2013 at 3:30 PM, bharat gupta bharat.85.m...@gmail.com wrote: thank you informing about g_rdf... Is it possible to dump the structure with those average water molecules interacting with the residues. I generated the hbond.log file which gives the details but I need to generate a figure for this ?? On Mon, Nov 11, 2013 at 10:40 AM, Justin Lemkul jalem...@vt.edu wrote: On 11/10/13 8:38 PM, bharat gupta wrote: But trjorder can be used to calculate the hydration layer or shell around residues ... Right ?? Yes, but I also tend to think that integrating an RDF is also a more straightforward way of doing that. With trjorder, you set some arbitrary cutoff that may or may not be an informed decision - with an RDF it is clear where the hydration layers are. -Justin On Mon, Nov 11, 2013 at 10:35 AM, Justin Lemkul jalem...@vt.edu wrote: On 11/10/13 8:30 PM, bharat gupta wrote: Thanks for your reply. I was missing the scientific notation part. Now everything is fine. Regarding trjorder, it doesn't measure h-bonds but gives the water nearest to protein. I wouldn't try to draw any sort of comparison between the output of trjorder and g_hbond. If you want to measure H-bonds, there's only one tool for that. -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users *
Re: [gmx-users] Calculating diffusion coefficient in three dimension
Thank you sir for the prompt reply. *g_msd -f traj.xtc -s topol.tpr -n msd.ndx -lateral z -o msd.xvg -tu ns* Here I am giving -lateral z (like for membrane simulations). Is it fine for spherical systems also? On Tue, Nov 12, 2013 at 5:32 PM, Dr. Vitaly Chaban vvcha...@gmail.comwrote: MSD is 3D by default. Dr. Vitaly V. Chaban On Tue, Nov 12, 2013 at 6:01 AM, Venkat Reddy venkat...@gmail.com wrote: Dear all, I am simulating a spherical lipid vesicle. I want to calculate the diffusion coefficient for each lipid component in 3D. How to calculate it using g_msd (or any other tool like g_velacc)? Thank you for your concern -- With Best Wishes Venkat Reddy Chirasani -- gmx-users mailing listgmx-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 -- gmx-users mailing listgmx-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 -- With Best Wishes Venkat Reddy Chirasani PhD student Laboratory of Computational Biophysics Department of Biotechnology IIT Madras Chennai INDIA-600036 -- gmx-users mailing listgmx-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
Re: [gmx-users] Calculating diffusion coefficient in three dimension
On 11/12/13 8:55 AM, Venkat Reddy wrote: Thank you sir for the prompt reply. *g_msd -f traj.xtc -s topol.tpr -n msd.ndx -lateral z -o msd.xvg -tu ns* Here I am giving -lateral z (like for membrane simulations). Is it fine for spherical systems also? No. The system is a sphere, so what use is it to calculate motion perpendicular to z when you have lipids moving in all three spatial dimensions? A vesicle is very different from a membrane, in which the lipids move in a plane, thus making -lateral z useful. -Justin On Tue, Nov 12, 2013 at 5:32 PM, Dr. Vitaly Chaban vvcha...@gmail.comwrote: MSD is 3D by default. Dr. Vitaly V. Chaban On Tue, Nov 12, 2013 at 6:01 AM, Venkat Reddy venkat...@gmail.com wrote: Dear all, I am simulating a spherical lipid vesicle. I want to calculate the diffusion coefficient for each lipid component in 3D. How to calculate it using g_msd (or any other tool like g_velacc)? Thank you for your concern -- With Best Wishes Venkat Reddy Chirasani -- gmx-users mailing listgmx-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 -- gmx-users mailing listgmx-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 -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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
[gmx-users] Re: ok, thank you
Send gmx-users mailing list submissions to gmx-users@gromacs.org To subscribe or unsubscribe via the World Wide Web, visit http://lists.gromacs.org/mailman/listinfo/gmx-users or, via email, send a message with subject or body 'help' to gmx-users-requ...@gromacs.org You can reach the person managing the list at gmx-users-ow...@gromacs.org When replying, please edit your Subject line so it is more specific than Re: Contents of gmx-users digest... Today's Topics: 1. Restarting a simulation after replacing an empty md.trr file (arun kumar) 2. Re: Re: Reaction field zero and ions (Justin Lemkul) 3. Re: Calculating diffusion coefficient in three dimension (Dr. Vitaly Chaban) 4. Re: Re: Reaction field zero and ions (Dr. Vitaly Chaban) 5. Re: Re: g_analyze (Justin Lemkul) 6. Re: Re: g_analyze (Justin Lemkul) -- Message: 1 Date: Tue, 12 Nov 2013 11:45:05 +0530 From: arun kumar arunjones.kuma...@gmail.com Subject: [gmx-users] Restarting a simulation after replacing an empty md.trr file To: gmx-users@gromacs.org Message-ID: cagm9vj8dn+fqb-cigjdvv+i1mq7mc2cxvkde3gsp5+lwhds...@mail.gmail.com Content-Type: text/plain; charset=ISO-8859-1 Dear Gromacs users, I am running a 50ns simulation of a protein having nearly 700 residues on 60 threads (Gromacs 4.6.3). At one point i got a disk space problem, so i have deleted the md.trr file and created an empty md.trr file. when i tried to restart the simulation from check point file on 100 threads, [ mdrun -v -deffnm md -cpi md.cpt -nt 100 ] i am getting a note and an error as fallows Reading checkpoint file md.cpt generated: #PME-nodes mismatch, current program: 100 checkpoint file: 60 Gromacs binary or parallel settings not identical to previous run. Continuation is exact, but is not guaranteed to be binary identical. ... Source code file: checkpoint.c, line: 1767 Fatal error: Can't read 1048576 bytes of 'md.trr' to compute checksum. The file has been replaced or its contents has been modified. please help me in overcoming this problem. Thanking you. -- Arun Kumar Somavarapu Project-JRF Dr. Pawan Gupta's lab Protein Science and Engineering Dept, Institute of Microbial Tecnology, Sec 39-A, Chandigarh - 160036. -- Message: 2 Date: Mon, 11 Nov 2013 13:06:32 -0500 From: Justin Lemkul jalem...@vt.edu Subject: Re: [gmx-users] Re: Reaction field zero and ions To: Discussion list for GROMACS users gmx-users@gromacs.org Message-ID: 52811ca8.5030...@vt.edu Content-Type: text/plain; charset=ISO-8859-1; format=flowed On 11/11/13 12:08 PM, Williams Ernesto Miranda Delgado wrote: Hello If I did the MD simulation using PME and neutralized with ions, and I want to rerun this time with reaction field zero, is there any problem if I keep the ions? This is for LIE calculation. I am using AMBER99SB. Why do you think it necessary to delete them? -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- Message: 3 Date: Tue, 12 Nov 2013 13:02:54 +0100 From: Dr. Vitaly Chaban vvcha...@gmail.com Subject: Re: [gmx-users] Calculating diffusion coefficient in three dimension To: Discussion list for GROMACS users gmx-users@gromacs.org Message-ID: capxdd+abay6mj_dkn6_k+mkbuty4eyzspdvdxj+m-m-2zae...@mail.gmail.com Content-Type: text/plain; charset=ISO-8859-1 MSD is 3D by default. Dr. Vitaly V. Chaban On Tue, Nov 12, 2013 at 6:01 AM, Venkat Reddy venkat...@gmail.com wrote: Dear all, I am simulating a spherical lipid vesicle. I want to calculate the diffusion coefficient for each lipid component in 3D. How to calculate it using g_msd (or any other tool like g_velacc)? Thank you for your concern -- With Best Wishes Venkat Reddy Chirasani -- gmx-users mailing listgmx-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 -- Message: 4 Date: Tue, 12 Nov 2013 13:04:10 +0100 From: Dr. Vitaly Chaban vvcha...@gmail.com Subject: Re: [gmx-users] Re: Reaction field zero and ions To: Discussion list for GROMACS users gmx-users@gromacs.org Message-ID:
[gmx-users] Invalid order for directive defaults
Dear all My system contains protein + cnt + water molecules. I have summarized what I did below: --- 1) By pdb2gmx and charmm27 force field, I obtained pr.top for protein then I converted it into pr.itp by deleting ; Include forcefield parameters #include charmm27.ff/forcefield.itp from begining of file and by deleting ; Include water topology #include charmm27.ff/tip3p.itp #ifdef POSRES_WATER ; Position restraint for each water oxygen [ position_restraints ] ; i funct fcxfcyfcz 11 1000 1000 1000 #endif ; Include topology for ions #include charmm27.ff/ions.itp [ system ] ; Name Protein [ molecules ] ; Compound#mols Protein_chain_A 1 from ending of file. 2) I used bonded and nonbonded parameters for cnt from paper: J. Phys. Chem. B 2001, 105, 9980-9987 (Carbon Nanotubes in Water: Structural Characteristics and Energetics). I created cnt.ff folder containing following files: ffcnt.atp / ffcnt.n2t / ffcnt.rtp / ffcntbon.itp / ffcntnonbon.itp / forcefield.itp Then, I put cnt.ff folder in GMXLIB directory. 3) By g_x2top and cnt.ff created in previous step, I obtained cnt.top for cnt then I converted it into cnt.itp by deleting ; Include forcefield parameters #include cnt.ff/forcefield.itp from begining of file and by deleting [ system ] ; Name CNT [ molecules ] ; Compound#mols CNT 1 from ending of file. 4) I combined cnt and protein to create one coordination file (system.gro). Order of molecules are as follows: 1) cnt 2) protein 3) water. 5) I wrote a topol.top file given below, ; Include forcefield parameters #include cnt.ff/forcefield.itp #include cnt.itp #include charmm27.ff/forcefield.itp #include pr.itp ; Include Position restraint file #ifdef POSRES #include posre.itp #endif ; Include water topology #include charmm27.ff/tip3p.itp #ifdef POSRES_WATER ; Position restraint for each water oxygen [ position_restraints ] ; i funct fcxfcyfcz 11 1000 1000 1000 #endif ; Include topology for ions #include charmm27.ff/ions.itp [ system ] ; Name CNT/Protein/SOL [ molecules ] ; Compound#mols CNT 1 Protein 1 SOL 1359 When I used grompp -f minim.mdp -c system.gro -p topol.top -o minim.tpr, I encountered with following error: Fatal error: Syntax error - File forcefield.itp, line 11 Last line read: '[ defaults ]' Invalid order for directive defaults Content of forcefield.itp in cnt.ff directory is as follows: *** *CHARMM port writted by * *Par Bjelkmar, Per Larsson, Michel Cuendet, * *Berk Hess and Erik Lindahl. * * Correspondance: * *bjelk...@cbr.su.se or lind...@cbr.su.se * *** #define _FF_CNT [ defaults ] ; nbfunccomb-rulegen-pairsfudgeLJfudgeQQ 12yes1.01.0 #include ffcntnonbon.itp #include ffcntbon.itp and Content of forcefield.itp in Charmm27.ff directory is as follows: *** *CHARMM port writted by * *Par Bjelkmar, Per Larsson, Michel Cuendet, * *Berk Hess and Erik Lindahl. * * Correspondance: * *bjelk...@cbr.su.se or lind...@cbr.su.se * *** #define _FF_CHARMM [ defaults ] ; nbfunccomb-rulegen-pairsfudgeLJfudgeQQ 12yes1.01.0 #include ffnonbonded.itp #include ffbonded.itp #include gb.itp #include cmap.itp ; Nucleic acids nonbonded and bonded parameters #include ffnanonbonded.itp #include ffnabonded.itp In both of forcefield.itp files, line 11 is [ defaults ]. I changed name of forcefield.itp file in cnt.ff directory to cntff.itp, but there is same error, again. That is all what I exactly did. Is anything wrong or missing? How to solve this error? Any help will highly appreciated --
Re: mdrun on 8-core AMD + GTX TITAN (was: Re: [gmx-users] Re: Gromacs-4.6 on two Titans GPUs)
As Mark said, please share the *entire* log file. Among other important things, the result of PP-PME tuning is not included above. However, I suspect that in this case scaling is strongly affected or by the small size of the system you are simulating. -- Szilárd On Sun, Nov 10, 2013 at 5:28 AM, Dwey Kauffman mpi...@gmail.com wrote: Hi Szilard, Thank you very much for your suggestions. Actually, I was jumping to conclusions too early, as you mentioned AMD cluster, I assumed you must have 12-16-core Opteron CPUs. If you have an 8-core (desktop?) AMD CPU, than you may not need to run more than one rank per GPU. Yes, we do have independent clusters of AMD, AMD opteron, Intel Corei7. All nodes of three clusters are installed with (at least) 1 GPU card. I have run the same test on these three clusters. Let's focus on a basic scaling issue: One GPU v.s Two GPUs within the same node of 8-core AMD cpu. Using 1 GPU, we can have a performance of ~32 ns/day. Using two GPU, we gain not much more ( ~38.5 ns/day ). It is about ~20% more performance. However, this is not really true because in some tests, I also saw only 2-5% more, which really surprised me. As you can see, this test was made on the same node regardless of networking. Can the performance be improved say 50% more when 2 GPUs are used on a general task ? If yes, how ? Indeed, as Richard pointed out, I was asking for *full* logs, these summaries can't tell much, the table above the summary entitled R E A L C Y C L E A N D T I M E A C C O U N T I N G as well as other reported information across the log file is what I need to make an assessment of your simulations' performance. Please see below. However, in your case I suspect that the bottleneck is multi-threaded scaling on the AMD CPUs and you should probably decrease the number of threads per MPI rank and share GPUs between 2-4 ranks. After I test all three clusters, I found it may NOT be an issue of AMD cpus. Intel cpus has the SAME scaling issue. However, I am curious as to how you justify the setup of 2-4 ranks sharing GPUs ? Can you please explain it a bit more ? You could try running mpirun -np 4 mdrun -ntomp 2 -gpu_id 0011 but I suspect this won't help because your scaling issue Your guess is correct but why is that ? it is worse. The more nodes are involved in a task, the performance is worse. in my experience even reaction field runs don't scale across nodes with 10G ethernet if you have more than 4-6 ranks per node trying to communicate (let alone with PME). What dose it mean let alone with PME ? how to do so ? by mdrun ? I do know mdrun -npme to specify PME process. Thank you. Dwey ### One GPU R E A L C Y C L E A N D T I M E A C C O U N T I N G Computing: Nodes Th. Count Wall t (s) G-Cycles % - Neighbor search18 11 431.81713863.390 1.6 Launch GPU ops.18501 472.90615182.556 1.7 Force 185011328.61142654.785 4.9 PME mesh 18501 11561.327 371174.09042.8 Wait GPU local 185016888.008 221138.11125.5 NB X/F buffer ops. 189911216.49939055.455 4.5 Write traj.18 1030 12.741 409.039 0.0 Update 185011696.35854461.226 6.3 Constraints185011969.72663237.647 7.3 Rest 11458.82046835.133 5.4 - Total 1 27036.812 868011.431 100.0 - - PME spread/gather 18 10026975.086 223933.73925.8 PME 3D-FFT 18 10023928.259 126115.97614.5 PME solve 18501 636.48820434.327 2.4 - GPU timings - Computing: Count Wall t (s) ms/step % - Pair list H2D 11 43.4350.434 0.2 X / q H2D501 567.1680.113 2.8 Nonbonded F kernel 400 14174.3163.54470.8 Nonbonded F+ene k.904314.4384.79421.5 Nonbonded F+ene+prune k. 11 572.3705.724 2.9 F D2H
Re: [gmx-users] Invalid order for directive defaults
On 11/12/13 10:07 AM, Atila Petrosian wrote: Dear all My system contains protein + cnt + water molecules. I have summarized what I did below: --- 1) By pdb2gmx and charmm27 force field, I obtained pr.top for protein then I converted it into pr.itp by deleting ; Include forcefield parameters #include charmm27.ff/forcefield.itp from begining of file and by deleting ; Include water topology #include charmm27.ff/tip3p.itp #ifdef POSRES_WATER ; Position restraint for each water oxygen [ position_restraints ] ; i funct fcxfcyfcz 11 1000 1000 1000 #endif ; Include topology for ions #include charmm27.ff/ions.itp [ system ] ; Name Protein [ molecules ] ; Compound#mols Protein_chain_A 1 from ending of file. 2) I used bonded and nonbonded parameters for cnt from paper: J. Phys. Chem. B 2001, 105, 9980-9987 (Carbon Nanotubes in Water: Structural Characteristics and Energetics). I created cnt.ff folder containing following files: ffcnt.atp / ffcnt.n2t / ffcnt.rtp / ffcntbon.itp / ffcntnonbon.itp / forcefield.itp Then, I put cnt.ff folder in GMXLIB directory. 3) By g_x2top and cnt.ff created in previous step, I obtained cnt.top for cnt then I converted it into cnt.itp by deleting ; Include forcefield parameters #include cnt.ff/forcefield.itp from begining of file and by deleting [ system ] ; Name CNT [ molecules ] ; Compound#mols CNT 1 from ending of file. 4) I combined cnt and protein to create one coordination file (system.gro). Order of molecules are as follows: 1) cnt 2) protein 3) water. 5) I wrote a topol.top file given below, ; Include forcefield parameters #include cnt.ff/forcefield.itp #include cnt.itp #include charmm27.ff/forcefield.itp #include pr.itp ; Include Position restraint file #ifdef POSRES #include posre.itp #endif ; Include water topology #include charmm27.ff/tip3p.itp #ifdef POSRES_WATER ; Position restraint for each water oxygen [ position_restraints ] ; i funct fcxfcyfcz 11 1000 1000 1000 #endif ; Include topology for ions #include charmm27.ff/ions.itp [ system ] ; Name CNT/Protein/SOL [ molecules ] ; Compound#mols CNT 1 Protein 1 SOL 1359 When I used grompp -f minim.mdp -c system.gro -p topol.top -o minim.tpr, I encountered with following error: Fatal error: Syntax error - File forcefield.itp, line 11 Last line read: '[ defaults ]' Invalid order for directive defaults Content of forcefield.itp in cnt.ff directory is as follows: *** *CHARMM port writted by * *Par Bjelkmar, Per Larsson, Michel Cuendet, * *Berk Hess and Erik Lindahl. * * Correspondance: * *bjelk...@cbr.su.se or lind...@cbr.su.se * *** #define _FF_CNT [ defaults ] ; nbfunccomb-rulegen-pairsfudgeLJfudgeQQ 12yes1.01.0 #include ffcntnonbon.itp #include ffcntbon.itp and Content of forcefield.itp in Charmm27.ff directory is as follows: *** *CHARMM port writted by * *Par Bjelkmar, Per Larsson, Michel Cuendet, * *Berk Hess and Erik Lindahl. * * Correspondance: * *bjelk...@cbr.su.se or lind...@cbr.su.se * *** #define _FF_CHARMM [ defaults ] ; nbfunccomb-rulegen-pairsfudgeLJfudgeQQ 12yes1.01.0 #include ffnonbonded.itp #include ffbonded.itp #include gb.itp #include cmap.itp ; Nucleic acids nonbonded and bonded parameters #include ffnanonbonded.itp #include ffnabonded.itp In both of forcefield.itp files, line 11 is [ defaults ]. I changed name of forcefield.itp file in cnt.ff directory to cntff.itp, but there is same error, again. That is all what I exactly did. Is anything wrong or missing? How to
[gmx-users] Re: segmentation fault on gromacs 4.5.5 after mdrun
I run grompp -f nvt.mdp -c em.gro -p topol.top -n index.ndx -o nvt.tpr and everything looks fine. I check the nvt.tpr, and temperature is ok. the real problem is with the mdrun function. could be a problem of the software? Thanks Javier Justin Lemkul wrote On 11/11/13 11:24 AM, Carlos Javier Almeciga Diaz wrote: Hello evryone, I doing a simulation of a ligand-protein interaction with gromacs 4.5.5. Everything looks fine after I equilibrate the protein-ligand complex. I'm running these commands: grompp -f nvt.mdp -c em.gro -p topol.top -n index.ndx -o nvt.tpr mdrun -deffnm nvt Nevertheless, I got this error: Reading file nvt.tpr, VERSION 4.5.5 (double precision) Segmentation fault What should I do? Instantaneous failure typically indicates that the forces are nonsensically high and the constraint algorithm immediately fails. Likely the previous energy minimization did not adequately complete. -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalemkul@.umaryland | (410) 706-7441 == -- gmx-users mailing list gmx-users@ 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-request@ . * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists -- View this message in context: http://gromacs.5086.x6.nabble.com/segmentation-fault-on-gromacs-4-5-5-after-mdrun-tp5012431p5012458.html Sent from the GROMACS Users Forum mailing list archive at Nabble.com. -- gmx-users mailing listgmx-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
[gmx-users] Re: Restarting a simulation after replacing an empty md.trr file
Hello Sir, Thanks for the reply. now is it fine if i use 100 threads in my restart. is there any impact on the over all simulation? -- View this message in context: http://gromacs.5086.x6.nabble.com/Restarting-a-simulation-after-replacing-an-empty-md-trr-file-tp5012443p5012459.html Sent from the GROMACS Users Forum mailing list archive at Nabble.com. -- gmx-users mailing listgmx-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
Re: [gmx-users] Re: segmentation fault on gromacs 4.5.5 after mdrun
On 11/12/13 10:58 AM, cjalmeciga wrote: I run grompp -f nvt.mdp -c em.gro -p topol.top -n index.ndx -o nvt.tpr and everything looks fine. I check the nvt.tpr, and temperature is ok. The fact that grompp completes indicates there is nothing syntactically wrong with the input files. Whether or not the content of the .mdp is physically sensible or the input configuration is plausible is an entirely different matter. Please tell us what the exact outcome of the previous energy minimization was (potential energy, maximum force, copied and pasted from screen output or .log file). the real problem is with the mdrun function. could be a problem of the software? You have presented no evidence that would lead anyone to believe the problem is with mdrun. In the vast majority of cases, user input is the problem. -Justin Thanks Javier Justin Lemkul wrote On 11/11/13 11:24 AM, Carlos Javier Almeciga Diaz wrote: Hello evryone, I doing a simulation of a ligand-protein interaction with gromacs 4.5.5. Everything looks fine after I equilibrate the protein-ligand complex. I'm running these commands: grompp -f nvt.mdp -c em.gro -p topol.top -n index.ndx -o nvt.tpr mdrun -deffnm nvt Nevertheless, I got this error: Reading file nvt.tpr, VERSION 4.5.5 (double precision) Segmentation fault What should I do? Instantaneous failure typically indicates that the forces are nonsensically high and the constraint algorithm immediately fails. Likely the previous energy minimization did not adequately complete. -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalemkul@.umaryland | (410) 706-7441 == -- gmx-users mailing list gmx-users@ 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-request@ . * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists -- View this message in context: http://gromacs.5086.x6.nabble.com/segmentation-fault-on-gromacs-4-5-5-after-mdrun-tp5012431p5012458.html Sent from the GROMACS Users Forum mailing list archive at Nabble.com. -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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
Re: [gmx-users] Re: Restarting a simulation after replacing an empty md.trr file
On 11/12/13 11:10 AM, arunjones wrote: Hello Sir, Thanks for the reply. now is it fine if i use 100 threads in my restart. is there any impact on the over all simulation? Only if that is the number of threads originally used in the run. If not, there will be a mismatch between the DD grid setup, the .cpt file will complain, and the run will fail. Rule of thumb: don't change settings or alter files mid-run ;) -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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
Re: [gmx-users] Calculating diffusion coefficient in three dimension
Then, how to mention the direction for spherical particles Sir? On Tue, Nov 12, 2013 at 7:28 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/12/13 8:55 AM, Venkat Reddy wrote: Thank you sir for the prompt reply. *g_msd -f traj.xtc -s topol.tpr -n msd.ndx -lateral z -o msd.xvg -tu ns* Here I am giving -lateral z (like for membrane simulations). Is it fine for spherical systems also? No. The system is a sphere, so what use is it to calculate motion perpendicular to z when you have lipids moving in all three spatial dimensions? A vesicle is very different from a membrane, in which the lipids move in a plane, thus making -lateral z useful. -Justin On Tue, Nov 12, 2013 at 5:32 PM, Dr. Vitaly Chaban vvcha...@gmail.com wrote: MSD is 3D by default. Dr. Vitaly V. Chaban On Tue, Nov 12, 2013 at 6:01 AM, Venkat Reddy venkat...@gmail.com wrote: Dear all, I am simulating a spherical lipid vesicle. I want to calculate the diffusion coefficient for each lipid component in 3D. How to calculate it using g_msd (or any other tool like g_velacc)? Thank you for your concern -- With Best Wishes Venkat Reddy Chirasani -- gmx-users mailing listgmx-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 -- gmx-users mailing listgmx-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 -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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 -- With Best Wishes Venkat Reddy Chirasani PhD student Laboratory of Computational Biophysics Department of Biotechnology IIT Madras Chennai INDIA-600036 -- gmx-users mailing listgmx-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
Re: [gmx-users] Calculating diffusion coefficient in three dimension
On 11/12/13 11:25 AM, Venkat Reddy wrote: Then, how to mention the direction for spherical particles Sir? Read g_msd -h again, paying specific attention to the -type flag. -Justin On Tue, Nov 12, 2013 at 7:28 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/12/13 8:55 AM, Venkat Reddy wrote: Thank you sir for the prompt reply. *g_msd -f traj.xtc -s topol.tpr -n msd.ndx -lateral z -o msd.xvg -tu ns* Here I am giving -lateral z (like for membrane simulations). Is it fine for spherical systems also? No. The system is a sphere, so what use is it to calculate motion perpendicular to z when you have lipids moving in all three spatial dimensions? A vesicle is very different from a membrane, in which the lipids move in a plane, thus making -lateral z useful. -Justin On Tue, Nov 12, 2013 at 5:32 PM, Dr. Vitaly Chaban vvcha...@gmail.com wrote: MSD is 3D by default. Dr. Vitaly V. Chaban On Tue, Nov 12, 2013 at 6:01 AM, Venkat Reddy venkat...@gmail.com wrote: Dear all, I am simulating a spherical lipid vesicle. I want to calculate the diffusion coefficient for each lipid component in 3D. How to calculate it using g_msd (or any other tool like g_velacc)? Thank you for your concern -- With Best Wishes Venkat Reddy Chirasani -- gmx-users mailing listgmx-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 -- gmx-users mailing listgmx-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 -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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 -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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
[gmx-users] Re: Restarting a simulation after replacing an empty md.trr file
Thank you Sir, initially i was running on 60 threads, now i changed it to 100. simulation is running with out any error, but i found a note in the log file as fallows #nodes mismatch, current program: 100 checkpoint file: 60 #PME-nodes mismatch, current program: -1 checkpoint file: 12 Gromacs binary or parallel settings not identical to previous run. Continuation is exact, but is not guaranteed to be binary identical. Initializing Domain Decomposition on 100 nodes is it a good idea to continue or shall i stick with 60 threads. -- View this message in context: http://gromacs.5086.x6.nabble.com/Restarting-a-simulation-after-replacing-an-empty-md-trr-file-tp5012443p5012463.html Sent from the GROMACS Users Forum mailing list archive at Nabble.com. -- gmx-users mailing listgmx-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
Re: [gmx-users] Re: Restarting a simulation after replacing an empty md.trr file
On 11/12/13 12:07 PM, arunjones wrote: Thank you Sir, initially i was running on 60 threads, now i changed it to 100. simulation is running with out any error, but i found a note in the log file as fallows #nodes mismatch, current program: 100 checkpoint file: 60 #PME-nodes mismatch, current program: -1 checkpoint file: 12 Gromacs binary or parallel settings not identical to previous run. Continuation is exact, but is not guaranteed to be binary identical. Initializing Domain Decomposition on 100 nodes is it a good idea to continue or shall i stick with 60 threads. Like I said, I think it is a bad idea to switch settings haphazardly during the run. As the note indicates, the continuation is exact, but not binary identical. Check the website for what all that means if you're not sure. -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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
[gmx-users] Re: segmentation fault on gromacs 4.5.5 after mdrun
The output of energy minimization was Potential Energy = -1.42173622068236e+06 Maximum force = 9.00312066109319e+02 on atom 148 Norm of force = 2.06087515037187e+01 Thanks Javier Justin Lemkul wrote On 11/12/13 10:58 AM, cjalmeciga wrote: I run grompp -f nvt.mdp -c em.gro -p topol.top -n index.ndx -o nvt.tpr and everything looks fine. I check the nvt.tpr, and temperature is ok. The fact that grompp completes indicates there is nothing syntactically wrong with the input files. Whether or not the content of the .mdp is physically sensible or the input configuration is plausible is an entirely different matter. Please tell us what the exact outcome of the previous energy minimization was (potential energy, maximum force, copied and pasted from screen output or .log file). the real problem is with the mdrun function. could be a problem of the software? You have presented no evidence that would lead anyone to believe the problem is with mdrun. In the vast majority of cases, user input is the problem. -Justin Thanks Javier Justin Lemkul wrote On 11/11/13 11:24 AM, Carlos Javier Almeciga Diaz wrote: Hello evryone, I doing a simulation of a ligand-protein interaction with gromacs 4.5.5. Everything looks fine after I equilibrate the protein-ligand complex. I'm running these commands: grompp -f nvt.mdp -c em.gro -p topol.top -n index.ndx -o nvt.tpr mdrun -deffnm nvt Nevertheless, I got this error: Reading file nvt.tpr, VERSION 4.5.5 (double precision) Segmentation fault What should I do? Instantaneous failure typically indicates that the forces are nonsensically high and the constraint algorithm immediately fails. Likely the previous energy minimization did not adequately complete. -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalemkul@.umaryland | (410) 706-7441 == -- gmx-users mailing list gmx-users@ 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-request@ . * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists -- View this message in context: http://gromacs.5086.x6.nabble.com/segmentation-fault-on-gromacs-4-5-5-after-mdrun-tp5012431p5012458.html Sent from the GROMACS Users Forum mailing list archive at Nabble.com. -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalemkul@.umaryland | (410) 706-7441 == -- gmx-users mailing list gmx-users@ 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-request@ . * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists -- View this message in context: http://gromacs.5086.x6.nabble.com/segmentation-fault-on-gromacs-4-5-5-after-mdrun-tp5012431p5012464.html Sent from the GROMACS Users Forum mailing list archive at Nabble.com. -- gmx-users mailing listgmx-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
Re: [gmx-users] Re: segmentation fault on gromacs 4.5.5 after mdrun
On 11/12/13 12:14 PM, cjalmeciga wrote: The output of energy minimization was Potential Energy = -1.42173622068236e+06 Maximum force = 9.00312066109319e+02 on atom 148 Norm of force = 2.06087515037187e+01 OK, reasonable enough. How about a description of what the system is, which force field you chose, how you derived the ligand topology, and the full contents of your .mdp file? -Justin Thanks Javier Justin Lemkul wrote On 11/12/13 10:58 AM, cjalmeciga wrote: I run grompp -f nvt.mdp -c em.gro -p topol.top -n index.ndx -o nvt.tpr and everything looks fine. I check the nvt.tpr, and temperature is ok. The fact that grompp completes indicates there is nothing syntactically wrong with the input files. Whether or not the content of the .mdp is physically sensible or the input configuration is plausible is an entirely different matter. Please tell us what the exact outcome of the previous energy minimization was (potential energy, maximum force, copied and pasted from screen output or .log file). the real problem is with the mdrun function. could be a problem of the software? You have presented no evidence that would lead anyone to believe the problem is with mdrun. In the vast majority of cases, user input is the problem. -Justin Thanks Javier Justin Lemkul wrote On 11/11/13 11:24 AM, Carlos Javier Almeciga Diaz wrote: Hello evryone, I doing a simulation of a ligand-protein interaction with gromacs 4.5.5. Everything looks fine after I equilibrate the protein-ligand complex. I'm running these commands: grompp -f nvt.mdp -c em.gro -p topol.top -n index.ndx -o nvt.tpr mdrun -deffnm nvt Nevertheless, I got this error: Reading file nvt.tpr, VERSION 4.5.5 (double precision) Segmentation fault What should I do? Instantaneous failure typically indicates that the forces are nonsensically high and the constraint algorithm immediately fails. Likely the previous energy minimization did not adequately complete. -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalemkul@.umaryland | (410) 706-7441 == -- gmx-users mailing list gmx-users@ 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-request@ . * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists -- View this message in context: http://gromacs.5086.x6.nabble.com/segmentation-fault-on-gromacs-4-5-5-after-mdrun-tp5012431p5012458.html Sent from the GROMACS Users Forum mailing list archive at Nabble.com. -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalemkul@.umaryland | (410) 706-7441 == -- gmx-users mailing list gmx-users@ 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-request@ . * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists -- View this message in context: http://gromacs.5086.x6.nabble.com/segmentation-fault-on-gromacs-4-5-5-after-mdrun-tp5012431p5012464.html Sent from the GROMACS Users Forum mailing list archive at Nabble.com. -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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
Re: [gmx-users] Calculating diffusion coefficient in three dimension
Dear Sir, Thanks for the quick reply. So, I have to declare -type no flag. Isn't it?? and I have recently gone through Justin's membrane protein tutorial, where he has calculated diffusion coefficient for lipids in a membrane by creating an index group for a particular atom. So, here also shall I do the same thing? Moreover, mine is a coarse-grained system. On Tue, Nov 12, 2013 at 9:57 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/12/13 11:25 AM, Venkat Reddy wrote: Then, how to mention the direction for spherical particles Sir? Read g_msd -h again, paying specific attention to the -type flag. -Justin On Tue, Nov 12, 2013 at 7:28 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/12/13 8:55 AM, Venkat Reddy wrote: Thank you sir for the prompt reply. *g_msd -f traj.xtc -s topol.tpr -n msd.ndx -lateral z -o msd.xvg -tu ns* Here I am giving -lateral z (like for membrane simulations). Is it fine for spherical systems also? No. The system is a sphere, so what use is it to calculate motion perpendicular to z when you have lipids moving in all three spatial dimensions? A vesicle is very different from a membrane, in which the lipids move in a plane, thus making -lateral z useful. -Justin On Tue, Nov 12, 2013 at 5:32 PM, Dr. Vitaly Chaban vvcha...@gmail.com wrote: MSD is 3D by default. Dr. Vitaly V. Chaban On Tue, Nov 12, 2013 at 6:01 AM, Venkat Reddy venkat...@gmail.com wrote: Dear all, I am simulating a spherical lipid vesicle. I want to calculate the diffusion coefficient for each lipid component in 3D. How to calculate it using g_msd (or any other tool like g_velacc)? Thank you for your concern -- With Best Wishes Venkat Reddy Chirasani -- gmx-users mailing listgmx-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 -- gmx-users mailing listgmx-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 -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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 -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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 -- With Best Wishes Venkat Reddy Chirasani PhD student Laboratory of Computational Biophysics Department of Biotechnology IIT Madras Chennai INDIA-600036 -- gmx-users mailing listgmx-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
[gmx-users] Change in the positon of structural Zinc and calcium ions during MD
Hi Gromacs users, I'm doing protein-Bilayer MD simulations. Enzyme contains structural zinc and calcium ions during Energy minimization, NVT and NPT stage, ions are changing there position even though I applied position restraints for the atoms and ions. Anyone could help me out. Thanks -- View this message in context: http://gromacs.5086.x6.nabble.com/Change-in-the-positon-of-structural-Zinc-and-calcium-ions-during-MD-tp5012467.html Sent from the GROMACS Users Forum mailing list archive at Nabble.com. -- gmx-users mailing listgmx-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
Re: [gmx-users] Calculating diffusion coefficient in three dimension
On 11/12/13 12:30 PM, Venkat Reddy wrote: Dear Sir, Thanks for the quick reply. So, I have to declare -type no flag. Isn't it?? The options for the -type flag are x, y, or z. You said you wanted the diffusion coefficient in each spatial dimension. That is precisely what this option will do. and I have recently gone through Justin's membrane protein tutorial, where You mean my tutorial :) he has calculated diffusion coefficient for lipids in a membrane by creating an index group for a particular atom. So, here also shall I do the same thing? Moreover, mine is a coarse-grained system. Yes, a representative atom is usually what is passed to g_msd. -Justin On Tue, Nov 12, 2013 at 9:57 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/12/13 11:25 AM, Venkat Reddy wrote: Then, how to mention the direction for spherical particles Sir? Read g_msd -h again, paying specific attention to the -type flag. -Justin On Tue, Nov 12, 2013 at 7:28 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/12/13 8:55 AM, Venkat Reddy wrote: Thank you sir for the prompt reply. *g_msd -f traj.xtc -s topol.tpr -n msd.ndx -lateral z -o msd.xvg -tu ns* Here I am giving -lateral z (like for membrane simulations). Is it fine for spherical systems also? No. The system is a sphere, so what use is it to calculate motion perpendicular to z when you have lipids moving in all three spatial dimensions? A vesicle is very different from a membrane, in which the lipids move in a plane, thus making -lateral z useful. -Justin On Tue, Nov 12, 2013 at 5:32 PM, Dr. Vitaly Chaban vvcha...@gmail.com wrote: MSD is 3D by default. Dr. Vitaly V. Chaban On Tue, Nov 12, 2013 at 6:01 AM, Venkat Reddy venkat...@gmail.com wrote: Dear all, I am simulating a spherical lipid vesicle. I want to calculate the diffusion coefficient for each lipid component in 3D. How to calculate it using g_msd (or any other tool like g_velacc)? Thank you for your concern -- With Best Wishes Venkat Reddy Chirasani -- gmx-users mailing listgmx-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 -- gmx-users mailing listgmx-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 -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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 -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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 -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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
Re: [gmx-users] Change in the positon of structural Zinc and calcium ions during MD
On 11/12/13 12:35 PM, Rama wrote: Hi Gromacs users, I'm doing protein-Bilayer MD simulations. Enzyme contains structural zinc and calcium ions during Energy minimization, NVT and NPT stage, ions are changing there position even though I applied position restraints for the atoms and ions. Anyone could help me out. If the atoms are moving drastically, then the position restraints aren't taking effect, either because they haven't been invoked or they are insufficient. Without your full .mdp file and relevant topology snippets, there's no real way to know. You can check the .log files, too - if there is not an energy component corresponding to position restraints, they haven't been activated. -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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
Re: [gmx-users] Calculating diffusion coefficient in three dimension
Thanks Justin. So, I have to calculate diffusion coefficient three times (x,y,z) and finally add-up together to get in 3D??? On Tue, Nov 12, 2013 at 11:25 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/12/13 12:30 PM, Venkat Reddy wrote: Dear Sir, Thanks for the quick reply. So, I have to declare -type no flag. Isn't it?? The options for the -type flag are x, y, or z. You said you wanted the diffusion coefficient in each spatial dimension. That is precisely what this option will do. and I have recently gone through Justin's membrane protein tutorial, where You mean my tutorial :) he has calculated diffusion coefficient for lipids in a membrane by creating an index group for a particular atom. So, here also shall I do the same thing? Moreover, mine is a coarse-grained system. Yes, a representative atom is usually what is passed to g_msd. -Justin On Tue, Nov 12, 2013 at 9:57 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/12/13 11:25 AM, Venkat Reddy wrote: Then, how to mention the direction for spherical particles Sir? Read g_msd -h again, paying specific attention to the -type flag. -Justin On Tue, Nov 12, 2013 at 7:28 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/12/13 8:55 AM, Venkat Reddy wrote: Thank you sir for the prompt reply. *g_msd -f traj.xtc -s topol.tpr -n msd.ndx -lateral z -o msd.xvg -tu ns* Here I am giving -lateral z (like for membrane simulations). Is it fine for spherical systems also? No. The system is a sphere, so what use is it to calculate motion perpendicular to z when you have lipids moving in all three spatial dimensions? A vesicle is very different from a membrane, in which the lipids move in a plane, thus making -lateral z useful. -Justin On Tue, Nov 12, 2013 at 5:32 PM, Dr. Vitaly Chaban vvcha...@gmail.com wrote: MSD is 3D by default. Dr. Vitaly V. Chaban On Tue, Nov 12, 2013 at 6:01 AM, Venkat Reddy venkat...@gmail.com wrote: Dear all, I am simulating a spherical lipid vesicle. I want to calculate the diffusion coefficient for each lipid component in 3D. How to calculate it using g_msd (or any other tool like g_velacc)? Thank you for your concern -- With Best Wishes Venkat Reddy Chirasani -- gmx-users mailing listgmx-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 -- gmx-users mailing listgmx-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 -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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 -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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 -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441
Re: [gmx-users] Calculating diffusion coefficient in three dimension
On 11/12/13 1:33 PM, Venkat Reddy wrote: Thanks Justin. So, I have to calculate diffusion coefficient three times (x,y,z) and finally add-up together to get in 3D??? If you just want the overall diffusion constant, that's what g_msd does without any additional options. -Justin On Tue, Nov 12, 2013 at 11:25 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/12/13 12:30 PM, Venkat Reddy wrote: Dear Sir, Thanks for the quick reply. So, I have to declare -type no flag. Isn't it?? The options for the -type flag are x, y, or z. You said you wanted the diffusion coefficient in each spatial dimension. That is precisely what this option will do. and I have recently gone through Justin's membrane protein tutorial, where You mean my tutorial :) he has calculated diffusion coefficient for lipids in a membrane by creating an index group for a particular atom. So, here also shall I do the same thing? Moreover, mine is a coarse-grained system. Yes, a representative atom is usually what is passed to g_msd. -Justin On Tue, Nov 12, 2013 at 9:57 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/12/13 11:25 AM, Venkat Reddy wrote: Then, how to mention the direction for spherical particles Sir? Read g_msd -h again, paying specific attention to the -type flag. -Justin On Tue, Nov 12, 2013 at 7:28 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/12/13 8:55 AM, Venkat Reddy wrote: Thank you sir for the prompt reply. *g_msd -f traj.xtc -s topol.tpr -n msd.ndx -lateral z -o msd.xvg -tu ns* Here I am giving -lateral z (like for membrane simulations). Is it fine for spherical systems also? No. The system is a sphere, so what use is it to calculate motion perpendicular to z when you have lipids moving in all three spatial dimensions? A vesicle is very different from a membrane, in which the lipids move in a plane, thus making -lateral z useful. -Justin On Tue, Nov 12, 2013 at 5:32 PM, Dr. Vitaly Chaban vvcha...@gmail.com wrote: MSD is 3D by default. Dr. Vitaly V. Chaban On Tue, Nov 12, 2013 at 6:01 AM, Venkat Reddy venkat...@gmail.com wrote: Dear all, I am simulating a spherical lipid vesicle. I want to calculate the diffusion coefficient for each lipid component in 3D. How to calculate it using g_msd (or any other tool like g_velacc)? Thank you for your concern -- With Best Wishes Venkat Reddy Chirasani -- gmx-users mailing listgmx-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 -- gmx-users mailing listgmx-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 -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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 -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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 -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20
[gmx-users] Re: Change in the position of structural Zinc and calcium ions during MD
Hi Justin, Below I pasted .mdp file and topology. In .log file I could see energy term for position restraints. .mdp file--- title = NPT Equilibration define = -DPOSRES ; position restraints for protein ; Run parameters integrator = md; leap-frog integrator nsteps = 50; 2 * 50 = 1000 ps (1 ns) dt = 0.002 ; 2 fs ; Output control nstxout = 500 ; save coordinates every 1 ps nstvout = 500 ; save velocities every 1 ps nstenergy = 500 ; save energies every 1 ps nstlog = 500 ; update log file every 1 ps ; Bond parameters continuation= yes ; Restarting after NVT constraint_algorithm = lincs; holonomic constraints constraints = all-bonds ; all bonds (even heavy atom-H bonds) constrained lincs_iter = 1 ; accuracy of LINCS lincs_order = 4 ; also related to accuracy ; Neighborsearching ns_type = grid ; search neighboring grid cels nstlist = 5 ; 10 fs rlist = 1.2 ; short-range neighborlist cutoff (in nm) rcoulomb= 1.2 ; short-range electrostatic cutoff (in nm) rvdw= 1.2 ; short-range van der Waals cutoff (in nm) ; Electrostatics coulombtype = PME ; Particle Mesh Ewald for long-range electrostatics pme_order = 4 ; cubic interpolation fourierspacing = 0.16 ; grid spacing for FFT ; Temperature coupling is on tcoupl = Nose-Hoover ; More accurate thermostat tc-grps = Protein_CA_ZN DMPC SOL_CL; three coupling groups - more accurate tau_t = 0.50.5 0.5 ; time constant, in ps ref_t = 298298 310 ; reference temperature, one for each group, in K ; Pressure coupling is on pcoupl = Parrinello-Rahman ; Pressure coupling on in NPT pcoupltype = semiisotropic ; uniform scaling of x-y box vectors, independent z tau_p = 5.0 ; time constant, in ps ref_p = 1.0 1.0 ; reference pressure, x-y, z (in bar) compressibility = 4.5e-54.5e-5 ; isothermal compressibility, bar^-1 ; Periodic boundary conditions pbc = xyz ; 3-D PBC ; Dispersion correction DispCorr= EnerPres ; account for cut-off vdW scheme ; Velocity generation gen_vel = no; Velocity generation is off ; COM motion removal ; These options remove motion of the protein/bilayer relative to the solvent/ions nstcomm = 1 comm-mode = Linear comm-grps = Protein_DMPC SOL_CL ; Scale COM of reference coordinates refcoord_scaling = com topol.top ; Include Position restraint file #ifdef POSRES #include posre.itp #endif ; Strong position restraints for InflateGRO #ifdef STRONG_POSRES #include strong_posre.itp #endif ; Include DMPC topology #include rama4LJ.ff/dmpcLJ.itp ; Include water topology #include rama4LJ.ff/spc.itp #ifdef POSRES_WATER ; Position restraint for each water oxygen [ position_restraints ] ; i funct fcxfcyfcz 11 1000 1000 1000 #endif ; Include topology for ions #include rama4LJ.ff/ions.itp ---.log file Energies (kJ/mol) AngleProper Dih. Ryckaert-Bell. Improper Dih. LJ-14 1.77761e+043.10548e+037.97673e+034.40586e+028.14131e+03 Coulomb-14LJ (SR) Disper. corr. Coulomb (SR) Coul. recip. 2.59758e+042.74092e+04 -2.56846e+03 -4.68637e+05 -1.67418e+05 Position Rest. PotentialKinetic En. Total EnergyTemperature 7.09403e+02 -5.47088e+058.83115e+04 -4.58777e+053.07118e+02 Pres. DC (bar) Pressure (bar) Constr. rmsd -2.00493e+021.00080e+000.0e+00 Thanks -- View this message in context: http://gromacs.5086.x6.nabble.com/Change-in-the-positon-of-structural-Zinc-and-calcium-ions-during-MD-tp5012467p5012474.html Sent from the GROMACS Users Forum mailing list archive at Nabble.com. -- gmx-users mailing listgmx-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
Re: [gmx-users] Calculating diffusion coefficient in three dimension
Hi Venkat, Can you make it a bit clearer what you actually want? If it is the diffusion of the lipids along the curved surface of the vesicle, rather than simply the overall 3D diffusion, this is not trivial to calculate as I don't believe g_msd will do this for you. This property has been studied before though, so I suggest you search the literature for papers simulating vesicles to see how the lipid diffusion was calculated. Cheers Tom On 11/12/2013 06:35 PM, Justin Lemkul wrote: On 11/12/13 1:33 PM, Venkat Reddy wrote: Thanks Justin. So, I have to calculate diffusion coefficient three times (x,y,z) and finally add-up together to get in 3D??? If you just want the overall diffusion constant, that's what g_msd does without any additional options. -Justin On Tue, Nov 12, 2013 at 11:25 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/12/13 12:30 PM, Venkat Reddy wrote: Dear Sir, Thanks for the quick reply. So, I have to declare -type no flag. Isn't it?? The options for the -type flag are x, y, or z. You said you wanted the diffusion coefficient in each spatial dimension. That is precisely what this option will do. and I have recently gone through Justin's membrane protein tutorial, where You mean my tutorial :) he has calculated diffusion coefficient for lipids in a membrane by creating an index group for a particular atom. So, here also shall I do the same thing? Moreover, mine is a coarse-grained system. Yes, a representative atom is usually what is passed to g_msd. -Justin On Tue, Nov 12, 2013 at 9:57 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/12/13 11:25 AM, Venkat Reddy wrote: Then, how to mention the direction for spherical particles Sir? Read g_msd -h again, paying specific attention to the -type flag. -Justin On Tue, Nov 12, 2013 at 7:28 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/12/13 8:55 AM, Venkat Reddy wrote: Thank you sir for the prompt reply. *g_msd -f traj.xtc -s topol.tpr -n msd.ndx -lateral z -o msd.xvg -tu ns* Here I am giving -lateral z (like for membrane simulations). Is it fine for spherical systems also? No. The system is a sphere, so what use is it to calculate motion perpendicular to z when you have lipids moving in all three spatial dimensions? A vesicle is very different from a membrane, in which the lipids move in a plane, thus making -lateral z useful. -Justin On Tue, Nov 12, 2013 at 5:32 PM, Dr. Vitaly Chaban vvcha...@gmail.com wrote: MSD is 3D by default. Dr. Vitaly V. Chaban On Tue, Nov 12, 2013 at 6:01 AM, Venkat Reddy venkat...@gmail.com wrote: Dear all, I am simulating a spherical lipid vesicle. I want to calculate the diffusion coefficient for each lipid component in 3D. How to calculate it using g_msd (or any other tool like g_velacc)? Thank you for your concern -- With Best Wishes Venkat Reddy Chirasani -- gmx-users mailing listgmx-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 -- gmx-users mailing listgmx-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 -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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 -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users * Please
Re: mdrun on 8-core AMD + GTX TITAN (was: Re: [gmx-users] Re: Gromacs-4.6 on two Titans GPUs)
Hi Mark and Szilard Thanks for your both suggestions. They are very helpful. Neither run had a PP-PME work distribution suitable for the hardware it was running on (and fixing that for each run requires opposite changes). Adding a GPU and hoping to see scaling requires that there be proportionately more GPU work available to do, *and* enough absolute work to do. mdrun tries to do this, and reports early in the log file, which is one of the reasons Szilard asked to see whole log files - please use a file sharing service to do that. This task involves GPU calculation. We would not see PP-PME work distribution. This is a good hint from the angle of PP-PME work distribution. And I guessed that two GPUs' calculations are fast / or no enough work for GPU calculation, which is aligned with your explanation. Please see logs below again. ONE GPU## http://pastebin.com/B6bRUVSa TWO GPUs## http://pastebin.com/SLAYnejP As you can see, this test was made on the same node regardless of networking. Can the performance be improved say 50% more when 2 GPUs are used on a general task ? If yes, how ? Indeed, as Richard pointed out, I was asking for *full* logs, these summaries can't tell much, the table above the summary entitled R E A L C Y C L E A N D T I M E A C C O U N T I N G as well as other reported information across the log file is what I need to make an assessment of your simulations' performance. Please see below. However, in your case I suspect that the bottleneck is multi-threaded scaling on the AMD CPUs and you should probably decrease the number of threads per MPI rank and share GPUs between 2-4 ranks. After I test all three clusters, I found it may NOT be an issue of AMD cpus. Intel cpus has the SAME scaling issue. However, I am curious as to how you justify the setup of 2-4 ranks sharing GPUs ? Can you please explain it a bit more ? NUMA effects on multi-socket AMD processors are particularly severe; the way GROMACS uses OpenMP is not well suited to them. Using a rank (or two) per socket will greatly reduce those effects, but introduces different algorithmic overhead from the need to do DD and explicitly communicate between ranks. (You can see the latter in your .log file snippets below.) Also, that means the parcel of PP work available from a rank to give to the GPU is smaller, which is the opposite of what you'd like for GPU performance and/or scaling. We are working on a general solution for this and lots of related issues in the post-5.0 space, but there is a very hard limitation imposed by the need to amortize the cost of CPU-GPU transfer by having lots of PP work available to do. Is this reason why the scaling of two GPUs won't happen because of smaller PP workload ? From the implication, I am wondering if we can increase PP workload through parameters in a mdp file. The question is what parameters are mostly related to PP workload ? Would you please give more specific suggestions ? NOTE: The GPU has 20% more load than the CPU. This imbalance causes performance loss, consider using a shorter cut-off and a finer PME grid. This note needs to be addressed before maximum throughput is achieved and the question of scaling is worth considering. Ideally, Wait GPU local should be nearly zero, achieved as suggested above. Note that launch+force+mesh+wait is the sum of gpu total! But much of the information needed is higher up the log file, and the whole question is constrained by things like rvdw. From the note, it clearly suggested a shorter cut-off and a finer PME grid. I am not sure how to set up a finer PME grid but I am able to set up shorter cut-offs . However, it is risky to do so based on others' reports. Indeed, I see differences among tests for 1 GPU. Here cutoffs refer to rlist, rvdw and rcoulomb. I found that the smaller values of cutoffs, the faster computations. The question is how small they can go because it is interesting to know if these different cutoffs generate equally good simulations. As to two GPUs, when I set up larger cut-offs, these two GPUs in the same node had been very busy. However, the outcome in such a configuration is worse in terms of ns/day and time. So what dose a finer PME grid mean, with respect to GPU workload ? You mention the sum of GPU total is launch + force + mesh + wait.I thought PME mesh is carried out by CPU instead of GPU. Do I miss something here ? I thought GPU is responsible for the calculation of short-ranged non-bonded force whereas CPU is responsible for that of bonded and PME long-ranged force. Can you clarify it here ? Also, would rvdw play an important role in improving the performance of GPU calculation ? Unfortunately you didn't copy the GPU timing stuff here! Roughly, all the performance gain you are seeing here is eliminating most of the single-GPU wait gpu term by throwing
Re: [gmx-users] Re: Change in the position of structural Zinc and calcium ions during MD
On 11/12/13 1:47 PM, Rama wrote: Hi Justin, Below I pasted .mdp file and topology. In .log file I could see energy term for position restraints. .mdp file--- title = NPT Equilibration define = -DPOSRES ; position restraints for protein ; Run parameters integrator = md; leap-frog integrator nsteps = 50; 2 * 50 = 1000 ps (1 ns) dt = 0.002 ; 2 fs ; Output control nstxout = 500 ; save coordinates every 1 ps nstvout = 500 ; save velocities every 1 ps nstenergy = 500 ; save energies every 1 ps nstlog = 500 ; update log file every 1 ps ; Bond parameters continuation= yes ; Restarting after NVT constraint_algorithm = lincs; holonomic constraints constraints = all-bonds ; all bonds (even heavy atom-H bonds) constrained lincs_iter = 1 ; accuracy of LINCS lincs_order = 4 ; also related to accuracy ; Neighborsearching ns_type = grid ; search neighboring grid cels nstlist = 5 ; 10 fs rlist = 1.2 ; short-range neighborlist cutoff (in nm) rcoulomb= 1.2 ; short-range electrostatic cutoff (in nm) rvdw= 1.2 ; short-range van der Waals cutoff (in nm) ; Electrostatics coulombtype = PME ; Particle Mesh Ewald for long-range electrostatics pme_order = 4 ; cubic interpolation fourierspacing = 0.16 ; grid spacing for FFT ; Temperature coupling is on tcoupl = Nose-Hoover ; More accurate thermostat tc-grps = Protein_CA_ZN DMPC SOL_CL; three coupling groups - more accurate tau_t = 0.50.5 0.5 ; time constant, in ps ref_t = 298298 310 ; reference temperature, one for each group, in K ; Pressure coupling is on pcoupl = Parrinello-Rahman ; Pressure coupling on in NPT pcoupltype = semiisotropic ; uniform scaling of x-y box vectors, independent z tau_p = 5.0 ; time constant, in ps ref_p = 1.0 1.0 ; reference pressure, x-y, z (in bar) compressibility = 4.5e-54.5e-5 ; isothermal compressibility, bar^-1 ; Periodic boundary conditions pbc = xyz ; 3-D PBC ; Dispersion correction DispCorr= EnerPres ; account for cut-off vdW scheme ; Velocity generation gen_vel = no; Velocity generation is off ; COM motion removal ; These options remove motion of the protein/bilayer relative to the solvent/ions nstcomm = 1 comm-mode = Linear comm-grps = Protein_DMPC SOL_CL ; Scale COM of reference coordinates refcoord_scaling = com topol.top ; Include Position restraint file #ifdef POSRES #include posre.itp #endif ; Strong position restraints for InflateGRO #ifdef STRONG_POSRES #include strong_posre.itp #endif ; Include DMPC topology #include rama4LJ.ff/dmpcLJ.itp ; Include water topology #include rama4LJ.ff/spc.itp #ifdef POSRES_WATER ; Position restraint for each water oxygen [ position_restraints ] ; i funct fcxfcyfcz 11 1000 1000 1000 #endif ; Include topology for ions #include rama4LJ.ff/ions.itp Do you have appropriate [position_restraints] assigned in this topology? None of the above, as shown, pertains to the ions, and the only relevant #ifdef block that would be triggered by -DPOSRES is for the protein. -Justin ---.log file Energies (kJ/mol) AngleProper Dih. Ryckaert-Bell. Improper Dih. LJ-14 1.77761e+043.10548e+037.97673e+034.40586e+028.14131e+03 Coulomb-14LJ (SR) Disper. corr. Coulomb (SR) Coul. recip. 2.59758e+042.74092e+04 -2.56846e+03 -4.68637e+05 -1.67418e+05 Position Rest. PotentialKinetic En. Total EnergyTemperature 7.09403e+02 -5.47088e+058.83115e+04 -4.58777e+053.07118e+02 Pres. DC (bar) Pressure (bar) Constr. rmsd -2.00493e+021.00080e+000.0e+00 Thanks -- View this message in context: http://gromacs.5086.x6.nabble.com/Change-in-the-positon-of-structural-Zinc-and-calcium-ions-during-MD-tp5012467p5012474.html Sent from the GROMACS Users Forum mailing list archive at Nabble.com. -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing
[gmx-users] help about ibi
Hi: When I am running the ibi procedure, I get the following error message: A coordinate in file conf.gro does not contain a '.' Additionally, I check the coordinate file of confout.gro in step_001. It showed that 'nan' symbol appeared in confout.gro. What is wrong with this? How can I fix it? I am very appreciating for anyone's help. Best Wishes! Zhicheng Guo-- gmx-users mailing listgmx-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
Re: [gmx-users] Calculating diffusion coefficient in three dimension
Dear Justin and Piggot, Thanks for the suggestions. Actually, I have constructed a CG lipid vesicle by placing lipids in random conformation in a simulation box. My lipid system is heterogeneous, i.e., it has different types of lipids (POPC,CHOLESTEROL,PPC...etc). Some of them occupy the outer layer of vesicle (POPC), and some stay in the intermediate region (CHOL). So, I want to calculate the diffusion rates of these lipids. Since POPC forms the surface (polar heads interacting with water and their tails points to the core), I suppose we have to calculate 2D diffusion for POPC. For the lipids in the core, they can diffuse in 3-dimension. So, it requires a 3D diffusion coefficient for these core lipids. How to calculate 2D and 3D diffusion coeff.? Hope I am clear. On Wed, Nov 13, 2013 at 12:58 AM, Thomas Piggot t.pig...@soton.ac.ukwrote: Hi Venkat, Can you make it a bit clearer what you actually want? If it is the diffusion of the lipids along the curved surface of the vesicle, rather than simply the overall 3D diffusion, this is not trivial to calculate as I don't believe g_msd will do this for you. This property has been studied before though, so I suggest you search the literature for papers simulating vesicles to see how the lipid diffusion was calculated. Cheers Tom On 11/12/2013 06:35 PM, Justin Lemkul wrote: On 11/12/13 1:33 PM, Venkat Reddy wrote: Thanks Justin. So, I have to calculate diffusion coefficient three times (x,y,z) and finally add-up together to get in 3D??? If you just want the overall diffusion constant, that's what g_msd does without any additional options. -Justin On Tue, Nov 12, 2013 at 11:25 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/12/13 12:30 PM, Venkat Reddy wrote: Dear Sir, Thanks for the quick reply. So, I have to declare -type no flag. Isn't it?? The options for the -type flag are x, y, or z. You said you wanted the diffusion coefficient in each spatial dimension. That is precisely what this option will do. and I have recently gone through Justin's membrane protein tutorial, where You mean my tutorial :) he has calculated diffusion coefficient for lipids in a membrane by creating an index group for a particular atom. So, here also shall I do the same thing? Moreover, mine is a coarse-grained system. Yes, a representative atom is usually what is passed to g_msd. -Justin On Tue, Nov 12, 2013 at 9:57 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/12/13 11:25 AM, Venkat Reddy wrote: Then, how to mention the direction for spherical particles Sir? Read g_msd -h again, paying specific attention to the -type flag. -Justin On Tue, Nov 12, 2013 at 7:28 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/12/13 8:55 AM, Venkat Reddy wrote: Thank you sir for the prompt reply. *g_msd -f traj.xtc -s topol.tpr -n msd.ndx -lateral z -o msd.xvg -tu ns* Here I am giving -lateral z (like for membrane simulations). Is it fine for spherical systems also? No. The system is a sphere, so what use is it to calculate motion perpendicular to z when you have lipids moving in all three spatial dimensions? A vesicle is very different from a membrane, in which the lipids move in a plane, thus making -lateral z useful. -Justin On Tue, Nov 12, 2013 at 5:32 PM, Dr. Vitaly Chaban vvcha...@gmail.com wrote: MSD is 3D by default. Dr. Vitaly V. Chaban On Tue, Nov 12, 2013 at 6:01 AM, Venkat Reddy venkat...@gmail.com wrote: Dear all, I am simulating a spherical lipid vesicle. I want to calculate the diffusion coefficient for each lipid component in 3D. How to calculate it using g_msd (or any other tool like g_velacc)? Thank you for your concern -- With Best Wishes Venkat Reddy Chirasani -- gmx-users mailing listgmx-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 -- gmx-users mailing listgmx-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 -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441
[gmx-users] Re: Bilayer COM removal issue: Large VCM
Hi All, Any suggestions? Thanks, On Mon, Nov 11, 2013 at 12:38 AM, rajat desikan rajatdesi...@gmail.comwrote: Hi All, I am experiencing a few problems in membrane simulations wrt COM removal. I downloaded a 400 ns pre-equilibrated Slipid-DMPC membrane with all the accompanying files. I then carried out the following steps: 1) energy minimization 2) NVT Eq - 100 ps 3) NPT Eq - 250 ps (Berendsen temp, Pres coupling) Then I used g_select to select the upper and lower DMPC leaflets. The then carried out a 250 ps NPT eq again. The only change was: comm-grps= SOL DMPC == comm-grps= SOL upper lower On every step in log file, I get the following message: *Step Time Lambda 124000 248.00.0 Large VCM(group lower): -0.00051, -0.00515, -0.00652, Temp-cm: 8.11828e+29 Energies (kJ/mol)U-BProper Dih. Improper Dih. LJ-14 Coulomb-147.23818e+044.19778e+046.46641e+024.54801e+03 -1.45245e+05 LJ (SR)LJ (LR) Disper. corr. Coulomb (SR) Coul. recip.2.79689e+04 -3.78407e+03 -2.10679e+03 -5.84134e+05 -8.87497e+04 PotentialKinetic En. Total EnergyTemperature Pres. DC (bar)-6.76497e+051.76468e+05 -5.00029e+05 3.10424e+02 -1.05704e+02 Pressure (bar) Constr. rmsd -1.85927e+02 6.42934e-06* *Large VCM(group lower): -0.00187, -0.00369, 0.00032, Temp-cm: 2.02076e+29 Large VCM(group lower): -0.00725, -0.00278, -0.00549, Temp-cm: 1.05988e+30Large VCM(group lower): 0.00020, 0.00308, -0.00176, Temp-cm: 1.48126e+29Large VCM(group lower): -0.00541, 0.00546, -0.00166, Temp-cm: 7.24656e+29 Large VCM(group lower): -0.00220, 0.00362, -0.00741, Temp-cm: 8.53812e+29Large VCM(group lower): 0.00140, -0.00160, 0.00029, Temp-cm: 5.39679e+28Large VCM(group lower): -0.00056, -0.00293, -0.00364, Temp-cm: 2.59422e+29 Large VCM(group lower): -0.00172, -0.00260, 0.00494, Temp-cm: 3.99945e+29Large VCM(group lower): 0.00252, 0.00594, 0.00068, Temp-cm: 4.93342e+29* *DD step 124999 vol min/aver 0.702 load imb.: force 1.3% pme mesh/force 0.636* I do not know what to make of it. There are no issues when I remove COM for the entire system. I have seen this issue come up a few times in the archives too, but I didn't find a satisfactory solution since the bilayer was very well equilibrated. I would appreciate any suggestions. Thank you. -- Rajat Desikan (Ph.D Scholar) Prof. K. Ganapathy Ayappa's Lab (no 13), Dept. of Chemical Engineering, Indian Institute of Science, Bangalore -- Rajat Desikan (Ph.D Scholar) Prof. K. Ganapathy Ayappa's Lab (no 13), Dept. of Chemical Engineering, Indian Institute of Science, Bangalore -- gmx-users mailing listgmx-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
Re: [gmx-users] Re: g_analyze
In addition to my previous question, I have another question about g_analyze. When I used the hbond.xvg file to get the average and plotted the average.xvg file I found that the average value is round 4 to 5 according to the graph. But g_analyze in its final calculation gives 7.150 as the average values... Here's the link for the graph and result of average value calculated by g_analyze :- std. dev.relative deviation of standard - cumulants from those of set average deviation sqrt(n-1) a Gaussian distribition cum. 3 cum. 4 SS1 * 7.150740e+00 * 8.803173e-01 1.760635e-02 0.0620.163 SS2 1.490604e+00 1.164761e+00 2.329523e-02 0.4950.153 https://www.dropbox.com/s/1vqixenyerha7qq/115-water.png Here's the link hbond.xvg file and its averaged file https://www.dropbox.com/s/4n0m47o3mrjn3o8/hbond_115-water.xvg https://www.dropbox.com/s/4n0m47o3mrjn3o8/hbond_115-water.xvg On Mon, Nov 11, 2013 at 3:30 PM, bharat gupta bharat.85.m...@gmail.comwrote: thank you informing about g_rdf... Is it possible to dump the structure with those average water molecules interacting with the residues. I generated the hbond.log file which gives the details but I need to generate a figure for this ?? On Mon, Nov 11, 2013 at 10:40 AM, Justin Lemkul jalem...@vt.edu wrote: On 11/10/13 8:38 PM, bharat gupta wrote: But trjorder can be used to calculate the hydration layer or shell around residues ... Right ?? Yes, but I also tend to think that integrating an RDF is also a more straightforward way of doing that. With trjorder, you set some arbitrary cutoff that may or may not be an informed decision - with an RDF it is clear where the hydration layers are. -Justin On Mon, Nov 11, 2013 at 10:35 AM, Justin Lemkul jalem...@vt.edu wrote: On 11/10/13 8:30 PM, bharat gupta wrote: Thanks for your reply. I was missing the scientific notation part. Now everything is fine. Regarding trjorder, it doesn't measure h-bonds but gives the water nearest to protein. I wouldn't try to draw any sort of comparison between the output of trjorder and g_hbond. If you want to measure H-bonds, there's only one tool for that. -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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 -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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 -- gmx-users mailing listgmx-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
[gmx-users] help about logfile
Hi I am confusing with the output file (.log) about the sample frequency (frames) in my simulation. The average information in .log file is 'Statistics over 31 steps using 3001 frames' where nstxout =4000 and nstlog =4000. While, 'Statistics over 31 steps using 20001 frames', appeared in .log file, where nstxout =15 and nstlog =15. In my opinion, the former sample frequency should be 75 frames rather than 3001 frames and the latter is right. I want to know how can I control the sample frequency, arbitrarily.-- gmx-users mailing listgmx-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
Re: [gmx-users] Re: g_analyze
On 11/11/13 1:30 AM, bharat gupta wrote: thank you informing about g_rdf... Is it possible to dump the structure with those average water molecules interacting with the residues. I generated the hbond.log file which gives the details but I need to generate a figure for this ?? g_select -Justin On Mon, Nov 11, 2013 at 10:40 AM, Justin Lemkul jalem...@vt.edu wrote: On 11/10/13 8:38 PM, bharat gupta wrote: But trjorder can be used to calculate the hydration layer or shell around residues ... Right ?? Yes, but I also tend to think that integrating an RDF is also a more straightforward way of doing that. With trjorder, you set some arbitrary cutoff that may or may not be an informed decision - with an RDF it is clear where the hydration layers are. -Justin On Mon, Nov 11, 2013 at 10:35 AM, Justin Lemkul jalem...@vt.edu wrote: On 11/10/13 8:30 PM, bharat gupta wrote: Thanks for your reply. I was missing the scientific notation part. Now everything is fine. Regarding trjorder, it doesn't measure h-bonds but gives the water nearest to protein. I wouldn't try to draw any sort of comparison between the output of trjorder and g_hbond. If you want to measure H-bonds, there's only one tool for that. -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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 -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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 -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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
Re: [gmx-users] Re: g_analyze
On 11/11/13 4:06 AM, bharat gupta wrote: In addition to my previous question, I have another question about g_analyze. When I used the hbond.xvg file to get the average and plotted the average.xvg file I found that the average value is round 4 to 5 according to the graph. But g_analyze in its final calculation gives 7.150 as the average values... Here's the link for the graph and result of average value calculated by g_analyze :- std. dev.relative deviation of standard - cumulants from those of set average deviation sqrt(n-1) a Gaussian distribition cum. 3 cum. 4 SS1 * 7.150740e+00 * 8.803173e-01 1.760635e-02 0.0620.163 SS2 1.490604e+00 1.164761e+00 2.329523e-02 0.4950.153 https://www.dropbox.com/s/1vqixenyerha7qq/115-water.png Here's the link hbond.xvg file and its averaged file https://www.dropbox.com/s/4n0m47o3mrjn3o8/hbond_115-water.xvg https://www.dropbox.com/s/4n0m47o3mrjn3o8/hbond_115-water.xvg Neither of these files produce output that corresponds to the PNG image above. Both files have values in 6-9 H-bond range and thus agree with the g_analyze output, which I can reproduce. I suspect you're somehow getting your files mixed up. -Justin On Mon, Nov 11, 2013 at 3:30 PM, bharat gupta bharat.85.m...@gmail.comwrote: thank you informing about g_rdf... Is it possible to dump the structure with those average water molecules interacting with the residues. I generated the hbond.log file which gives the details but I need to generate a figure for this ?? On Mon, Nov 11, 2013 at 10:40 AM, Justin Lemkul jalem...@vt.edu wrote: On 11/10/13 8:38 PM, bharat gupta wrote: But trjorder can be used to calculate the hydration layer or shell around residues ... Right ?? Yes, but I also tend to think that integrating an RDF is also a more straightforward way of doing that. With trjorder, you set some arbitrary cutoff that may or may not be an informed decision - with an RDF it is clear where the hydration layers are. -Justin On Mon, Nov 11, 2013 at 10:35 AM, Justin Lemkul jalem...@vt.edu wrote: On 11/10/13 8:30 PM, bharat gupta wrote: Thanks for your reply. I was missing the scientific notation part. Now everything is fine. Regarding trjorder, it doesn't measure h-bonds but gives the water nearest to protein. I wouldn't try to draw any sort of comparison between the output of trjorder and g_hbond. If you want to measure H-bonds, there's only one tool for that. -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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 -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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 -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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
Re: [gmx-users] help about logfile
On 11/11/13 5:04 AM, guozhicheng222 wrote: Hi I am confusing with the output file (.log) about the sample frequency (frames) in my simulation. The average information in .log file is 'Statistics over 31 steps using 3001 frames' where nstxout =4000 and nstlog =4000. While, 'Statistics over 31 steps using 20001 frames', appeared in .log file, where nstxout =15 and nstlog =15. In my opinion, the former sample frequency should be 75 frames rather than 3001 frames and the latter is right. I want to know how can I control the sample frequency, arbitrarily. Check to make sure you set nstlog correctly in the first run. Look for the value in the .log file itself. It is highly unlikely that something so fundamental and simple is being executed incorrectly in one run, but correctly in another. -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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
[gmx-users] installation error under openSuse 12.2
Dear gmx-users, I am trying to install gromacs-4.6.3 on an older Dell that is running openSuse 12.2 Using DGMX_BUILD_OWN_FFTW=ON failed for me so to get through cmake I had to compile fftw from scratch and I followed the recommendation of going directly to their website. Now I am at the make stage and get the following message: [ 67%] Built target gmxfftw make[2]: *** No rule to make target `//home/koln/bin/gromacs-4.6.3/build/src/contrib/fftw/gmxfftw-prefix/lib/libfftw3f.a', needed by `src/gmxlib/libgmx.so.8'. Stop. make[1]: *** [src/gmxlib/CMakeFiles/gmx.dir/all] Error 2 make: *** [all] Error 2 When I check for libfftw, I get: koln@linux-5bim:~/bin/gromacs-4.6.3/build2 rpm -qf /usr/lib/libfftw3f.* libfftw3-3-32bit-3.3.2-1.1.2.x86_64 Any suggestions for how to get past this obstacle would be greatly appreciated! Koln -- View this message in context: http://gromacs.5086.x6.nabble.com/installation-error-under-openSuse-12-2-tp5012430.html Sent from the GROMACS Users Forum mailing list archive at Nabble.com. -- gmx-users mailing listgmx-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
Re: [gmx-users] segmentation fault on gromacs 4.5.5 after mdrun
On 11/11/13 11:24 AM, Carlos Javier Almeciga Diaz wrote: Hello evryone, I doing a simulation of a ligand-protein interaction with gromacs 4.5.5. Everything looks fine after I equilibrate the protein-ligand complex. I'm running these commands: grompp -f nvt.mdp -c em.gro -p topol.top -n index.ndx -o nvt.tpr mdrun -deffnm nvt Nevertheless, I got this error: Reading file nvt.tpr, VERSION 4.5.5 (double precision) Segmentation fault What should I do? Instantaneous failure typically indicates that the forces are nonsensically high and the constraint algorithm immediately fails. Likely the previous energy minimization did not adequately complete. -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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
[gmx-users] segmentation fault on gromacs 4.5.5 after mdrun
Hello evryone, I doing a simulation of a ligand-protein interaction with gromacs 4.5.5. Everything looks fine after I equilibrate the protein-ligand complex. I'm running these commands: grompp -f nvt.mdp -c em.gro -p topol.top -n index.ndx -o nvt.tpr mdrun -deffnm nvt Nevertheless, I got this error: Reading file nvt.tpr, VERSION 4.5.5 (double precision) Segmentation fault What should I do? Carlos Javier Alméciga Díaz, QF., PhD. Profesor Asistente Pontificia Universidad Javeriana Facultad de Ciencias Instituto de Errores Innatos del Metabolismo Tel: 57-1-3208320 Ext 4140-4099 Fax: 57-1-3208320 Ext 4099 Bogotá. D.C. - COLOMBIA cjalmec...@javeriana.edu.co http://www.javeriana.edu.co/ieim AVISO LEGAL: El presente correo electronico no representa la opinion o el consentimiento oficial de la PONTIFICIA UNIVERSIDAD JAVERIANA. Este mensaje es confidencial y puede contener informacion privilegiada la cual no puede ser usada ni divulgada a personas distintas de su destinatario. Esta prohibida la retencion, grabacion, utilizacion, aprovechamiento o divulgacion con cualquier proposito. Si por error recibe este mensaje, por favor destruya su contenido y avise a su remitente. En este aviso legal se omiten intencionalmente las tildes. Este mensaje ha sido revisado por un sistema antivirus, por lo que su contenido esta libre de virus. This e-mail has been scanned by an antivirus system, so its contents is virus free. -- gmx-users mailing listgmx-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
Re: [gmx-users] installation error under openSuse 12.2
On 11/11/13 11:27 AM, kolnkempff wrote: Dear gmx-users, I am trying to install gromacs-4.6.3 on an older Dell that is running openSuse 12.2 Using DGMX_BUILD_OWN_FFTW=ON failed for me so to get through cmake I had to compile fftw from scratch and I followed the recommendation of going directly to their website. Now I am at the make stage and get the following message: [ 67%] Built target gmxfftw make[2]: *** No rule to make target `//home/koln/bin/gromacs-4.6.3/build/src/contrib/fftw/gmxfftw-prefix/lib/libfftw3f.a', needed by `src/gmxlib/libgmx.so.8'. Stop. make[1]: *** [src/gmxlib/CMakeFiles/gmx.dir/all] Error 2 make: *** [all] Error 2 When I check for libfftw, I get: koln@linux-5bim:~/bin/gromacs-4.6.3/build2 rpm -qf /usr/lib/libfftw3f.* libfftw3-3-32bit-3.3.2-1.1.2.x86_64 Any suggestions for how to get past this obstacle would be greatly appreciated! Are you building in a clean installation directory? It looks to me that cmake is still trying to build the internal FFTW code, rather than looking for your manually installed FFTW. Providing your full cmake command would also be very useful (hint: always provide it). -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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
[gmx-users] Re: Reaction field zero and ions
Hello If I did the MD simulation using PME and neutralized with ions, and I want to rerun this time with reaction field zero, is there any problem if I keep the ions? This is for LIE calculation. I am using AMBER99SB. Thanks Williams -- gmx-users mailing listgmx-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
[gmx-users] hydrogen bond calculation problem
Dear gromacs users, I have NaCl-water system of finite concentration .I want to calculate water-water hydrogen bond in the first solvation shell around ion..So I used g_hbond -shell option.For single ion in water,its easy to calculate..But say,for water-water hydrogen bonds around two Na+ ions,it is saying to Select one atom for shell. Although I have more than 1 Na+ ions and I want to see water-water hydrogen bond around all Na+ ions simultaneously.Is it possible to take solvation shell of each ion at a same time? Thanks in advance. Regards, *Sushma* -- gmx-users mailing listgmx-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
Re: [gmx-users] About Compiler Compatibility for Gromacs 4.6.2 Compilation
On Nov 10, 2013 10:04 AM, Mark Abraham mark.j.abra...@gmail.com wrote: Yes (unless you are using AMD cpus), per the installation instructions, although you will probably do slightly better with GCC 4.7, and should not do a new install of 4.6.2 after 4.6.3 is released. In particular, 4.6.2 has an affinity-related performance regression when using external MPI libraries. Mark On Nov 10, 2013 7:55 AM, vidhya sankar scvsankar_...@yahoo.com wrote: Dear Justin and Mark Thank you for your Previous reply Can i Use the Following Intel Compiler for grmacs 4.6.2 in centos Linux OS ? Intel® C++ Composer XE 2013 for Linux it Includes Intel® C++ Compiler, Intel® Integrated Performance Primitives 7.1, Intel® Math Kernel Library 11.0, Intel Cilk™ Plus, the Intel® Threading Building Blocks (Intel® TBB)” -- gmx-users mailing listgmx-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 -- gmx-users mailing listgmx-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
Re: mdrun on 8-core AMD + GTX TITAN (was: Re: [gmx-users] Re: Gromacs-4.6 on two Titans GPUs)
On Sun, Nov 10, 2013 at 5:28 AM, Dwey Kauffman mpi...@gmail.com wrote: Hi Szilard, Thank you very much for your suggestions. Actually, I was jumping to conclusions too early, as you mentioned AMD cluster, I assumed you must have 12-16-core Opteron CPUs. If you have an 8-core (desktop?) AMD CPU, than you may not need to run more than one rank per GPU. Yes, we do have independent clusters of AMD, AMD opteron, Intel Corei7. All nodes of three clusters are installed with (at least) 1 GPU card. I have run the same test on these three clusters. Let's focus on a basic scaling issue: One GPU v.s Two GPUs within the same node of 8-core AMD cpu. Using 1 GPU, we can have a performance of ~32 ns/day. Using two GPU, we gain not much more ( ~38.5 ns/day ). It is about ~20% more performance. However, this is not really true because in some tests, I also saw only 2-5% more, which really surprised me. Neither run had a PP-PME work distribution suitable for the hardware it was running on (and fixing that for each run requires opposite changes). Adding a GPU and hoping to see scaling requires that there be proportionately more GPU work available to do, *and* enough absolute work to do. mdrun tries to do this, and reports early in the log file, which is one of the reasons Szilard asked to see whole log files - please use a file sharing service to do that. As you can see, this test was made on the same node regardless of networking. Can the performance be improved say 50% more when 2 GPUs are used on a general task ? If yes, how ? Indeed, as Richard pointed out, I was asking for *full* logs, these summaries can't tell much, the table above the summary entitled R E A L C Y C L E A N D T I M E A C C O U N T I N G as well as other reported information across the log file is what I need to make an assessment of your simulations' performance. Please see below. However, in your case I suspect that the bottleneck is multi-threaded scaling on the AMD CPUs and you should probably decrease the number of threads per MPI rank and share GPUs between 2-4 ranks. After I test all three clusters, I found it may NOT be an issue of AMD cpus. Intel cpus has the SAME scaling issue. However, I am curious as to how you justify the setup of 2-4 ranks sharing GPUs ? Can you please explain it a bit more ? NUMA effects on multi-socket AMD processors are particularly severe; the way GROMACS uses OpenMP is not well suited to them. Using a rank (or two) per socket will greatly reduce those effects, but introduces different algorithmic overhead from the need to do DD and explicitly communicate between ranks. (You can see the latter in your .log file snippets below.) Also, that means the parcel of PP work available from a rank to give to the GPU is smaller, which is the opposite of what you'd like for GPU performance and/or scaling. We are working on a general solution for this and lots of related issues in the post-5.0 space, but there is a very hard limitation imposed by the need to amortize the cost of CPU-GPU transfer by having lots of PP work available to do. You could try running mpirun -np 4 mdrun -ntomp 2 -gpu_id 0011 but I suspect this won't help because your scaling issue Your guess is correct but why is that ? it is worse. The more nodes are involved in a task, the performance is worse. in my experience even reaction field runs don't scale across nodes with 10G ethernet if you have more than 4-6 ranks per node trying to communicate (let alone with PME). What dose it mean let alone with PME ? how to do so ? by mdrun ? I do know mdrun -npme to specify PME process. If using PME (rather than RF), network demands are more severe. Thank you. Dwey ### One GPU R E A L C Y C L E A N D T I M E A C C O U N T I N G Computing: Nodes Th. Count Wall t (s) G-Cycles % - Neighbor search18 11 431.81713863.390 1.6 Launch GPU ops.18501 472.90615182.556 1.7 Force 185011328.61142654.785 4.9 PME mesh 18501 11561.327 371174.09042.8 Wait GPU local 185016888.008 221138.11125.5 NB X/F buffer ops. 189911216.49939055.455 4.5 Write traj.18 1030 12.741 409.039 0.0 Update 185011696.35854461.226 6.3 Constraints185011969.72663237.647 7.3 Rest 11458.82046835.133 5.4 - Total 1 27036.812 868011.431 100.0
[gmx-users] Umbrella Sampling tutorial
Dear Justin Thanks for your reply. You are right. I should not extrapolate too literally from your tutorial to my system. But, I have a general question: There is 2 groups in COM pulling method (reference group + pull group). If I want to use pull_geometry = distance, so, I should fix reference group to be immobile. Is it true? On the other hand, I want to know exactly using position restraining on reference group is optional or mandatory in COM pulling method? Best wishes -- gmx-users mailing listgmx-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
[gmx-users] Re: choosing force field
Thank you Justin for your kind help. The simple reason for considering only gromos parameter sets is that the parameters for the metal ions (in my protein) are not defined in other force fields. On Sat, Nov 9, 2013 at 7:18 PM, Justin Lemkul [via GROMACS] ml-node+s5086n5012376...@n6.nabble.com wrote: On 11/9/13 12:48 AM, pratibha wrote: Sorry for the previous mistake. Instead of 53a7, the force field which I used was 53a6. 53A6 is known to under-stabilize helices, so if a helix did not appear in a simulation using this force field, it is not definitive proof that the structure does not populate helical structures. I generally see mixed opinions in the literature in terms of which Gromos parameter set is the most reliable. As was asked by someone else, is there a reason you are only considering Gromos parameter sets? Others may be better suited to your study. -Justin On Fri, Nov 8, 2013 at 12:10 AM, Justin Lemkul [via GROMACS] [hidden email] http://user/SendEmail.jtp?type=nodenode=5012376i=0 wrote: On 11/7/13 12:14 PM, pratibha wrote: My protein contains metal ions which are parameterized only in gromos force field. Since I am a newbie to MD simulations, it would be difficult for me to parameterize those myself. Can you please guide me as per my previous mail which out of the two simulations should I consider more reliable-43a1 or 53a7? AFAIK, there is no such thing as 53A7, and your original message was full of similar typos, making it nearly impossible to figure out what you were actually doing. Can you indicate the actual force field(s) that you have been using in case someone has any ideas? The difference between 53A6 and 54A7 should be quite pronounced, in my experience, thus any guesses as to what 53A7 should be doing are not productive because I don't know what that is. -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 [hidden email] http://user/SendEmail.jtp?type=nodenode=5012325i=0 | (410) 706-7441 == -- gmx-users mailing list[hidden email] http://user/SendEmail.jtp?type=nodenode=5012325i=1 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 [hidden email] http://user/SendEmail.jtp?type=nodenode=5012325i=2. * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists -- If you reply to this email, your message will be added to the discussion below: . NAML http://gromacs.5086.x6.nabble.com/template/NamlServlet.jtp?macro=macro_viewerid=instant_html%21nabble%3Aemail.namlbase=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespacebreadcrumbs=notify_subscribers%21nabble%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_instant_email%21nabble%3Aemail.naml -- View this message in context: http://gromacs.5086.x6.nabble.com/choosing-force-field-tp5012242p5012370.html Sent from the GROMACS Users Forum mailing list archive at Nabble.com. -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 [hidden email] http://user/SendEmail.jtp?type=nodenode=5012376i=1 | (410) 706-7441 == -- gmx-users mailing list[hidden email]http://user/SendEmail.jtp?type=nodenode=5012376i=2 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 [hidden email]http://user/SendEmail.jtp?type=nodenode=5012376i=3. * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists -- If you reply to this email, your message will be added to the discussion below: http://gromacs.5086.x6.nabble.com/choosing-force-field-tp5012242p5012376.html To unsubscribe from choosing force field, click herehttp://gromacs.5086.x6.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_codenode=5012242code=a2Fwb29ycHJhdGliaGE3QGdtYWlsLmNvbXw1MDEyMjQyfC02NjkwNjY5MjU= .
[gmx-users] problem in running mdrun command
Dear all I encounter a problem while running command mdrun_mpi -v -deffnm em in gromacs. I am new to the gromacs. i just ran test calculation, *simulation of lyzozyme in water*. i am able to generate gro, tpr files. But in the final step i got following error. Thanks in advance. [localhost.localdomain:23122] mca: base: component_find: paffinity mca_paffinity_linux uses an MCA interface that is not recognized (component MCA v1.0.0 != supported MCA v2.0.0) -- ignored [localhost.localdomain:23123] mca: base: component_find: paffinity mca_paffinity_linux uses an MCA interface that is not recognized (component MCA v1.0.0 != supported MCA v2.0.0) -- ignored [localhost.localdomain:23123] mca: base: component_find: ras mca_ras_dash_host uses an MCA interface that is not recognized (component MCA v1.0.0 != supported MCA v2.0.0) -- ignored [localhost.localdomain:23123] mca: base: component_find: ras mca_ras_gridengine uses an MCA interface that is not recognized (component MCA v1.0.0 != supported MCA v2.0.0) -- ignored [localhost.localdomain:23123] mca: base: component_find: ras mca_ras_localhost uses an MCA interface that is not recognized (component MCA v1.0.0 != supported MCA v2.0.0) -- ignored [localhost.localdomain:23123] mca: base: component_find: errmgr mca_errmgr_hnp uses an MCA interface that is not recognized (component MCA v1.0.0 != supported MCA v2.0.0) -- ignored [localhost.localdomain:23123] mca: base: component_find: errmgr mca_errmgr_orted uses an MCA interface that is not recognized (component MCA v1.0.0 != supported MCA v2.0.0) -- ignored [localhost.localdomain:23123] mca: base: component_find: errmgr mca_errmgr_proxy uses an MCA interface that is not recognized (component MCA v1.0.0 != supported MCA v2.0.0) -- ignored [localhost.localdomain:23123] mca: base: component_find: iof mca_iof_proxy uses an MCA interface that is not recognized (component MCA v1.0.0 != supported MCA v2.0.0) -- ignored [localhost.localdomain:23123] mca: base: component_find: iof mca_iof_svc uses an MCA interface that is not recognized (component MCA v1.0.0 != supported MCA v2.0.0) -- ignored [localhost.localdomain:23122] mca: base: component_find: rcache mca_rcache_rb uses an MCA interface that is not recognized (component MCA v1.0.0 != supported MCA v2.0.0) -- ignored [localhost:23122] *** Process received signal *** [localhost:23122] Signal: Segmentation fault (11) [localhost:23122] Signal code: Address not mapped (1) [localhost:23122] Failing at address: 0x4498 [localhost:23122] [ 0] /lib64/libpthread.so.0() [0x3fcee0f500] [localhost:23122] [ 1] /usr/local/lib/libmpi.so.1(PMPI_Comm_size+0x4e) [0x2acc6d93727e] [localhost:23122] [ 2] /usr/local/gromacs/bin/../lib/libgmx_mpi.so.8(gmx_setup+0x32) [0x2acc6d195e02] [localhost:23122] [ 3] /usr/local/gromacs/bin/../lib/libgmx_mpi.so.8(init_par+0x51) [0x2acc6d234251] [localhost:23122] [ 4] mdrun_mpi(cmain+0x11f9) [0x435799] [localhost:23122] [ 5] /lib64/libc.so.6(__libc_start_main+0xfd) [0x3fcea1ecdd] [localhost:23122] [ 6] mdrun_mpi() [0x406ee9] [localhost:23122] *** End of error message *** Segmentation fault (core dumped) -- gmx-users mailing listgmx-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
Re: [gmx-users] Re: g_analyze
On 11/10/13 12:20 AM, bharat gupta wrote: Hi, I used the command g_hbond to find h-bond between residues 115-118 and water. Then I used g_analyze to find out the average and it gives the value for the hbonds like this :- std. dev.relative deviation of standard - cumulants from those of set average deviation sqrt(n-1) a Gaussian distribition cum. 3 cum. 4 SS1 6.877249e-02 2.546419e-01 5.092839e-03 2.1813.495 SS2 6.997201e-02 2.673450e-01 5.346901e-03 2.4215.001 When I calculated the average manually, by taking the average of numbers in second column of hbnum.xvg file, I got a value of around 13.5.. What is the reason for such a large difference. Hard to say, but I've never known g_analyze to be wrong, so I'd suspect something is amiss in your manual calculation. The difference between 13.5 and 0.0069 is huge; you should be able to scan through the data file to see what the expected value should be. In another case, g_analyze gives avg values of aroun 6.9 for hbond between two residues and when I calculated it maually I got the avg values as 6.8 .. Whats the meaning of SS1 and SS2,?? Does it mean that SS1 refers to time and SS2 refers to hbond numbers in the hbnum.xvg obtained from g_hbond analysis ?? Data sets 1 and 2. You will note that there are two columns of data in the -hbnum output produced by g_hbond, with titles explaining both. -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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
Re: [gmx-users] Umbrella Sampling tutorial
On 11/10/13 5:31 AM, shahab shariati wrote: Dear Justin Thanks for your reply. You are right. I should not extrapolate too literally from your tutorial to my system. But, I have a general question: There is 2 groups in COM pulling method (reference group + pull group). If I want to use pull_geometry = distance, so, I should fix reference group to be immobile. Is it true? What you described earlier should not be attempted with distance geometry. It won't work very well. The use of restraints is almost NEVER necessary, especially in the case where the reference group is much more massive than the pulling group. On the other hand, I want to know exactly using position restraining on reference group is optional or mandatory in COM pulling method? Almost never used. More explicitly - you do not need position restraints. -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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
Re: [gmx-users] problem in running mdrun command
Hi, There's nothing GROMACS-specific here - something about your MPI installation, configuration or use is pretty wrong, but we can't help work out what. Mark On Sun, Nov 10, 2013 at 12:31 PM, S.Chandra Shekar chandrashe...@iisertvm.ac.in wrote: Dear all I encounter a problem while running command mdrun_mpi -v -deffnm em in gromacs. I am new to the gromacs. i just ran test calculation, *simulation of lyzozyme in water*. i am able to generate gro, tpr files. But in the final step i got following error. Thanks in advance. [localhost.localdomain:23122] mca: base: component_find: paffinity mca_paffinity_linux uses an MCA interface that is not recognized (component MCA v1.0.0 != supported MCA v2.0.0) -- ignored [localhost.localdomain:23123] mca: base: component_find: paffinity mca_paffinity_linux uses an MCA interface that is not recognized (component MCA v1.0.0 != supported MCA v2.0.0) -- ignored [localhost.localdomain:23123] mca: base: component_find: ras mca_ras_dash_host uses an MCA interface that is not recognized (component MCA v1.0.0 != supported MCA v2.0.0) -- ignored [localhost.localdomain:23123] mca: base: component_find: ras mca_ras_gridengine uses an MCA interface that is not recognized (component MCA v1.0.0 != supported MCA v2.0.0) -- ignored [localhost.localdomain:23123] mca: base: component_find: ras mca_ras_localhost uses an MCA interface that is not recognized (component MCA v1.0.0 != supported MCA v2.0.0) -- ignored [localhost.localdomain:23123] mca: base: component_find: errmgr mca_errmgr_hnp uses an MCA interface that is not recognized (component MCA v1.0.0 != supported MCA v2.0.0) -- ignored [localhost.localdomain:23123] mca: base: component_find: errmgr mca_errmgr_orted uses an MCA interface that is not recognized (component MCA v1.0.0 != supported MCA v2.0.0) -- ignored [localhost.localdomain:23123] mca: base: component_find: errmgr mca_errmgr_proxy uses an MCA interface that is not recognized (component MCA v1.0.0 != supported MCA v2.0.0) -- ignored [localhost.localdomain:23123] mca: base: component_find: iof mca_iof_proxy uses an MCA interface that is not recognized (component MCA v1.0.0 != supported MCA v2.0.0) -- ignored [localhost.localdomain:23123] mca: base: component_find: iof mca_iof_svc uses an MCA interface that is not recognized (component MCA v1.0.0 != supported MCA v2.0.0) -- ignored [localhost.localdomain:23122] mca: base: component_find: rcache mca_rcache_rb uses an MCA interface that is not recognized (component MCA v1.0.0 != supported MCA v2.0.0) -- ignored [localhost:23122] *** Process received signal *** [localhost:23122] Signal: Segmentation fault (11) [localhost:23122] Signal code: Address not mapped (1) [localhost:23122] Failing at address: 0x4498 [localhost:23122] [ 0] /lib64/libpthread.so.0() [0x3fcee0f500] [localhost:23122] [ 1] /usr/local/lib/libmpi.so.1(PMPI_Comm_size+0x4e) [0x2acc6d93727e] [localhost:23122] [ 2] /usr/local/gromacs/bin/../lib/libgmx_mpi.so.8(gmx_setup+0x32) [0x2acc6d195e02] [localhost:23122] [ 3] /usr/local/gromacs/bin/../lib/libgmx_mpi.so.8(init_par+0x51) [0x2acc6d234251] [localhost:23122] [ 4] mdrun_mpi(cmain+0x11f9) [0x435799] [localhost:23122] [ 5] /lib64/libc.so.6(__libc_start_main+0xfd) [0x3fcea1ecdd] [localhost:23122] [ 6] mdrun_mpi() [0x406ee9] [localhost:23122] *** End of error message *** Segmentation fault (core dumped) -- gmx-users mailing listgmx-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 -- gmx-users mailing listgmx-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
[gmx-users] Umbrella Sampling tutorial
Dear Justin Very thanks for your reply. What you described earlier should not be attempted with distance geometry. It won't work very well. The use of restraints is almost NEVER necessary, especially in the case where the reference group is much more massive than the pulling group. I want to calculate Potential of mean force as a function of the distance between the centers of mass of drug and the lipid bilayer. You said distance geometry won't work very well in my case. What is your better suggestion about my case? Best wishes -- gmx-users mailing listgmx-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
Re: [gmx-users] Umbrella Sampling tutorial
On 11/10/13 9:14 AM, shahab shariati wrote: Dear Justin Very thanks for your reply. What you described earlier should not be attempted with distance geometry. It won't work very well. The use of restraints is almost NEVER necessary, especially in the case where the reference group is much more massive than the pulling group. I want to calculate Potential of mean force as a function of the distance between the centers of mass of drug and the lipid bilayer. You said distance geometry won't work very well in my case. What is your better suggestion about my case? http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin/gmx-tutorials/umbrella/05a_pull_tips.html pull_geometry = position -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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
[gmx-users] Bilayer COM removal issue: Large VCM
Hi All, I am experiencing a few problems in membrane simulations wrt COM removal. I downloaded a 400 ns pre-equilibrated Slipid-DMPC membrane with all the accompanying files. I then carried out the following steps: 1) energy minimization 2) NVT Eq - 100 ps 3) NPT Eq - 250 ps (Berendsen temp, Pres coupling) Then I used g_select to select the upper and lower DMPC leaflets. The then carried out a 250 ps NPT eq again. The only change was: comm-grps= SOL DMPC ==comm-grps = SOL upper lower On every step in log file, I get the following message: *Step Time Lambda 124000 248.00.0Large VCM(group lower): -0.00051, -0.00515, -0.00652, Temp-cm: 8.11828e+29 Energies (kJ/mol)U-BProper Dih. Improper Dih. LJ-14 Coulomb-147.23818e+044.19778e+046.46641e+024.54801e+03 -1.45245e+05LJ (SR)LJ (LR) Disper. corr. Coulomb (SR) Coul. recip.2.79689e+04 -3.78407e+03 -2.10679e+03 -5.84134e+05 -8.87497e+04 PotentialKinetic En. Total EnergyTemperature Pres. DC (bar) -6.76497e+051.76468e+05 -5.00029e+05 3.10424e+02 -1.05704e+02 Pressure (bar) Constr. rmsd -1.85927e+02 6.42934e-06* *Large VCM(group lower): -0.00187, -0.00369, 0.00032, Temp-cm: 2.02076e+29Large VCM(group lower): -0.00725, -0.00278, -0.00549, Temp-cm: 1.05988e+30Large VCM(group lower): 0.00020, 0.00308, -0.00176, Temp-cm: 1.48126e+29Large VCM(group lower): -0.00541, 0.00546, -0.00166, Temp-cm: 7.24656e+29Large VCM(group lower): -0.00220, 0.00362, -0.00741, Temp-cm: 8.53812e+29Large VCM(group lower): 0.00140, -0.00160, 0.00029, Temp-cm: 5.39679e+28Large VCM(group lower): -0.00056, -0.00293, -0.00364, Temp-cm: 2.59422e+29Large VCM(group lower): -0.00172, -0.00260, 0.00494, Temp-cm: 3.99945e+29Large VCM(group lower): 0.00252, 0.00594, 0.00068, Temp-cm: 4.93342e+29* *DD step 124999 vol min/aver 0.702 load imb.: force 1.3% pme mesh/force 0.636* I do not know what to make of it. There are no issues when I remove COM for the entire system. I have seen this issue come up a few times in the archives too, but I didn't find a satisfactory solution since the bilayer was very well equilibrated. I would appreciate any suggestions. Thank you. -- Rajat Desikan (Ph.D Scholar) Prof. K. Ganapathy Ayappa's Lab (no 13), Dept. of Chemical Engineering, Indian Institute of Science, Bangalore -- gmx-users mailing listgmx-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
Re: [gmx-users] Re: g_analyze
I checked the file hbnum.xvg file and it contains three columns - time, hbonds, hbonds that donot follow the angle criteria. In that case SS1 is the average of actual hbonds (2nd column ) and SS2 is the average of 3rd column. Am I right here or not ?? I tried to calculate the h-bond for residues 115-118 individually, and then checked the average for each residue. For single residue calculation, the g_analyze average value is correct. But when I calculate the h-bond as a range 115-118, I get the g_analyze value as 1.62 . I calculated the average manually in excel, got the average values as 16.2 [which is (g_analyze avg value)/10]. I then added up the average values of h-bonds of individual residues and the final comes around 16.2, same as that of the 115-118 range h-bonds. This means that my calculation is correct. I also used trjorder to calculate h-bond at distance 0.34 for residues 115-118. I got the average value around 2.51 from g_analyze, where as manual calculation gives 25.1. I don't knw why for the range the g_analyze give avg as (actual avg value)/10 ?? Why does trjorder and g_hbond gives different number of hydrogen bonds for the same residue set?? Thanks --- BHARAT On Sun, Nov 10, 2013 at 10:01 PM, Justin Lemkul jalem...@vt.edu wrote: On 11/10/13 12:20 AM, bharat gupta wrote: Hi, I used the command g_hbond to find h-bond between residues 115-118 and water. Then I used g_analyze to find out the average and it gives the value for the hbonds like this :- std. dev.relative deviation of standard - cumulants from those of set average deviation sqrt(n-1) a Gaussian distribition cum. 3 cum. 4 SS1 6.877249e-02 2.546419e-01 5.092839e-03 2.1813.495 SS2 6.997201e-02 2.673450e-01 5.346901e-03 2.4215.001 When I calculated the average manually, by taking the average of numbers in second column of hbnum.xvg file, I got a value of around 13.5.. What is the reason for such a large difference. Hard to say, but I've never known g_analyze to be wrong, so I'd suspect something is amiss in your manual calculation. The difference between 13.5 and 0.0069 is huge; you should be able to scan through the data file to see what the expected value should be. In another case, g_analyze gives avg values of aroun 6.9 for hbond between two residues and when I calculated it maually I got the avg values as 6.8 .. Whats the meaning of SS1 and SS2,?? Does it mean that SS1 refers to time and SS2 refers to hbond numbers in the hbnum.xvg obtained from g_hbond analysis ?? Data sets 1 and 2. You will note that there are two columns of data in the -hbnum output produced by g_hbond, with titles explaining both. -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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 -- gmx-users mailing listgmx-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
[gmx-users] Re: Thankful
Justin, thank you very much for your kind help about LIE and PME Williams -- gmx-users mailing listgmx-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
Re: [gmx-users] Re: g_analyze
On 11/10/13 7:18 PM, bharat gupta wrote: I checked the file hbnum.xvg file and it contains three columns - time, hbonds, hbonds that donot follow the angle criteria. In that case SS1 is The third column is not actually H-bonds, then ;) the average of actual hbonds (2nd column ) and SS2 is the average of 3rd column. Am I right here or not ?? Yes. I tried to calculate the h-bond for residues 115-118 individually, and then checked the average for each residue. For single residue calculation, the g_analyze average value is correct. But when I calculate the h-bond as a range 115-118, I get the g_analyze value as 1.62 . I calculated the average manually in excel, got the average values as 16.2 [which is (g_analyze avg value)/10]. That is impossible. You cannot get a different average by examining the same numbers. Read the g_analyze output again - I am willing to bet that you're not seeing the exponent of the scientific notation. I then added up the average values of h-bonds of individual residues and the final comes around 16.2, same as that of the 115-118 range h-bonds. This means that my calculation is correct. I also used trjorder to calculate h-bond at distance 0.34 for residues 115-118. I got the average value around 2.51 from g_analyze, where as manual calculation gives 25.1. I don't knw why for the range the g_analyze give avg as (actual avg value)/10 ?? Why does trjorder and g_hbond gives different number of hydrogen bonds for the same residue set?? All of this comes down to correctly reading the screen output. I have no idea what you're doing with trjorder, though. It doesn't measure H-bonds. -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- gmx-users mailing listgmx-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