<br><br>----- Original Message -----<br>From: Nilesh Dhumal &lt;ndhumal@andrew.cmu.edu&gt;<br>Date: Tuesday, October 26, 2010 10:50<br>Subject: Re: [gmx-users] (no subject)<br>To: Discussion list for GROMACS users &lt;gmx-users@gromacs.org&gt;<br><br>&gt; I run a test simulation for -rerun. I didn't change the topology file.<br>&gt; <br>&gt; grompp -f md.mdp&nbsp; -c&nbsp; solvent-bmi-pf6-128.pdb&nbsp; -<br>&gt; p&nbsp; solvent-bmi-pf6-128.top<br>&gt; &nbsp;-o 3.tpr<br>&gt; mpirun -machinefile cp -np 8 mdrun -s 3.tpr -o 3.trr&nbsp;&nbsp; <br>&gt; -c<br>&gt; solvent-bmi-pf6-128.pdb -e 3.edr -g 3.log<br>&gt; <br>&gt; with -rerun<br>&gt; grompp -f md.mdp&nbsp; -c&nbsp; solvent-bmi-pf6-128.pdb&nbsp; -<br>&gt; p&nbsp; solvent-bmi-pf6-128.top<br>&gt; &nbsp;-o 6.tpr<br>&gt; mpirun -machinefile cp -np 8 mdrun -s 6.tpr -o 6.trr -rerun <br>&gt; 3.trr&nbsp; -c<br>&gt; solvent-bmi-pf6-128.pdb -e 6.edr -g 6.log<br>&gt; <br>&gt; I calculate the total energy by<br>&gt; g_energy -f 3.edr -o 3.xvg<br>&gt; g_energy -f 6.edr -o 6.xvg<br>&gt; <br>&gt; The total energy varies between +- 30.00 KJ/mol.<br>&gt; It should be constant since I using same topology file and trajectory.<br>&gt; Why the total energy is not constant.<br><br>Those .tpr should be identical - but you can check that with gmxcheck. Reruns do neighbour-searching every step, whereas your normal simulation followed the nstlist setting. That's part of why my earlier advice suggested doing reruns for each simulation you wish to compare. You should be able to get good/better agreement for steps where nstlist directed neighbour-searching in the original run. Also, whether or not constraints have been applied (and when!) could influence the energies to about this degree. I don't recall the details here.<br><br>Even once you've removed all algorithm-specific sources of difference, there are other sources of non-reproducibility, such as the assignment of particles to DD cells. Your original mdrun probably used dynamic load-balancing, and that cannot be reproduced in the DD used by the rerun. (Or indeed by a repeat of your original mdrun!) Setting -dlb no in the original simulation might be enough to get agreement here, or maybe mdrun -reprod will be required.<br><br>Mark<br><br>&gt; NIlesh<br>&gt; <br>&gt; <br>&gt; <br>&gt; <br>&gt; On Thu, October 21, 2010 10:24 am, Mark Abraham wrote:<br>&gt; &gt; On 22/10/2010 1:17 AM, Nilesh Dhumal wrote:<br>&gt; &gt;<br>&gt; &gt;&gt; I am doing solvation dynamics for my system.<br>&gt; &gt;&gt; I have system with diatomic (PA---NE)solute surrounded by water<br>&gt; &gt;&gt; molecules.<br>&gt; &gt;&gt;<br>&gt; &gt;&gt; I want to run simulation with two differcent cases.<br>&gt; &gt;&gt; 1. PA charge=0 and NE charge=0 : No charge on solute<br>&gt; &gt;&gt; 2. PA charge=+1 and NE charge=-1 : Charge on solute<br>&gt; &gt;&gt;<br>&gt; &gt;&gt;<br>&gt; &gt;&gt; I want to calculate the energy at each step keeping the solvent<br>&gt; &gt;&gt; configration same.<br>&gt; &gt;&gt;<br>&gt; &gt;&gt; IF I start a simulation with no charge on solute (case1), I <br>&gt; have the<br>&gt; &gt;&gt; energy for 1 step. I want to calculate the energy with charge <br>&gt; on solute<br>&gt; &gt;&gt; (case 2) with same configration water molecules.<br>&gt; &gt;&gt; Each step&nbsp; I want to calculate the energy with and <br>&gt; without charge on<br>&gt; &gt;&gt; solute since the configration of solvent will be same for <br>&gt; that step.<br>&gt; &gt;&gt;<br>&gt; &gt;&gt; I was thinking two make two topologies file with charge and <br>&gt; with out<br>&gt; &gt;&gt; charge on solute. I don't know how to use them simultaneously <br>&gt; during the<br>&gt; &gt;&gt;&nbsp; simulation.<br>&gt; &gt;<br>&gt; &gt; Well, you don't use them simultaneously. You run a simulation on<br>&gt; &gt; whatever you think will generate a relevant conformational <br>&gt; ensemble. Then<br>&gt; &gt; you want to use mdrun -rerun twice on the resulting <br>&gt; trajectory, using .tpr<br>&gt; &gt; files based on .top files corresponding to the two cases in <br>&gt; order to<br>&gt; &gt; create your comparison.<br>&gt; &gt;<br>&gt; &gt; Mark<br>&gt; &gt; --<br>&gt; &gt; gmx-users mailing list&nbsp;&nbsp;&nbsp; gmx-users@gromacs.org<br>&gt; &gt; http://lists.gromacs.org/mailman/listinfo/gmx-users<br>&gt; &gt; Please search the archive at<br>&gt; &gt; http://www.gromacs.org/Support/Mailing_Lists/Search before posting!<br>&gt; &gt; Please don't post (un)subscribe requests to the list. Use the<br>&gt; &gt; www interface or send it to gmx-users-request@gromacs.org. <br>&gt; Can't post? Read<br>&gt; &gt; http://www.gromacs.org/Support/Mailing_Lists<br>&gt; &gt;<br>&gt; &gt;<br>&gt; &gt;<br>&gt; <br>&gt; <br>&gt; <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