Dear Justin,<br />
<br />
Thanks a lot for your reply.<br />
Here I am providing the data and explaining the problem in detail:<br />
<br />
I hope that it will help you in understanding in a better way. 
Considering the tetramer as 'abcd'.  'ab' as one dimer and 'cd' as 
another, I want to separate 'ab' and 'cd'. Length of each dimer is 8 nm,
 so the tetramer (initial configuration) is approximately 16 nm and I 
have used a cubic box of 30 nm. Two groups (reference group and pull 
group) are decided on the basis of residue numbers forming 'ab' as 
reference grp and 'bc' as pull group by making an index file. <br />
<br />
I have verified that the initial configuration does not split across 
periodic boundaries, it is always inside the box throughout the 
simulation. Moreover, the initial distance is not equal to box vector. <br />
<br />
By <span style="font-weight: bold;">g_dist</span> , I have tried to calculate the distance between the two groups and it came out to be approx. 15 nm (Please have a look at the <span style="font-weight: bold;">dist.xvg</span>
 file at the end). Do you think that in the first ps time frame, should 
it get separated by 15 nm?? Please correct me if I have done anything 
wrong.<br />
<br />
By mistake, in my last mail I wrote that I need to pull along Y axis. In
 actual, it is Z-axis and I have used pull_dim N N Y. So it is just the 
mistake in writing the mail.<br />
<br />
Please let me know do you need any other data or parameters that I have used during this simulation.<br />
An early reply would me much appreciated.<br />
<br />
<br />
-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------<br />
-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------<br />
dist.xvg<br />
------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------<br />
<br />
# This file was created Tue Sep 13 11:09:23 2011<br />
# by the following command:<br />
# g_dist -f tub_pull.xtc -s tub_pull.tpr -n index.ndx <br />
#<br />
# g_dist is part of G R O M A C S:<br />
#<br />
# Giving Russians Opium May Alter Current Situation<br />
#<br />
@    title &quot;Distance&quot;<br />
@    xaxis  label &quot;Time (ps)&quot;<br />
@    yaxis  label &quot;Distance (nm)&quot;<br />
@TYPE xy<br />
@ view 0.15, 0.15, 0.75, 0.85<br />
@ legend on<br />
@ legend box on<br />
@ legend loctype view<br />
@ legend 0.78, 0.8<br />
@ legend length 2<br />
@ s0 legend &quot;|d|&quot;<br />
@ s1 legend &quot;d\sx\N&quot;<br />
@ s2 legend &quot;d\sy\N&quot;<br />
@ s3 legend &quot;d\sz\N&quot;<br />
   0.0000000   15.0342455   -0.7714367   -0.9718266   14.9829559<br />
   0.5000000   15.0345469   -0.7726345   -0.9766102  -14.9828863<br />
   1.0000000   15.0345631   -0.7665710   -0.9862967  -14.9825792<br />
   1.5000000   15.0343399   -0.7575226   -0.9916229   14.9824638<br />
   2.0000000   15.0335970   -0.7490749   -0.9957008  -14.9818726<br />
   2.5000000   15.0322971   -0.7446661   -0.9996347  -14.9805260<br />
   3.0000000   15.0341187   -0.7364883   -1.0029345  -14.9825373<br />
   3.5000000   15.0331717   -0.7328205   -1.0078001   14.9814405<br />
   4.0000000   15.0311899   -0.7334194   -1.0118599   14.9791489<br />
   4.5000000   15.0334759   -0.7351160   -1.0128736   14.9812908<br />
   5.0000000   15.0336533   -0.7327585   -1.0102243   14.9817638<br />
   5.5000000   15.0327044   -0.7342167   -1.0076447   14.9809132<br />
   6.0000000   15.0335188   -0.7365904   -1.0090990  -14.9815168<br />
   6.5000000   15.0296631   -0.7375288   -1.0111809  -14.9774609<br />
   7.0000000   15.0294476   -0.7318525   -1.0155020  -14.9772310<br />
   7.5000000   15.0322533   -0.7323456   -1.0173597  -14.9798965<br />
   8.0000000   15.0353861   -0.7333517   -1.0172930   14.9829950<br />
   8.5000000   15.0323915   -0.7285290   -1.0122719   14.9805660<br />
   9.0000000   15.0284719   -0.7241497   -1.0072327   14.9771843<br />
   9.5000000   15.0285101   -0.7232113   -1.0051012   14.9774113<br />
  10.0000000   15.0306396   -0.7257242   -1.0020084   14.9796333<br />
  10.5000000   15.0324621   -0.7242069   -0.9959688   14.9819384<br />
  11.0000000   15.0324860   -0.7183657   -0.9881392  -14.9827623<br />
  11.5000000   15.0308609   -0.7121477   -0.9848318  -14.9816465<br />
  12.0000000   15.0318832   -0.7091122   -0.9842033   14.9828577<br />
  12.5000000   15.0316772   -0.7063818   -0.9834919  -14.9828262<br />
  13.0000000   15.0291939   -0.7024584   -0.9846754  -14.9804420<br />
  13.5000000   15.0307999   -0.7052488   -0.9864941   14.9818020<br />
  14.0000000   15.0286407   -0.7049980   -0.9902172   14.9794016<br />
  14.5000000   15.0301456   -0.7010813   -0.9916439  -14.9810019<br />
  15.0000000   15.0305882   -0.6967783   -0.9922857  -14.9816036<br />
  15.5000000   15.0319242   -0.6941795   -0.9950562   14.9828815<br />
  16.0000000   15.0315027   -0.6892128   -0.9928131   14.9828358<br />
  16.5000000   15.0314274   -0.6889334   -0.9883585  -14.9830685<br />
  17.0000000   15.0296593   -0.6906672   -0.9852924   14.9814167<br />
  17.5000000   15.0277662   -0.6964712   -0.9844646  -14.9793034<br />
  18.0000000   15.0301647   -0.7006731   -0.9863977  -14.9813862<br />
  18.5000000   15.0285692   -0.7027130   -0.9882107   14.9795704<br />
  19.0000000   15.0304384   -0.7040758   -0.9874868  -14.9814291<br />
  19.5000000   15.0310335   -0.7008915   -0.9824295  -14.9825077<br />
  20.0000000   15.0304193   -0.6977320   -0.9794283   14.9822359<br />
  20.5000000   15.0306473   -0.6955061   -0.9762669   14.9827747<br />
  21.0000000   15.0293531   -0.6917162   -0.9746027   14.9817610<br />
  21.5000000   15.0295181   -0.6870613   -0.9746943  -14.9821339<br />
  22.0000000   15.0302134   -0.6855507   -0.9771519   14.9827404<br />
  22.5000000   15.0306492   -0.6869965   -0.9785366   14.9830208<br />
  23.0000000   15.0285063   -0.6840801   -0.9742985  -14.9812803<br />
  23.5000000   15.0265427   -0.6824503   -0.9738894   14.9794121<br />
  24.0000000   15.0262260   -0.6803436   -0.9716139   14.9793377<br />
  24.5000000   15.0266447   -0.6778135   -0.9723339  -14.9798260<br />
  25.0000000   15.0236931   -0.6710443   -0.9711208  -14.9772491<br />
