19.02.2013 Views

Reviews in Computational Chemistry Volume 18

Reviews in Computational Chemistry Volume 18

Reviews in Computational Chemistry Volume 18

SHOW MORE
SHOW LESS
  • No tags were found...

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

108 Polarizability <strong>in</strong> Computer Simulations<br />

U(q)−E 0 (eV)<br />

15<br />

10<br />

5<br />

0<br />

−5<br />

−1 0<br />

q (e)<br />

Figure 3 Energy versus charge for chlor<strong>in</strong>e (solid l<strong>in</strong>e) and sodium (dashed l<strong>in</strong>e). The<br />

l<strong>in</strong>es are a quadratic fit through the energies of the ions relative to the neutral atom.<br />

from the Coulomb <strong>in</strong>teraction between delocalized charge distributions rðrÞ,<br />

rather than po<strong>in</strong>t charges,<br />

ð<br />

riðriÞrjðrjÞ JijðrÞ ¼<br />

jri rj rj dri drj ½37Š<br />

The charge distributions are frequently assumed to be spherical, for simplicity.<br />

125–128 Directional <strong>in</strong>teractions can be <strong>in</strong>corporated with nonspherical<br />

charge distributions, at some added computational expense. 129–131<br />

The partial charges on each atom of the molecule are found by m<strong>in</strong>imiz<strong>in</strong>g<br />

the energy, subject to a constra<strong>in</strong>t that the total charge is conserved.<br />

X<br />

½38Š<br />

i<br />

qi ¼ qtot<br />

The charge conservation constra<strong>in</strong>t can be enforced us<strong>in</strong>g an undeterm<strong>in</strong>ed<br />

multiplier,<br />

UðqÞ ¼UðqÞ l X<br />

!<br />

½39Š<br />

i<br />

qi qtot<br />

M<strong>in</strong>imiz<strong>in</strong>g this expression for the energy with respect to each of the qi<br />

under the assumption that the molecule <strong>in</strong> question is neutral ðqtot ¼ 0Þ gives<br />

1

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!