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.7 Theoretical and computing modeling 467<br />

8.7.1.2 Perturbation theories<br />

In perturbation theories 44,58 one relates the properties (e.g., the distribution functions or free<br />

energy) <strong>of</strong> the real system, for which the intermolecular potential energy is E, to those <strong>of</strong> a<br />

reference system where the potential is E0, usually by an expansion in powers <strong>of</strong> the perturbation<br />

potential E1=E-E0. The methods used are conveniently classified according to whether the reference system<br />

potential is spherically symmetric or anisotropic. Theories <strong>of</strong> the first type are most appropriate<br />

when the anisotropy in the full potential is weak and long ranged; those <strong>of</strong> the<br />

second type have a greater physical appeal and a wider range <strong>of</strong> possible application, but<br />

they are more difficult to implement because the calculation <strong>of</strong> the reference-system properties<br />

poses greater problems.<br />

In considering perturbation theory for liquids it is convenient first to discuss the historical<br />

development for atomic liquids.<br />

In many cases the intermolecular pair potential can be separated in a natural way into a<br />

sharp, short-range repulsion and a smoothly varying, long range attraction. A separation <strong>of</strong><br />

this type is an explicit ingredient <strong>of</strong> many empirical representations <strong>of</strong> the intermolecular<br />

forces including, for example, the Lennard-Jones potential. It is now generally accepted that<br />

the structure <strong>of</strong> simple liquids, at least <strong>of</strong> high density, is largely determined by geometric<br />

factors associated with the packing <strong>of</strong> the molecular hard cores. By contrast, the attractive<br />

interactions may, in the first approximation, be regarded as giving rise to a uniform background<br />

potential that provides the cohesive energy <strong>of</strong> the liquid but has little effect on its<br />

structure. A further plausible approximation consists <strong>of</strong> modeling the short-range forces by<br />

the infinitely steep repulsion <strong>of</strong> the hard-sphere potential. In this way, the properties <strong>of</strong> a<br />

given liquid can be related to those <strong>of</strong> a hard-sphere reference system, the attractive part <strong>of</strong><br />

the potential being treated as a perturbation. The choice <strong>of</strong> the hard-sphere fluid as a reference<br />

system is an obvious one, since its thermodynamic and structural properties are well<br />

known.<br />

The idea <strong>of</strong> representing a liquid as a system <strong>of</strong> hard spheres moving in a uniform, attractive<br />

potential well is an old one; suffice here to recall the van der Waals equation.<br />

Roughly one can thus regard perturbation methods as attempts to improve the theory <strong>of</strong> van<br />

der Waals in a systematic fashion.<br />

The basis <strong>of</strong> all the perturbation theories we shall consider is a division <strong>of</strong> the pair potential<br />

<strong>of</strong> the form<br />

( 12) ( 12) ( 12)<br />

v = v + w<br />

[8.83]<br />

0<br />

where v 0 is the pair potential <strong>of</strong> the reference system and w(12) is the perturbation. The following<br />

step is to compute the effect <strong>of</strong> the perturbation on the thermodynamic properties<br />

and pair distribution function <strong>of</strong> the reference system. This can be done systematically via<br />

an expansion in powers either <strong>of</strong> inverse temperature (the “λ expansion”) or <strong>of</strong> a parameter<br />

that measures the range <strong>of</strong> the perturbation (the “γ expansion”).<br />

In spite <strong>of</strong> the fact that hard spheres are a natural choice <strong>of</strong> the reference system, for the<br />

reasons discussed above, realistic intermolecular potentials do not have an infinitely steep<br />

hard core, and there is no natural separation into a hard-sphere part and a weak perturbation.<br />

A possible improvement is that to take proper account <strong>of</strong> the “s<strong>of</strong>tness”’ <strong>of</strong> the repulsive<br />

part <strong>of</strong> the intermolecular potential. A method <strong>of</strong> doing this was proposed by Rowlinson, 59

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

Saved successfully!

Ooh no, something went wrong!