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 439<br />

We have here explicitly written the monomer wave functions, with a generic indica-<br />

A<br />

tion <strong>of</strong> their electronic definition: Ψ0 is the starting wave function <strong>of</strong> A, already indicated<br />

A<br />

with ΨA, ΨK is another configuration for A with one or two occupied one-electron orbitals<br />

replaced by virtual MOs.<br />

The first order contribution exactly corresponds to the definition we have done <strong>of</strong> ES<br />

with equation [8.16]: it is the Coulomb interaction between the charge distributions <strong>of</strong> A<br />

and B.<br />

The second term <strong>of</strong> the expansion, i.e., the first element <strong>of</strong> the second order contributions,<br />

corresponds to the polarization <strong>of</strong> A, due to the fixed unperturbed charge distribution<br />

<strong>of</strong> B; the next term gives the polarization <strong>of</strong> B, due to the fixed charge distribution <strong>of</strong> A. The<br />

two terms, summed together, approximate IND. We have already commented that perturbation<br />

theory in a standard formulation cannot give IND with a unique term: further refinements<br />

regarding mutual polarization effects have to be searched at higher order <strong>of</strong> the PT<br />

expansion.<br />

The last term <strong>of</strong> the second order contribution is interpreted as a dispersion energy<br />

contribution. The manifold <strong>of</strong> excited monomer states is limited here at single MO replacements<br />

within each monomer; the resulting energy contribution should correspond to a preliminary<br />

evaluation <strong>of</strong> DIS, with refinements coming from higher orders in the PT<br />

expansion.<br />

Remark that in PT methods, the final value <strong>of</strong> ΔE is not available. It is not possible here<br />

to get a numerical appraisal <strong>of</strong> correction to the values obtained at a low level <strong>of</strong> the expansion.<br />

All contributions are computed separately and added together to give ΔE.<br />

In conclusion, this PT formulation, which has different names, among which “standard”<br />

PT and RS (Rayleigh-Schrödinger) PT, gives us the same ES as in the variational<br />

methods, a uncompleted value <strong>of</strong> IND and a uncompleted appraisal <strong>of</strong> DIS: higher order PT<br />

contributions should refine both terms. One advantage with respect to the variational approach<br />

is evident: DIS appears in PT as one <strong>of</strong> the leading terms, while in variational treatments<br />

one has to do ad hoc additional calculations.<br />

Conversely, in the RS-PT formulation CT and EX terms are not present. The absence<br />

<strong>of</strong> CT contributions is quickly explained: RS-PT works on separated monomers and CT<br />

contributions should be described by replacements <strong>of</strong> an occupied MO <strong>of</strong> A with an empty<br />

MO <strong>of</strong> B (or by an occupied MO <strong>of</strong> B with an empty MO <strong>of</strong> A), and these electronic configurations<br />

do not belong to the set <strong>of</strong> state on which the theory is based. For many years the<br />

lack <strong>of</strong> CT terms has not been considered important. The attention focused on the examination<br />

<strong>of</strong> the interaction energy <strong>of</strong> very simple systems, such as two rare gas atoms, in which<br />

CT effects are in fact <strong>of</strong> very limited importance.<br />

The absence <strong>of</strong> EX terms, on the contrary, indicates a serious deficiency <strong>of</strong> the RS formulation,<br />

to which we have to pay more attention.<br />

The wave function ΦK used in the RS formulation does not fully reflect the electron<br />

permutational symmetry <strong>of</strong> the dimer: the permutations among electrons <strong>of</strong> A and B are neglected.<br />

This leads to severe inconsistencies and large errors when RS PT is applied over the<br />

whole range <strong>of</strong> distances. One has to rework the perturbation theory in the search <strong>of</strong> other<br />

0 0<br />

approaches. The simplest way would simply replace Φ0 =| ΨAΨB > with |AABΦ0 > where<br />

AAB is the additional antisymmetry operator we have already introduced. Unfortunately<br />

|AABΦ0 > is not an eigenfunction <strong>of</strong> H 0 as asked by the PT. Two ways <strong>of</strong> overcoming this difficulty<br />

are possible. One could abandon the natural partitioning [8.36] and define another

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

Saved successfully!

Ooh no, something went wrong!