Potential modifications

hard_core

Add a hard core of radius \(r_\mathrm{core}\) to point particles by shifting the potential radially outwards:

\[\tilde U(r) = U(r - r_\mathrm{core}), \qquad r > r_\mathrm{core}.\]

Example:

potential = potential:modify({"hard_core", radius = 0.5})