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.

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

proportional to the electric field at that site, Ei. The proportionality constant is<br />

the polarizability tensor, ai. The dipole feels an electric field both from the permanent<br />

charges of the system and from the other <strong>in</strong>duced dipoles. The expression<br />

for the l i is<br />

l i ¼ ai Ei ¼ ai E 0 i<br />

" #<br />

X<br />

j 6¼ i<br />

Tijl j<br />

where E0 is the field from the permanent charges. (There also may be permanent<br />

dipoles or other multipoles present contribut<strong>in</strong>g to E0 .) The <strong>in</strong>duced<br />

dipoles <strong>in</strong>teract through the dipole field tensor, Tij,<br />

0<br />

1<br />

Tij ¼ 1<br />

I<br />

r3 ½3Š<br />

3<br />

r5 x2 xy xz<br />

yx y2 yz<br />

zx zy z2 B<br />

C<br />

@<br />

A ½4Š<br />

where I is the identity matrix, r is the distance between i and j, and x, y, and z<br />

are the Cartesian components of the vector between i and j.<br />

The energy of the <strong>in</strong>duced dipoles, U <strong>in</strong>d, can be split <strong>in</strong>to three terms,<br />

U <strong>in</strong>d ¼ Ustat þ Umm þ U pol<br />

The energy Ustat is the <strong>in</strong>teraction energy of the N <strong>in</strong>duced dipoles with the<br />

permanent, or static, field<br />

Ustat ¼ XN<br />

i ¼ 1<br />

½5Š<br />

l i E 0 i ½6Š<br />

the energy Umm is the <strong>in</strong>duced dipole–<strong>in</strong>duced dipole <strong>in</strong>teraction<br />

Umm ¼ 1<br />

2<br />

and the polarization energy, U pol,<br />

X N<br />

X<br />

li Tij lj ½7Š<br />

i ¼ 1 j 6¼ i<br />

Upol ¼ 1<br />

2<br />

X N<br />

i ¼ 1<br />

l i Ei ½8Š<br />

is that required to distort the electron distribution to create the dipoles. 4,14<br />

Any polarizable model <strong>in</strong> which dipole moments, charges, or other multipoles

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

Saved successfully!

Ooh no, something went wrong!