sarkas.potentials.coulomb.coulomb_force

sarkas.potentials.coulomb.coulomb_force(r_in, pot_matrix)

Numba’d function to calculate the bare coulomb potential and force between two particles.

Parameters

• r_in (float) – Distance between two particles.

• pot_matrix (numpy.ndarray) – It contains potential dependent variables.

  Shape = (sarkas.core.Parameters.num_species, sarkas.core.Parameters.num_species, 3) .

Returns

• u_r (float) – Potential value.

• f_r (float) – Force between two particles.

Examples

>>> import numpy as np
>>> r = 2.0
>>> pot_matrix = np.array([ 1.0, 0.0, 0.0])
>>> coulomb_force(r, pot_matrix)
(0.5, 0.25)