This keyword specifies a Welch potential[69,70] of the form
$\displaystyle E$   $\displaystyle = \sum_{i<j}\Bigg[ {\displaystyle q_i q_j\over\displaystyle r_{ij...
-{\displaystyle q_j({\bf\mu}_i\cdot{\bf r}_{ji})\over\displaystyle r_{ij}^3}$  
    $\displaystyle -3{\displaystyle ({\bf\mu}_i\cdot{\bf r}_{ij})
...le r_{ij}^3}\Bigg]
+\sum_i {\displaystyle \mu_i^2\over\displaystyle 2\alpha_i},$  


$\displaystyle {\bf r}_{ij}^{\rm eff}={\bf r}_{ij}+{{\bf\mu}_i\over Q_i}-{{\bf\m...
...r Q_j},
\qquad r_{ij}^{\rm eff}=\left\vert{\bf r}_{ij}^{\rm eff}\right\vert. $

Welch parameters in atomic units should be entered after the Welch keyword in the order $ A_{++}$, $ A_{--}$, $ A_{+-}$, $ \rho$, $ Q_+$, $ Q_-$, $ \alpha_+$, $ \alpha_-$. The ions are then specified using atom types PL and MI for plus and minus, respectively, and can be in any order. There need not be equal numbers of positive and negative ions.

David Wales 2017-09-21