htmd.adaptive.adaptiverun module#

class htmd.adaptive.adaptiverun.AdaptiveMD#

Bases: AdaptiveBase

Adaptive sampling using Markov state models for conformation selection.

Projects completed simulations according to the specified projection, clusters the projected data, and builds a Markov state model. From the model it selects conformations from the various states based on the chosen criteria to spawn the next epoch of simulations.

Parameters:

  • app (SimQueue object, default None) – A SimQueue class object used to retrieve and submit simulations

  • project (str, default 'adaptive') – The name of the project

  • nmin (int, default 0) – Minimum number of running simulations

  • nmax (int, default 1) – Maximum number of running simulations

  • nepochs (int, default 1000) – Stop adaptive once we have reached this number of epochs

  • nframes (int, default 0) – Stop adaptive once we have simulated this number of aggregate simulation frames.

  • inputpath (str, default 'input') – The directory used to store input folders

  • generatorspath (str, default 'generators') – The directory containing the generators

  • dryrun (boolean, default False) – A dry run means that the adaptive will retrieve and generate a new epoch but not submit the simulations

  • updateperiod (float, default 0) – When set to a value other than 0, the adaptive will run synchronously every updateperiod seconds

  • coorname (str, default 'input.coor') – Name of the file containing the starting coordinates for the new simulations

  • boxname (str, default 'input.xsc') – Name of the file containing the starting box dimensions for the new simulations. Set to ‘none’ to disable box writing.

  • lock (bool, default False) – Lock the folder while adaptive is ongoing

  • mps (int, default 0) – If mps > 0, it will run simulations using the Multi-Process Service (MPS) with the number of processes specified. If set to 0, mps is disabled

  • datapath (str, default 'data') – The directory in which the completed simulations are stored

  • filter (bool, default True) – Enable or disable filtering of trajectories.

  • filtersel (str, default 'not water') – Atom selection string for filtering. See more here

  • filteredpath (str, default 'filtered') – The directory in which the filtered simulations will be stored

  • projection (Projection object, default None) – A Projection class object or a list of objects which will be used to project the simulation data before constructing a Markov model

  • truncation (str, default None) – Method for truncating the prob distribution (None, ‘cumsum’, ‘statecut’

  • statetype ((``’micro’, ``'cluster', 'macro'), str, default 'micro') – What states (cluster, micro, macro) to use for calculations.

  • macronum (int, default 8) – The number of macrostates to produce

  • skip (int, default 1) – Allows skipping of simulation frames to reduce data. i.e. skip=3 will only keep every third frame

  • lag (int, default 1) – The lagtime used to create the Markov model. The units are in frames.

  • clustmethod (ClusterMixin class, default <class 'htmd.clustering.kcenters.KCenter'>) – Clustering algorithm used to cluster the contacts or distances

  • method (str, default '1/Mc') – Criteria used for choosing from which state to respawn from

  • ticalag (int, default 20) – Lagtime to use for TICA in frames. When using skip remember to change this accordinly.

  • ticadim (int, default 3) – Number of TICA dimensions to use. When set to 0 it disables TICA

  • contactsym (str, default None) – Contact symmetry

  • save (bool, default False) – Save the model generated

Examples

>>> adapt = AdaptiveMD()
>>> adapt.nmin = 2
>>> adapt.nmax = 3
>>> adapt.nepochs = 2
>>> adapt.ticadim = 3
>>> adapt.projection = [MetricDistance('name CA', 'resname MOL', periodic='selections'), MetricDihedral()]
>>> adapt.generatorspath = htmd.home()+'/data/dhfr'
>>> adapt.app = LocalGPUQueue()
>>> adapt.run()