<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta http-equiv="Content-Type" content="text/html;charset=ISO-8859-1">
<title></title>
</head>
<body text="#000000" bgcolor="#ffffff">
Ryan-<br>
A quick look at your errors suggests that the energy evaluator is
looking for the Lennard-Jones parameters. I could be wrong, but check
out the argon entry in Database/Atoms/ar for comparison. <br>
Good luck.<br>
Eric<br>
<br>
Ryan Gutenkunst wrote:<br>
<blockquote type="cite" cite="mid40FD6F2C.5090607@cornell.edu">Hi all,
<br>
<br>
I'm trying to simulate a protein along with a magnesium ion present in
the pdb file, and I'm getting errors in the energy evaluation.
<br>
<br>
Below I've attached the traceback, and my Database files Atoms/mg and
Molecules/mg
<br>
<br>
Does anyone know what's wrong?
<br>
<br>
Thanks for any help,
<br>
<br>
The traceback is here:
<br>
<br>
Traceback (most recent call last):
<br>
File "run.py", line 78, in ?
<br>
minimizer(convergence = 1.e-3, steps = minSteps)
<br>
File
"/afs/msc.cornell.edu/home/jsethna/ryang/installed/lib/python2.3/site-packages/MMTK/Minimization.py",
line 124, in __call__
<br>
mpi_communicator=comm)
<br>
File
"/afs/msc.cornell.edu/home/jsethna/ryang/installed/lib/python2.3/site-packages/MMTK/Universe.py",
line 615, in energyEvaluator
<br>
threads, mpi_communicator)
<br>
File
"/afs/msc.cornell.edu/home/jsethna/ryang/installed/lib/python2.3/site-packages/MMTK/ForceFields/ForceField.py",
line 183, in __init__
<br>
self.global_data)
<br>
File
"/afs/msc.cornell.edu/home/jsethna/ryang/installed/lib/python2.3/site-packages/MMTK/ForceFields/MMForceField.py",
line 320, in evaluatorTerms
<br>
global_data)
<br>
File
"/afs/msc.cornell.edu/home/jsethna/ryang/installed/lib/python2.3/site-packages/MMTK/ForceFields/ForceField.py",
line 93, in evaluatorTerms
<br>
global_data)
<br>
File
"/afs/msc.cornell.edu/home/jsethna/ryang/installed/lib/python2.3/site-packages/MMTK/ForceFields/MMForceField.py",
line 301, in evaluatorTerms
<br>
subset2, global_data)
<br>
File
"/afs/msc.cornell.edu/home/jsethna/ryang/installed/lib/python2.3/site-packages/MMTK/ForceFields/ForceField.py",
line 93, in evaluatorTerms
<br>
global_data)
<br>
File
"/afs/msc.cornell.edu/home/jsethna/ryang/installed/lib/python2.3/site-packages/MMTK/ForceFields/NonBondedInteractions.py",
line 132, in evaluatorTerms
<br>
eps, sigma, mix = self._ljParameters(t, global_data)
<br>
File
"/afs/msc.cornell.edu/home/jsethna/ryang/installed/lib/python2.3/site-packages/MMTK/ForceFields/MMForceField.py",
line 68, in _ljParameters
<br>
return self.dataset.ljParameters(type)
<br>
File
"/afs/msc.cornell.edu/home/jsethna/ryang/installed/lib/python2.3/site-packages/MMTK/ForceFields/Amber/AmberData.py",
line 296, in ljParameters
<br>
p = ljp.getEntry(name)
<br>
File
"/afs/msc.cornell.edu/home/jsethna/ryang/installed/lib/python2.3/site-packages/MMTK/ForceFields/Amber/AmberData.py",
line 426, in getEntry
<br>
return self.entries[name]
<br>
KeyError: 'MG'
<br>
<br>
<pre wrap="">
<hr width="90%" size="4">
name = 'MG'
MG = Atom('Mg')
pdbmap = [('MG', {'MG': MG})]
amber_atom_type = {MG: 'MG'}
amber_charge = {MG: 2}
configurations = {'default': Cartesian({MG: (0, 0, 0)})}
</pre>
<pre wrap="">
<hr width="90%" size="4">
#info from <a class="moz-txt-link-freetext" href="http://www.webelements.com/webelements/elements/text/Mg/radii.html">http://www.webelements.com/webelements/elements/text/Mg/radii.html</a>
name = 'magnesium'
symbol = 'Mg'
mass = [(23.9850423, 78.99), (24.9858374, 10.00), (25.9825937, 11.01)]
color = 'gray'
vdW_radius = 0.173
</pre>
<pre wrap="">
<hr width="90%" size="4">
_______________________________________________
mmtk maillist - <a class="moz-txt-link-abbreviated" href="mailto:mmtk@starship.python.net">mmtk@starship.python.net</a>
<a class="moz-txt-link-freetext" href="http://starship.python.net/mailman/listinfo/mmtk">http://starship.python.net/mailman/listinfo/mmtk</a>
</pre>
</blockquote>
</body>
</html>