An normal mode analysis plugin based on basics objects of the OpenMM Simulation Module
This module uses classical way to analyze normal modes of a minimized protein structure. After minimization, the protein is assumed harmonically aprroximated. To obtain the hessian, structure corrdinates are tweak each steps and then the first and second derivatives of potential energy can be calculated. The algorithms are designed so that the sizes of tweaked positions would produce a small enough, and at a similar scale for different lengths of protein, forces such that the harmonic approximation of potential energy will be automatically guaranteed.
More details would be commited later...
- Two step (CPU, followed by GPU accelerated) energy minimization of protein structure before calculating normal modes
- Calculating demass-weighted hessian
- Obtaining normal modes and their corresponding vibrational frequencies
- Plotting vibrational spectrum to be compared with IR spectrum data
- Calculation of classical or semi-classical (QM formula) vibrational entropy
The OpenMM topology object
The OpenMM system object
The OpenMM integrator object
The N*3 array of positions (N is the number of atoms)
The CPU platform-specific properties used for generating simulation object in OpenMM
The CUDA platform-specific properties used for generating simulation object in OpenMM
To be added later...