<br><br>----- Original Message -----<br>From: Sascha Rehm &lt;Sascha.Rehm@itb.uni-stuttgart.de&gt;<br>Date: Tuesday, September 21, 2010 22:34<br>Subject: [gmx-users] Problems with coupling groups in gas phase simulations<br>To: gmx-users@gromacs.org<br><br>&gt; Dear Gromacs-User (and Developers),<br>&gt; <br>&gt; I tried to set up and simulate some systems with 1000 Argon <br>&gt; atoms and roughly a thousand water molecules that were randomly <br>&gt; distributed (like gaseous water in argon gas).<br>&gt; The starting systems are set up with OPLS/AA and 100nm box side <br>&gt; length, so the systems are much too large and are expected to <br>&gt; shrink in the first few picoseconds.<br><br>Generating velocities at 5K, annealing to 318K in 2ps at NPT at a density that's far too low is being pretty rough on the system. I'm surprised the integration didn't break completely. I've no idea what your objective is, but I wouldn't expect to be able to make equilibrium observations for many nanoseconds :-)<br><br>&gt; Two different coupling groups were used (Ar and SOL) with two <br>&gt; different compressibilities (1.27e-2 for Argon and 5.3e-5 for Water).<br><br>T-coupling works best in the limit of large groups of atoms of similar heat capacity. Standard MD advice is not to couple (say) ions in water to separate groups, because the resulting coupling leads to artifacts. This would seem to apply even more strongly to a mixture of components such as you have.<br><br>&gt; The systems are heated up in the first ns to 318K and then <br>&gt; continued for another 100 ns. During that time the boxes shrink <br>&gt; only to a box size of 99.5 nm, leading to a totally wrong density.<br>&gt; <br>&gt; When I set up the systems with only one couling group ("System") <br>&gt; the systems shrink to roughly 38nm and show a correct density in <br>&gt; a few hundred picoseconds but then the density starts to <br>&gt; oszilate every few ns. I guess this happens because the <br>&gt; compressibility is set to 1.2e-2 and the water starts to form <br>&gt; droplets (which should have a compressibility of 5.3e-5)? <br>&gt; Increasing tau_p from 0.6 to 2.0 does not show big effects on <br>&gt; the oszilation.<br>&gt; <br>&gt; 1. Is there an explanation why the systems behave strange when <br>&gt; using 2 coupling groups? My search on the mailing list did not <br>&gt; lead to an answer (also, there are only a few topics on <br>&gt; simulating gas or condensing water droplets)<br><br>Yes, see above.<br><br>&gt; 2. Are there reasons, why I should not use only one coupling <br>&gt; group for the whole system? E.g. in protein-solvent simulations, <br>&gt; the manual says that using only one coulpling group leads to an <br>&gt; increase of temperature in the protein and a decrease in the solvent.<br><br>There, the issue is a difference in heat capacity between two groups of notable size. Search literature for "hot solute cold solvent" for discussion.<br><br>Since you're simulating a system that's intended to start mixed and then do something, this situation cannot arise until (perhaps) your phases separate.<br><br>Mark<br><br>&gt; Hopefully someone can help me out of that problem...<br>&gt; <br>&gt; Best Regards,<br>&gt; <br>&gt; Sascha Rehm<br>&gt; <br>&gt; ---------------------------------------------------------<br>&gt; The used MDP file looks like the following:<br>&gt; <br>&gt; integrator&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;        = md<br>&gt; dt&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.002<br>&gt; nsteps&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 5000000<br>&gt; <br>&gt; annealing = single<br>&gt; annealing_npoints = 3<br>&gt; annealing_time = 0 100 5000000<br>&gt; annealing_temp = 5 318.15 318.15<br>&gt; <br>&gt; nstxout&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&gt; nstvout&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&gt; <br>&gt; nstenergy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1000<br>&gt; nstxtcout&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1000<br>&gt; <br>&gt; ;group(s) to write to energy file<br>&gt; energygrps&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = system<br>&gt; <br>&gt; nstlist&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 10<br>&gt; ns_type&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = grid<br>&gt; <br>&gt; ;cut-off distance for the short-range neighbor list<br>&gt; rlist&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1.5<br>&gt; <br>&gt; ;treatment of electrostatic interactions<br>&gt; coulombtype&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = Cut-off<br>&gt; rcoulomb&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1.5<br>&gt; <br>&gt; ;treatment of van der waals interactions<br>&gt; vdwtype&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = cut-off<br>&gt; rvdw&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1.5<br>&gt; <br>&gt; ; Periodic boudary conditions in all the directions<br>&gt; pbc&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = xyz<br>&gt; <br>&gt; ;Temperature coupling<br>&gt; tcoupl&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = berendsen<br>&gt; tc-<br>&gt; grps&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = system<br>&gt; tau_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.1<br>&gt; ref_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 318.15<br>&gt; <br>&gt; ;Pressure coupling<br>&gt; Pcoupl&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = berendsen<br>&gt; Pcoupltype&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = isotropic<br>&gt; tau_p&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.6<br>&gt; compressibility&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1.27e-2<br>&gt; ref_p&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1.0<br>&gt; <br>&gt; ;Velocity generation<br>&gt; gen_vel&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = yes<br>&gt; gen_temp&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 5<br>&gt; gen_seed&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 173529<br>&gt; <br>&gt; <br>&gt; <br>&gt; <br>&gt; <br>&gt; -- <br>&gt; Dipl. Bioinf. Sascha Rehm<br>&gt; Institut fuer Technische Biochemie<br>&gt; Allmandring 31<br>&gt; D-70569 Stuttgart<br>&gt; http://www.itb.uni-stuttgart.de/<br>&gt; -- <br>&gt; gmx-users mailing list&nbsp;&nbsp;&nbsp; gmx-users@gromacs.org<br>&gt; http://lists.gromacs.org/mailman/listinfo/gmx-users<br>&gt; Please search the archive at <br>&gt; http://www.gromacs.org/Support/Mailing_Lists/Search before posting!<br>&gt; Please don't post (un)subscribe requests to the list. Use the <br>&gt; www interface or send it to gmx-users-request@gromacs.org.<br>&gt; Can't post? Read http://www.gromacs.org/Support/Mailing_Lists