[gmx-users] x2top - wrong mass in .top file

Jon jon.ellis at utoronto.ca
Tue Mar 22 19:37:48 CET 2005


Hello Users

I am trying to make a topology file of a quartz surface, from a 
coordinate file (.gro) that I generated myself. I started with a single 
unit cell of quartz (SiO2)_3, and am using the GMX force field. After 
modifying the .n2t file so that x2top would run, it generated a .top 
file. The mass of Silicon is fine (28.080), but the mass of the oxygen 
is set at 12, not the 15.99 that is set in ffgmx.atp.

I am using standard atomtypes (OS for oxygen, SI for silicon). I had to 
play around with the n2t file, such that some of the bonding is not 
realistic, and there are no hydrogens, but this is just a first attempt, 
so its okay for now.

Here are excerpts from my gro (input) and top (output) files:

.gro
    1QRTZ    SI    1   0.000   0.246   0.078
    1QRTZ   OS1    2   0.098   0.237   0.208
    1QRTZ   OS2    3   0.026   0.387   0.000
    2QRTZ    SI    4   0.139   0.143   0.334
    2QRTZ   OS1    5   0.008   0.116   0.428
    2QRTZ   OS2    6   0.196   0.000   0.279
    3QRTZ    SI    7   0.348   0.200   0.555
    3QRTZ   OS1    8   0.254   0.217   0.422
    3QRTZ   OS2    9   0.443   0.331   0.571

.top
    [ atoms ]
     1         SI      1   QRTZ     SI      0          0      28.08
     2         OS      1   QRTZ    OS1      1          0         12   
     3         OS      1   QRTZ    OS2      2          0         12    
     4         SI      2   QRTZ     SI      3          0      28.08
     5          O      2   QRTZ    OS1      4          0         12     
     6          O      2   QRTZ    OS2      5          0         12    
     7         SI      3   QRTZ     SI      6          0      28.08    
     8         OS      3   QRTZ    OS1      7          0         12    
     9         OS      3   QRTZ    OS2      8          0         12

Why am I getting the wrong mass (12 instead of 15.99)? As well, I told 
it to make an rtp file, but I cannot find it anywhere. Is this normal? 
My command was "x2top -f X.gro -o Y.top -r Z.rtp -param"

Cheers
Jon


-- 
----------
Jon Ellis
Thompson Group
Institute of Biomaterials and Biomedical Engineering/Dept of Chemistry
University of Toronto
LM147 - (416) 978-6568
jon.ellis at utoronto.ca




More information about the gromacs.org_gmx-users mailing list