<br />
<br />
<br />
<br />
Best regards,<br />
Shilpi Chaurasia<br />
<br />
<br />
<br />
<br />
<br />
<br /><br /><span>On 09/09/11, <b class="name">&quot;Justin A. Lemkul&quot; </b> &lt;jalemkul@vt.edu&gt; wrote:</span><blockquote cite="mid:4E6A4D6E.2090505@vt.edu" class="iwcQuote" style="border-left: 1px solid rgb(0, 0, 255); padding-left: 13px; margin-left: 0pt;" type="cite"><div class="mimepart text plain"><br /><br />Shilpi Chaurasia wrote:<br />&gt; I am trying to pull apart two dimers of tubulin protein joined together, to form two separate dimers by steered molecular dynamics using Pull code. The objective is to separate the dimers by pulling along one axis only (Y axis). I ran the simulation for 25 ps and it accomplished well but while visualizing the trajectory in VMD, in the very beginning (first frame), dimers are at distance which is much greater than expected. I expected that in first frame they should not have separated too much (i.e. in just 1 ps, they are at greater distance) I guess this is unusual. Please let me know where I am doing wrong? I have used the following pull code:<br />&gt;<br /><br />Actual numbers would help.  Is your initial configuration simply split across periodic boundaries?  If the initial distance is roughly equivalent to one of the box vectors, then that's the case.  Otherwise, please provide real data and check to make sure your pull groups are what you think they are.<br /><br />&gt;pull            = umbrella<br />&gt;pull_geometry   = distance<br />&gt;pull_dim        = N N Y<br /><br />Note that you're not actually pulling along Y as you state above.  Here, you're pulling along Z.  Perhaps that's one issue.<br /><br />-Justin<br /><br />&gt;pull_start      = yes         pull_init1      = 0.0<br />&gt;pull_ngroups    = 1<br />&gt;pull_group0     = dimer_ref               ; reference group (based on index file)<br />&gt;pull_group1     = dimer_pull           ; pull group<br />&gt;pull_rate1      = 0.1          pull_k1         = 1000         pull_nstfout    = 10          pull_nstxout    = 10<br />&gt;<br />&gt;best regards,<br />&gt;Shilpi Chaurasia<br />&gt;<br /><br />-- <br />========================================<br /><br />Justin A. Lemkul<br />Ph.D. Candidate<br />ICTAS Doctoral Scholar<br />MILES-IGERT Trainee<br />Department of Biochemistry<br />Virginia Tech<br />Blacksburg, VA<br />jalemkul[at]vt.edu | (540) 231-9080<br /><a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="l">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br /><br />========================================<br />-- <br />gmx-users mailing list    gmx-users@gromacs.org<br /><a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="l">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br />Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="l">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br />Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-request@gromacs.org.<br />Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="l">http://www.gromacs.org/Support/Mailing_Lists</a><br /></div></blockquote>