Simulation¶

The simulation finally connects a Topology with an EntityHandle. While applying minimization or md tasks, the current positions of the simulation object can be mapped back to the attached structure at any time.

class Simulation(entity, settings)¶

First constructor that takes an ost.mol.EntityHandle as an input and automatically constructs an internal Topology according to the provided settings. Be aware, that the entity will be altered depending on hydrogen/termini definitions in the Forcefield attached to the settings

Parameters:	entity (EntityHandle) – Entity from which a default Topology will be created according to the settings settings (Settings) – Controls the parametrization of the Simulation and the buildup of the internal Topology
Raises:	RuntimeError when construction of Topology fails

class Simulation(topology, entity, settings)

Second constructor that takes a Topology, a consistent ost.mol.EntityHandle and a Settings as input.

Parameters:	topology (Topology) – Topology to initialize the Simulation entity (EntityHandle) – Entity that is consistent with the provided topology settings (Settings) – Controls the parametrization of the Simulation
Raises:	RuntimeException when number of atoms in entity is not equal the number of particles in topology

Save(filename, settings)¶

Dumps current status of the simulation with its connected entity onto disk. All informations from the entity except plain chains, residues, atoms and bonds are lost.

Parameters:	filename (str) – Filename

static Load(filename)¶

Loads a dumped simulation from disk. You have to make sure, that the provided settings are consistent with those from the saved simulation. Undefined behaviour must be expected otherwise.

Parameters:	filename (str) – Filename settings (Settings) – Settings to reparametrize loaded simulation

static IsPlatformAvailable(settings)¶

Returns:	True, if platform defined in settings is available. Otherwise, construction of a simulation object will fail with these settings.
Return type:	bool
Parameters:	settings (Settings) – Controls the parametrization of this class. Only Settings.openmm_plugin_directory, Settings.custom_plugin_directory and Settings.platform are relevant.

ApplyLBFGS([tolerance=1.0, max_iterations=1000])¶

Run minimization using the Limited-memory Broyden–Fletcher–Goldfarb–Shanno (L-BFGS) method.

Parameters:	tolerance (float) – Tolerance describing abort criteria of minimization. The tolerance is interpreted as the root mean square value of all force components. max_iterations (int) – Maximum number of minimization iterations

ApplySD([tolerance=100, max_iterations=1000])¶

Run minimization using the steepest descend method.

Parameters:	tolerance (float) – Tolerance describing abort criteria of minimization. The tolerance is interpreted as the norm of the maximum force acting onto the particles. max_iterations (int) – Maximum number of minimization iterations
Returns:	bool whether tolerance has been reached

Steps(steps)¶

Runs the simulation for steps number of steps.

Parameters:	steps (int) – Number of steps to perform

GetPositions([enforce_periodic_box = false, in_angstrom = true])¶

Parameters:	enforce_periodic_box (bool) – True if periodic box has to be enforced in_angstrom (bool) – True if positions should be in Angstrom, they’re in nm otherwise
Returns:	A ost.geom.Vec3List of all the positions of the atoms in the Entity.

SetPositions(positions[, in_angstrom = true])¶

Parameters:	positions (ost.geom.Vec3List) – Positions to be set in_angstrom (bool) – True if positions are in Angstrom, they’re expected to be in nm otherwise

UpdatePositions([enforce_periodic_box=false])¶

Sets the positions of the atoms in the EntityHandle to the current positions of the simulation.

Parameters:	enforce_periodic_box (bool) – True if periodic box should be enforced

GetVelocities()¶

Returns:	A Vec3List of all the velocities of the atoms in the Entity.

GetForces()¶

Returns:	A Vec3List of all the resulting forces for each atom in the Entity.

SetVelocities(velocities)¶

Parameters:	velocities (ost.geom.Vec3List) – Velocities to be set

GetEntity()¶

Returns:	EntityHandle attached to the simulation

GetEnergy()¶

Returns:	float Potential + Kinetic energy in kJ/mol of current simulation state

GetPotentialEnergy()¶

Returns:	float Potential energy in kJ/mol of current simulation state

GetKineticEnergy()¶

Returns:	float Kinetic energy in kJ/mol of current simulation state

GetTopology()¶

Returns:	Topology attached to the simulation

GetPeriodicBoxExtents()¶

Returns:	Vec3 Extents of the periodic box

SetPeriodicBoxExtents(vec)¶

Parameters:	vec (ost.geom.Vec3) – Periodic box extents in A

Register(observer)¶

Adds an observer to the simulation object. This can be used for example to write out the trajectory during the simulation.

Parameters:	observer (Observer) – Observer to be attached

AddPositionConstraint(index)¶

Fixes the position of the particle with index given in the argument. This requires to reinitialize the internal openmm Context (this is expensive!). Positions, velocities, forces, energies etc. are properly preserved but e.g. states of random number generators etc. might be lost.

Parameters:	index (int) – Particle to be fixed
Raises:	RuntimeError when index exceeds number of particles

AddPositionConstraints(indices)¶

Fixes the position of the atoms with the indices given in in the argument. This requires to reinitialize the internal openmm Context (this is expensive!). Positions, velocities, forces, energies etc. are properly preserved but e.g. states of random number generators etc. might be lost.

Parameters:	indices (list) – Particles to be fixed
Raises:	RuntimeError when an index in indices exceeds number of particles

ResetPositionConstraints()¶

