Let's keep going with what we did last time. It will get even more exciting than it already was.
Well I think it's nice to see what's behind some of those interactions and if I write it down here I'll find it again. But I admit that I usually read blogs with pictures rather than lots of weird signs.
The first expression was
Let us assume that 2 or more orbitals are different between the Slater determinants (pn-2). Then one of the terms in every product has to be 0 and the expression vanishes.
If p=n-1, then we have to claim the following:
This of course means that π=σ. There are now (n-1)! possible permutations that work with this. And the expression reduces to: (the n! cancels against the n! we started out with)
If p=n, we also need π=σ to have all the factors unequal to 0 (and equal to 1). The expression reduces to
The squared sign is of course 1. For every i there are (n-1)! peruations π with π(1)=i and ui=vi. Then we have
The second expression was
This time it turns out that one of the overlap integrals (and hence the whole product) equals zero if pn-3.
I'll only talk about the case p=n-2. If you want to have no overlap integral equal to 0, you first need
Aside from that you have to make sure that the mutually different orbitals are in the last integral. 1 and 2 have to be mapped onto n-1 and n.
As a next step you can acknowledge that π and σ have to either be the same or different by one transposition.
In the second case their signs are opposite. We get the following expression.
There are appropriately (n-2)!*2 permutations π that fulfill the requirements above and since the signs are opposite, this reduces to
The derivations seem to be similar for p=n-1 and p=n. As Levine tells me (who was in turn told by Parr) we get for p=n-1
and for p=n
(the n! is only there because I started with an extra n!)
The first integral is the Coulomb integral (with g12=1/r12). It corresponds to the Coulomb repulsion of two blurred electron clouds corresponding to the MOs. The second integral, called exchange integral, comes in because of the Slater determinant and is a consequence of the Pauli principle.
The next question is using the spin functions from two posts ago and see if some of the integrals vanish. And compare this to what you expected.
Boosting Molecular Dynamics with Socket-Based Communication
-
MD simulations can be 10x faster by replacing files with socket
communication. In brief: Our latest publication in the Journal of Physical
Chemistry Letter...
2 weeks ago
No comments:
Post a Comment