Reviews in Computational Chemistry Volume 18
Reviews in Computational Chemistry Volume 18
Reviews in Computational Chemistry Volume 18
- 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