Struct Yukawa 

pub struct Yukawa { /* private fields */ }

Screened Coulomb (Yukawa / Debye–Hückel) interaction U = k q_i q_j e^{−r/λ} / r.

Implementations

impl Yukawa

pub fn new( charge_i: Fix128, charge_j: Fix128, screening_length: Fix128, ) -> Result<Self, PairPotentialError>

Charges in coulombs with COULOMB_CONSTANT and screening length λ > 0.

Errors

PairPotentialError::NonPositiveScreeningLength.

pub fn with_constant( coulomb_constant: Fix128, charge_i: Fix128, charge_j: Fix128, screening_length: Fix128, ) -> Result<Self, PairPotentialError>

Charges with an explicit constant k > 0 and screening length λ > 0.

Errors

PairPotentialError::NonPositiveCoulombConstant, PairPotentialError::NonPositiveScreeningLength.

Trait Implementations

impl Clone for Yukawa

fn clone(&self) -> Self

Returns a duplicate of the value.

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

Performs copy-assignment from source.

impl Copy for Yukawa

impl Debug for Yukawa

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

Formats the value using the given formatter.

impl Eq for Yukawa

impl PairPotential for Yukawa

fn energy(&self, r: Fix128) -> Result<Fix128, PairPotentialError>

Potential energy U(r) at distance r > 0.

fn force(&self, r: Fix128) -> Result<Fix128, PairPotentialError>

Radial force F(r) = −dU/dr at r > 0 (positive = repulsive).

fn force_on_first(&self, d: Vec3Fix) -> Result<Vec3Fix, PairPotentialError>

Force on the first particle of a pair separated by d = x_i − x_j: F(|d|) d / |d|. The second particle receives the negative.

impl PartialEq for Yukawa

fn eq(&self, other: &Self) -> bool

Equality operator ==.

fn ne(&self, other: &Rhs) -> bool

Inequality operator !=.

impl StructuralPartialEq for Yukawa