paratool for MG

From: Tianjiao Wang (tjwang_at_iastate.edu)
Date: Fri Nov 30 2007 - 12:48:05 CST

Hello,

  I am trying to use paratool from VMD 1.8.6 to generate topology and
parameter files for malachite green (MG) simulation in NAMD. I first loaded
MG.pdb as base molecules. Then I eidted atom properties. I loaded
top_all27_prot_na.rtf as topology file and Par_all27_prot_na.prm as parameter
files. Optimized geometry and single point calc (Hessian + charges) were done
In GAMESS. When I tried to import geometry log file from GAMESS into paratool
main, I got error message as followed:

list index out of range
list index out of range
    while executing
"lset zmat 0 4 $ndihed "
    (procedure "find_improper_dihedrals" line 22)
    invoked from within
"find_improper_dihedrals $zmat"
    (procedure "read_gaussian_log" line 191)
    invoked from within
"read_gaussian_log $file $molid"
    (procedure "load_gaussian_log" line 17)
    invoked from within
"load_gaussian_log $file"
    ("log" arm line 1)
    invoked from within
"switch $type {
      int { read_intcoords $file; }
      zmt { set molid [read_zmtfile $file] }
      com { set molid [read_gaussian_input $file] }
  ..."
    (procedure "::QMtool::load_file" line 6)
    invoked from within
"::QMtool::load_file $file [string range [file extension $file] 1 end"
    (procedure "load_molecule" line 14)
    invoked from within
"load_molecule OPT $file"
    (procedure "::Paratool::opendialog" line 128)
    invoked from within
"::Paratool::opendialog loadoptgeom "[file rootname ${::Paratool::molnamebase}]
_opt""
    (menu invoke)

Could anybody give me some advice how to solve this problem? Alternatively,
are there other ways to generate par_ and top_ file for malachite green, an
organic dye (MW:329 dalton).

Thanks a lot!

Tianjiao
 

Tianjiao Wang
Rm 3288 MBB
Iowa State Univ
Ames, IA 50010
e-mail: tjwang_at_iastate.edu
tel: 515-441-6417(cell)
     515-294-5682(O)

This archive was generated by hypermail 2.1.6 : Wed Feb 29 2012 - 15:45:38 CST