[MMTK] Creating arbitrary molecules
stoppeml at purdue.edu
Wed Jul 30 16:57:25 EDT 2003
I believe you need to add bonds between the atoms. You will have to dig
through the codebase a little. There is bond code in the Peptide chain
class that will give you a lead.
mike at stoppelman.com
> After a number of years of subscribing to the MMTK list, I've finally
> started to use MMTK for my day to day research. I'm tripping over a
> couple of minor bugs, and I was hoping that someone on the list could
> point out what I'm doing wrong.
> I can run several of the simple examples, particularly the ones that
> use molecules already in the database. However, when I try and create
> my own molecule from coordinates I run into trouble.
> Here's a trivial example:
> universe = InfiniteUniverse(Amber94ForceField())
> temperature = 298*Units.K
> integrator = VelocityVerletIntegrator(universe)
> If I replace the addObject calls with universe.molecule('water'), the
> example runs without problem. So I'm doing something wrong when I add
> the objects. Can anyone point out what I'm doing wrong?
> Rick Muller
> rmuller at sandia.gov
> mmtk maillist - mmtk at starship.python.net
More information about the mmtk