class Jagla : public feasst::ModelTwoBody

The Jagla model as described in https://doi.org/10.1103/PhysRevE.98.012138 but with the variables renamed as follows: \(\epsilon_1 \rightarrow \gamma\), \(\epsilon_2 \rightarrow \epsilon\), \(\lambda_1 \rightarrow \lambda\), \(\lambda_2 \rightarrow r_c\).

\( U(r) = \left\{ \begin{array}{ll} \infty & : r < \sigma, \\ \frac{\gamma(\lambda - r)-\epsilon(r-\sigma)}{\lambda-\sigma} & : \sigma \le r \le \lambda, \\ -\frac{\epsilon(r_c - r)}{r_c-\lambda} & : \lambda \le r \le r_c, \\ 0 & : r \ge r_c. \end{array} \right. \)

Public Functions

double num_discretized_steps() const

Return the argument num_discretized_steps in read-only fashion.

void precompute(const ModelParams &existing)

Precompute model parameters based on existing model parameters.

void serialize(std::ostream &ostr) const

Output a serialized version of the existing model.


  • num_discretized_steps: convert the continuous ramp potential into a number of dicontinuous, discretized steps. See Fig. 1(d) of https://doi.org/10.1103/PhysRevE.98.012138 . If 0, then use the linear ramp potential (default: 0).