Differences between revisions 12 and 13
Revision 12 as of 2010-10-20 09:11:16
Size: 2593
Editor: localhost
Comment: converted to 1.6 markup
Revision 13 as of 2010-11-28 13:32:18
Size: 2595
Deletions are marked like this. Additions are marked like this.
Line 76: Line 76:
.. _addhydrogen : inline : addhydrogen.py .. _addhydrogen : attachment:addhydrogen.py

ASE interface to the Molecular Modelling Toolkit (MMTK)

1   Introduction

The Molecular Modelling Toolkit (MMTK) by Konrad Hinsen is a powerfull program library for molecular simulation.

MMTK's functionality includes construction of molecular systems, including support for proteins and nucleic acids, dynamics, normal mode analysis and much more.

The ASE calculator interface enables the use of ASE dynamics/filters etc, with the MMTK force-fields, in this way the ASE functionality, like the Nudged Elastic Band method and Langevin dynamics, can be added to the MMTK functionality.

For general information on MMTK see the Molecular Modelling Toolkit homepage.

2   MMTK Calculator

The MMTK ASE calculator can not be attached to a general ListOfAtoms, but must use a ListOfAtom constructed from the MMTK atoms. The reason is that in constructing the force-fields additional information for each atoms must be present; a Carbon atom is not just a Carbon atom. This information is usually supplied in a PDB file.

Setting up the force-field starting from PDB files, the MMTK universe must be handle by MMTK commands.

The examples in the next section shows examples of setting up a MMTK universe.

Once the MMTK universe is constructed a ASE ListOfAtoms can be constructed and attached to a MMTK calculator:

atoms = MMTKListOfAtoms(universe)

3   Examples

3.1   Adding Hydrogen

In this example MMTK will be used to add missing Hydrogen for a experimental pdb file.

The example will used the dipeptide.pdb pdb file as an example.

The following MMTK python script will add the missing Hydrogen atoms for the dipeptide:


The added Hydrogen atoms are shown here (orange):


3.2   Minimizing using quasi-Newton

In this example the ASE quasi-Newton minimizer will be used to find the ground state for the dipeptide molecule:


ase2: MMTK (last edited 2010-11-28 13:32:18 by OleHolmNielsen)