52 throw std::domain_error(
"Parameter 'kappa' must be >= 0");
55 throw std::domain_error(
"Parameter 'r_cut' must be >= 0");
78 double const dist)
const {
83 auto fac =
prefactor * q1q2 / Utils::int_pow<3>(dist);
86 fac *= std::exp(-
kappa * dist) * (1. +
kappa * dist);
95 double pair_energy(
double const q1q2,
double const dist)
const {
103 energy *= std::exp(-
kappa * dist);
Vector implementation and trait types for boost qvm interoperability.
void set_prefactor(double new_prefactor)
double prefactor
Electrostatics prefactor.
void sanity_checks_charge_neutrality() const
This file contains the defaults for ESPResSo.
double kappa
Ionic strength.
void on_node_grid_change() const
void on_activation() const
void on_periodicity_change() const
double r_cut
Interaction cutoff.
double pair_energy(double const q1q2, double const dist) const
Compute the pair energy.
void on_boxl_change() const
DebyeHueckel(double prefactor, double kappa, double r_cut)
void on_cell_structure_change() const
Utils::Vector3d pair_force(double const q1q2, Utils::Vector3d const &d, double const dist) const
Compute the pair force.
void sanity_checks() const