VMD-L Mailing List
From: Mayne, Christopher G (cmayne2_at_illinois.edu)
Date: Thu Sep 17 2015 - 16:21:12 CDT
- Next message: Gianluca Interlandi: "Plugin to read DCD files"
- Previous message: Amy Rice: "ffTK: Charge optimization"
- In reply to: Amy Rice: "ffTK: Charge optimization"
- Next in thread: Amy Rice: "Re: ffTK: Charge optimization"
- Reply: Amy Rice: "Re: ffTK: Charge optimization"
- Messages sorted by: [ date ] [ thread ] [ subject ] [ author ] [ attachment ]
Amy,
I'm not quite sure what's going on here. It seems like it's looking for an atom type that is a concatenation of several atoms, and I'm not sure why it would do that since you're PSF seems OK. I can come back to this problem, but I see a couple of other issues here.
Based on the contents of your initial PAR file I can tell that you didn't include the par_all36_prot.prm file as an "associated parameter file". Many of the parameters that you have in your PAR file are defined by the standard CHARMM force field for proteins. Including this file when generating the initial parameter file will significantly reduce the number of missing parameters that you will have to work on.
You should check out the CGenFF topology file, specifically how residue GUA changes to MGUA. Comparing these two entries will be helpful in assessing the changes associated with methylating the nitrogens of guanidinium, e.g., how the partial atomic charges change, and it will give you proper LJ parameters for the terminal methyl. Many of the missing parameters should be defined, albeit in CGenFF terms (i.e., using CGenFF atom types rather than CHARMM protein atom types). It may be a little tedious, but you should be able to translate these.
Anything that remains missing won't required parameterizing the full molecule that you've constructed. At worst you can truncate the system to an N-methyl and N-propyl guanidinium.
Regards,
Christopher Mayne
On Sep 17, 2015, at 2:39 PM, Amy Rice wrote:
Hi all,
I am attempting to parameterize a modified arginine residue following the ffTK screencasts, using VMD 1.9.2. I've been able to complete the first 4 steps without issue, but I'm running into some errors during the charge optimization, in both downhill mode and simulated annealing. The error message reads:
"can't read "ljPars(HA2 HA2 HA2 H)": no such element in array"
I noticed a similar problem was discussed previously (http://www.ks.uiuc.edu/Research/vmd/mailing_list/vmd-l/22158.html Thank you for the help,
- Amy
--
Amy Rice
Ph.D. Student
Physics Department
Illinois Institute of Technology
<m_arginine_files.zip>