[gmx-users] fix a group in truncated octahedron
Hi GMX users, I want to fix a group in a truncate octahedron. How can I dealt with the question below, There is One spherical rigid body consists of face-center cubic lattices fixed in the center of the box. I do not want to calculate the force and energy between the paritcles of this rigid body, so that no matter how large force between them shoud not blow up the rigid body. How can I do this? Thanks in advance. -- wende -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/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/mailing_lists/users.php
Re: [gmx-users] fix a group in truncated octahedron
Wende, This is the third time you post this exact message. Apparently, nobody felt equipped to reply to your post, either because nobody knows the answer, or your message and intent were not clear enough. In either case, reposting your message over and over is completely useless and annoying. Please read http://catb.org/~esr/faqs/smart-questions.html Tsjerk On Tue, Dec 29, 2009 at 9:07 AM, lammps lammps lammp2fo...@gmail.com wrote: Hi GMX users, I want to fix a group in a truncate octahedron. How can I dealt with the question below, There is One spherical rigid body consists of face-center cubic lattices fixed in the center of the box. I do not want to calculate the force and energy between the paritcles of this rigid body, so that no matter how large force between them shoud not blow up the rigid body. How can I do this? Thanks in advance. -- wende -- gmx-users mailing list gmx-us...@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/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/mailing_lists/users.php -- Tsjerk A. Wassenaar, Ph.D. Computational Chemist Medicinal Chemist Neuropharmacologist -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/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/mailing_lists/users.php
[gmx-users] linc warnings
Hi My protein has a Fe4S4 cluster. I have been able to minimize the structure using Steepest descent. The log file shows following message: Steepest Descents converged to Fmax 2000 in 372 steps Potential Energy = -4.1580605e+06 Maximum force = 1.7402682e+03 on atom 820 Norm of force = 3.0794989e+01 - But during the position restraining step, following Linc warning is coming before stopping all-together Step 228, time 0.456 (ps) LINCS WARNING relative constraint deviation after LINCS: rms 51159085.131913, max 1100381184.00 (between atoms 5437 and 5444) bonds that rotated more than 30 degrees: atom 1 atom 2 angle previous, current, constraint length 1262 1264 38.80.1009 0.1054 0.1000 1262 1263 83.60.1035 0.0960 0.1000 1261 1265 42.10.1349 0.1366 0.1340 1259 1261 42.90.1359 0.1372 0.1340 1237 1241 33.00.1480 0.1788 0.1470 1235 1237 75.60.1360 0.4797 0.1330 1235 1236 40.90.1254 0.3422 0.1230 1232 1233 71.30.1834 0.3421 0.1830 1231 1232 43.10.1559 0.3294 0.1530 1230 1235 106.00.1589 0.2070 0.1530 1230 1231 111.00.1606 0.3819 0.1530 1228 1230 81.90.1535 8.3688 0.1470 1228 1229 100.80.1040 7.2523 0.1000 1226 1228 108.70.1480 29.6051 0.1330 1226 1227 112.00.1271 25.3002 0.1230 1224 1226 67.70.2210 50.5775 0.1530 1224 1225 73.60.2217 43.8606 0.1530 1222 1224 146.53.0118 273.3403 0.1470 1222 1223 173.15.0766 350.6593 0.1000 1220 1222 140.6 50.1464 1239.3472 0.1330 1220 1221 142.5 50.5481 1210.8007 0.1230 1218 1219 46.50.1089 0.0687 0.1090 1216 1218 106.30.1532 2.4558 0.1390 1216 1217 57.30.1434 2.4980 0.1090 1214 1218 106.40.1028 2.6673 0.1390 1214 1215 166.40.2538 2.1798 0.1090 1212 1216 119.00.5844 15.3761 0.1390 1212 1213 50.70.5406 17..8588 0.1090 1210 1214 116.20.9490 16.1642 0.1390 1210 1211 150.40.7387 23.8303 0.1090 1209 1212 57.44.3532 165.6521 0.1390 1209 1210 53.54.3615 166.9433 0.1390 1208 1209 130.9 34.7441 858.7278 0.1530 1207 1220 174.9 199.1050 3552.6362 0.1530 1207 1208 172.1 195.9491 3367.8894 0.1530 1205 1207 176.1 490.9904 5776.2954 0.1470 1205 1206 177.8 631.7326 5326.1118 0.1000 1203 1205 179.0 1242.0747 9168.7764 0.1330 1203 1204 178.9 1057.3107 7388.7915 0.1230 1202 1203 179.3 1221.3862 8532.1836 0.1530 1200 1202 177.8 455.7037 3702.9392 0.1470 1200 1201 172.7 158.7020 1978.1310 0.1000 1198 1200 174.2 160.3432 1536.1411 0.1330 1198 1199 92.0 16.9060 423.7516 0.1230 1197 1198 94.1 16.9010 421.8028 0.1530 1195 1197 99.1 10.3514 9.3388 0.1470 1195 1196 96.29.0920 8.6023 0.1000 1193 1195 91.2 21.5753 58.1620 0..1330 1193 1194 89.2 21.1892 62.7649 0.1230 1191 1192 90.1 50.3829 12.4455 0.1330 1190 1191 89.4 13.1419 10.9133 0.1830 1189 1193 114.3 75.7546 37.4387 0.1530 1189 1190 91.1 24.9488 56.3341 0.1530 1187 1189 146.9 53.8226 42.4679 0.1470 1187 1188 139.0 22.0031 100.4932 0.1000 1185 1187 163.3 21.1472 184.6258 0.1330 1185 1186 170.3 61.0351 169.5860 0.1230 1182 1184 153.34.0050 30.0847 0.1530 1182 1183 109.92.7482 31.3314 0.1530 1181 1185 170.1 69.8200 192.5774 0.1530 1181 1182 129.0 21.4261 109.4787 0.1530 1179 1181 132.8 19.8073 116.5793 0.1470 1179 1180 159.42.3810 48.3401 0..1000 1177 1179 175.21.4786 53.0837 0.1330 1177 1178 166.40.7559 28.5291 0.1230 1174 1176 90.00.1801 80.6845 0.1530 1174 1175 99.70.2091 19.2129 0.1530 1173 1177 149.10.2721 28.0611 0.1530 1173 1174 73.60.7968 36.5323 0.1530 1171 1173 87.80.8744 115.1496 0.1470 1171 1172 91.20.1603 111.8982 0.1000 1169 1171 89.70.2456 106.5120 0.1330 1169 1170 90.00.1499 17.6986 0.1230 1166 1168 76.10.1256 0.4912 0.1250 1166 1167 77.80.1256 0.4782 0.1250 1165 1166 71.90.1546 2.7698 0.1530 1164 1169 83.50.1781 17.7449 0.1530 1164 1165 89.30.1594 8.0829 0..1530 1162 1164 88.00.1531 8.2075 0.1470 1162 1163 75.00.1015 0.7203 0.1000 1160 1162 77.50.1343 1.0789
Re: [gmx-users] Problems with umbrella sampling
Thanks Justin, it seems to help, i have another general question that you might help me with: should i define the distances for the umbrella sampling run in some file and then run it?(i saw something about *.ppa file but couldn't find any further instructions in the mannual..) will this file apply for the WHAM analysis after the run? thanks, Amir On Mon, Dec 28, 2009 at 6:14 PM, Justin A. Lemkul jalem...@vt.edu wrote: Amir Marcovitz wrote: Hi, my system is constructed of 2 parallel plates in a box with solvent. each plate is made of 36 atoms which are positively charged on one of them (SRP) and negatively charged on the other plate (SRN). i want to perform a simulation with umbrella sampling between the 2 so i defined the pulling section parameters in the .mdp parameter file as following: ; COM PULLING ; Pull type: no, umbrella, constraint or constant_force pull = umbrella ; Pull geometry: distance, direction, cylinder or position pull_geometry= distance ; Select components for the pull vector. default: Y Y Y pull_dim = Y Y Y ; Cylinder radius for dynamic reaction force groups (nm) pull_r1 = 1 ; Switch from r1 to r0 in case of dynamic reaction force pull_r0 = 1.5 pull_constr_tol = 1e-06 pull_start = no pull_nstxout = 10 pull_nstfout = 1 ; Number of pull groups pull_ngroups = 2 ; Group name, weight (default all 1), vector, init, rate (nm/ps), kJ/(mol*nm^2) pull_group0 = SRP pull_weights0= 1 pull_pbcatom0= 0 pull_group1 = SRN pull_weights1= 1 pull_pbcatom1= 0 pull_vec1= 0.0 1.0 0.0 pull_init1 = 1.5 pull_rate1 = 0 pull_k1 = 1 pull_kB1 = 0 when proccesing the file with grompp i get the following error: *Fatal error: Number of weights (1) for pull group 0 'SRP' does not match the number of atoms (36)* is someone recognizing my mistake? Please refer to the manual (manual.gromacs.org is quite handy), you will find: Optional relative weights which are multiplied with the masses of the atoms to give the total weight for the COM. The number should be 0, meaning all 1, or the number of atoms in the pull group. I also think your value for pull_ngroup is wrong. It appears you are pulling SRN with respect to SRP, so you only have one pull group, not two. -Justin does someone has an experience with umbrella sampling in GROMACS? thanks' amir * * -- Justin A. Lemkul Ph.D. Candidate ICTAS Doctoral Scholar MILES-IGERT Trainee Department of Biochemistry Virginia Tech Blacksburg, VA jalemkul[at]vt.edu | (540) 231-9080 http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/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/mailing_lists/users.php -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/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/mailing_lists/users.php
Re: [gmx-users] Problems with umbrella sampling
Amir Marcovitz wrote: Thanks Justin, it seems to help, i have another general question that you might help me with: should i define the distances for the umbrella sampling run in some file and then run it?(i saw something about *.ppa file but couldn't find any further instructions in the mannual..) will this file apply for the WHAM analysis after the run? All settings are defined in the .mdp file (see, i.e. pull_start and pull_init). Using .ppa files is obsolete (version 3.x and before). The WHAM analysis later reads these distances from the .tpr files passed to g_wham. There are many discussions (some quite recent) in the list archive about proper settings. I suggest you have a look at the manual and then go searching. -Justin thanks, Amir On Mon, Dec 28, 2009 at 6:14 PM, Justin A. Lemkul jalem...@vt.edu mailto:jalem...@vt.edu wrote: Amir Marcovitz wrote: Hi, my system is constructed of 2 parallel plates in a box with solvent. each plate is made of 36 atoms which are positively charged on one of them (SRP) and negatively charged on the other plate (SRN). i want to perform a simulation with umbrella sampling between the 2 so i defined the pulling section parameters in the .mdp parameter file as following: ; COM PULLING ; Pull type: no, umbrella, constraint or constant_force pull = umbrella ; Pull geometry: distance, direction, cylinder or position pull_geometry= distance ; Select components for the pull vector. default: Y Y Y pull_dim = Y Y Y ; Cylinder radius for dynamic reaction force groups (nm) pull_r1 = 1 ; Switch from r1 to r0 in case of dynamic reaction force pull_r0 = 1.5 pull_constr_tol = 1e-06 pull_start = no pull_nstxout = 10 pull_nstfout = 1 ; Number of pull groups pull_ngroups = 2 ; Group name, weight (default all 1), vector, init, rate (nm/ps), kJ/(mol*nm^2) pull_group0 = SRP pull_weights0= 1 pull_pbcatom0= 0 pull_group1 = SRN pull_weights1= 1 pull_pbcatom1= 0 pull_vec1= 0.0 1.0 0.0 pull_init1 = 1.5 pull_rate1 = 0 pull_k1 = 1 pull_kB1 = 0 when proccesing the file with grompp i get the following error: *Fatal error: Number of weights (1) for pull group 0 'SRP' does not match the number of atoms (36)* is someone recognizing my mistake? Please refer to the manual (manual.gromacs.org http://manual.gromacs.org/ is quite handy), you will find: Optional relative weights which are multiplied with the masses of the atoms to give the total weight for the COM. The number should be 0, meaning all 1, or the number of atoms in the pull group. I also think your value for pull_ngroup is wrong. It appears you are pulling SRN with respect to SRP, so you only have one pull group, not two. -Justin does someone has an experience with umbrella sampling in GROMACS? thanks' amir * * -- Justin A. Lemkul Ph.D. Candidate ICTAS Doctoral Scholar MILES-IGERT Trainee Department of Biochemistry Virginia Tech Blacksburg, VA jalemkul[at]vt.edu http://vt.edu/ | (540) 231-9080 http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin -- gmx-users mailing listgmx-users@gromacs.org mailto:gmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/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 mailto:gmx-users-requ...@gromacs.org. Can't post? Read http://www.gromacs.org/mailing_lists/users.php -- Justin A. Lemkul Ph.D. Candidate ICTAS Doctoral Scholar MILES-IGERT Trainee Department of Biochemistry Virginia Tech Blacksburg, VA jalemkul[at]vt.edu | (540) 231-9080 http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/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
Re: [gmx-users] linc warnings
subarna thakur wrote: Hi My protein has a Fe4S4 cluster. I have been able to minimize the structure using Steepest descent. The log file shows following message: Steepest Descents converged to Fmax 2000 in 372 steps Potential Energy = -4.1580605e+06 Maximum force = 1.7402682e+03 on atom 820 Norm of force = 3.0794989e+01 - But during the position restraining step, following Linc warning is coming before stopping all-together If EM worked nicely, then something else you are doing is incorrect. snip the atoms 5437,5438,5439 are the FE atoms of the Fe4S4 cluster. these atoms are rotating hugely and showing large rms. Please suggest how to tackle the FE atoms during simulation run ? Without the metal cluster, the same protein is working fine under simulation and there is no such linc warinings? Without seeing your .mdp file and perhaps even snippets of your topology pertaining to the Fe4S4 cluster, there's nothing else to say except please read the standard advice: http://www.gromacs.org/Documentation/Errors#LINCS.2fSETTLE.2fSHAKE_warnings -Justin Subarna The INTERNET now has a personality. YOURS! See your Yahoo! Homepage http://in.rd.yahoo.com/tagline_yyi_1/*http://in.yahoo.com/. -- Justin A. Lemkul Ph.D. Candidate ICTAS Doctoral Scholar MILES-IGERT Trainee Department of Biochemistry Virginia Tech Blacksburg, VA jalemkul[at]vt.edu | (540) 231-9080 http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/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/mailing_lists/users.php
[gmx-users] mpi code in the gromacs utilities
Hi, I wonder to know which of the gromacs programs except mdrun support parallel execution? Occasionally I tried g_energy and it seems not to support MPI, so it looks like the each node performs the same job... -- Vitaly V. Chaban, Ph.D. http://www-rmn.univer.kharkov.ua/chaban.html -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/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/mailing_lists/users.php
Re: [gmx-users] mpi code in the gromacs utilities
Vitaly V. Chaban wrote: Hi, I wonder to know which of the gromacs programs except mdrun support parallel execution? Occasionally I tried g_energy and it seems not to support MPI, so it looks like the each node performs the same job... That's because only mdrun is MPI-enabled. I believe some framework is in place in some of the memory-intensive tools, but none of them can actually be executed using MPI. -Justin -- Vitaly V. Chaban, Ph.D. http://www-rmn.univer.kharkov.ua/chaban.html -- Justin A. Lemkul Ph.D. Candidate ICTAS Doctoral Scholar MILES-IGERT Trainee Department of Biochemistry Virginia Tech Blacksburg, VA jalemkul[at]vt.edu | (540) 231-9080 http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/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/mailing_lists/users.php
Re: [gmx-users] linc warnings
Hi, I am using the following keywords in pr.mdp file for the simulation for protein with fe4s4 cluster; title = Yo cpp = /usr/bin/cpp define = -DPOSRES constraints = all-bonds integrator = md dt = 0.002; ps ! nsteps = 150; total 3000 ps. nstcomm = 1 nstxout = 500 nstvout = 1000 nstfout = 0 nstlog = 10 nstenergy = 10 nstlist = 10 ns_type = grid rlist = 1.0 rcoulomb= 1.0 rvdw= 1.0 ; Berendsen temperature coupling is on in two groups Tcoupl = berendsen tc-grps= ProteinNon-Protein tau_t = 0.10.1 ref_t = 300300 ; Energy monitoring energygrps= Protein Non-Protein ; Pressure coupling is not on Pcoupl = no tau_p = 0.5 compressibility = 4.5e-5 ref_p = 1.0 ; Generate velocites is on at 300 K. gen_vel = yes gen_temp= 300.0 gen_seed= 173529 But the Fe4s4 cluster atoms are showing large rms and rotating hugely in position restraining step and Linc warnings are comming. In some published papers it is mentioned that during the molecular dynamics simulations, the proteins' heavy atoms and the counter ions were harmonically restrained to their initial positions of pdb, How can I do this for my fe4s4 cluster. What keywords should I use for that and in which step I can do? subarna From: Justin A. Lemkul jalem...@vt.edu To: Discussion list for GROMACS users gmx-users@gromacs.org Sent: Tue, 29 December, 2009 5:44:03 PM Subject: Re: [gmx-users] linc warnings subarna thakur wrote: Hi My protein has a Fe4S4 cluster. I have been able to minimize the structure using Steepest descent. The log file shows following message: Steepest Descents converged to Fmax 2000 in 372 steps Potential Energy = -4.1580605e+06 Maximum force = 1.7402682e+03 on atom 820 Norm of force = 3.0794989e+01 - But during the position restraining step, following Linc warning is coming before stopping all-together If EM worked nicely, then something else you are doing is incorrect. snip the atoms 5437,5438,5439 are the FE atoms of the Fe4S4 cluster. these atoms are rotating hugely and showing large rms. Please suggest how to tackle the FE atoms during simulation run ? Without the metal cluster, the same protein is working fine under simulation and there is no such linc warinings? Without seeing your .mdp file and perhaps even snippets of your topology pertaining to the Fe4S4 cluster, there's nothing else to say except please read the standard advice: http://www.gromacs.org/Documentation/Errors#LINCS.2fSETTLE.2fSHAKE_warnings -Justin Subarna The INTERNET now has a personality. YOURS! See your Yahoo! Homepage http://in.rd.yahoo.com/tagline_yyi_1/*http://in.yahoo.com/. -- Justin A. Lemkul Ph.D. Candidate ICTAS Doctoral Scholar MILES-IGERT Trainee Department of Biochemistry Virginia Tech Blacksburg, VA jalemkul[at]vt.edu | (540) 231-9080 http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/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/mailing_lists/users.php The INTERNET now has a personality. YOURS! See your Yahoo! Homepage. http://in.yahoo.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/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/mailing_lists/users.php