Chirality potential

namespace mdk
class Chirality : public mdk::Force
#include <Chirality.hpp>

Chirality potential for the models with a native structure (in particular the positions of the residues). It’s supposed to be a more replacement for bond and dihedral potentials, if the native structure is provided.

Public Functions

virtual void bind(Simulation &simulation) override

Bind the force field to the simulation.

Parameters

simulationSimulation to bind to.

virtual void asyncPart(Dynamics &dynamics) override

Asynchronous part of the computation.

Parameters

dynamicsDynamics object to add potential energy and forces to.

Public Members

double e_chi = 1.0 * eps

The amplitude of the potential.

Private Members

Scalars d0_cube_inv

A list of cube inverses of $d_0$, where $d_0$ is the length of $v_i = r_{i+1} - r_i$ in the native structure.

Scalars C_nat

A list of the values of C_i in the native structure.

Bytes inRange

inRange[i] = 1 if a quadruple (i-2, i+1, i, i+1) is connected, i.e. in a single chain.