Re: [gmx-users] Fatal error: Atomtype F not found

2012-10-04 Thread David van der Spoel

On 2012-10-04 04:39, Nur Syafiqah Abdul Ghani wrote:

Dear Users,

Right now i already done for creating the a gro file from antechamber
to gromacs format of my molecule which is hexafluoroisopropanol.
But when i want to minimize it in vacuum it show atomtype F not found.
Im using oplsaa force field and i already change the atom type
according to the force field.


Do not combine force fields. Use antechamber and the gromacs conversion 
script to make a gromacs topology from it (amb2gmx.pl, search on 
google), then use amber force field in combination with this if you want 
to dissolve biomolecules or something like that.


More organic molecules topologies can be found at 
http://virtualchemistry.org


--
David van der Spoel, Ph.D., Professor of Biology
Dept. of Cell  Molec. Biol., Uppsala University.
Box 596, 75124 Uppsala, Sweden. Phone:  +46184714205.
sp...@xray.bmc.uu.sehttp://folding.bmc.uu.se
--
gmx-users mailing listgmx-users@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-requ...@gromacs.org.

* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


Re: [gmx-users] Fatal error: Atomtype F not found

2012-10-03 Thread Justin Lemkul



On 10/3/12 10:43 PM, Nur Syafiqah Abdul Ghani wrote:

Dear Users,

Right now i already done for creating the a gro file from antechamber
to gromacs format of my molecule which is hexafluoroisopropanol.
But when i want to minimize it in vacuum it show atomtype F not found.
Im using oplsaa force field and i already change the atom type
according to the force field.

Below are my files that i need to used for my minimization :

PDB file :

HETATM1  opls_164  HFI A  101 -1.433  -1.291  -0.584  1.00  -0.219790
HETATM2  opls_161  HFI A  101 -1.274  -0.111  -0.039  1.00   0.600750
HETATM3  opls_164  HFI A  101 -1.331  -0.244   1.257  1.00  -0.208510
HETATM4  opls_164  HFI A  101 -2.292   0.639  -0.413  1.00  -0.235630
HETATM5  opls_158  HFI A  101  0.013   0.562  -0.520  1.00   0.064010
HETATM6  opls_078  HFI A  101  0.068   1.864  -0.080  1.00  -0.549590
HETATM7  opls_079  HFI A  101 -0.560   2.404  -0.527  1.00   0.438940
HETATM8  opls_140  HFI A  101  0.014   0.485  -1.602  1.00   0.103910
HETATM9  opls_161  HFI A  101  1.289  -0.122  -0.036  1.00   0.652070
HETATM   10  opls_164  HFI A  101  2.323   0.430  -0.620  1.00  -0.216720
HETATM   11  opls_164  HFI A  101  1.282  -1.395  -0.353  1.00  -0.223970
HETATM   12  opls_164  HFI A  101  1.450  -0.022   1.254  1.00  -0.205480
CONECT12
CONECT21345
CONECT32
CONECT42
CONECT52678
CONECT65
CONECT75
CONECT859   10   11
CONECT98
CONECT   108
CONECT   118   12
CONECT   12   11
END

also the topology ;

; Include forcefield parameters
#include oplsaa.ff/forcefield.itp
; Include organic solvent
#include hfi.itp
; Include water topology
#include oplsaa.ff/spc.itp
#ifdef POSRER_WATER
; Position restraint for each water oxygen
[position_restraints]
; i funct  fcx   fcy   fcz
   11   1000  1000  1000
#endif
[system]
;Name
solute
[molecules]
;Compound  #mols
HFI 1


gro file;

solvent_HFI.gro.gro created by rdparm2gmx.pl Sat May 28 10:25:18 MYT 2005
  12
 1HFI   opls_1641   0.354   0.142   0.000
 1HFI   opls_1612   0.369   0.260   0.055
 1HFI   opls_1643   0.364   0.247   0.184
 1HFI   opls_1644   0.267   0.335   0.017
 1HFI   opls_1585   0.498   0.328   0.006
 1HFI   opls_0786   0.503   0.458   0.050
 1HFI   opls_0797   0.440   0.512   0.006
 1HFI   opls_1408   0.498   0.320  -0.102
 1HFI   opls_1619   0.626   0.260   0.055
 1HFI   opls_164   10   0.729   0.315  -0.004
 1HFI   opls_164   11   0.625   0.132   0.023
 1HFI   opls_164   12   0.642   0.270   0.184
 0.625250.132330.02309

is there something wrong of my file?


Your .pdb and .gro files call opls_* atom names, which are actually atom types. 
 Do not confuse the two.



Why it cant recognize the F atom??


Without seeing the contents of hfi.itp, it's hard to say aside from the general 
statement that you're using an atom type that does not exist.  I suspect you've 
confused atom types and names, given the coordinate files shown.


-Justin

--


Justin A. Lemkul, Ph.D.
Research Scientist
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 listgmx-users@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-requ...@gromacs.org.

* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists