[gmx-users] Question about topologies

Thanasis Koukoulas koukoulas_th at yahoo.gr
Fri Jun 18 11:54:27 CEST 2010


Hello
[ atomtypes ];   type      mass    charge    ptype       c6            c12    CH3     15.035     0.000    A       9.0638e-3      2.5206e-5    CH2     14.027     0.000    A       5.8108e-3      2.2071e-5    CS2     14.027     0.000    A       5.8108e-3      2.2071e-5    OS      15.999     0.000    A       8.8147e-4      4.2477e-7     OA      15.999     0.000    A       2.3465e-3      1.7802e-6    HO       1.008     0.000    A       0.0000e00      0.0000e00 
at atomtypes i can give c6 and c12 values and gromacs will generate the c6ij and c12ij values from the formulas i have chosen with selecting combination rule Am i right????
in my case i gave c6 and c12 as c6i=4*ei*σi^6 and c12i=4*ei*σi^12 and i chose combination rule 1. so gromacs will generate the c6ij and c12ij as c6ij=sqrt(C6iC6j) and c12ij=sqrt(c6i*c12j)
So if in my topology I put [nonbond_params] what happens then? the values from [nonbond_params] c6 and c12 are being replaced by the values gromacs finds from the equations???
for example
if i have 
[nonbond_params] 
CH3 CH3 5e-3 6e-5
and gromacs calculates from the [atomtypes] are different lets say (i didnt do the calculation it is an assumption  value) 3e-2 and 3e-5 which values gromacs accepts?i have read the manual but it isnd clear how this works. 
please someone enlighten me.
thanks

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20100618/07c93caf/attachment.html>


More information about the gromacs.org_gmx-users mailing list