[gmx-users] smooth trajectory

ana johari ana_j0000 at yahoo.com
Sat Apr 2 21:23:50 CEST 2011


dear Justin,
tanks for your kindness.I would like to follow up what you say,but I want to 
have simulation without wraping.
could these command be helpful?
trjconv -s first1ns.tpr -f whatever_trajectory_you_want(xtc after 3ns 
simulation).xtc -pbc mol -center -o fit.gro
 grompp -f md300.mdp -c fit.gro -p n.top -o md300.tpr -n n.ndx
 nohup mpirun -np 8 mdrun -s md300.tpr  -o md300.trr -c md300.gro -g md300.log 
-e md300.edr&
best regards


ana johari wrote:
> dear user,
> 
> Tanks for your advice,
> 
> But the mdrun is know on the step of 33562880  but the protein was split at the 
>first 1ns of MD. is <http://md.is/> it necessary to use “tpbconv” command and 
>extract the *.tpr file just after first 1 ns and  then use trjconv –ur compact 
>as command.
> 

No.  If the starting configuration had the protein properly centered in the box, 
then a simple:

trjconv -s first.tpr -f whatever_trajectory_you_want.xtc -pbc mol -center

should work to make molecules whole and center the protein in the box.

-Justin

> Best regards,
> 
> 
> ------------------------------------------------------------------------
> *From:* Justin A. Lemkul <jalemkul at vt.edu>
> *To:* Discussion list for GROMACS users <gmx-users at gromacs.org>
> *Sent:* Fri, April 1, 2011 2:52:53 PM
> *Subject:* Re: [gmx-users] FW: protein split over boundary
> 
> 
> 
> ana johari wrote:
>  > Dear user,
>  >
>  > Tanks for your attention ,I read about trjconv,but I want to know is it 
>necessary  to back to the exact fram befor molecule split happens and then 
>center the molecule by trjconv command?
>  >
> 
> If the protein starts in the center of the box, just use trjconv in conjunction 
>with your original .tpr file.  A suggested trjconv workflow is on the page 
>Tsjerk pointed you to.
> 
>  > The other point,if you attention to my value”edit conf  d=0.9 and during MD 
>rvdw=1..4” is it broke the rule of periodic boundary condition or not?
>  >
> 
> In principle, no, as long as your box does not significantly deform.  Check 
>with g_mindist -pi.
> 
> -Justin
> 
>  > tanks
>  >
>  >
>  >
>  > ------------------------------------------------------------------------
>  > *From:* Tsjerk Wassenaar <tsjerkw at gmail.com <mailto:tsjerkw at gmail.com>>
>  > *To:* Discussion list for GROMACS users <gmx-users at gromacs.org 
><mailto:gmx-users at gromacs.org>>
>  > *Sent:* Fri, April 1, 2011 1:29:19 PM
>  > *Subject:* Re: [gmx-users] FW: protein split over boundary
>  >
>  > 
>http://www.gromacs.org/Documentation/Terminology/Periodic_Boundary_Conditions
>  >
>  > -TAW
>  >
>  > On Fri, Apr 1, 2011 at 10:53 AM, anahita <ana_j0000 at yahoo.com 
><mailto:ana_j0000 at yahoo.com> <mailto:ana_j0000 at yahoo.com 
><mailto:ana_j0000 at yahoo.com>>> wrote:
>  >  >
>  >  >
>  >  >
>  >  >
>  >  > From: anahita [mailto:ana_j0000 at yahoo.com <mailto:ana_j0000 at yahoo.com> 
><mailto:ana_j0000 at yahoo.com <mailto:ana_j0000 at yahoo.com>>]
>  >  > Sent: Friday, April 01, 2011 1:13 PM
>  >  > To: 'gmx-users-request at gromacs.org <mailto:gmx-users-request at gromacs.org> 
><mailto:gmx-users-request at gromacs.org <mailto:gmx-users-request at gromacs.org>>'
>  >  > Subject: protein split over boundary
>  >  >
>  >  >
>  >  >
>  >  > Dear user,
>  >  >
>  >  > Hello, I appreciate if somebody help me. During the simulation of my 
>protein
>  >  > I didn’t get any error but after 1ns of MD
>  >  >
>  >  > my molecule was split over boundary in cubic box. it means some part of 
>it
>  >  > enter the other side of the box.
>  >  >
>  >  > At first,For decreasing the cost of simulation, in “editconf” command I 
>set
>  >  > the number “d” on 0.9  to decease the box size. I want to mention that
>  >  > during the calculation my rvdw is 1.4.
>  >  >
>  >  > I want to know instead of problem of visualization, the other things is
>  >  > fine?
>  >  >
>  >  > Best regards.
>  >  >
>  >  > A. johari
>  >  >
>  >  > --
>  >  > gmx-users mailing list    gmx-users at gromacs.org 
><mailto:gmx-users at gromacs.org> <mailto:gmx-users at gromacs.org 
><mailto:gmx-users at gromacs.org>>
>  >  > http://lists.gromacs.org/mailman/listinfo/gmx-users
>  >  > Please search the archive at
>  >  > http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
>  >  > Please don't post (un)subscribe requests to the list. Use the
>  >  > www interface or send it to gmx-users-request at gromacs.org 
><mailto:gmx-users-request at gromacs.org> <mailto:gmx-users-request at gromacs.org 
><mailto:gmx-users-request at gromacs.org>>.
>  >  > Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>  >  >
>  >
>  >
>  >
>  > -- Tsjerk A. Wassenaar, Ph.D.
>  >
>  > post-doctoral researcher
>  > Molecular Dynamics Group
>  > * Groningen Institute for Biomolecular Research and Biotechnology
>  > * Zernike Institute for Advanced Materials
>  > University of Groningen
>  > The Netherlands
>  > --
>  > gmx-users mailing list    gmx-users at gromacs.org 
><mailto:gmx-users at gromacs.org> <mailto:gmx-users at gromacs.org 
><mailto:gmx-users at gromacs.org>>
>  > http://lists.gromacs.org/mailman/listinfo/gmx-users
>  > Please search the archive at 
>http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
>  > Please don't post (un)subscribe requests to the list. Use the
>  > www interface or send it to gmx-users-request at gromacs.org 
><mailto:gmx-users-request at gromacs.org> <mailto:gmx-users-request at gromacs.org 
><mailto:gmx-users-request at gromacs.org>>.
>  > Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>  >
> 
> -- ========================================
> 
> 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 list    gmx-users at gromacs.org 
><mailto:gmx-users at gromacs.org>
> http://lists.gromacs.org/mailman/listinfo/gmx-users
> Please search the archive at 
>http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
> Please don't post (un)subscribe requests to the list. Use the www interface or 
>send it to gmx-users-request at gromacs.org <mailto:gmx-users-request at gromacs.org>.
> Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

-- ========================================

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 list    gmx-users at gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/Search 
before posting!
Please don't post (un)subscribe requests to the list. Use the www interface or 
send it to gmx-users-request at gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20110402/97b5e62f/attachment.html>


More information about the gromacs.org_gmx-users mailing list