Struct Torsion 

pub struct Torsion;

Periodic torsion potential for dihedral angles.

Physics

Models the rotational barrier around a bond axis using a periodic cosine function.

• Formula: $$ E = \frac{1}{2} V [1 - \cos(n(\phi - \phi_0))] $$
• Derivative (diff): $$ T = \frac{dE}{d\phi} = \frac{1}{2} V \cdot n \cdot \sin(n(\phi - \phi_0)) $$

Parameters

• v_half: Half barrier height $V_{half} = V/2$.
• n: Periodicity/multiplicity.
• cos_n_phi0: $\cos(n\phi_0)$, pre-computed phase cosine.
• sin_n_phi0: $\sin(n\phi_0)$, pre-computed phase sine.

Inputs

• cos_phi: Cosine of the dihedral angle $\cos\phi$.
• sin_phi: Sine of the dihedral angle $\sin\phi$.

Implementation Notes

• Uses optimized closed-form formulas for common periodicities ($n = 1, 2, 3$).
• Falls back to Chebyshev recurrence for higher periodicities.
• All intermediate calculations are shared between energy and torque computations.
• Branchless and panic-free.

Trait Implementations

impl Clone for Torsion

fn clone(&self) -> Torsion

Returns a duplicate of the value.

fn clone_from(&mut self, source: &Self)

Performs copy-assignment from source.

impl Debug for Torsion

fn fmt(&self, f: &mut Formatter<'_>) -> Result

Formats the value using the given formatter.

impl Default for Torsion

fn default() -> Torsion

Returns the “default value” for a type.

impl<T: Real> TorsionKernel<T> for Torsion

fn energy( cos_phi: T, sin_phi: T, (v_half, n, cos_n_phi0, sin_n_phi0): Self::Params, ) -> T

Computes only the potential energy.

Formula

$$ E = V_{half} [1 - \cos(n(\phi - \phi_0))] $$

fn diff( cos_phi: T, sin_phi: T, (v_half, n, cos_n_phi0, sin_n_phi0): Self::Params, ) -> T

Computes only the torque $T$.

Formula

$$ T = V_{half} \cdot n \cdot \sin(n(\phi - \phi_0)) $$

This factor allows computing forces via the chain rule: $$ \vec{F} = -T \cdot \nabla \phi $$

fn compute( cos_phi: T, sin_phi: T, (v_half, n, cos_n_phi0, sin_n_phi0): Self::Params, ) -> EnergyDiff<T>

Computes both energy and torque efficiently.

This method reuses intermediate calculations to minimize operations.

type Params = (T, u8, T, T)

Associated constants/parameters required by the potential (e.g., $V, n$).

impl Copy for Torsion