sp;&nbsp;&nbsp;&nbsp; <code><strong>P</strong> = V<sup>-1</sup> <strong>&sum;</strong> q<sub>i</sub> <strong>r</strong><sub>i</sub></code></p>

<p>où V est le volume de la boîte. Par défaut la polarisation totale est censée être en e/&Aring;<sup>2</sup>, mais là encore cela dépend des unités utilisées dans les fichiers d'entrée. Notez que si le système a une charge totale nulle alors le vecteur de polarisation est indépendant du point d'observation (qui est arbitrairement défini en (0,0,0) ici)