\[ C_\alpha(r) = 3\frac{\mathrm{erfc}(\alpha r)}{r^5} + \frac{2\alpha}{\sqrt{\pi}}\left(2\alpha^2 + \frac{3}{r^2}\right) \frac{e^{-\alpha^2 r^2}}{r^2} \]
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=DP), | intent(in) | :: | alpha | |||
| real(kind=DP), | intent(in) | :: | r |