Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Interpolation mechanics

  • Explain what is interpolation mechanics

  • Deriving the Force-Field

  • Running a dynamics

Force field generation

We define first a VeloxChem Molecule for which the interpolation mechaics force field (IM-MM) will be constructed.

Loading...

We then initialize the QM-Driver determining the underlying quality of the interpolation data.

Define dyanmics and interpolation settings which are crucial for the quality in timing of the database. Important variables:

.desired_point_density -> how many datapoints should the database contain for each conformer.
.converged_cycle -> after how many cycles the construction for 1 conformer terminates (1 cycle: 3000 steps without database expansion).
.energy_thrsh -> allowed devaition of the IM energy and QM reference calculations.

As this notebook is intended to be a practival example on how to use the IMFroceFieldGenerator but keeping the computational effort at a reasoable level, the shown settings have been used. For production run we recommend using the following settings:

.nsteps = 50000
.snapshots = 10000
.desired_point_density = 80
.converged_cycle = 5
.energy_threshold = 2.0 kcal/mol
IMPORTANT: IM ForceFieldFile is initalized as im_database.h5

The compute function starts the construction of a new or given database. The construction run ends by confirming the qulity of the database and the result are stored.

Note that here we provide the constructed database (‘im_database.h5’) as the computational cost of the second derivatives are very demanding.

To use the constructed IM force field one has to initialize the InterpolationDriver() with user-defined settings (recommended). The current implementation is restricted to well seperated IM and MM regions and is not yet included into a linking atom framework.