28.02.2013 Views

Handbook of Solvents - George Wypych - ChemTech - Ventech!

Handbook of Solvents - George Wypych - ChemTech - Ventech!

Handbook of Solvents - George Wypych - ChemTech - Ventech!

SHOW MORE
SHOW LESS

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

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

8.4 Two-body interaction energy 447<br />

Limitations due to the non-convergence <strong>of</strong> these expansion series are reduced by the<br />

use <strong>of</strong> multicenter expansions: in such a way it is possible to reach a satisfactory representation<br />

even for polyatomic molecules at short distances using the first expansion term only.<br />

In general, the expansion sites are the (heavy) atoms <strong>of</strong> the dimer, using a simpler formulation<br />

due to London. The London formulation uses the averaged (isotropic) value <strong>of</strong> the<br />

static dipole polarizabilities, α j , <strong>of</strong> the atoms and their first ionization potentials I i and I j:<br />

with<br />

6 ( ) ∑ ∑ 6(<br />

)<br />

− −<br />

DIS( 6) ≈ C6 AB R = C i, j Rij 3<br />

C6( i, j)<br />

=<br />

2<br />

α α II<br />

I + I<br />

i j i j<br />

i j<br />

A<br />

i<br />

B<br />

j<br />

6<br />

[8.56]<br />

[8.57]<br />

This formula derives from the standard RS-PT with the expansion over the excited state<br />

truncated at the first term. For isolated atoms i and j, the coefficients C6(i,j) can be drawn<br />

from experimental data, while for atomic fragments <strong>of</strong> molecules only from computations.<br />

The exchange (or repulsion) term<br />

The exchange terms are related to the introduction <strong>of</strong> the intermonomer antisymmetrizer<br />

AAB in the expression giving the electrostatic contribution. The correction factors introduced<br />

in the SAPT methods 15 may be related to a renormalization factor that may be written<br />

in the following form when AAB is replaced with 1+P for simplicity:<br />

1 1 1<br />

= =<br />

Φ| PΦΨΨ| PΨΨ P<br />

0 0<br />

Α<br />

B A B<br />

( 1+<br />

)<br />

[8.58]<br />

This factor is different from 1, because contains all the multiple overlap values between<br />

MOs <strong>of</strong> A and B. may be then expanded into terms containing an increasing number <strong>of</strong><br />

multiple overlaps (or exchanges <strong>of</strong> orbitals):<br />

P = P + P + P + P + � [8.59]<br />

1 2 3 4<br />

The exchange energy may be expanded into increasing powers <strong>of</strong> MO overlap<br />

integrals S:<br />

with<br />

( ) ( ) ( )<br />

EX = EX S + EX S + EX S + � [8.60]<br />

( )<br />

2 3 4<br />

EX S = VP − P<br />

[8.61]<br />

2 2 2<br />

This first term is sufficient for modeling intermolecular potentials. The overlaps depend<br />

on the monomer separation roughly as exp(-αR) so the following terms give small<br />

contribution at large-medium distances. At short distances the positive exchange term rapidly<br />

increases. Reasons <strong>of</strong> uniformity with the other terms <strong>of</strong> the interaction potential sug-

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

Saved successfully!

Ooh no, something went wrong!