[gmx-users] topology file for methane.

Justin A. Lemkul jalemkul at vt.edu
Mon Oct 27 23:30:14 CET 2008



xianghong qi wrote:
> Hi, Justin:
> I am reading chapter 5 now. Still feel confused. I wrote methane.top 
> file like this:
> 
> ; Methane topology file
> ;
> ; The force field files to be included
> #include "ffoplsaa.itp"
> 
> [ moleculetype ]
> ; Name            nrexcl
>  METH             3
> 
> [ atoms ]
> ;   nr       type  resnr residue  atom   cgnr     charge       mass  
> typeB    chargeB      massB
>      1   opls_138      1    METH     C      1      -0.24     12.011   ; 
> qtot 0.61
>      2   opls_140      1    METH     H      1       0.06      1.008   ; 
> qtot 0.67
>      3   opls_140      1    METH     H      1       0.06      1.008   ; 
> qtot 0.67
>      4   opls_140      1    METH     H      1       0.06      1.008   ; 
> qtot 0.73
>      5   opls_140      1    METH     H      1       0.06      1.008   ; 
> qtot 0.79
> 
> [ bonds ]
> ;  ai    aj funct            c0            c1            c2            c3
>     1     2     1
>     1     3     1
>     1     4     1
>     1     5     1
> 
> ; Include Position restraint file
> #ifdef POSRES
> #include "posre.itp"
> #endif
> 
> ; Include water topology
> #include "spc.itp"
> 
> [ system ]
> ; Name
> Methane
> 
> [ molecules ]
> ; Compound        #mols
> METH             1
> ; ;;;;;;;
> 
> Then I should go to ffoplsaa.rtp file to add a new molecule like :
> 
> [METH]
>  [ atoms ]
>      C    opls_135   -0.240     1
>    HC1    opls_140    0.060     1
>    HC2    opls_140    0.060     1
>    HC3    opls_140    0.060     1
>    HC4    opls_140    0.60      1
>  [ bonds ]
>      C  HC1
>      C  HC2
>      C  HC3
>      C  HC4
> ;;;;;;;
> 
> Is that reasonable?  Thanks for your help.
> 
> 

Seems fine to me.

-Justin

> On Mon, Oct 27, 2008 at 1:54 PM, Justin A. Lemkul <jalemkul at vt.edu 
> <mailto:jalemkul at vt.edu>> wrote:
> 
> 
> 
>     xianghong qi wrote:
> 
>         Hello,
> 
>         Could anyone give me some suggestions about how to write a
>         toplogy file for methane with OPLS-AA force field?
>         Thanks in advance.
> 
> 
>     For a very simple molecule like methane, it should be very
>     straightforward to do from reading Chapter 5.
> 
>     -Justin
> 
>         -Xianghong Qi
> 
>         -- 
>         Some people make the world more special just by being in it.
> 
> 
>         ------------------------------------------------------------------------
> 
>         _______________________________________________
>         gmx-users mailing list    gmx-users at gromacs.org
>         <mailto:gmx-users at gromacs.org>
>         http://www.gromacs.org/mailman/listinfo/gmx-users
>         Please search the archive at http://www.gromacs.org/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/mailing_lists/users.php
> 
> 
>     -- 
>     ========================================
> 
>     Justin A. Lemkul
>     Graduate Research Assistant
>     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://www.gromacs.org/mailman/listinfo/gmx-users
>     Please search the archive at http://www.gromacs.org/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/mailing_lists/users.php
> 
> 
> 
> 
> -- 
> Some people make the world more special just by being in it.

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

Justin A. Lemkul
Graduate Research Assistant
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin

========================================



More information about the gromacs.org_gmx-users mailing list