Removes all position constraints. This requires to reinitialize the internal openmm Context (this is expensive!). Positions, velocities, forces, energies etc. are properly preserved but e.g. states of random number generators etc. might be lost.

ResetHarmonicBond(index, bond_length, force_constant)¶

Update of the harmonic bond parameters in the simulation and in the attached Topology

Parameters:	index (int) – Harmonic Bond to be reset bond_length (float) – New bond length in nm force_constant (float) – New force constant in kJ/mol/nm^2
Raises:	RuntimeError when index exceeds number of harmonic bonds

ResetHarmonicAngle()¶

Update of the harmonic angle parameters in the simulation and in the attached Topology

Parameters:	index (int) – Harmonic Angle to be reset angle (float) – New Angle in radians force_constant (float) – New force constant in kJ/mol/radian^2
Raises:	RuntimeError when index exceeds number of harmonic angles

ResetUreyBradleyAngle(index, angle, angle_force_constant, bond_length, bond_force_constant)¶

Update of the Urey-Bradley angle parameters in the simulation and in the attached Topology

Parameters:	index (int) – Urey-Bradley angle that has to be reset angle (float) – New angle in radians angle_force_constant (float) – New angle force constant kJ/mol/radian^2 bond_length (float) – New bond length in nm bond_force_constant (float) – New bond force constant in kJ/mol/nm^2
Raises:	RuntimeError when index exceeds number of urey-bradley angles

ResetPeriodicDihedral(index, multiplicity, phase, force_constant)¶

Update of the periodic dihedral parameters in the simulation and in the attached Topology

Parameters:	index (int) – Periodic dihedral to be reset multiplicity (int) – New Multiplicity phase (float) – New phase in radians force_constant (float) – New force constant in kJ/mol/radian^2
Raises:	RuntimeError when index exceeds number of periodic dihedrals

ResetPeriodicImproper(index, multiplicity, phase, force_constant)¶

Update of the periodic improper parameters in the simulation and in the attached Topology

Parameters:	index (int) – Periodic improper to be reset multiplicity (int) – New multiplicity phase (float) – New phase in radians force_constant (float) – New force constant kJ/mol/radian^2
Raises:	RuntimeError when index exceeds number of periodic impropers

ResetHarmonicImproper(index, phase, force_constant)¶

Update of the harmonic improper parameters in the simulation and in the attached Topology

Parameters:	index (int) – Harmonic improper to be reset phase (float) – New phase in radians force_constant (float) – New force constant kJ/mol/radian^2
Raises:	RuntimeError when index exceeds number of harmonic impropers

ResetLJPair(index, sigma, epsilon)¶

Update of the LJ pair parameters in the simulation and in the attached Topology

Parameters:	index (int) – LJ pair to be reset sigma (float) – New sigma in nm epsilon (float) – New epsilon in kJ/mol
Raises:	RuntimeError when index exceeds number of LJ pairs

ResetDistanceConstraint(index, constraint_length)¶

Update of the distance constraint parameters in the simulation and in the attached Topology. This requires to reinitialize the internal openmm Context (this is expensive!). Positions, velocities, forces, energies etc. are properly preserved but e.g. states of random number generators etc. might be lost.

Parameters:	index (int) – Distance constraint to be reset constraint_length (float) – New constraint length in nm
Raises:	RuntimeError when index exceeds number of distance constraints

ResetHarmonicPositionRestraint(index, ref_position, force_constant, x_scale, y_scale, z_scale)¶

Update of the harmonic position restraint parameters in the simulation and in the attached Topology

Parameters:	index (int) – Harmonic position restraint to be reset ref_position (ost.geom.Vec3) – New reference position in nm force_constant (float) – New force constant in kJ/mol/nm^2 x_scale (float) – New x-scale y_scale (float) – New y-scale z_scale (float) – New z-scale
Raises:	RuntimeError when index exceeds number of harmonic position restraints

ResetHarmonicDistanceRestraint(index, length, force_constant)¶

Update of the harmonic distance restraint parameters in the simulation and in the attached Topology

Parameters:	index (int) – Harmonic distance restraint to be reset length (float) – New restraint length in nm force_constant (float) – New force constant kJ/mol/nm^2
Raises:	RuntimeError when index exceeds number of harmonic distance restraints

ResetLJ(index, sigma, epsilon)¶

Update of the lj parameters in the simulation and in the attached Topology

Parameters:	index (int) – LJ to be reset sigma (float) – New sigma in nm epsilon (float) – New epsilon in kJ/mol
Raises:	RuntimeError when index exceeds number of particles

ResetGBSA(index, radius, scaling)¶

Update of the GBSA parameters in the simulation and in the attached Topology

Parameters:	index (int) – : GB to be reset radius (float) – New radius in nm scaling (float) – New OBC scaling
Raises:	RuntimeError when index exceeds number of particles

ResetCharge(index, charge)¶

Update of the charge in the simulation and in the attached Topology

Parameters:	index (int) – Charge to be reset charge (float) – New Charge
Raises:	RuntimeError when index exceeds number of particles

ResetMass(index, mass)¶

Update of the mass in the simulation and in the attached Topology This requires to reinitialize the internal openmm Context (this is expensive!). Positions, velocities, forces, energies etc. are properly preserved but e.g. states of random number generators etc. might be lost.

Parameters:	index (int) – Mass to be reset mass (float) – New mass
Raises:	RuntimeError when index exceeds number of particles