<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style></head>
<body class='hmmessage'><div dir='ltr'>
Hi,<br>the center of mass of my channel is at the middle of the ion channel, and it is a symmetric system, so I suppose these results would be OK. Anyway, I will check the options you propose.<br>Thanks a lot for all your comments!!<br>Best wishes,<br>Rebeca.<br><br><div><div id="SkyDrivePlaceholder"></div>&gt; Date: Thu, 31 May 2012 20:08:26 +0200<br>&gt; From: schlesi@uni-mainz.de<br>&gt; To: gmx-users@gromacs.org<br>&gt; Subject: [gmx-users] Re: boundaries in PMF<br>&gt; <br>&gt; Where is the center of mass of reference group (MOL) located?<br>&gt; <br>&gt; It seems that the COM is near the middle of the ion channel. Since you <br>&gt; use 'pull_geometry=distance', g_wham will look only for the distance <br>&gt; between 'MOL' and 'Na' and that leads to problem.<br>&gt; If the com of 'MOL' sits in the center of the channel (which is around <br>&gt; 5nm long), one has 2.5nm in both directions. Since g_wham uses only the <br>&gt; distance you get the PMF for half of the channel, but with the data of <br>&gt; both parts.<br>&gt; If the channel would be symmetric and the com of 'MOL' would lie exactly <br>&gt; in the middle of the channel, this could be ok. But if even one of both <br>&gt; assumptions is wrong, the results would be errorous.<br>&gt; <br>&gt; A better approach would be to use 'pull_geometry=direction', since the <br>&gt; you define a vector along which the windows lie and do not have the <br>&gt; problem that the distance is in both directions (along positive and <br>&gt; negative vector) the same.<br>&gt; Only problem could be that g_wham supports 'pull_geometry=direction' <br>&gt; only from gromacs 4.5.x (don't know this, since instead of umbrella <br>&gt; smapling i use thermodynamik integration, where one uses constraints <br>&gt; (instead of restraints) and integrates the constraint-force; but the <br>&gt; conceptual problem with 'distance/direction' is the same).<br>&gt; <br>&gt; Another approach (with 'pull_geometry=distance') would be to use a <br>&gt; reference group which is just outside of the channel (like going some <br>&gt; steps away from the channel, along the vector which goes through the <br>&gt; channel). If there is only water, it would be bad, because then the <br>&gt; reference group would be move away.<br>&gt; But then on could use the entry and exit of the channel as a reference <br>&gt; point for two simulations. In the case that the entry is the reference <br>&gt; group, the PMF would be ill defined near the entry, but from the <br>&gt; simulation with the exit as reference you would get the right PMF for <br>&gt; the entry region and vice versa.<br>&gt; <br>&gt; Greetings<br>&gt; Thomas<br>&gt; <br>&gt; <br>&gt; Am 31.05.2012 19:39, schrieb gmx-users-request@gromacs.org:<br>&gt; &gt; Thanks a lot for your quick answer. The mdp file I have used is copied<br>&gt; &gt; below. What is strange is that when I look at the *gro files for the<br>&gt; &gt; different windows (100 windows in total), i. e: window 1: 8704Na Na56458<br>&gt; &gt; 5.134 5.085 5.824 window 50: 8704Na Na56458 5.053 5.081 3.459 window<br>&gt; &gt; 100: 8704Na Na56458 4.990 5.042 0.951 you can see that the z-coordinate<br>&gt; &gt; goes from 0.951 to 5.824 nm I should have a total distance in the x-axis<br>&gt; &gt; of about 5 nm, and however, the boundaries calculated by g_wham are<br>&gt; &gt; "Determined boundaries to 0.000035 to 2.603290 " Can you see anything in<br>&gt; &gt; the mdp file which is causing this behaviour...? Thanks again for your<br>&gt; &gt; help! MDP FILE USED: title = Umbrella pulling simulation define =<br>&gt; &gt; -DPOSRES_B define = ; Run parameters integrator = md dt = 0.002 tinit =<br>&gt; &gt; 0 nsteps = 500000 ; 1 ns nstcomm = 10 ; Output parameters nstxout = 5000<br>&gt; &gt; ; every 10 ps nstvout = 5000 nstfout = 5000 nstxtcout = 5000 ; every 10<br>&gt; &gt; ps nstenergy = 5000 ; Bond parameters constraint_algorithm = lincs<br>&gt; &gt; constraints = all-bonds continuation = yes ; Single-range cutoff scheme<br>&gt; &gt; nstlist = 5 ns_type = grid rlist = 1.4 rcoulomb = 1.4 rvdw = 1.4 ; PME<br>&gt; &gt; electrostatics parameters coulombtype = PME fourierspacing = 0.12<br>&gt; &gt; fourier_nx = 0 fourier_ny = 0 fourier_nz = 0 pme_order = 4 ewald_rtol =<br>&gt; &gt; 1e-5 optimize_fft = yes ;Berendsen temperature coupling is on Tcoupl =<br>&gt; &gt; V-rescale tau_t = 0.1 0.1 0.1 tc-grps = MOL dop WAT_Cl_Na ref_t = 300<br>&gt; &gt; 300 300 ; Pressure coupling Pcoupl = Parrinello-Rahman Pcoupltype =<br>&gt; &gt; Semiisotropic ; Time constant (ps), compressibility (1/bar) and<br>&gt; &gt; reference P (bar) tau-p = 1.0 1.0 compressibility = 4.6E-5 4.6E-5 ref-p<br>&gt; &gt; = 1.0 1.0 ; Generate velocities is off gen_vel = no ; Periodic boundary<br>&gt; &gt; conditions are on in all directions pbc = xyz ; Long-range dispersion<br>&gt; &gt; correction DispCorr = EnerPres ; Pull code pull = umbrella pull_geometry<br>&gt; &gt; = distance pull_dim = N N Y pull_start = yes pull_ngroups = 1<br>&gt; &gt; pull_group0 = MOL pull_group1 = Na pull_init1 = 0 pull_rate1 = 0.0<br>&gt; &gt; pull_k1 = 1000 ; kJ mol^-1 nm^-2 pull_nstxout = 1000 ; every 2 ps<br>&gt; &gt; pull_nstfout = 1000 ; every 2 ps<br>&gt; &gt;&gt; &gt;  Date: Thu, 31 May 2012 13:25:26 -0400<br>&gt; &gt;&gt; &gt;  From:jalemkul@vt.edu<br>&gt; &gt;&gt; &gt;  To:gmx-users@gromacs.org<br>&gt; &gt;&gt; &gt;  Subject: Re: [gmx-users] boundaries in PMF<br>&gt; &gt;&gt; &gt;<br>&gt; &gt;&gt; &gt;<br>&gt; &gt;&gt; &gt;<br>&gt; &gt;&gt; &gt;  On 5/31/12 1:20 PM, Rebeca Garc?a Fandi?o wrote:<br>&gt; &gt;&gt;&gt; &gt;  &gt;  Hi,<br>&gt; &gt;&gt;&gt; &gt;  &gt;  I am trying to calculate a PMF for an ion along a channel. Everything went OK,<br>&gt; &gt;&gt;&gt; &gt;  &gt;  but when I used g_wham I got a profile with strange dimensions for the x-axis.<br>&gt; &gt;&gt;&gt; &gt;  &gt;  What are the boundaries g_wham is using for calculating the units of x-axis?<br>&gt; &gt;&gt;&gt; &gt;  &gt;<br>&gt; &gt;&gt; &gt;<br>&gt; &gt;&gt; &gt;  The values are based on the restraint distances along the reaction coordinate.<br>&gt; &gt;&gt; &gt;<br>&gt; &gt;&gt;&gt; &gt;  &gt;  I have used:<br>&gt; &gt;&gt;&gt; &gt;  &gt;  g_wham -it tpr-files.dat -if pullf-files.dat -o file_output.xvg -hist<br>&gt; &gt;&gt;&gt; &gt;  &gt;  file_histo_output.xvg -unit kCal -cycl yes<br>&gt; &gt;&gt;&gt; &gt;  &gt;<br>&gt; &gt;&gt;&gt; &gt;  &gt;  (version 4.0.7)<br>&gt; &gt;&gt;&gt; &gt;  &gt;<br>&gt; &gt;&gt;&gt; &gt;  &gt;  and the units I got in the x-axis are determined by the boundaries:<br>&gt; &gt;&gt;&gt; &gt;  &gt;<br>&gt; &gt;&gt;&gt; &gt;  &gt;  "Determined boundaries to 0.000035 to 2.603290"<br>&gt; &gt;&gt;&gt; &gt;  &gt;<br>&gt; &gt;&gt;&gt; &gt;  &gt;  Which are these units? nm?<br>&gt; &gt;&gt;&gt; &gt;  &gt;<br>&gt; &gt;&gt; &gt;<br>&gt; &gt;&gt; &gt;  Yes.<br>&gt; &gt;&gt; &gt;<br>&gt; &gt;&gt;&gt; &gt;  &gt;  The z coordinate for my ion explores at least 5 nm!!<br>&gt; &gt;&gt;&gt; &gt;  &gt;<br>&gt; &gt;&gt;&gt; &gt;  &gt;  I am a bit confused about that.<br>&gt; &gt;&gt;&gt; &gt;  &gt;<br>&gt; &gt;&gt; &gt;<br>&gt; &gt;&gt; &gt;  The exact output depends on how you set up the umbrella sampling (in the .mdp<br>&gt; &gt;&gt; &gt;  file).  The range of values corresponds to whatever the distances are that are<br>&gt; &gt;&gt; &gt;  sampled in the various windows.  Perhaps there is a sign issue here?  Do you<br>&gt; &gt;&gt; &gt;  have some restraints that are at negative displacement and others at positive?<br>&gt; &gt;&gt; &gt;  Did you set up the .mdp files properly to account for this behavior?<br>&gt; &gt;&gt; &gt;<br>&gt; &gt;&gt; &gt;  -Justin<br>&gt; <br>&gt; -- <br>&gt; gmx-users mailing list    gmx-users@gromacs.org<br>&gt; http://lists.gromacs.org/mailman/listinfo/gmx-users<br>&gt; Please search the archive at 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<br></div>                                               </div></body>
</html>